Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South Western Uganda

*Denis Nseka, Vincent Kakembio, Frank Mugagga, Henry Semakula, Hosea Opedes, Hannington Wasswa and Patience Ayesiga*

#### **Abstract**

Generally, soil characteristics have a significant influence on landslide occurrence. This issue has, however, not yet been adequately analysed in Kigezi highlands of South Western Uganda. In this study, soil properties such as dispersion, grain size distribution, Atterberg limits, shear strength and clay mineralogy were analysed to establish their contribution to the spatial distribution of landslides in Kigezi highlands. The results demonstrate that deep soil profiles ranging between 2.5 and 7 meters were dominated by clay-pans at a depth between 0.75 and 3 meters. Although the uppermost surface horizons of the soil profile are loamy sand, the clay content is more than 35% especially in the sub soil. This suggests that the soil materials are Vertic in nature. In addition, the upper soil layers predominantly contain quartz, while subsurface horizons have considerable amounts of illite as the dominant clay minerals, ranging from 43–47%. The average liquid limit and plasticity index was 58.43% and 33.3% respectively. Besides, high average computed weighted plasticity index (28.4%) and expansiveness (38.6%) were obtained. These soil characteristics have great implication on the timing and nature of landslide processes in the study area. A change in soil material due to varying moisture content is thought to be a major trigger of landslides in Kigezi highlands of South Western Uganda. This understanding of soil characteristics is a key step in mitigating landslide hazards in the area.

**Keywords:** soil properties, landslide occurrence, Kigezi highlands

#### **1. Introduction**

Landslides are among the major life threatening global natural disasters that cause great environmental and developmental challenges [1]. Landslides normally occur on terrains with steep slopes [1, 2]. Slope failures are attributed to factors like topography, geotechnical properties of the material as well as the existence of discontinuities [3, 4]. Most of these factors are specific to particular areas and thus, site-specific studies are important [2]. In Ref. [5] noted that soil types and soil

texture are primary-level factors, while elevation, land cover types, and drainage density are next in terms of landslide inducement. In [6] the low permeability of fine textured clayey soils exacerbates the vulnerability to landslides. This is due to increased saturation and pore water pressure which reduces soil shear strength [1]. In [7] it is also noted that landslide occurrence is common in regions where sandy clay loams are underlain by sandy clay soils. The stability of slope materials is also influenced by the presence of swelling clay minerals [8]. Soils with high clay content are considered to be the leading cause of landslides in most East African highlands [7]. Soil materials with clay content exceeding 30% and classified as vertic soils, have swell-shrink characteristics [9] which make them highly vulnerable to landslides [10].

In Ref. [11] it is indicated that clay mineralogy affects the shear and frictional resistance of the soils. Among the clay minerals, smectite, particularly montmorillonite, and illite decrease the soil residual strength, owing to their peculiar colloidchemical characteristic [12] and contribute to landslide occurrence [7]. Montmorillonite a swelling clay mineral, has a negative behaviour, very strong attraction for water and can induce soil collapse due to susceptibility to volume change [13]. The dispersive, collapsible and expansive nature of clay soils constitutes what is referred as problem soils [14]. In [15] it is noted that problem soils are widespread around the world, notwithstanding the little attention they receive in many landslide studies. At various moisture contents, landslides may be induced by problem soils due to their distinct shrink-swell properties [14].

The stability of any slope is affected by specific soil parameters including bulk density, shear strength, clay mineralogy and particle size distribution [1, 7, 16]. Such soil properties vary significantly in space and require site-specific investigations to understand their contribution to landslide occurrence [1, 17]. Whereas previous landslide studies have focused on volcanic soils of Mount Elgon region in Eastern Uganda [1, 7, 18, 19], the present study examines the influence of selected morphological properties of non-volcanic soils on landslide occurrence in Kigezi highlands. The study results provide a comparative analysis of the correlation between soil properties and landslide occurrence in both mountainous as well as highland areas of Uganda.

#### **2. Materials and methods**

#### **2.1 Study area**

The study was carried out in Rukiga uplands located within Kigezi highlands of South Western Uganda. The study area is situated between 01°21<sup>0</sup> <sup>25</sup>″ and 0°58<sup>0</sup> 08″ South, and 29°43<sup>0</sup> <sup>30</sup>″ and 30°05<sup>0</sup> <sup>51</sup>″ East (**Figure 1**). Rukiga uplands an area covering 427 km<sup>2</sup> [20], span the attitudinally heterogeneous landscapes of Rwamucucu, Maziba, Kashambya, Bubaare, Bukinda, Kamwezi and Kaharo [21] was selected for this study on the basis of its unique topography which is synonymous with reported high incidences of landslides, with visible scars unlike other parts of the highlands where landslide scars have been concealed owing to rapid regeneration [22]. The topography is very rugged, with narrow steep convex slopes and high valleys between hills. Most of these valleys have drainage lines connecting to the main valley [23]. The topography has substantial flat-topped ridges and hills, with short, steep-sided deep valleys fragmented by fluted spurs [24]. Landslide scars were, therefore, identified following an inventory with the help of local communities (**Figure 2**). The landslide scar zones provided soil sampling sites for analysis of soil characteristics and their influence on landslides (**Figure 3**).

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

**Figure 1.** *Location of the study area and soil sampling sites.*

**Figure 2.**

*Field investigations to identify landslide scars and soil sampling. Photo credit, Nseka. November, 2019.*

The highland's geology is sedimentary in nature of the Precambrian rock system [25], categorised by [26], as phyllite, shale, sandstone, quartzite, granite and gneisses of granitic composition. Other rock types of the study area include grades of schists like quartz-schists and fine textured mica-schists which belong to both the Ankole-Karagwe rock system as well as the Achaean basement complex [25].

Since rock types have influence on slope factors like slope angle and stability, slopes with phyllites and shales underneath are more affected by instability processes compared to those covered by quartzite and micaceous sandstones [25]. Additionally, slope sections with relatively weaker rocks beneath like shales have deep soil profiles attributed to high weathering rates [23]. Besides, slopes underlain

#### **Figure 3.**

