4.1. Soil Geomechanical Properties and Their Influence on Slope Stability
The bivariate correlation of soil geomechanical properties was analyzed by computing the Pearson correlation coefficient and using the textural and Casagrande diagrams. Geomechanical properties determined at the Sol Solution Afrique Centrale Laboratory, Yaoundé Cameroon (PT01, PT02, PT03, PT04, and PT05 samples) and the Engineering Geology Department of the Technische Universität Berlin, Germany (PT06, PT07, PT08, PT09, PT10, and PT11 samples) are presented in
Table 2, with PT meaning a pit from which soils samples were removed.
These results, obtained in both laboratories, are homogeneous with minor variations among each group of values. Soil sampling depths in this study fluctuated between 0.25–7.3 m, corresponding to the average depth of landslide slip surfaces recorded at the western flank of Mount Oku. The degree of linear relationship between them pairwise was determined by calculating the Pearson correlation coefficient. The Pearson correlation matrix displayed in
Table 3 shows that these properties are all linearly interconnected. In this table, there are one strong negative correlation ≤ 0.7, 16 moderate negative correlations ≤ 0.4, 10 negligible/weak negative correlations (between −0.40), 14–negligible/weak positive correlations (between 0–0.4), 10 moderate positive correlations > 0.4, and two strong positive correlations > 7, based on the classification proposed by [
52]. It is also noted that relationships between bulk and particle densities, sand, silt and clay contents, plasticity index, friction angle, and cohesion are stronger with many properties compared to porosity, water content, and methylene blue value, as shown by their high correlation coefficients.
Additionally, bulk and particle densities, porosity, fine particle content (clay and silt), cohesion, and friction angle have positive correlation coefficient signs, supposing that they increase or decrease together. Contrarily, plasticity index, methylene blue values, water content, and sand percentages are supposed to vary in the opposite direction with cohesion, regarding their correlation coefficient signs, which are negative. When they are increasing, cohesion is supposed to decrease and the liquefaction potential increases, leading to landslides. Additionally, sandy soils exhibit no cohesion [
66], but they display a high liquefaction potential. Furthermore, Ref. [
67] in their study concerning the influence of ants on soil and water losses in eastern Spain, they found a reduction in soil bulk density and an increase in macropore flow in ant-affected soils, which makes them prone to landslides. This supports our observations of densities being interrelated with soil grain sizes, water content, porosity, consistency limits, absorption capacity, cohesion, and friction angle. Therefore, it can be concluded that landslides at Mount Oku are closely influenced by their soil geomechanical properties.
Although the linear correlation is strong between most of these properties, there could also be a positive or negative nonlinear, monotonic relationship [
67,
68]. Therefore, other correlation approaches are needed to confirm these soil properties correlations before introducing them in landslide-causative factors systems. In addition to the Pearson correlation coefficient, many authors investigated the relationships among soil properties, as presented below.
Moreover, soil samples from Mount Oku exhibited water amounts that exceed their plasticity index values, as already noticed above. In this state, water exerts pressure on soil pores, decreasing the friction just as the shear strength responsible for the material stability is also reduced. This points out a possible initiation mechanism of slope instabilities that can trigger landslides, if they are combined with other factors, such as rainfall and steepness of the slope [
69,
70,
71]. Additionally, water causes a decrease in shear strength either by reducing the apparent soil cohesion or by creating or extending cracks, which represent potential slip surfaces when moistened. This is directly related to intense or long-lasting rainfall events, as also shown by [
72].
The proportions of clay, sand, and silt particles presented in
Table 2 were plotted on the texture diagram, as shown in
Figure 6. It can be seen that nonfailure sites correspond to PT01, PT02, PT06, PT07, PT08, PT09, PT10, and PT11 samples and failure or landslide sites are locations of PT03, PT04, and PT05 samples.
The samples from landslide sites are rather in the sandy regime, while the samples from the nonlandslide sites have higher fine particle contents, since failure phenomena depend greatly on the grain size, as concluded by [
73] after their investigation of the pore-pressure generation and movement of rainfall-induced landslides in laboratory flume tests.
Figure 6.
Texture triangle showing proportions of clay, sand, and silt particles of samples at nonfailures (PT01, PT02, PT01, PT02, PT06, PT07, PT08, PT09, PT10, and PT11) and failures, or landslides, sites (PT03, PT04, and PT05); modified from [
74,
75].
Figure 6.
Texture triangle showing proportions of clay, sand, and silt particles of samples at nonfailures (PT01, PT02, PT01, PT02, PT06, PT07, PT08, PT09, PT10, and PT11) and failures, or landslides, sites (PT03, PT04, and PT05); modified from [
74,
75].
Furthermore, liquid limit and plasticity index values of Mount Oku soils reported on the Casagrande diagram describe a cloud of points with low dispersion (
Figure 7). This diagram shows that the samples PT03, PT06, PT07, and PT10 are highly plastic with a great swelling potential. PT01, PT02, PT04, PT05, PT08, PT09, and PT11 are very plastic with excessive swelling potential.
In view of all this, soil samples were finally classified as highly plastic clays with a high to very high swelling behavior, showing that they can experience swelling and significant shrinkage in the presence or absence of water [
76,
77]. Swelling and shrinkage alternation can create weaknesses in these soil textures, easing their failure. Similar observations have been made by [
9,
77] in their investigations concerning slope shear and strength deterioration through drying–wetting succession and influence of synthetic wick fibers on clayey soils’ behavior, respectively.
It, therefore, follows that the interdependence of soil geomechanical properties is directly related to landslide occurrences. This is supported by the statistical correlation of these properties with landslides as displayed below.
4.2. Spatial Analysis of Soil Geomechanical Properties with the Fuzzy Membership Approach
The average values of soil geomechanical properties were computed for each rock type, as can be seen in
Table 4. These soils on basalt, trachy-rhyolite, and migmatite displayed almost similar average values of geomechanical properties. Particle and bulk densities varied between 2.59–2.66 g/cm
3 and 1.15–1.58 g/cm
3, respectively. Water content and porosity varied from 39.3–44.1% and 43–45%, respectively. The differences in sand, silt, and clay content among soils are less than 25%. Soils on basalt and trachy-rhyolite display the most similar geomechanical properties. Soil on migmatite displays the highest bulk and particle densities, porosity, silt, clay, friction angle, and cohesion mean values. Moreover, they show the lowest water content, sand, and plasticity limit mean values. Rhyolites were assumed, based on literature, to have a bulk density of 26.5 g/cm
3, cohesion of 1000 kPa, and 46° as friction angle [
78]. They display no porosity, individual grain sizes, plasticity limits, water content, and methylene blue values.
Many authors have shown that the increasing or decreasing tendency of soil geomechanical properties can create weaknesses in these soil textures, easing their failure or sliding [
7,
69,
74]. In other words, some of these geomechanical properties should be decreasing to make the site more susceptible to landslides. This is the case for bulk and particle densities, internal angle of friction, and cohesion (
Table 5). The fuzzy membership values of bulk density classes vary from 0–0.99, while those of particle density, cohesion, and friction angle fluctuate between 0–1.
However, in other cases, rising of soil geomechanical properties increases the probability of landslide occurrence. This is the case for porosity, water content, MBV, PI, sand, and clay content (
Table 6). The fuzzy membership values of porosity, water content, MBV, PI, sand, and clay contents classes vary between 0–1.
4.3. Bivariate Correlations of Geo-Environmental Factors with Landslides
Slope angle and aspect, land use, elevation, distances from the main road, lithology, and distances from the major stream maps are presented in
Figure 8.
(a) Slope Angle
Slope angle, which is the inclined surface of the land, is among the most determining factors of landslide types and velocities [
18]. In this area, very gentle slopes (0–15°) and gentle slopes (15–25°) occupy the largest surface of this area, but they have less influence on landslide occurrences, while slope angles ranging from 24–31° and >31° (
Figure 8a), respectively, present the highest positive weights of 0.18 and 0.73, correspondingly (
Table 7a and
Figure 9). Their highest positive weights can be due to interactions between slope steepness and other factors, such as shear stress produced by the material weight and road and river cuttings, which modify the natural slopes and destabilize the block material [
10,
11].
(b) Slope Aspect
Slope aspect (
Figure 8b) is the direction toward which the surface of the slope faces. At Mount Oku, most slopes are directed toward north (0–22.5°), southwest, west, northwest, and north (337.5–360°), as shown in
Table 7 and
Figure 9, and are more influential with, respectively, 0.05, 0.12, 0.26, 0.23, and 0.20 as weights, while others are negative. Indeed, meteorological events in Mount Oku are likely to be more intense mainly on northern-, southwestern-, western-, and northwestern-oriented slopes. This can be attributed to the fact that slope aspect, vegetation cover, soil water retention, rainfall depending on the wind direction, and sunshine intensity are interrelated, as shown by [
79] in their investigation of the effect of elevation and aspect on wind temperature and humidity. The low intensity of sunshine, intensified soil moisture, and weathering on these slopes is leading to landslides events [
11,
80,
81].
(c) Elevation
The elevation classes (
Figure 8c) between 1574–1731 m and 1731–1880 m display the highest weight values of 0.24 and 0.25 each, as shown in
Table 7a and
Figure 9. Elevation is an indirect landslide factor as it has an influence on other factors, such as rainfall, temperature, soil development, and vegetation [
82,
83].
(d) Land Cover
Land cover of the study area (
Figure 8d) has been categorized into four classes, namely, barren land, shrub land, shrub with emergent trees, and forest. The forest areas have the highest weight values of 0.22 (
Figure 9 and
Table 7a). The positive weight displayed by the forest cover can be explained by the fact that in Mount Oku the population used to stabilize landslide scars by planting trees, such as eucalyptus, which later grew into forest. Moreover, in this area, the natural forest plays a role as a stabilizer for steep slopes through promotion of infiltration and drainage. Similar conclusions were made by [
81] concerning landslide-susceptibility mapping with information value and logistic regression methods in the Bailongjiang watershed in China. Similarly, Ref. [
72] established that the amount of water entering a slope depends on several geomorphological factors, anthropogenic activities, and atmospheric conditions including vegetation type, drainage, soil type, and rock structure.
(e) Proximity to Major Rivers
On the western flank of Mount Oku, most landslides are found at distances between 400–600 m, 600–800 m, 800–1000 m, 1200–1400 m, and >1400 m from rivers (
Figure 8e). These represent the highest information values of 0.14, 0.12, 0.20, 0.54, and 1.87, respectively, while others are negative (
Figure 9 and
Table 7b). Deep, incised river channels can modify the natural slope at the toe of mountains through water erosion. It also influences groundwater level fluctuation, which is related to intensive soil wetting and drying phenomenon that lead to slope instability [
81]. Here, however, there is no clear relationship between proximity to rivers and landslide occurrence and the classes with high weight values seem relatively random. Especially the high weight values in the classes >1200 m could be related to the small class size, an effect that could be considered to be an artifact.
(f) Proximity to Main Roads
Major landslide scars are found at 1000–1200 m and >1400 m from the roads, with weight values of 0.13 and 0.45, respectively, and highly negative values closer to the roads (
Table 7b,
Figure 8f and
Figure 9). Usually, the presence of roads is believed to trigger landslides through undercutting of slopes. As a result of an increase in stress on the back of the slope, due to changes in topography and decrease of load on the toe, tension cracks may develop [
59,
80]. Here, however, the high weight values observed in the classes far away from the roads, at 1000–1200 m and >1400 m, and negative values close to roads seem, again, quite random.
This effect can be explained by the fact that these landslides close to rivers and roads are of rather small sizes. Therefore, even if they are many landslides, they could be regarded as insignificant at the scale of the study area.
(g) Curvature
Curvature allows identifying the shape of the slope. It plays an important role in erosion and deposition processes by defining the convergence and divergence of water flow. Curvature classes ranging from −314 to −6°, −6 to −3°, 2 to 4°, 4 to 7°, 7 to 13°, and >13° have the highest information values of 0.82, 0.20, 0.05, 0.42, 0.73, and 1.08, each (
Table 7b,
Figure 8g and
Figure 9). As presented by [
84,
85], the landslide movement direction together with driving and resisting stresses along the failure slope are influenced by curvature, since it controls the speed and convergence or divergence of landslide-displaced material and water flowing down the slope.
(h) Lithology
The study area is covered by volcanic rocks, precisely basalt (highly weathered with thick, residual soil), rhyolite (slightly weathered, no residual soil, and indicated by an arrow on the map), trachy-rhyolite (moderately weathered with very steep slopes), and migmatites (highly weathered with moderately steep slopes). Trachy-rhyolite covers 37%, while basalt covers 54% of the study area (
Figure 8h). Trachy-rhyolite displays the highest information values of 0.27. Migmatite, basalt, and rhyolite present negative information values of −1.66, −0.23, and −3, respectively (
Figure 9 and
Table 7b). This is probably due to variations in thickness, steepness of slopes, and strength of soils developed from the weathering of these rocks. This has also been noted by [
18,
86] in their research on soft rock mass-weathering effect on slope stability and debris flow susceptibility assessment in Subao river valley, respectively.
Geo-environmental factor classes, with their corresponding area percentages and weights, calculated using the information value technique, are presented in
Table 7a,b and
Figure 9. Some classes of these features have the highest positive weight values, demonstrating their higher landslide prediction ability in this zone, as also noticed by [
23] in their work on landslide-susceptibility mapping on the Bamenda mountain. The sum of each geo-environmental factor-positive class weights show that lithology is the dominant landslide factor in the western flank of Mount Oku, followed, respectively, by slope angle, curvature, land use, aspect, proximity to road, proximity to rivers, and elevation. However, it was difficult to determine which geomechanical properties were the most important, since their spatial distribution was identical to rock type.
4.4. Landslide-Susceptibility Model Results and Discussion
In order to thoroughly evaluate the significance of geo-environmental factors and soil geomechanical properties, three landslide-susceptibility models were established with the first one merging only soil properties, geo-environmental factors for the second, and the last one combining all geo-environmental factors and soil properties. These models present areas of identified landslides and areas with similar predisposing conditions, where landslides have not yet been experienced, as also noticed by [
18,
19], among others. The resulting landslide-susceptibility map indexes have been normalized and classified using the equal interval method, so the results could easily be compared with each other. The landslide-susceptibility maps were classified into four susceptibility areas: High, moderate, low, and very low (
Figure 10). The very-low-susceptibility class ranged between 0–0.25, the low between 0.25–0.50, the moderate from 0.5–0.75, and the high one from 0.75–1. For the model with only geo-environmental factors, 6% of the study area was highly susceptible, 69% was moderately susceptible, 22% was lowly susceptible, and 3% was very lowly susceptible (
Figure 10 and
Figure 11). The landslide-susceptibility model results presented here are justified based on established methods and existing literature on geomechanical properties, relationships among the geomechanical properties, and geo-environmental and geomechanical factors in landslide-susceptibility assessment.
The model based only on soil properties indicates 99% to be highly susceptible, 0% moderately susceptibility, 0% lowly susceptible, and 1% to be very lowly susceptible. The combined model shows most of the study area, 67%, to be highly susceptible, 30% moderately susceptible, 2% lowly susceptible, and 1% to be very lowly susceptible (
Figure 11). Basalt and trachyte are highly weathered with thick, residual, clayey soils that are highly plastic and, hence, sensitive to the variations of water content. Soils on trachyte and basalt may rapidly pass from liquid to plastic or solid state, predisposing this area to landslides, when combined with steep slopes. This makes these soils more susceptible to landslides than shallow, residual soils developed on migmatites, which is reflected by the model. Furthermore, soils on migmatites display relatively gentle slopes, moderately weathered with shallower, residual soil, low porosity, and high cohesion values, making them less susceptible to landslides. Besides this, Ref. [
86,
87] stated that when the water content values vary between 30 and 40%, the degree of swelling is medium, as it was the case in this work with soils on migmatite. This soil also presents a high cohesion value, which means that the grains are strongly cemented.
A slight variation was observed between the susceptibility class percentages of the model with only geo-environmental factors and those of the model merging geo-environmental factors and soil properties. This is mostly because the weights of geomechanical properties obtained with the fuzzy membership method are small (from 0–0.74). These models allow observing some variations between the susceptibility class percentages. A distinctive maximum of the high susceptibility class percentage was observed in the model combining only soil geomechanical properties, while the minimum was displayed by the geo-environmental factors model. Moreover, the lowest percentages of very low and low susceptibilities were observed in the soil properties and the soil-geo-environment models. In other words, the susceptibility class percentages of the model with only geo-environmental factors and the one with only soil properties presented either the lowest or the highest susceptibility class percentages (
Figure 11). The combined model of soil properties and geo-environmental factors tended to be the most reasonable and stable. When soil properties or geo-environmental factors were used alone, the resulting landslide-susceptibility model tended to misjudge the susceptibility degree. The soil, which is the material concerned by slides, allows a more precise examination of landslide predisposition conditions when combined to other landslide-predisposing factors, such as geo-environmental factors in this case. The soil of an area classified with low landslide predisposition regarding geo-environmental conditions can hide some characteristics that situate it at the stability limit or out. These characteristics can include tension cracking, spring lines mostly due to high swelling, and shrinkage capacities, as previously stated by [
3,
4,
5,
6,
9].
One limitation of the soil geomechanical properties model is clearly the lack of spatial variability of the three different soil units that follow the spatial distribution of the geological units. This explains also the large size of the high susceptibility class that is also affecting the result of the combined models. The combined model, however, showed a much better spatial distinction of areas of different landslide susceptibility. Moreover, also the geo-environmental factors model showed a very large medium- and high-susceptibility class. Such large areas of high landslide susceptibility may lead to high prediction rates, but they are not very useful in terms of efficient targeting of measures to reduce landslide damage. A more adequate method for the classification of the landslide susceptibility could help to improve the effectiveness of the resulting landslide-susceptibility maps.
The prediction accuracy of the landslide-susceptibility models was evaluated using the receiver operator characteristic (ROC) curve and the area under the curve (AUC). Comparison of ROC curves and corresponding AUC values are presented in
Figure 12. The ROC curve of the model combining only soil properties was not computed because of its low spatial variability (this map could not be classified into more than three pixel categories). The ROC curves of the geo-environmental and combined soil-geo-environmental models display for each, an experimental and a model-fitting curves. The experimental curves of these models describe an exponential profile according to equation (Equation (6)), generated in the Origin 6.1 software (OriginLab Corporation, Northampton, Massachusetts, USA) using the “Nonlinear Curve Fit” option and presented below:
The characteristics of Equation (6): true positive rate (
TPR0), coefficient (
A1,) false positive rate (FPR) and the threshold
t1 are given in
Table 8.
The curves of the models with only geo-environmental factors and the combined model present AUC values of 0.80 and 0.93, respectively (
Figure 12). These AUC values traduce the high efficiency of these models in landslide-susceptibility prediction and show that the combined model with soil properties and geo-environmental factors is the most efficient in the identification of future landslide events of the western flank of Mount Oku. Thus, it can be noted that soil properties increase the predictive power of the model with only geo-environmental factors. Therefore, to take efficient measures in order to reduce landslide damages, geo-environmental factors and geomechanical properties should be combined in landslide-susceptibility assessment.