*In-situ soil property analysis within the landslide scars. Photo credit, Nseka. November, 2019.*

by quartzite and granitic intrusions are less prone to landslides because they are covered by shallow soils and in some cases by bare rocks [25]. Most of the water courses within the highlands contain varied patches of alluvia sand and clay [26]. For example, coarse sands have been oberved in the warped valleys with gneiss and granite intrusions beneath, while shale and phyllite are associated with clay deposits [25]. By implication, the local geology has an important influence on landslide occurrence in the highlands. The soil types in the study area include Luvisols, Histosols, Acric Ferralsols and Dystric Regosols [23].

The climate of Kigezi highlands is warm to cool humid, characterised by a bimodal rainfall pattern with annual rainfall of 1092 mm [20] classified as moderate. Rainfall, however, increases from 1250 to 1540 mm or more in high altitude areas of greater than 2000 m above sea level [27, 28]. The main rainfall seasons span from mid-February to May with a peak in March/April, and September to December with a peak in October/November [29]. Since the highland receives significant amounts of rainfall, the study area has various highlands streams which drain valleys incised within the ridges and hills [30]. Until about a century ago, the highlands' vegetation cover was characterised by montane forests [23]. Depletion of vegetation cover has occurred in the highlands due to increased human interference [31]. The highlands of the study area are currently characterised by *Eucalyptus globulus*, *Pinus leiophylla* and farmlands [28].

#### **2.2 Soil morphological analysis**

To evaluate the soil morphology-landslide relationships, field investigations were conducted and soil samples taken at different depths and points along the slope profile and positions. In this study, 120 soil samples were collected and used in the analysis of soil morphological characteristics. Soil description was done according to FAO guidelines [32]. Visible landslide scars were identified and categorised into 10 groups based on their morphological characteristics.

Profiles in the upper slope sections were dug to a depth of 1 to 1.5 m (classified as shallow). Profiles in the middle slopes ranged between 2 to 4 m (medium) while those in the lower sections were greater than 5 m (deep). Soil samples were obtained from within and outside the landside scars. Onsite analysis of physical soil properties was conducted within landslide scars, min pits, auger holes, and full profile representative sites (**Figure 3**). Site analysis sought to characterise among

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

other things, soil matrices and bedrock within each landslide. Soil horizons were analysed in detail with a specific focus on depth, colour, texture, presence and location of claypans, and structure.

#### *2.2.1 Soil-water infiltration*

Characterisation of soil infiltration levels is an importance aspect when analysing any landslide occurrence in a given region [26]. Quantification of infiltration rates is important to understand the mechanism of slope failure [25]. In situ infiltration tests were performed purposely to measure soil-water infiltration rates. Infiltration tests were performed at various slope positions including lower-bottom, lower-middle, upper-middle and uppermost. This was done owing to variations between topographic configurations and soil characteristics. A total of 32 experiments were performed on 4 landslide sites with different land use and cover categories. Field measurements of soil-water infiltration were done using the double ring infiltrometer method [33] consisting of two concentric metal rings [34, 35].

#### *2.2.2 Rainfall data*

The frequent heavy rainfall received within Kigezi highlands have been observed to be a main trigger of landslides because, the rains increase pore water pressure in voids [28]. Given the importance of pore water pressure and antecedent moisture in triggering landslides in the study area [23, 29], rainfall data are important in their analysis [24]. The rainfall data used in this study were obtained from Kabale Meteorology Station, weather data 2015: WMO No. 63726, National No. 91290000, station name Kabale, at Elevation 1867 m, Latitude 01°15<sup>0</sup> , Longitude 29° <sup>59</sup>″. Rainfall data were compared with landslide occurrence periods to ascertain their relationships.

#### **2.3 Laboratory analyses**

The scope of laboratory testing comprised shear box, Atterberg limits, sieve and hydrometer analysis, specific gravity using test standards and XRD analysis. For purposes of defining the dominant soil mechanical and physical properties, several laboratory tests were undertaken. The grain size distribution, unit weight, natural water content, degree of saturation and Atterberg limits were among the major physical and mechanical soil properties examined. Soil property tests were undertaken in accordance with the British Standard BS 1377 procedures [36]. These tests were used in the estimation of porosity, dispersion, saturated specific weight, and saturation humidity. These soil properties were considered relevant in understanding and characterising site-specific landslide controls ("For example, see [1, 7, 17]").

#### *2.3.1 Determination of particle size*

Particle size was determined using mechanical analysis. In this method, the size range of particles within a soil sample was determined and the results were expressed as a percentage of the total dry weight. To determine the soil particle-size distribution, *Sieve analysis* and *Hydrometer analysis* methods were used on particle sizes greater and less than 0.075 mm in diameter respectively. Sieve analysis involved shaking of the soil sample through a set of sieves (ranging from 75.00 mm to 0.075 mm). The percentage of the finest size and that of the weights retained on each of a series of standard sieves of decreasing size were used to infer particle size distribution.

#### *2.3.2 Atterberg limits*

A number of soil samples obtained from the landslide sites were used for Atterberg limit tests undertaken in the geotechnical laboratory. The tests were conducted to determine among others the clay plasticity and measure the threshold water contents of a fine soil. Atterberg limits were used to determine soil plasticity which provides a clue to the type of mass movement that would characterise a given area ("For example, see [37]"). The Atterberg limits were important in determining among others, soil strength, type, stability, behaviour, and state of consolidation. Using the procedures in ASTM Standard D 4318 and D431/84 and also CNR UNI 100141 ("For example, see [38]") the limits were determined as Liquid Limit (LL) and Plastic Limit (PL). The tests gave an indication of the levels of saturation and response of a soil material to landslides. The results of grain-size distribution and Atterberg limits tests were used to classify the colluvium soils according to the Unified Soil Classification System which enabled further classification of fine materials.

#### *2.3.3 Soil dispersion tests*

A double hydrometer test was used based on Stoke's law of settling velocity as an indicative laboratory test for identification of dispersive soils ("For example, see [11, 14]").

#### *2.3.4 Soil expansiveness*

To determine whether the soils in the study area are susceptible to slope failure, investigations were done to determine the expansiveness of the soil. Three procedures were done to determine the expansiveness of the soil in the study area. During such field surveys, a description of the crack patterns on the soil surface was done at different sites to determine the swelling-shrinking characteristics of particularly clay soils. The second procedure entailed calculating the weighted Plasticity Index (PIw) on the fraction <425 μm and weighted for the sample's actual content of particles <425 μm using the formula:

$$\text{PI}\_{\text{w}} = \text{PI} \ast (\text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree}) / \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\textdegree} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\textdegree}} \text{\text{\text$$

Where PI is the Plasticity index.

The third procedure involved calculating the expansiveness (ɛex) using the formula:

$$
\varepsilon\_{\rm ex} = 2.4\mathbf{w}\_{\rm p} - 3.9\mathbf{w}\_{\rm s} + 32.5,\tag{2}
$$

Where; wp = (Plastic Limit) \* (%passing 425 mm)/100 and ws = (Shrinkage Limit) \* (%passing 425 mm)/100.

#### *2.3.5 Clay mineralogy analyses*

Clay mineralogical composition was analysed using X-ray diffraction. Soil samples were prepared for XRD analysis using the back loading preparation method ("For example, see [39]"). A PANalytical Aeris diffractometer with PIXcel detector in combination with fixed relative phase amounts (weight %) were estimated by X'Pert Highscore plus software and Rietveld method respectively.

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

#### *2.3.6 Shear strength tests*

To determine material strength, shear strength tests were executed on samples obtained from the landslide sites. These tests included Shear Box and Unconsolidated – Undrained (UU) tests. The tests were conducted to determine the shear strength parameters of soil cohesion (c) and the angle of internal friction (ϕ) in accordance to British Standard (BS) 1377: Part 7: 1990 [36].

#### **3. Results**

#### **3.1 Soil morphological properties**

#### *3.1.1 Soil profile characteristics*

At the depth less than 0.85 m, shallow soil groups were observed and these occured on very steep (i.e. between 35° and 45°), and precipitous (i.e. >45°) slopes. Medium and deep soil groups in the range of 1.5–4 m and greater than 6 m, were found occupying midslopes along topographic hollows and lower slopes respectively. Surface layers were covered by deposited black soils with an average depth of 0.5 m to 1 m. The B and C horizons had depths ranging from 0.65 m to 3.44 m and 0.88 m and 5.7 m, with a reddish-brown colour and coarse-textured materials concentrated in the lowest layers, respectively. At a 1.3 m depth along the spur slopes, unweathered materials were observed. Other unweathered materials were observed at a depth of 4.5 m along topographic hollows and 7 m deep in the valley bottoms. Along the soil profiles, visible colours and textural gradation within 3 distinct horizons were seen. An abrupt change in colour from yellow (5Y8/2) to brown (7.5YR5/2) clay with Gs 2.9 g/ cm3 and γ 2.06 g/cm3 was exhibited in the upper slope elements for most of the analysed profiles. Dark reddish brown (5YR 3/6) and orange (7YR 5/6) colour in moist and dry conditions respectively characterised most of the surface horizons. Surface horizons with brown (7.5YR 4/4) colour in both moist and dry conditions and a dark to dull reddish-brown colour (7.5YR 4/6) in underlying horizons dominated most of lower slope element profiles. Soil colour in combination with other physical properties, were used to differentiate horizon types of the same and different soil profiles.

A relationship was noticed between the regions topographic characteristics and the depth of soil. Basing on the soil profile description, it was noted that the depth of soil reduced with an increase in slope angle. Whereas very deep soils (>5 m) were dominant on slope sections with slope angles less than 10°, moderately deep soils (2 to 5 m) were noticed on slope angles between 15° and 35°. Slope elements with slope angles greater than 35°, were characterised with shallow thin soils of less than 1 m (Appendix 1). A relationship was also established between slope curvature and the depth of the soil. Very deep, deep and moderately deep soil profiles (>2.5 m) were dominant along concave profile elements common in topographic hollows. Convex profile elements along spur slopes and hilltops were dominated by shallow thin soils of less than 1 m (Appendix 1).

#### *3.1.2 Location of clay pans*

Field observations within dug profiles indicated the presence of claypans at a depth of 0.75 m to 3 m (**Figure 4**). Variations in the location of claypans is noticeable along the soil profile. In some profiles, the claypans were found close to the surface while in some horizons they were identified at greater depths. Whereas the upper slope soil profiles had claypans close to the surface at less than 0.85 m depth,

**Figure 4.** *Soil profiles dominated by claypans in the sub-surface soil horizons. Photo credit, Nseka, may 2019.*

lower slope soil profiles had pans at greater than 2 m depth. Claypans in soil profiles along the middle slope sections existed at depth ranging between 1 and 2 m. Characterising variations in depth of claypans horizon is an important step in explaining other soil properties, as expounded in the discussion section.

#### **3.2 Effect of particle size on landslide occurrence**

The study area is dominated by clay soils in the subsurface horizons. The particle size determination (**Figure 5**) shows that soils are dominated by clay presence, except for the uppermost surface horizons.

From the materials analysed, fine grained materials of either clay or silt were predominant (**Figure 5**). The texture of the soil varied with profile depth as

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

revealed by sieve analysis. In the subsurface horizons, finer clay and silt materials dominates, while sand particles are the most dominant in surface soils. It was observed that all surface horizons had sand (33–55%), silt (22–40% and clay (10–30%). In the deeper horizons, sand was observed to reduce drastically to less than 23%, while clay increased to greater than 50% (Appendix 1). The clay content for all the samples in the sub soil and deep soil horizons is well above the 32% threshold for vertic soils. The swelling and shrinkage characteristics of vertic soils are very important in the localisation of landslides, as explored in the subsequent section.

#### **3.3 Clay mineralogy and landslide occurrence**

The results of the X-ray diffraction tests are provided in **Figure 6** and **Table 1**. Following the X-ray diffraction patterns, it was established that quartz, illite/ muscovite and kaolinite are the dominant minerals in the soil materials. From the soil mineralogical analysis conducted, it was established that quartz was the dominant soil constituent within the top layers. By implication, these high amounts of quartz within surface soil horizons affect the behaviour of the soil to incoming rainfall as will be unravelled in the subsequent section. Considerable amounts of illite/muscovite as the dominant clay minerals, ranging between 43% and 47% (**Table 1** and **Figure 6**) were present in the subsurface soil horizons. Notwithstanding the absence of smectite clays, the soils contain large amounts of moderately expansive clays, particularly illite/muscovite.

#### **3.4 Soil dispersion and landslide occurrence**

The LL and Plasticity Index for all the tested samples are greater than 50% and 30% respectively (**Figure 7** and **Table 2**).

The fine-grained soil materials tested had LL with an average value of 58.43%, ranging between 50.43% and 66.43%, which is considered to be of high plasticity. Whereas the plasticity index ranged between 22.4% and 44.2%, the plastic value ranged between 21.3% and 28.9% (**Figure 7**). High plasticity index of more than 30

**Figure 6.** *XRD patterns of clay minerals.*


#### **Table 1.**

*XRD mineral distribution in percentages.*

was detected in the clay materials, thus rendering them as highly plastic and expansive. The computed Linear Shrinkage (LS) ranged between 10.53 and 20.76. Whereas soil expansiveness (ɛex) ranged between 10.7 and 54.8 averaging 38.6, the computed weighted Plasticity Index (PIw) ranged between 17.92% and 34.92% averaging 28.4% from all the analysed soil samples (**Figure 7**). As regards to soil


#### *Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

**Table 2.** *Plasticity index and shear strength parameters.* expansiveness (ɛex), all the tested samples were in the range between 20 and 50, and this indicated medium expansive soils which were susceptible to landslides.

From the plasticity chart (**Figure 8**), all the soil samples tested were inorganic clay of high plasticity belonging to the CH group. Such soils can easily move when saturated, leading to a high incidence of landslides.

From **Figure 8**, it is evident that soils in the study area have high plasticity as already revealed by XRD results. Such soils can easily move when saturated, leading to high incidence of landslide occurrence. The presence of many cracks, observed on the soil surface during field investigations confirms the high plasticity nature of the soil materials (**Figure 9**). The presence of such cracks on the surface is also a characteristic of Vertisols with high expansive potential.

From the double hydrometer test, it can be vividly observed that most of the samples have dispersion values greater than 30%. By implication, such soil materials are susceptible to landslide occurrence. From the soil samples examined, critical dispersion values greater than 50% were established from more than 90% of the samples (**Figure 10**). Such high dispersion values imply greater susceptibility to landslides in the region. This study established that, highly dispersive soils are particularly dominant in the surface soil layers associated with greater percentages of illite/muscovite clay minerals.

**Figure 8.** *The distribution of samples on the plasticity chart for the USCS.*

**Figure 9.** *Highly expansive materials with numerous surface cracks. Photo credit, Nseka, march 2019.*

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

**Figure 10.** *Double hydrometer test results.*

The linear shrinkage computed from the soil samples ranged between 10.53 and 20.76 (**Table 2**). The computed weighted plasticity index (PIw) from the analysed soil materials ranged between 17.92% and 34.92%, averaging 28.4%.

The computed weighted plasticity index for most of the tested soils is above the 20% threshold for expansive soils, signifying highly unstable soils. All the tested samples had expansiveness (ɛex) between 20 and 50 indicating medium expansive soils which are highly vulnerable to slope failure. More than 90% of the analysed samples have expansiveness above the 20% threshold. More than 80% of the samples have values between 20 and 50% showing medium expansive soils (**Table 2**).

#### **3.5 Shear strength parameters**

Shear strength parameter test results show that the soils have low cohesion (C). From the samples analysed, soil cohesion ranged between 5.2 kPa and 11.1 kPa averaging 8.2 kPa. The angle of internal friction obtained from the analysed soil samples ranged between 2.6° and 8.1° averaging 5.4° (**Table 2**). From the soils completely saturated with water, a lower cohesion of 5.2 kPa was detected. A minimal cohesion value of 8.2 kPa was considered as the critical state equilibrium all over the area. This is due to the fact that slope failure is expected to occur when soil materials are saturated. The present study established that most of tested soil samples had a very low internal friction angle (<8.5°). Such soil materials with very low internal friction angle were considered weak with high vulnerable to landslides.

#### **3.6 Soil water infiltration**

High infiltration rates were observed in the top soils with depth ranging from 0.3 to 0.8 m, but drastically reduced in the sub soils. Steady infiltration rates greatly varied with topographic characteristics and land use/cover types. Higher infiltration values were noticed in the lower slope elements than the upper sections (**Figure 11**). From the infiltration experiments conducted, it was revealed that the average soilwater infiltration was 24 cm/h<sup>1</sup> in uppermost slope sections. On the upper-middle, lower-middle slope sections and bottom valleys, soil-water infiltration rates of 30 cm/h<sup>1</sup> , 70 cm/h<sup>1</sup> and greater than 80 cm/h<sup>1</sup> respectively were achieved from the experiments conducted. Along the hollows within the upper slope sections, the

**Figure 11.** *Infiltration rates along slope positions, topographic configurations and land uses/covers.*

observed infiltration rates were greater than 30 cm/h<sup>1</sup> and less than 12 cm/h<sup>1</sup> along the spur slope sections. Whereas the infiltration rates within topographic hollows along the middle slope sections was greater than 70 cm/h<sup>1</sup> , it was less than 45 cm/h<sup>1</sup> on the spur slope elements (**Figure 11**). Within the top soil layers, high infiltration rates are explained by the dominancy of loamy sandy soils. The predominance of clay materials in the subsoil, with claypans distinctly underlying the top soils, limits water infiltration in the study area. Such soil characteristics affect their response to incoming rainfall and consequently the timing of landslides, as explored in the discussion section.

The land use and cover characteristics influence soil water infiltration in the study area. Within the agricultural land uses particularly along cultivated zones, infiltration rates were noticed to be higher than those observed on natural land cover types. The infiltration rates along the agricultural land uses were generally greater than 65 cm/h<sup>1</sup> with the exception of beans covered areas where infiltration rates were noticed to be lower than 42 cm/h<sup>1</sup> . On the natural land cover types including grasslands, thickets and shrubs, infiltration values lower than 30 cm/h<sup>1</sup> were generally observed (**Figure 11**). Areas covered by forests, however, tended to have infiltration values higher than 45 cm/h<sup>1</sup> (**Figure 11**). Basing on these results, it was deduced that infiltration values in the study area vary greatly between rapid and very rapid. They also vary with slope characteristics including gradient and curvature as well as land use and cover properties. Following the experiments conducted, it was concluded that the steady state water infiltration values in the study area range from 12.2 cm h<sup>1</sup> to 88.5 cm h<sup>1</sup> .

#### **3.7 Rainfall distribution and soil behaviour**

The behaviour of soil materials and its susceptibility to landslide occurrence greatly depends on rainfall amounts and distribution in the region. During the

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

#### **Figure 12.**

*Long term average monthly rainfall distribution for 1980 to 2020.*

analysis of landslide occurrence, considering the influence of soil antecedent moisture and its implication on soil pore water pressure is very important. Monthly rainfall distribution shows that March, April, October and November are the wettest months in the study area (**Figure 12**). Rainfall amounts and distribution have great implications for soil behaviour and hence landslide occurrence in the study area. Seasonal rainfall distribution shows more rains are received during the MAM (352.5 mm) and SON (327.8 mm) seasons, while DJF (238.4 mm) and JJA (104.72 mm) seasons receive less rainfall (Appendix 2). It is noteworthy, however, that landslide occurrence in the study area is not linked to individual rainfall events, but correspond with seasonal rainfall distribution.

Following an interaction with the local communities as well as local government reports, it was established that most of the landslides are experienced during months of May and November. These are, however, not the wettest months in the region (Appendix 2). In 2010, for example, more landslides occurred during the month of May which, however, received lower rainfall amounts (97.7 mm) compared to the preceding months of March (149 mm) and April (133 mm). During 2013, landslides were similarly experienced during the month of November with lower rainfall amounts (122. 2 mm) compared to the preceding months of September (134.1 mm) and October (154 mm) (Appendices 2 and 3). It is therefore, noteworthy that landslides in the study area do not necessarily occur in the wettest months of the year. The implications of this phenomenon is unravelled in the subsequent section.

#### **4. Discussion**

#### **4.1 Soil profile characteristics**

Deep soil profiles ranging between 2.5 and 7 m are a major characteristic of the study area. Deeper soil profiles are more pronounced along topographic hollows and valley bottoms. Soil depth forms one of the conditions for assessing the stability of the soil materials and landslide susceptibility of the landscape [40]. Soil depth and its moisture content determine how water can be stored in the soil before saturation is reached [41]. Although most of the ridges in the study are characterised by deep to very deep soils profiles (greater than 6 m on most slopes), majority of the landslide features are shallow, concentrating within 1 to 3 m of the profile. Shallow landslide scars are not uncommon in slope sections covered by deep and very deep soil materials. As has been illustrated by [7], most sections covered by deep soils on the slopes of Mt. Elgon in Eastern Uganda experience deep seated landslides. Shallow landslides on deep soil covered slopes in the Kigezi highlands was, therefore, considered as an anomaly. Soil profile analysis was undertaken to establish the cause of this irregularity. The profile characterisation indicated the existence of 0.9 to 3 m dense claypans within the profiles. Following the infiltration experiments conducted, it was established that claypans decreased the rate of water infiltration within the soil materials. The restraint of perpendicular water flow within the soil materials by claypans has also been confirmed by [42]. The accumulating water leads to saturation of claypans sandwiched between more stable materials. Saturated claypans can act as a sliding surface for the overlying materials, consequently inducing landslides. Other studies elsewhere also confirm that variations of claypan profile properties over the landscape greatly influence soil-water holding capacity [35, 43]. It can, therefore, be inferred that the occurrence and characteristics of landslides in Kigezi highlands is highly influenced by the presence and position of claypan horizons within the soil profile.

#### **4.2 Particle size distribution**

Fine-grained silt and clay soils predominate the study area. From the mechanical analysis conducted to determine soil particle sizes, it was established that clay is the dominant material, with greater than 40%. The distribution of sand and silt in the soil materials was less than 35% and 25% respectively. The soil materials in Kigezi highlands can be classified as Vertisols due to the high clay content of more than than 35% on average and plasticity index (PI) greater than 33%. Such soils are known for inducing landslides [1]. Vertic soils characteristically expand when wet and shrink in dry conditions due to their high clay content [9]. The presence of large amounts of clay in soils of the study area is a major factor in landslide occurrence, since it affects the stability of the soils when wet. Likewise [7], also observed that the susceptibility to landslide occurrence on the slopes of Mount Elgon is due to abundance of fine-grained materials in the subsurface. The study results are, therefore, consistent with several studies elsewhere which have demonstrated the influence of high clay content on landslide occurrence [1, 18, 44].

#### **4.3 Clay mineralogy**

XRD Clay mineral analyses indicated the presence of moderately expansive clays, particularly illite/muscovite. Previous studies indicate that the presence of illite clays can lead to landslide occurrence due to their swelling potential and low shear strength [10, 45]. In the same vein [7], also confirmed that the occurrence of landslides on Mount Elgon slopes in Eastern Uganda is associated with the existence of greater amounts of kaolinite and illite clay minerals in the soil materials. The existence of significant amounts of illite/muscovite clay minerals in the study area further confirmed that Vertisols with high shrink-swell properties can result into landslides. The behaviour of the soil materials greatly determine the timing of landslide occurrence during rainfall seasons. There is fast flow of water through the surface soil materials with greater amounts of quartz into the deep soil profiles with clay abundance at the start of the rainfall season. Greater amounts of

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

illite/muscovite clay minerals in the subsurface absorb the arriving water resulting into moisture build-up. The amassing water in the soil materials leads to soil behavioural change which can swell and loose cohesion [46]. This phenomenon, therefore, explains why landslides in the study area are not experienced at the beginning of the rainfall season or immediately after extreme rainfall events, as is the case with Mt. Elgon region in Eastern Uganda.

#### **4.4 Soil dispersion**

The study area is characterised by soil materials associated with high plasticity and are, therefore, inorganic in nature (CH), signifying weak soils of high saturation (**Figure 7**). Results from soil dispersion tests revealed high plasticity index of greater than 30% signifying Vertic soils. During continuous rainfall events, such soils with high plasticity can easily slide. The LLs for all samples analysed was above 50%, signifying high plasticity. Following soil expansiveness analysis, it was established that the average weighted plasticity index (PIw) and expansiveness were 26.4% and 32.8% respectively. By implication, such soils are highly dispersive and excessively susceptible to landslide occurrence [10]. Landslide occurrence in the study area is, therefore, associated with expansive soils which shrink and swell leading to loss of soil strength. The expansive potential of the soils influenced by high clay content and type especially illite/muscovite is one of the major factors promoting landslides in Kigezi highlands.

#### **4.5 Soil water infiltration**

Soil water infiltration was noted to vary across slope position and topographic configuration. Whereas upper slope sections and spur slopes with shallow soils experience low infiltration rates, lower slope sections and topographic hollows with deeper soils are associated with high infiltration rates. This variation in infiltration rates also signifies differences in soil saturation levels ("For example, see [47]"). The lower slope sections are associated with greater saturation rates which result into saturation overland flow processes. Along the topographic slope elements, the saturation overland flow incrementally moves upslope from the slope base [48]. Materials along topographic hollows remain saturated most times due to greater clay content dominated by illite/muscovite minerals. This phenomenon leads to decreased soil material strength along the topographic hollows, inducing landslide occurrence. Equally noticed in this study was a relationship between topographic characteristics and soil water infiltration. Greater soil-water infiltration values are experienced along topographic hollows and lower slope elements than the spurs and upper sections. Such variations in infiltration values in relation to topographic characteristics are associated with differences in soil depth along the slope configurations. This study, therefore, approves that using the soil-water infiltration experiments along different slope sections and gradients it is possible to predict landslide occurrence.

#### **4.6 Change in soil behaviour with rainfall distribution**

As opposed to the situation on the slopes of Mt. Elgon in Eastern Uganda, landslide occurrence in Kigezi highlands is not related to extreme rainfall events. On the slopes of Mt. Elgon in Eastern Uganda, majority of the landslides are commonly witnessed during or immediately after extreme rainfall events [1, 7, 18]. Local communities and regional government reports confirmed that landslide occurrence in the study area is usually not experienced during or immediately after peak rainfall seasons. Paradoxically, they occur during the less wet months of the rainfall season. Most landslides in the study area are experienced during the months of May and November, despite the preceding months of April and October receiving more rainfall amounts. This phenomenon can be explained by the unique infiltration dynamics through quartz dominated top soil layers and saturation of the claypans dominated by illite/muscovite clays in the lower soil horizons. This leads to antecedent moisture building up in the sub soil materials as more rainfall is received hence landslide occurrence in the region. Antecedent soil moisture condition prior to a rainfall event has been confirmed as the most significant factor in landslide occurrence [49].

#### **5. Conclusion**

Deep soil profiles are a major characteristic of the study area. Notwithstanding the deep soil profiles on most slope elements in the study area, majority of the landslide scars are shallow, occurring within less than half of the profile due to presence of claypans. The claypans act as a slipping zone for the overlying soil materials. The study area is dominated by fine silt and clay soil materials. In association with greater amounts of clay percentage of more than 35% on average and PI greater than 33%, the soil materials in the study area are classified as Vertisols, which are synonymous with landslide occurrence. The predominance of reasonable amounts of expansive clays, mainly illite in the study area influences the stability and vulnerability of slope materials to landslide occurrence. In Kigezi highlands, landslides are not normally experienced during or immediately after extreme rainfall events but occur later in the rainfall season due to initial infiltration through quartz dominated upper soil layers, before illite/muscovite clays in the lower soil horizons get saturated. This behavioural change in the soil material due to moisture content is, therefore, the major trigger of landslides in Kigezi highlands. An understanding of these soil characteristics is an important step in landslide hazard mitigation in Kigezi highlands.

#### **Acknowledgements**

The authors gratefully acknowledge the research grant from Makerere University – Swedish International Development Cooperation Agency (SIDA) Phase IV (2015/2020 Agreement) – Building Resilient Ecosystems and livelihoods to Climate Change and Disaster Risk (BREAD) project 331 research component, which funded travel and fieldwork for this study.

#### **Conflict of interest**

The authors have no conflict of interest.

**Appendices**



#### *Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*


#### *Landslides*

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*


#### **Appendix 2. Temporal distribution of rainfall in Kabale highlands**

*Source: Kabale meteorology station, weather data: WMO No. 63726, National No. 91290000, station name KABALE, Elevation 1867 m, Latitude 01°15*<sup>0</sup> *, Longitude 29°59*<sup>0</sup> *.*


#### **Appendix 3. Relationship between seasonal rainfall distribution and landslide occurrence**

## **Author details**

Denis Nseka1 \*, Vincent Kakembio<sup>2</sup> , Frank Mugagga<sup>1</sup> , Henry Semakula<sup>1</sup> , Hosea Opedes<sup>1</sup> , Hannington Wasswa<sup>1</sup> and Patience Ayesiga<sup>3</sup>

1 Department of Geography, Geo-Informatics and Climatic Sciences, School of Forestry, Environmental and Geographical Sciences, Makerere University, Kampala, Uganda

2 Department of Geosciences, School of Environmental Studies, Nelson Mandela University, Port Elizabeth, South Africa

3 Department of Geography, Faculty of Education, Bishop Stuart University, Mbarara, Uganda

\*Address all correspondence to: denisnseka1@gmail.com

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

#### **References**

[1] Mugagga, F. V., Kakembo, V., Buyinza, M., 2011. A characterisation of the physical properties of soil and the implications for landslide occurrence on the slopes of mount Elgon, eastern Uganda. Journal of the International Society for the Prevention and Mitigation of natural hazards ISSN 0921-030X. DOI:10.1007/s11069-011- 9896-3.

[2] Kirschbaum, D. B. and Zhou, Y. 2015. Spatial and temporal analysis of a global landslide catalog. Geomorphology DOI: 10.1016/Journal of Geomorphology. 2015.03.016.

[3] Bamutaze, Y. (2019). Morphometric Conditions Underpinning the Spatial and Temporal Dynamics of Landslide Hazards on the Volcanics of Mt. Elgon, Eastern Uganda. In Emmerging Voices on Natural Hazards. Elsevier Inc.

[4] López-Davalillo, B., Monod, M. I., Alvarez-Fernandez, G., Herrera Garcia, J., Darrozes, C., Gonzalez-Nicieza, Olivier, M. 2014. Morphology and causes of landslides in Portalet area (Spanish Pyrenees): Probabilistic analysis by means of numerical modelling. Engineering Failure Analysis 36.

[5] Hong, Y., Adler, R., Huffman, G., 2007. Use of satellite remote sensing data in the mapping of global landslide susceptibility. Natural Hazards 43, 245-256.

[6] Vagen, T.G., 2010. Africa soil information service: Hydrologically corrected/adjusted SRTM DEM (AfrHySRTM). Nairobi, Kenya and Palisades, NY: International Centre for Tropical Agriculture – Tropical Soil Biology and Fertility Institute (CIAT-TSBF), World Agroforestry Centre (ICRAF), Centre for International Earth Science Information Network (CIESIN), Columbia University.

[7] Kitutu, M. G., Muwanga, A., Poesen, J., Deckers, J. A. 2009. Influence of soil properties on landslide occurrence in Bududa district, eastern Uganda. African Journal of Agriculture Res 4(7), 611–620.

[8] Zinck, J.A., 2013. Geopedology: Elements of geomorphology for soil and geo-hazard studies. ITC Special Lecture Notes Series: ISBN: 978-90-6164-352-4, 13.

[9] Mukasa-Tebandeke, I. Z., Ssebuwufu, P. J. M., Nyanzi, S. A., Schumann, A., Nyakairu, G. W. A., Ntale, M., and Lugolobi, F., 2015. The elemental, mineralogical, IR, DTA and XRD analyses characterized clays and clay minerals of central and eastern Uganda. Advances in Materials Physics and Chemistry 5, 67-86.

[10] Yalcin, A. 2011. A geotechnical study on the landslides in the Trabzon Province, NE, Turkey. Applied Clay Science 42, 11–19.

[11] Fauziah, A., Yahaya, A, S. and Farooqi, M. A. 2006. Characterization and geotechnical properties of Penang residual soils with emphasis on landslides. American Journal of Environmental Sciences 2, 121-128.

[12] Bell, F.G. and Walker, D.J.H. 2000. A further examination of the nature of dispersive soils in Natal, South Africa. Quarterly Journal of Engineering Geology and Hydrogeology, 33, 197–199.

[13] Bell, F.G. 2004. Engineering Geology and Construction. CRC Press, Taylor and Francis. ISBN; 0415259398, 9780415259392

[14] Baynes, F.J. 2008. Anticipating problem soils on linear projects. Conference proceedings on problem soils in South Africa, 3-4 November 2008, 9-21.

[15] Williams, A.A.B., Pidgeon, J.T. and Day, P.W. 1985. Expansive soils. Transactions of the South African Institution of Civil Engineers, 27, 367-397.

[16] Zung, A.B., Sorensen, C.J., Winthers, E., 2009. Landslide soils and geomorphology in Bridger/Teton Forest Northwest Wyoming. Physical Geography, 30(6), 501–516.

[17] Roller, S., Wittmann, H., Kastowski, M., Hinderer, M., 2012. Erosion of the Rwenzori Mountains, east African rift, from in situ-produced cosmogenic 10Be. Journal of Geophysical Research, 117, F03003.

[18] Knapen, A., Kitutu, M.G., Poesen, J., Breugelmans, W., Deckers, J., Muwanga, A., 2006. Landslides in a densely populated county at the footsteps of mount Elgon (Uganda): Characteristics and causal factors. Geomorphology, 73, 149–165.

[19] Kitutu, M.G., Poesen, J.M., Deckers, J., 2011. Farmer's perception on landslide occurrences in Bududa District, eastern Uganda. African Journal of Agricultural Research, 6, 7-18.

[20] UBOS (Uganda Bureau of Statistics). 2017. Statistical Abstracts 2017. Ministry of Finance, Planning and Economic Development, Uganda. http://www.ubos.org

[21] Mugagga, F., N. Nakanjakko., B. Nakileza., D. Nseka 2020. Vulnerability of smallholder sorghum farmers to climate variability in a heterogeneous landscape of Kigezi highlands, South-Western Uganda. Journal of Disaster Risk Studies (Jàmbá), 12(1), a849. DOI: 10.4102/jamba.v12i1.849. https://jamba. org.za/index.php/jamba/article/view/ 849

[22] Nseka, D., Mugagga, F., Opedes, H., Ayesiga, P., Wasswa, H., Mugume, I., Nimusiima, A., & Nalwanga, F. 2021.

The damage caused by landslides in socio-economic spheres within the Kigezi highlands of South Western Uganda. Journal of environmental and socio-economic studies, DOI:10.2478/e nviron-2021-0003, 9, 1: 23-34.

[23] Bagoora, F.D.K., 1993. An assessment of some causes and effects of soil erosion hazard in Kabale Highland, South Western Uganda, and people's attitude towards conservation. In Abdellatif (ed.) resource use and conservation: Faculty of Social Sciences; Mohammed V. university, Rabat Morocco. Mountain Research and Development, Vol 8.

[24] Nseka, D., Kakembo, V., Bamutaze, Y. and Mugagga, F. 2019. Analysis of topographic parameters underpinning landslide occurrence in Kigezi highlands of South Western Uganda, Journal of Natural Hazards. Springer, DOI: 10.1007/s11069-019-03787-x

[25] Bagoora, F.D.K., 1989. A preliminary investigation into the consequences of inadequate conservation policies on steep slopes of the Rukiga highlands, South Western Uganda, In Thomas D. B, E. K Biamah, A. M Kilewe (eds.); soil conservation in Kenya, Dept. of Agri. Er, Univ. of Nairobi, Kenya.

[26] Bagoora FDK (1988) Soil erosion, mass wasting risk in the highland areas of Uganda: Mountain Research and Development, Vol 8.

[27] NEMA (National Environment Management Authority). 2018. State of Environment Report for Uganda for 2017/18. National Environment Management Authority, Kampala, Uganda. http://www.nemaug.org

[28] UBOS (Uganda Bureau of Statistics). 2016. Statistical Abstracts 2016. Ministry of Finance, Planning and Economic Development, Uganda. http://www.ubos.org

*Implications of Soil Properties on Landslide Occurrence in Kigezi Highlands of South… DOI: http://dx.doi.org/10.5772/intechopen.99865*

[29] NEMA (National Environment Management Authority). 2017. State of Environment Report for Uganda for 2016/17. National Environment Management Authority, Kampala, Uganda. http://www.nemaug.org

[30] NEMA (National Environment Management Authority). 2016. State of Environment Report for Uganda for 2015/16. National Environment Management Authority, Kampala, Uganda. http://www.nemaug.org

[31] NEMA (National Environment Management Authority). 2020. State of Environment Report for Uganda for 2019/20. National Environment Management Authority, Kampala, Uganda. http://www.nemaug.org

[32] FAO., 2006. Food and Agriculture Organization of the United Nations: Guidelines for soil description. Rome (Italy). ISBN: 92-5-105521-1.

[33] ASTM D3385., 2003. Standard test method for infiltration rate of soils in field using double-ring infiltrometer. Annual book of ASTM standards, vol. 4.02. ASTM International.

[34] Philips, C. E. and Kitch. W.A., 2011. A review of methods for characterization of site infiltration with design recommendations. Paper presented at 43rd Symposium on Engineering Geology and Geotechnical Engineering, University of Las Vegas, Las Vegas, NV, March 23-25, 2011.

[35] Ahmed, F., Gulliver, J.S. and Nieber, J.L., 2011. A new technique to measure infiltration rate for assessing infiltration of BMPs. Presented at 12th international conference on urban drainage, Porto Alegre/Brazil, 11-16 September 2011.

[36] British Standards Institution. British Standards 1377: 1995 and British Standard 1377: 1990. Methods of test for soils for civil engineering purposes, London, 1990.

[37] Moussadek, R., Laghrour, M., Mrabet, R., Van Ranst, E., Badraoui, M., Mekkaoui, M., 2017. Morocco's Vertisols characterization. Journal of Materials and Environmental Sciences ISSN; 2028–2508, Volume 8, Issue 11, pp. 3932-3942.

[38] Gao, J. and Maro, J., 2010. Topographic controls on evolution of shallow landslides in pastoral Wairarapa, New Zealand, 1979-2003. Geomorphology, 114(3), 373–381. DOI: 10.1016/j.geomorph.2009.08.002

[39] Okalebo, J.R., Gathua, K.W., Woomer, P.L., 1993. Laboratory Methods of Soil and Plant Analysis: A Working Manual, Soil Science Society of East Africa, EPZ (Kenya) Limited, Nairobi, Kenya.

[40] Liang, W. L. and Uchida, T., 2014. Effects of topography and soil depth on saturated-zone dynamics in steep hillslopes explored using the threedimensional Richards' equation. Journal of Hydrology, 510, 124–136. DOI: 10.1016/j.jhydrol.2013.12.029.

[41] Merino-Martín, L., Moreno-de las Heras, M., Espigares, T. and Nicolau, J. M., 2015. Overland flow directs soil moisture and ecosystem processes at patch scale in Mediterranean restored hillslopes. Catena, 133, 71–84. DOI: 10.1016/j.catena.2015.05.002

[42] Jiang, F.S., Huang, Y.H., Wang, M. K., Lin, J.S., Zhao, G., Ge, H.L., 2014. Effects of Rainfall Intensity and Slope Gradient on Steep Colluvial Deposit Erosion in Southeast China. Soil Science. Society of American Journal.

[43] Xiang, W., Cui, D., Liu, L., 2007. Experimental study on sliding soil of ionic soil stabilizer-reinforces. Earth Sciences 32(3), 397–402.

[44] Wati, S.E., Hastuti, T., Wijojo, S. and Pinem, F., 2010. Landslide susceptibility mapping with heuristic

#### *Landslides*

approach in mountainous area. A case study in Tawangmangu sub district, Central Java, Indonesia. International achieves of the photogrammetry, Remote Sensing and Spatial Information Science, 38, 248–253.

[45] Yalcin, A., 2007. The effects of clay on landslides: A case study. Applied Clay Science, 38, 78–85.

[46] Van Den Eeckhaut, M., Poesen, J., Hervas, J., 2013. Mass-movement causes: Overloading. In: Schroder, J.F. (Ed.). Treatise on Geomorphology, Vol 7, Mountain and Hillslope Geomorphology. Academic, San Diego, 200–206

[47] Rousseau, M., Cerdan, O., Ern, A., Le Maître, O. and Sochala, P., 2012. Study of overland flow with uncertain infiltration using stochastic tools. Advances in Water Resources, 38, 1–12. DOI:10.1016/j.advwatres.2011.12.004

[48] Reynolds, W.D., Elrick, D.E., Youngs, E.G., Amoozegar, A., Booltink, H.W.G. and Bouma, J., 2002. Saturated and field-saturated water flow parameters. In Dane J. H and G.C. Topp (Eds.) methods of soil analysis. Part 4. Physical methods. SSSA book Ser. 5. SSSA, Madison, WI, 797–878.

[49] Guan-Wei, L. and Hongey, C., 2012. The relationship of rainfall energy with landslides and sediment delivery. Engineering Geology, 125, 108-118.

#### **Chapter 4**
