British Geological Survey

Quick links

  • Climate change
  • Multiple hazards
  • Decarbonisation

Landslide case studies

BGS Research

The  landslides team  at the BGS has studied numerous landslides. This work informs our geological maps, memoirs and sheet explanations and provides data for our  National Landslide Database , which underpins much of our research. If you’ve seen a landslide or want to contact us about landslides, you can email us at  [email protected] .

View our landslide case studies

Aldbrough Coastal Landslide Observatory.

Aldbrough Coastal Landslide Observatory – landslide case study

View of the Berry Hill landslide and disused quarry slope taken from Bank End Close © BGS

Berry Hill quarry, Mansfield, Nottinghamshire

In 2019, a landslide occurred on a disused former sandstone quarry slope in Mansfield and affected over 35 properties.

Tynemouth Castle, Priory ruins, 2013. Photo taken across Short Sands beach and King Edward’s Bay looking south-south-east. Permian Upper and Lower Magnesian Limestone on Yellow Sands on Middle Coal Measures sandstone.

Tynemouth Headland, North Shields

Seaside Road in Aldbrough in2003 cut off by the cliff.

Landslides and coastal erosion at Aldbrough, East Riding of Yorkshire

Evans Cliff Landslide, a minor slip on the western flank of Stonebarrow Hill.

Stonebarrow Hills, Dorset

Derailed engine below the rail track.

Stob Coire Sgriodain, Highlands

Cracks in a building caused by the 1994 landslide.

St Dogmaels, Pembrokeshire

Terrestrial LiDAR in action at Sidestrand.

Sidestrand, Norfolk

The road has buckled up under compressional forces as the ground beneath it moved downslope.

Shortwood Farm, Pencombe, Herefordshire

Contact the landslide response team, you may also be interested in:.

landslides and rainfall

Landslides and rainfall

We combine rainfall statistics and landslide data to produce graphs demonstrating the relationship between landslides and rainfall.

Downloading meteorological data

Hazard monitoring

Hazard monitoring provides valuable information for process understanding and for the science required for developing early warning systems to protect people and infrastructure.

Was this page helpful?

How can we make this section better?*

  • How can we make this section better? *

Please select a reason*

  • Please select a reason * Select reason Didn't answer my question Too confusing Links are broken/Doesn't work

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • View all journals
  • Explore content
  • About the journal
  • Publish with us
  • Sign up for alerts
  • Open access
  • Published: 27 August 2019

Comprehensive analysis of landslide stability and related countermeasures: a case study of the Lanmuxi landslide in China

  • Zheng Han 1 , 2 ,
  • Yange Li 1 , 2 ,
  • Yangfan Ma 1 ,
  • Weidong Wang 1 , 3 &
  • Guangqi Chen 4  

Scientific Reports volume  9 , Article number:  12407 ( 2019 ) Cite this article

6622 Accesses

12 Citations

Metrics details

  • Geomorphology
  • Natural hazards

We report on a comprehensive method for analyzing landslide stability and the mitigation effect of countermeasures in this paper. The proposed method is a combination of theoretical method and numerical method. To address the uncertainties of the soil strength parameters, the rational values of these parameters are comprehensively determined by the back-analysis result of the reliability method and the result by the quantitative method, as well as the in-situ geological test. To evaluate the slope stability, the limit analysis using the 2D upper bound method and the FEM simulation using strength reduction method are performed, respectively. In order to illustrate the presented method, the so-called Lanmuxi landslide in China is selected as a case study. Results demonstrated that the stress and strain majorly concentrated at the toe and crown of the slope. According to the analysis results, countermeasures consisting of anchor lattice beams, landslide piles, and cracks filling, are suggested to reduce the failure risk of the landslide. Effect assessment based on the FEM analysis verifies the feasibility and effectiveness of the recommended countermeasures.

Similar content being viewed by others

case study on landslides

Failure probability assessment of landslides triggered by earthquakes and rainfall: a case study in Yadong County, Tibet, China

case study on landslides

The failure mechanism of the Baishi landslide in Beichuan County, Sichuan, China

case study on landslides

Assessment of landslide behaviour in colluvium deposit at Doi Chang, Thailand

Introduction.

Landslides are a common geological phenomenon in mountainous regions worldwide, posing a severe risk to local infrastructures. During the years from 2004 to 2010, 2620 fatal landslides in total were recorded, causing 32,322 fatalities 1 . To protect residents and infrastructures against landslides, a rational design of countermeasures is necessary. The commonly used countermeasures include supporting measures and drainage measures 2 . Current studies in landslide hazard mitigation mainly focus on the overall design principles, structure optimization, and the mitigation effect evaluation. However, it is widely accepted that an effective countermeasure against landslide depends on the deep understanding of landslide mechanisms and rational analysis of landslide stability, which remain a scientific challenge.

Up-to-date studies on landslides can be briefly summarized in several categories, e.g., landslide mechanisms 3 , 4 , 5 , failure behaviour simulation 6 , 7 , 8 , 9 , sensitivity analysis 10 , 11 , 12 , as well as countermeasure design and optimization 13 , 14 . These studies have underlaid a solid theoretical foundation for landslide mitigation work.

The studies on the landslide mechanisms provide a fundamental understanding of the landslide process. Previous studies 15 , 16 , 17 have long proposed various failure models of landslides. These remarkable studies also substantiated the key factors dominating the landslide process, e.g., topography, material strength 3 , 18 , as well as some other triggering factors, e.g., rainfall precipitation and earthquakes 19 , 20 . These studies explore important features of landslide process and benefit the practical work.

Most studies regarding landslide mechanisms are commonly empirically or semi-empirically based. However, a physical and quantitative analysis of the landslide process is required, because slope stability, as well as magnitude and possible runout extension are the key parameters to the countermeasure design. In this context, some attempts have been made to evaluate slope stability. In order to reproduce the landslides process, it is essential to determine the key parameters of landslides, cohesion strength c , and internal friction angle φ . The two parameters are mainly based on the in-situ test, engineering experience analogy, and back analysis 21 at present. Being an efficient solution to determine the key parameters, the back-analysis method can be briefly summarized into two categories, the deterministic method 22 , 23 and reliability method 24 , 25 , 26 , 27 , 28 . Previous studies 22 also support that the back-analysis provides reliable results approximating the expected shear strength parameters.

In order to estimate the potential risk of slope failure, the potential magnitude and runout extension should be estimated. It is a major issue and a complex task for landslide mitigation work. Numerical simulation provides an alternative solution for this purpose. Several models have been developed and applied to practical work. They use the discontinuous deformation analysis (DDA) method 29 , 30 , 31 , the finite element method (FEM) 32 , 33 , the discrete element method (DEM) 34 , 35 , the smoothed particle hydrodynamics method (SPH) 7 , 9 , to analyse the stability and simulate the failure process of the landslide. Some other studies also incorporate hydraulically based models to simulate landslide behaviour, such as DDA-SPH coupled method 36 and shallow water assumption-based model 37 . However, difficulties remain in the measurement of related parameters in these numerical models. Owing to the significant individual differences from case to case, as well as the temporal and spatial variation, some important parameters are difficult to determine, requiring trial-and-error adjusting during simulation 38 . In this sense, a comprehensive analysis of landslide stability and failure process simulation remains a major topic.

In this paper, we report on a comprehensive method for analysing landslide stability and its failure process. Two essential issues are discussed, i.e., comprehensive estimation of soil strength parameters using the back-analysis results using the reliability method, and the evaluation of the slope stability using the limit analysis and the FEM simulation. The so-called Lanmuxi landslide in China is selected as a case study to illustrate the presented method.

Back analysis of key parameters

In this paper, we first use reliability theory to back analyse the shearing strength parameters. Back analysis is an effective method for determining landslide parameters. For landslide mitigation work, the safety factor F s is commonly used as an indicator for the stability. The slope is likely to fail when F s  ≤ 1, while is under a limit state when F s  = 1. In the reliability method, the function equation is incorporated to represent the limit state of the slope. The function is as below,

where X 1 , X 2 , X 3 , …, X n denotes random variables regarding to slope stability. In our study, only two parameters, the cohesion c and friction angle φ are considered and assumed as random variables. Slope stability can be quantitated by the value of Z . Z  < 0 denotes that the slope is at risk of failure, while Z  = 0 represents that the slope stays in the limit state. When Z  > 0, the slope keeps stable. Another key parameter in reliability theory is the slope reliability index β , which refers to the probability of a slope that completes the pre-determined function under the specified condition and time. Based on previous studies 39 , the reliability index β can be expressed as:

where μ i denotes the averaged value, and σ i represents the standard deviation of the variable. The failure model is established based on the limit state (Z = 0) in the back-analysis of parameters using the reliability method. The final back-analysis result is determined according to the minimum reliability index β   40 , which means, the key parameters X i are iteratively solved until a minimum β is obtained.

The spreadsheet-based method 41 , 42 is used to solve the key parameters X i . The reliability index β is expressed as the forms of Hasofer-Lind’s reliability index 43 ,

where X denotes random variables, μ is the averaged value of the random variables, C x is the covariance matrix of variables, F is the failure domain. Equations ( 1 , 3 ) can be iteratively solved according to the optimization algorithm. Compared with the traditional deterministic back analysis method, reliability method generates better results because the effect of the uncertainties in the landslide parameters can be considered.

Stability analysis based on the 2D upper bound (UB) theory

Presently, the majority of slope stability analyses evaluate a safety factor using a 2D representation of the slope, while 3D analyses of slope stability are much less reported 44 , 45 . It can be explained in part by the fact that a 3D model is likely to introduce a much larger number of degrees-of-freedom, which demands significantly more computational time and effort than the 2D model 46 . Meanwhile, in the 3D model, the failure surface does not only cross weak soil layers, but also strong ones with uncertainties, consequently increases the calculated factor of safety, e.g., 14% to 18% increasement by Reyes and Parra 47 and 13.9% by Xie 48 , compared with the minimum factor of safety in the 2D model. In a word, the 2D analysis could lead to more conservative results than the 3D analysis in slope stability problems.

A primary concern in the practice of landslide mitigation is safety. For this consideration, we use the 2D upper bound (UB) limit analysis (LA) method to evaluate a conservative slope stability in this paper. The UB method is an energy-based method, using the principle of virtual work 49 , 50 . According to the virtual work equations, the internal work of the slope should be equal to its external work (i.e., D int  =  W ext ) in the limit state. A kinematically admissible velocity field 51 is presented for calculating the virtual work. The velocity formula can be described as follows,

where v i and v i −1 are the virtual velocities of a slice on the sliding surface. [ v ] i −1, i represents the relative velocity that defined as the vector difference from v i to v i −1 . It should be noticed that v i , v i −1 and [ v ] i −1, i satisfy the closure relations. φ i and φ i− 1 are the internal friction angles of the neighboured slice, respectively. [ φ ] i −1, i denotes the relative internal friction angle in the vertical direction of slices. α i and α i− 1 are the inclined angles. Therefore, the UB solution of the safety factor can be calculated by the following equation:

where φ f  = arctan(tan φ/F s ). l i is the sliding surface length of the slice i . P w and U represent the external pressure arisen by pore water. P w can be calculated by \({P}_{w}=\frac{1}{2}{\rho }_{w}{Z}_{w}\) , where P w denotes the density of water and the Z w represents the depth of pore water. U is obtained by constructing the equilibrium equation with W i and α i , at vertical direction. W i is the gravity of the slice i . [ h ] i−1, i is the length of the interface between the slice i and the slice i- 1. Equations ( 4 – 6 ) can be further explained by Fig.  1 .

figure 1

The 2D stability analysis based on the UB theory. ( a ) Virtual velocity field of slope divided into vertical slices; ( b ) Pore water pressure on the sliding layer.

Stability analysis using the FEM strength reduction method

Compared with the traditional LA method, the finite element method (FEM) performs better with the advantage of considering much complex boundary conditions, as well as the non-homogeneity of the soil and rock mass. Another advantage of the FEM method is that the stress and deformation field can be also obtained. The FEM strength reduction method has been proposed and applied in the analysis of slope stability 52 , 53 . It regarded the safety factor F S as the reduction degree of shear strength of the soil material when the slope reaches the limit state. The safety factor F S can be redefined as F S  =  c / c f or F S  = tan φ /tan φ f , where c and φ are the initially-input shear strength parameters, while c f and φ f are the output shear strength parameters in the limit state after reduction, respectively.

The definition of the safety factor mentioned above is consistent with the definition introduced by Bishop 54 , that the safety factor F S is expressed as bellow,

where, τ f represents the shear strength of the slope and can be calculated by the Mohr-Coulomb model with the cohesion and internal friction angle. τ denotes the actual shear stress of slope. Equation ( 7 ) can be transformed as follow,

Equation ( 8 ) means that the slope will reach the limit state with the c ′ =  c / F s and tan φ ′ = tan φ / F s , which is the same with the definition mentioned above. Therefore, in our study, the reduction equations of shear strength are expressed as bellow, deducing with the assumption of the constant external load,

where R denotes the reduction coefficient with an initiate value of R  = 1.0. In each iteration step, R increases with the stress and deformation analysis based on FEM implemented. The Mohr-Coulomb model is introduced to describe the constitutive relationship of the landslide mass. Commonly, the vertical boundary of the model is fixed in the horizontal direction, while the bottom boundary is constrained in both horizontal and vertical directions. The iteration of simulation breaks until the slope reaches the limit state. The value of the variable R under the limit state is regarded as the safety factor Fs for the selected landslide profile.

Background of the case study

The so-called Lanmuxi landslide is located in the northwest of Xikou village, Fenghuang town, Hunan province, China (as shown in Fig.  2a ). Evidence of previous landslide had been previously observed on July 2014. Owing to the impact of continuous heavy rainfall, the slope deformed, causing cracks in the crown and middle part of the slope body. Subsequently, on September 2014, two secondary landslides were triggered due to a heavy rainstorm, resulting in new cracks occurred in the crown region of the slope, with a maximum width up towards to 0.50 m. The total length of the landslide area is approximately 100.00 m. The maximum width and average thickness of the landslide reach 206.00 m and 4.40 m, respectively. The area of this landslide is estimated 1.37 × 10 4 m 2 , with a total mass volume of 6.08 × 10 4 m 3 . The slope body and the sliding layer mainly composed of silty clay, while the substrate consists of argillaceous siltstone. Geological investigation shows that there are only tiny amounts of pore water in the stratum and no groundwater is observed.

figure 2

Overview of the Lanmuxi landslide. ( a ) Location of the landslide. ( b ) Map of the Lanmuxi landslide (using Grapher® 10, https://www.goldensoftware.com/products/grapher ). ( c ) Photographic view (Photograph was taken by Z. Han).

The in-situ investigation after the secondary landslide event illustrates that the slope is still unstable, threatening 134 residents and 272 local buildings, with a total potential economic loss of 1.8 million dollars. The overview of the Lanmuxi landslide is shown in Fig.  2b,c .

Back analysis of the shearing strength parameters

In order to describe the state of the landslide under diverse condition reasonably, two different conditions are considered. In Condition 1, the mass of the slope is dry, while in Condition 2, it is supposed to be saturated under the impact of heavy rainfall. The saturated soil mass of the slope increases self-weight and consequently reduces shearing strength, leading to a greater risk for landslide failure.

According to the laboratory test on the three groups of soil material samples that obtained in-situ , the mean values of the parameters in this case are listed in Table  1 . In the back analysis using the UB theory and the reliability method, the standard deviations are required. The laboratory test shows that the standard deviations of the shear strength parameters in the sliding layer are ± 2.98 kPa and ± 1.04° in the Condition 1, while ± 3.02 kPa and ± 1.96° in the Condition 2. The back analysis based on the deterministic method and reliability method are conducted, respectively. Results are shown in Table  2 .

Laboratory test and back analysis are comprehensively considered to obtain the parameter of the sliding layer. We use the laboratory test results in Table  1 and the back-analysis results in Table  2 to attain the averaged values of the parameters in Table  3 , which are suggested as the best-fitting parameters for stability analysis and FEM simulation.

Stability analysis

As mentioned in the above section, it has been widely accepted that the 2D analysis could lead to more conservative results than the 3D analysis in slope stability problems. For the safety consideration, we evaluate a conservative slope stability using a 2D model. To simplify the actual 3D slope into a 2D model, the length along the plane direction, i.e., C1-C1′ and D1-D1′ in Fig.  2b is hereafter referred to as the two typical slope profiles. Most of the local buildings and infrastructures are distributed at the landslide toes along the directions of the selected profiles. The slopes along the both profiles are presumed to be infinitely wide in the 2D model, negating the 3D effects caused by the infinite width of the sliding mass.

Owing to the significant reduction of the safety factor when soil mass is saturated, we mainly focus on the stability analysis under Condition 2. Two analysis methods, the UB method and the FEM strength reduction method, are conducted, respectively. The results in Table  4 illustrate that the landslide along Profile 2 is unstable ( F s  < 1).

Stress and deformation analysis using the FEM method

In order to perform stress and deformation analysis using the FEM method, the both profiles are meshed into triangle blocks (Fig.  3 ). Profile 1 is meshed into 1494 blocks with 2515 nodes, while Profile 2 is separated into 1209 blocks with 2029 nodes. In each iteration step of the FEM simulation, stress and deformation are calculated and landslide strength is gradually reduced. Iteration breaks until the slope section comes to the limit state.

figure 3

Finite element mesh. ( a ) Profile 1 (C1-C1′). ( b ) Profile 2 (D1-D1′).

Figure  4 and 5 illustrate the stress and deformation distribution. Figure  4a–d demonstrate the distribution of shear stress, effective stress, strain, and shear strain of Profile 1. While the plastic strain, the deformation along both horizontal and vertical directions under the limit state are shown in Fig.  4e–g . In contrast, Fig.  5 only shows the limit state of Profile 2, because the slope along this profile is unstable when soil mass is saturated, limiting the FEM simulation converged. As such, we only use the Eq. ( 9 ) to calculate the plastic strain and deformations of this profile in the limit state. Different from the strength reduction method which increases the reduction coefficient R as illustrated in Eq. ( 9 ), in order to obtain a converged FEM simulation results, we use a decreased R in each iteration step in the analysis of Profile 2. The shear strength consequently increases in each step until the unstable slope along Profile 2 reaches a limit state. Figure  4 indicates that the stress and deformation of Profile 1 majorly concentrate at the toe, while Fig.  5 reveals that the stress and strain of Profile 2 in the limit state majorly concentrate at the crown.

figure 4

The FEM analysis results of Profile 1. ( a ) Shear stress. ( b ) Effective stress. ( c ) Strain. ( d ) Shear strain. ( e ) Plastic strain (the limit state). ( f ) Deformation in x direction (the limit state). ( g ) Deformation in z direction (the limit state).

figure 5

The FEM analysis results of Profile 2. ( a ) Plastic strain (the limit state) ( b ). Deformation in x direction (the limit state) ( c ). Deformation in z direction (the limit state).

Suggestion of countermeasures

The comprehensive results of stability analysis demonstrate that the slope along Profile 1 (C1-C1′) approximates the limit state, while Profile 2 (D1-D1′) is unstable. Therefore, in order to protect local buildings at the downslope area, countermeasures are required. As shown in Figs  4 and 5 , the FEM analysis indicates the weak parts of the slope along both profiles. The stress and strain majorly concentrate on the toe and main body along the Profile 1, as well as the landslide crown along the Profile 2. These weak parts are supposed to subject obvious surface deformation that up towards to 142 mm as shown in Fig.  5c . In this context, we suggest structural strengthening at these weak parts, using anchor lattice beams at the landslide body, anti-slide piles at the landslide toe, as well as intercepting drains and cracks filling at the landslide crown. The length of the anchorage section is 4 m, and the anchoring force exerted by each anchor is 500 kN. The overall configuration of the landslide countermeasures is shown in Figs  6 and 7 .

figure 6

The plane layout of landslide mitigation measures (the figure was generated by Grapher 11, version 11.7.825, https://www.goldensoftware.com/products/grapher ).

figure 7

The profile layout of landslide mitigation measures. ( a ) Profile 1. ( b ) Profile 2.

To evaluate the effect of the countermeasures, an effect assessment is conducted in this paper by using the 2D UB method and the FEM strength reduction method. The analysis results are shown in Fig.  8 and Table  5 . Figure  8 reveals an obvious control of the slope deformation after settling the structural strengthening. The maximum positive deformation in x direction of Profile 1 (C1-C1′) is reduced from 166.8 mm to 41.6 mm, while 162 mm to 95.2 mm of Profile 2 (D1-D1′). Table  5 shows that the safety factors of Profile 1 and 2 have been increased 43.6% and 43.8% respectively after the slope reinforcement. These results have demonstrated the feasibility and effectiveness of the recommended countermeasures.

figure 8

The FEM analysis results of slope deformation in x-direction (after the slope reinforcement). ( a ) Profile 1. ( b ) Profile 2.

Limitations

For the simplification of calculation, the reliability method introduced in this study is performed with the assumption that the shear strength parameters are independent to each other. Therefore, the negative correlation between these parameters as revealed in previous studies 55 , 56 , and the impact of spatial correlation of shear strength remains unknown and were not considered in this study. The above limitation should be considered in future works.

Another limitation is with respect to the FEM simulation. Presently, effects of tension cracks are not considered in our analysis. However, previous study 57 applied the kinematic approach of limit analysis to assess the stability of uniform cohesive friction slopes with cracks, indicating that failure mechanisms departing from the crack tip can lead to a significant overestimation of the stability of the slope. For this reason, the slope profile 1 in the case study may be instable in view of many observed tension cracks in the crown. Improvement on this issue is ongoing.

The third limitation in this paper relates to the 2D simplification of the actual complex 3D slope. The composition of the slope is usually heterogeneous and, in combination with the complicated soil layers, it is commonly difficult to select appropriate profiles to analyse its stability. In the 2D limit analysis of slope stability, the slope is presumed to be infinitely wide along the out-of-plane direction, negating the 3D effects caused by the infinite width of sliding mass. In consequence, the 2D slope stability analysis may be over-conservative and insufficient. In this sense, a large and realistic 3D model should be built and analysed in order to overcome such apparent instability found in the 2D analysis.

In this paper, we use a comprehensive analysis method to evaluate the stability of the Lanmuxi landslide and discuss the related countermeasures. The presented method has advantages in considering the uncertainties of the soil shear strength parameters, and generating more conservative results compared to the deterministic method.

The rational values of the sliding layer parameters are comprehensively determined by using the mean value of the back-analysis result of the reliability method, as well as the results by the deterministic method and the in-situ test. Following the determination of parameters, the 2D UB method of limit analysis and the FEM strength reduction method are performed for stability evaluation. The safety factors Fs along two typical profiles of the slope are calculated, indicating that the slopes along both profiles are unstable. The FEM-based analysis furthermore demonstrates the weak parts of the slope where stress and strain concentrate.

Structural countermeasures, using anchor lattice beams, landslide piles, and cracks filling, are suggested based on the comprehensive analysis results. Subsequently, an effect assessment based on UB method and FEM simulation is implemented. Results show notable decreasing of the deformation and about 40% increasing of the safety factors, which demonstrates the feasibility and effectiveness of the recommended countermeasures.

Petley, D. Global patterns of loss of life from landslides. Geology. 40 (10), 927–930 (2012).

Article   ADS   Google Scholar  

Ma, N., Ma, H. M. & Feng, W. Q. Interaction between countermeasure works and landsliding mass: a case study. The 10th Asian Regional Conference of IAEG. 1–6 (2015).

Leroueil, S. Natural slopes and cuts: movement and failure mechanisms. Geotechique. 51 (3), 197–243 (2001).

Article   Google Scholar  

Take, W. A., Bolton, M. D., Wong, P. C. P. & Yeung, F. J. Evaluation of landslide triggering mechanisms in model fill slopes. Landslides. 1 (3), 173–184 (2004).

Wang, F. W., Sassa, K. & Wang, G. Mechanism of a long-runout landslide triggered by the august 1998 heavy rainfall in Fukushima prefecture, Japan. Engineering Geology. 63 (1), 169–185 (2002).

Wu, J. H. Seismic landslide simulations in discontinuous deformation analysis. Computers & Geotechnics. 37 (5), 594–601 (2010).

Article   CAS   Google Scholar  

Pastor, M., Haddad, B., Sorbino, G., Cuomo, S. & Drempetic, V. A depth‐integrated, coupled SPH model for flow‐like landslides and related phenomena. International Journal for Numerical & Analytical Methods in Geomechanics. 33 (2), 143–172 (2010).

Article   ADS   MATH   Google Scholar  

Yuan, R. M., Tang, C. L., Hu, J. C. & Xu, X. W. Mechanism of the donghekou landslide triggered by the 2008 wenchuan earthquake revealed by discrete element modeling. Natural Hazards & Earth System Sciences. 14 (5), 3895 (2014).

Han, Z. et al . Numerical simulation of debris-flow behavior based on the SPH method incorporating the Herschel-Bulkley-Papanastasiou rheology model. Engineering Geology. 255 , 26–36 (2019).

Ashis, K. S., Ravi, P. G. & Irene, S. An approach for GIS-based statistical landslide susceptibility zonation—with a case study in the Himalayas. Landslides. 2 (1), 61–69 (2005).

Catani, F., Lagomarsino, D., Segoni, S. & Tofani, V. Landslide susceptibility estimation by random forests technique: sensitivity and scaling issues. Natural Hazards & Earth System Sciences. 13 (11), 2815–2831 (2013).

Han, Z. et al . An integrated method for rapid estimation of the valley incision by debris flows. Engineering Geology. 232 (8), 34–45 (2018).

Zaruba, Q. & Mencl, V. Landslides and their control. Engineering Geology. 7 (3), 281–282 (1982).

Google Scholar  

Wei, Z., Li, S., Wang, J. G. & Ling, W. A dynamic comprehensive method for landslide control. Engineering Geology. 84 (1), 1–11 (2006).

Terzaghi, K. Mechanism of landslides. Application of Geology to Engineering Practice. 83–123 (1950).

Wieczorek, G.F. Chapter 4-landslide triggering mechanisms. Landslides: investigation and mitigation. (1996).

Huang, R. Q. Large-scale lan 357 dslide and their sliding mechanisms in China since the 20th century. Chinese Journal of Rock Mechanics and Engineering . 26 (3),433–454 (in Chinese) (2007).

Zhou, S. H., Chen, G. Q. & Fang, L. G. Distribution pattern of landslides triggered by the 2014 Ludian earthquake of China: implications for regional threshold topography and the seismogenic fault identification. ISPRS International Journal of Geo-information. 5 (4), 46 (2016).

Iverson, R. M. Landslide triggering by rain infiltration. Water Resources Research. 36 (7), 1897–1910 (2000).

Malamud, B. D., Turcotte, D. L., Guzzetti, F. & Reichenbach, P. Landslide inventories and their statistical properties. Earth Surface Processes and Landforms. 29 (6), 687–711 (2004).

Sonmez, H., Ulusay, R. & Gokeceoglu, C. A practical procedure for the back analysis of slope failures in closely jointed rock masses. Int J Rock Mech Mining Sci. 35 (2), 219–233 (1998).

Tang, W. H., Stark, T. D. & Angulo, M. Reliability in back analysis of slope failures. Soils Found. 39 (5), 73–80 (1999).

Wang, L., Hwang, J. H., Luo, Z., Jiang, C. H. & Xiao, J. Probabilistic back analysis of slope failure—a case study in Taiwan. Comput Geotech. 51 , 12–23 (2013).

Low, B. K. & Tang, W. H. Efficient reliability evaluation using spreadsheet. J Eng Mech. 123 (7), 749–752 (1997).

Low, B. K. & Tang, W. H. Reliability analysis using object-oriented constrained optimization. Struct Saf. 26 (1), 69–89 (2004).

Low, B. K. FORM, SORM, and spatial modeling in geotechnical engineering. Struct Saf. 49 , 56–64 (2014).

Zhang, J., Tang, W. H. & Zhang, L. M. Efficient probabilistic back-analysis of slope stability model parameters. J Geotech Geoenviron. 136 (1), 99–109 (2009).

Article   MathSciNet   Google Scholar  

Zhang, L. L., Zhang, J., Zhang, L. M. & Tang, W. H. Back analysis of slope failure with Markov chain Monte Carlo simulation. Comput Geotech. 37 (7), 905–912 (2010).

Yin, K. L., Jiang, Q. H. & Wang, Y. Numerical simulation 387 on the movement process of Xintan landslide by DDA method. Chinese Journal of Rock Mechanics and Engineering . 21(7), 959–962 (in Chinese) (2002).

Zhang, Y. et al . DDA validation of the mobility of earthquake-induced landslides. Engineering Geology. 194 , 1069–1075 (2015).

Wu, J. H. & Chen, C. H. Application of DDA to simulate characteristics of the Tsaoling landslide. Computers and Geotechnics. 38 (5), 741–750 (2011).

Chen, H. & Lee, C. F. A dynamic model for rainfall-induced landslides on natural slopes. Geomorphology. 51 (4), 269–288 (2003).

Crosta, G. B., Chen, H. & Lee, C. F. Replay of the 1987 Val Pola Landslide, Italian Alps. Geomorphology. 60 (1), 127–146 (2004).

Boon, C. W., Houlsby, G. T. & Utili, S. New insights into the 1963 Vajont slide using 2D and 3D distinct-element method analyses. Geotechnique. 64 (10), 800–816 (2014).

Li, X. P., He, S. M., Luo, Y. & Wu, Y. Simulation of the sliding process of Donghekou landslide triggered by the Wenchuan earthquake using a distinct element method. Environmental Earth Sciences. 65 (4), 1049–1054 (2012).

Wang, W., Chen, G., Zhang, Y., Zheng, L. & Zhang, H. Dynamic simulation of landslide dam behavior considering kinematic characteristics using a coupled DDA-SPH method. Engineering Analysis with Boundary Elements. 80 , 172–183 (2017).

Article   MathSciNet   MATH   Google Scholar  

Han, Z. et al . Numerical Simulation of Debris-flow Behavior Incorporating a Dynamic Method for Estimating the Entrainment. Engineering Geology. 190 , 52–64 (2015).

Han, Z. et al . Numerical simulation for run-out extent of debris flows using an improved cellular automaton model. Bulletin of Engineering Geology & the Environment. 76 (3), 961–974 (2017).

Verma, A. K., Srividya, A. & Karanki, D. R. Reliability and safety engineering. Springer, London. (2010).

Sun, L. C., Wang, H. X. & Zhou, N. Q. Application of reliability theory to back analysis of rocky slope wedge failure. Chin J Rock. Mech Eng. 31 (A01), 2660–2667 (2012).

Low, B. K. Reliability analysis of rock slopes involving correlated nonnormals. Int J Rock Mech Mining Sci. 44 (6), 922–933 (2007).

Zhao, L. H., Zuo, S., Lin, Y. L. & Zhang, Y. B. Reliability back analysis of shear strength parameters of landslide with three-dimensional upper bound limit analysis theory. Landslides. 13 (4), 711–724 (2016).

Lü, Q. & Low, B. K. Probabilistic analysis of underground rock excavations using response surface method and SORM. Comput Geotech. 38 (8), 1008–1021 (2011).

Hicks, M. A., Nuttall, J. D. & Chen, J. Influence of heterogeneity on 3D slope reliability and failure consequence. Computers and Geotechnics. 61 , 198–208 (2014).

Xiao, T., Li, D. Q., Cao, Z. J. & Zhang, L. M. CPT-based probabilistic characterization of three-dimensional spatial variability using MLE. Journal of Geotechnical and Geoenvironmental Engineering. 144 (5), 04018023 (2018).

Liu, Y. et al . Probabilistic stability analyses of undrained slopes by 3D random fields and finite element methods. Geoscience Frontiers. 9 , 1657–1664 (2018).

Reyes, A. & Parra, D. 3D slope stability analysis by the using limit equilibrium method analysis of a mine waste dump. In: Proceedings Tailings and Mine Waste 2014, Keystone, Colorado, USA, October 5–8 (2014).

Xie, M. W., Wang, Z. F., Liu, X. Y. & Xu, B. Three-dimensional critical slip surface locating and slope stability assessment for lava lobe of Unzen volcano. Journal of Rock Mechanics and Geotechnical Engineering. 3 (1), 82–89 (2011).

Chen, W. F. Limit analysis and soil plasticity. Elsevier Scientific Pub. Co. (1975).

Chen, W. F. Limit analysis in soil mechanics. Elsevier Scientific Pub. Co. (1990).

Michalowski, R. L. Slope stability analysis: a kinematical approach. Geotechnique. 45 (2), 283–293 (1995).

Zienkiewicz, O. C., Humpheson, C. & Lewis, R. W. Associated and non-associated visco-plasticity in soil mechanics. Geotechnique. 25 (4), 671–689 (1975).

Duncan, J. M. State of the art: limit equilibrium and finite element analysis of slopes. discussion and closure. Journal of Geotechnical Engineering. 123 (9), 577–596 (1996).

Bishop, A. W. The use of the slip circle in the stability analysis of slopes. Geotechnique. 5 (1), 7–17 (1955).

Cherubini, C. Reliability evaluation of shallow foundation bearing capacity on c′, φ′ soils. Can Geotech J. 37 (1), 264–269 (2000).

Wolff, T. F. Analysis and design of embankment dam slope: a probabilistic approach. University Microfilms. (1985).

Utili, S. Investigation by limit analysis on the stability of slopes with cracks. Geotechnique. 63 (2), 140–154 (2013).

Download references

Acknowledgements

This study was financially supported by the National Key R&D Program of China (Grant No. 2018YFC1505401, Z. Han); the National Natural Science Foundation of China (Grant No. 41702310, Z. Han); the Natural Science Foundation of Hunan (Grant No. 2018JJ3644, Z. Han); the Innovation-Driven Project of Central South University (Grant No. 2019CX011); and the Fundamental Research Funds for the Central Universities of Central South University. These financial supports are gratefully acknowledged.

Author information

Authors and affiliations.

School of Civil Engineering, Central South University, Changsha, Hunan, 410075, China

Zheng Han, Bin Su, Yange Li, Yangfan Ma & Weidong Wang

State Key Laboratory of Geohazard Prevention and Geo-environment Protection, Chengdu University of Technology, Chengdu, Sichuan, 610059, China

Zheng Han & Yange Li

Key Laboratory of Heavy-haul Railway, Ministry of Education, Changsha, Hunan, 410075, China

Weidong Wang

Department of Civil and Structural Engineering, Kyushu University, 819-0395, Fukuoka, Japan

Guangqi Chen

You can also search for this author in PubMed   Google Scholar

Contributions

Z.H., Y.G.L. and B.S. designed the study. Y.F.M. and W.D.W. performed the in-situ investigation. G.Q.C. simulates the performance of the suggested countermeasures. Z.H. and B.S. wrote the manuscript. All authors discussed the results and commented on the manuscript.

Corresponding author

Correspondence to Yange Li .

Ethics declarations

Competing interests.

The authors declare no competing interests.

Additional information

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ .

Reprints and permissions

About this article

Cite this article.

Han, Z., Su, B., Li, Y. et al. Comprehensive analysis of landslide stability and related countermeasures: a case study of the Lanmuxi landslide in China. Sci Rep 9 , 12407 (2019). https://doi.org/10.1038/s41598-019-48934-3

Download citation

Received : 10 December 2018

Accepted : 16 August 2019

Published : 27 August 2019

DOI : https://doi.org/10.1038/s41598-019-48934-3

Share this article

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

This article is cited by

Assessing the effectiveness of landslide slope stability by analysing structural mitigation measures and community risk perception.

  • A. S. M. Maksud Kamal
  • Farhad Hossain
  • Peter Sammonds

Natural Hazards (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines . If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Quick links

  • Explore articles by subject
  • Guide to authors
  • Editorial policies

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

case study on landslides

  • Share full article

Advertisement

Supported by

Deadly Landslides in India Made Worse by Climate Change, Study Finds

Extreme rainfall made 10 percent heavier by human-caused climate change triggered landslides that killed hundreds, according to a new study.

People wearing red shirts walk through muddy area littered with debris

By Austyn Gaffney

A sudden burst of rainfall on July 30 caused a cascade of landslides that buried hundreds of people in the mountainous Kerala region of southern India.

That downpour was 10 percent heavier because of human-caused climate change, according to a study by World Weather Attribution, a group of scientists who quantify how climate change can influence extreme weather. Nearly six inches, or 150 millimeters, of rain fell on soils already highly saturated from two months of monsoon and marked the third highest single-day rain event on record for India.

“The devastation in northern Kerala is concerning not only because of the difficult humanitarian situation faced by thousands today, but also because this disaster occurred in a continually warming world,” said Maja Vahlberg, a climate risk consultant at the Red Cross Red Crescent Climate Centre. “The increase in climate-change-driven rainfall found in this study is likely to increase the number of landslides that could be triggered in the future.”

In a state that is highly prone to landslides, the Wayanad district is considered the riskiest part. As of Tuesday, at least 231 people had died and 100 remained missing.

The Kerala landslides were the second extreme landslide event in July, following one in Ethiopia that killed 257 people. July was the second-worst month on record, after July 2019, with 95 landslide events that caused 1,167 fatalities , according to data maintained by Dave Petley, the vice-chancellor of the University of Hull. Together, they caused roughly one-third of the more than 3,600 deaths resulting from some 429 fatal landslides recorded this year, Dr. Petley said in an email.

Already, 2024 is an outlier, Dr. Petley posted to The Landslide Blog on Tuesday . He wrote that he could “only speculate on the likely underlying reasons for this very high incidence of fatal landslides,” but “the most likely cause continues to be the exceptionally high global surface temperatures, and the resultant increase in high intensity rainfall events.”

We are having trouble retrieving the article content.

Please enable JavaScript in your browser settings.

Thank you for your patience while we verify access. If you are in Reader mode please exit and  log into  your Times account, or  subscribe  for all of The Times.

Thank you for your patience while we verify access.

Already a subscriber?  Log in .

Want all of The Times?  Subscribe .

Information

  • Author Services

Initiatives

You are accessing a machine-readable page. In order to be human-readable, please install an RSS reader.

All articles published by MDPI are made immediately available worldwide under an open access license. No special permission is required to reuse all or part of the article published by MDPI, including figures and tables. For articles published under an open access Creative Common CC BY license, any part of the article may be reused without permission provided that the original article is clearly cited. For more information, please refer to https://www.mdpi.com/openaccess .

Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for future research directions and describes possible research applications.

Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive positive feedback from the reviewers.

Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world. Editors select a small number of articles recently published in the journal that they believe will be particularly interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the most exciting work published in the various research areas of the journal.

Original Submission Date Received: .

  • Active Journals
  • Find a Journal
  • Proceedings Series
  • For Authors
  • For Reviewers
  • For Editors
  • For Librarians
  • For Publishers
  • For Societies
  • For Conference Organizers
  • Open Access Policy
  • Institutional Open Access Program
  • Special Issues Guidelines
  • Editorial Process
  • Research and Publication Ethics
  • Article Processing Charges
  • Testimonials
  • Preprints.org
  • SciProfiles
  • Encyclopedia

sustainability-logo

Article Menu

case study on landslides

  • Subscribe SciFeed
  • Recommended Articles
  • Google Scholar
  • on Google Scholar
  • Table of Contents

Find support for a specific problem in the support section of our website.

Please let us know what you think of our products and services.

Visit our dedicated information section to learn more about MDPI.

JSmol Viewer

Study on landslide susceptibility based on multi-model coupling: a case study of sichuan province, china.

case study on landslides

1. Introduction

2. study area and methodology, 2.1. study area, 2.2. mathematical and statistical methods.

  • Information Value Method

2.3. Machine Learning Models

  • Decision Tree C5.0 Model

3. Data and Impact Factors

3.1. data sources, 3.2. selection of impact factors, 3.3. covariance diagnostics, 4. results analysis and accuracy validation, 4.1. mathematical statistical models, 4.2. coupled models, 4.3. distribution of landslide conditions, 4.4. roc curve accuracy validation, 5. conclusions, author contributions, institutional review board statement, informed consent statement, data availability statement, acknowledgments, conflicts of interest.

YearNumber of Geological Disasters (Incidents)LandslidesDebris FlowsGround Subsidence
20059736706024
200617598284
2007295156466
200877074883235138
20096205339544525
201093458110817
20112161148234518
20121997141833015
20133149226746611
201427584424429
2015275818554429
20163491531210
2017227119702
201817564760
20195632381251
20207253022870
2021251317374472
2022403239853
202314854581
  • Ali, S.A.; Parvin, F.; Vojteková, J.; Costache, J.; Linh, N.T.T.; Pham, Q.B.; Vojtek, M.; Gigović, L.; Ahmad, A.; Ghorbani, M.A. GIS-Based Landslide Susceptibility Modeling: A Comparison between Fuzzy Multi-Criteria and Machine Learning Algorithms. Geosci. Front. 2020 , 12 , 857–876. [ Google Scholar ] [ CrossRef ]
  • Napoli, M.D.; Martire, D.D.; Bausilio, G.; Calcaterra, D.; Confuorto, P.; Firpo, M.; Pepe, G.; Cevasco, A. Rainfall-Induced Shallow Landslide Detachment, Transit and Runout Susceptibility Mapping by Integrating Machine Learning Techniques and GIS-Based Approaches. Water 2021 , 13 , 488. [ Google Scholar ] [ CrossRef ]
  • Huang, Y.; Zhao, L. Review on landslide susceptibility mapping using support vector machines. Catena 2018 , 165 , 520–529. [ Google Scholar ] [ CrossRef ]
  • Sun, D.; Gu, Q.; Wen, H.; Shi, S.; Mi, C.; Zhang, F. A Hybrid Landslide Warning Model Coupling Susceptibility Zoning and Precipitation. Forests 2022 , 13 , 827. [ Google Scholar ] [ CrossRef ]
  • Wang, B.; Li, S.; Xu, W.; Yang, Y.; Li, Y. A Comparative Study of Landslide Susceptibility Evaluation Based on Three Different Machine Learning Algorithms. A Comparative Study of Landslide Susceptibility Evaluation Based on Three Different Machine Learning Algorithms. Northwestern Geol. 2024 , 57 , 34–43. [ Google Scholar ]
  • Jia, J.; Mao, Y.; Meng, X.; Gao, B.; Gao, M.; Wu, W. Comparison of Landslide Susceptibility Evaluation by Deep Random Forest and Random Forest Model: A Case Study of Lueyang County, Hanzhong City. Northwestern Geol. 2023 , 56 , 239–249. [ Google Scholar ]
  • Li, P.; Ye, H.; Tan, S. Evaluation of Geological Hazards in Yongde County Based on Analytic Hierarchy Process. Res. Soil Water Conserv. 2021 , 28 , 394–406. [ Google Scholar ]
  • Yi, J.; Wang, F.; Cheng, Y.; Zhang, Y. Study on the risk assessment of geological disasters in alpine valley area: A case study in Aba County, Sichuan Province. Chin. J. Geol. Hazard Control 2022 , 33 , 134–142. [ Google Scholar ]
  • Moragues, S.; Lenzano, M.G.; Lanfri, M.; Moreiras, S.; Lannutti, E.; Lenzano, L. Analytic hierarchy process applied to landslide susceptibility mapping of the North Branch of Argentino Lake, Argentina. Nat. Hazards 2021 , 105 , 915–941. [ Google Scholar ] [ CrossRef ]
  • Chen, Y.; Yang, H.; Zhong, Y. Evaluation of landslide susceptibility based on IV-IOE-LR coupling model—A case study of Lushui city, Yunnan. Ind. Miner. Process. 2023 , 52 , 62–79. [ Google Scholar ]
  • Fan, H.; Lu, Y.; Hu, Y.; Fang, J.; Lv, C.; Xu, C.; Feng, X.; Liu, Y. A Landslide Susceptibility Evaluation of Highway Disasters Based on the Frequency Ratio Coupling Model. Sustainability 2022 , 14 , 7740. [ Google Scholar ] [ CrossRef ]
  • Zhang, S.; Tan, S.; Zhou, J.; Sun, Y.; Ding, D.; Li, J. Geological Disaster Susceptibility Evaluation of a Random-Forest-Weighted Deterministic Coefficient Model. Sustainability 2023 , 15 , 12691. [ Google Scholar ] [ CrossRef ]
  • Shirani, K.; Pasandi, M.; Arabameri, A. Landslide susceptibility assessment by Dempster–Shafer and Index of Entropy models, Sarkhoun basin, Southwestern Iran. Nat. Hazards 2018 , 93 , 1379–1418. [ Google Scholar ] [ CrossRef ]
  • Wen, H.; Wu, X.; Ling, S.; Sun, C.; Qiang, L.; Zhou, G. Characteristics and susceptibility assessment of the earthquake-triggered landslides in moderate-minor earthquake prone areas at southern margin of Sichuan Basin, China. Bull. Eng. Geol. Environ. 2022 , 81 , 346. [ Google Scholar ] [ CrossRef ]
  • Bhandari, B.P.; Dhakal, S.; Tsou, C.-Y. Assessing the Prediction Accuracy of Frequency Ratio, Weight of Evidence, Shannon Entropy, and Information Value Methods for Landslide Susceptibility in the Siwalik Hills of Nepal. Sustainability 2024 , 16 , 2092. [ Google Scholar ] [ CrossRef ]
  • Holloway, J.; Mengersen, K. Statistical Machine Learning Methods and Remote Sensing for Sustainable Development Goals: A Review. Remote Sens. 2018 , 10 , 1365. [ Google Scholar ] [ CrossRef ]
  • Liu, X.; Shao, S.; Shao, S. Landslide susceptibility prediction and mapping in Loess Plateau based on different machine learning algorithms by hybrid factors screening: Case study of Xunyi County, Shaanxi Province, China. Adv. Space Res. 2024 , 74 , 192–210. [ Google Scholar ] [ CrossRef ]
  • Agboola, G.; Beni, L.H.; Elbayoumi, T.; Thompson, G. Optimizing landslide susceptibility mapping using machine learning and geospatial techniques. Ecol. Inf. 2024 , 81 , 102583. [ Google Scholar ] [ CrossRef ]
  • Huang, F.; Cao, Z.; Jiang, S.; Zhou, C.; Huang, J.; Guo, Z. Landslide susceptibility prediction based on a semi-supervised multiple-layer perceptron model. Landslides 2020 , 17 , 2919–2930. [ Google Scholar ] [ CrossRef ]
  • Su, Y.; Chen, Y.; Lai, X.; Huang, S.; Lin, C.; Xie, X. Feature adaptation for landslide susceptibility assessment in ”no sample“ areas. Gondwana Res. 2024 , 131 , 1–17. [ Google Scholar ] [ CrossRef ]
  • Li, Y.; Wang, X.; Mao, H. Influence of human activity on landslide susceptibility development in the Three Gorges area. Nat. Hazards 2020 , 104 , 2115–2151. [ Google Scholar ] [ CrossRef ]
  • Orhan, O.; Bilgilioglu, S.S.; Kaya, Z.; Ozcan, A.K.; Bilgilioglu, H. Assessing and mapping landslide susceptibility using different machine learning methods. Geocarto Int. 2020 , 37 , 2795–2820. [ Google Scholar ] [ CrossRef ]
  • Yang, C.; Liu, L.; Huang, F.; Huang, L.; Wang, X. Machine learning-based landslide susceptibility assessment with optimized ratio of landslide to non-landslide samples. Gondwana Res. 2023 , 123 , 198–216. [ Google Scholar ] [ CrossRef ]
  • Fan, H.; Lu, Y.; Shao, S.; Li, L.; Wang, Y.; Lu, M.; Li, J.; Yao, K.; Sun, Y. Evaluation and analysis of statistical and coupling models for highway landslide susceptibility. Geomat. Nat. Hazards Risk 2023 , 14 , 2167612. [ Google Scholar ] [ CrossRef ]
  • Wang, H.; Xu, J.; Tan, S.; Zhou, J. Landslide Susceptibility Evaluation Based on a Coupled Informative–Logistic Regression Model—Shuangbai County as an Example. Sustainability 2023 , 15 , 12449. [ Google Scholar ] [ CrossRef ]
  • Chen, W.; Peng, J.; Hong, H.; Shahabi, H.; Pradhan, B.; Liu, J.; Zhu, A.; Pei, X.; Duan, Z. Landslide susceptibility modelling using GIS-based machine learning techniques for Chongren County, Jiangxi Province, China. Sci. Total Environ. 2018 , 626 , 1121–1135. [ Google Scholar ] [ CrossRef ]
  • Ngo, P.T.T.; Panahi, M.; Khosravi, K.; Ghorbanzadeh, O.; Kariminejad, N.; Cerda, A.; Lee, S. Evaluation of deep learning algorithms for national scale landslide susceptibility mapping of Iran. Geosci. Front. 2021 , 12 , 505–519. [ Google Scholar ]
  • Kavzoglu, T.; Teke, A. Predictive Performances of ensemble machine learning algorithms in landslide susceptibility mapping using random forest, extreme gradient boosting (XGBoost) and natural gradient boosting (NGBoost). Arab. J. Sci. Eng. 2022 , 47 , 7367–7385. [ Google Scholar ] [ CrossRef ]
  • Xue, Y.; He, S.; Wu, S. Assessment of Social Vulnerability to Natural Disasters in Sichuan Province. North China Earthq. Sci. 2018 , 36 , 33–40. [ Google Scholar ]
  • Wang, W.; He, Z.; Han, Z.; Li, Y.; Dou, J.; Huang, J. Mapping the susceptibility to landslides based on the deep belief network: A case study in Sichuan Province, China. Nat. Hazards 2020 , 103 , 3239–3261. [ Google Scholar ] [ CrossRef ]
  • Hou, F.; Wu, G.; Zhang, S. Evaluation of landslide susceptibility in Jinsha River Gorge section based on information method. Sichuan Arch. 2019 , 39 , 147–150. [ Google Scholar ]
  • Zhang, Q.; Ling, S.; Li, X.; Sun, C.; Xu, J.; Huang, T. Comparison of landslide susceptibility mapping rapid assessment models in Jiuzhaigou County, Sichuan province, China. Rock Mech. Rock Eng. 2020 , 39 , 1595–1610. [ Google Scholar ]
  • Xu, C.; Dai, F.; Yao, X.; Zhao, Z.; Xiao, J. GIS platform and certainty factor analysis method based Wen Chuan earthquake-induced landslide susceptibility evaluation. Eng. Geol. 2010 , 18 , 15–25. [ Google Scholar ]
  • Zhang, Z. A comparative study of landslide susceptibility in Wangmu County coupled with MaxEnt modeling. J. Chifeng Univ. 2022 , 38 , 108–113. [ Google Scholar ]
  • Kadavi, P.R.; Lee, C.; Lee, S. Landslide-susceptibility mapping in Gangwon-do, South Korea, using logistic regression and decision tree models. Environ. Earth Sci. 2019 , 78 , 116. [ Google Scholar ] [ CrossRef ]
  • Zhao, J.; Chen, H.; Yang, S.; Ma, Z. Landslide risk assessment based on decision tree arithmetic. J. Zhejiang Univ. 2004 , 31 , 465–470. [ Google Scholar ]
  • Sun, D.; Wen, H.; Xu, J.; Zhang, Y.; Wang, D.; Zhang, J. Improving Geospatial Agreement by Hybrid Optimization in Logistic Regression-Based Landslide Susceptibility Modelling. Front. Earth Sci. 2021 , 9 , 713803. [ Google Scholar ] [ CrossRef ]
  • Li, K.; Su, Q.; Zhang, X.; Fan, J.; Bai, D. Evaluation of Landslide Susceptibility in Xishan Coalfield Based on Information Value Method Coupled with Machine Learning Model. Radio Eng. 2024 , 54 , 390–401. [ Google Scholar ]
  • Che, W.; Qin, S.; Miao, Q.; Su, G.; Chen, Y.; Yao, J. Study on Factor Classification Method in the Evaluation of landslide sensitivity. In Proceedings of the 11th National Engineering Geology Conference, Wuhan, China, 16–18 October 2020. [ Google Scholar ]
  • Günther, A.; Reichenbach, P.; Malet, J.; Eeckhaut, M.V.D.; Hervás, J.; Dashwood, C.; Guzzetti, F. Tier-based approaches for landslide susceptibility assessment in Europe. Landslides 2013 , 10 , 529–546. [ Google Scholar ] [ CrossRef ]
  • Wu, X.; Wang, L.; Ding, Z.; Zhang, X.; Zhou, Y.; Wang, Y. The susceptibility evaluation of township landslide in HengDuan mountain area of Yunnan. Geol. Hazards Environ. Preserv. 2023 , 4 , 13–19. [ Google Scholar ]
  • Liu, M. Research on Sensitivity of Influence Factors on Soil Landslide Stability—Taking Qingchuan Earthquake-Stricken Areas as An Example. Master’s Thesis, Chengdu University of Technology, Chengdu, China, June 2011. [ Google Scholar ]
  • Liu, M.; Chen, X.; Yang, S. Zonation of landslide risk with logistic regression model and certainty factor. Eng. Geol. 2014 , 22 , 1250–1256. [ Google Scholar ]
  • Wang, J. Landslide Risk Assessment in Wanzhou Country, Three Gorges Reservoir. Ph.D. Thesis, China University of Geosciences, Wuhan, China, June 2015. [ Google Scholar ]
  • Wang, G.; Lei, X.; Chen, W.; Shahabi, H.; Shirzadi, A. Hybrid Computational Intelligence Methods for Landslide Susceptibility Mapping. Symmetry 2020 , 12 , 325. [ Google Scholar ] [ CrossRef ]
  • Kadirhodjaev, A.; Rezaie, F.; Lee, M.; Lee, S. Landslide Susceptibility Assessment Using an Optimized Group Method of Data Handling Mode. ISPRS Int. J. Geo-Inf. 2020 , 9 , 566. [ Google Scholar ] [ CrossRef ]
  • Mai, J.; Xi, Y.; Liu, G. Predicting Potential Rainfall- Triggered Landslides Sites in Guangdong Province (China) using MaxEnt Model under Climate Changes Scenarios. Int. J. Geog. Inf. Sci. 2021 , 23 , 2042–2054. [ Google Scholar ]
  • Moore, I.D.; Gearson, R.B.; Ladson, A.R. Digital terrain modelling: A review of hydrological, geomorphological, and biological applications. Hydrol. Proc. 1991 , 5 , 3–30. [ Google Scholar ] [ CrossRef ]
  • He, Q.; Shahabi, H.; Ataollah, S.; Li, S.; Chen, W.; Wang, N.; Chai, H.; Bian, H.; Ma, J.; Chen, Y.; et al. Landslide spatial modelling using novel bivariate statistical based Naïve Bayes, RBF Classifier, and RBF Network machine learning algorithms. Sci. Total Environ. 2019 , 663 , 1–15. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Zeng, T.; Wang, L.; Zhang, Y.; Cheng, P.; Wu, F. Landslide susceptibility modeling and interpretability based on CatBoost-SHAP model. Chin. J. Geol. Hazard Control 2024 , 35 , 37–50. [ Google Scholar ]
  • Zhao, Z. Landslide Susceptibility Evaluation Based on Spatial Database and Machine Learning: A Case Study of Xinyuan County. Master’s Thesis, Xinjiang University, Ürümqi, China, June 2020. [ Google Scholar ]
  • Qiu, H.; Regmi, A.D.; Cui, P.; Cao, M.; Lee, J.; Zhu, X. Size distribution of loess slides in relation to local slope height within different slope morphologies. Catena 2016 , 145 , 155–163. [ Google Scholar ] [ CrossRef ]
  • Cantarino, I.; Carrion, M.A.; Goerlich, F.; Ibañez, V.M. A ROC analysis-based classification method for landslide susceptibility maps. Landslides 2019 , 16 , 265–282. [ Google Scholar ] [ CrossRef ]

Click here to enlarge figure

FactorToleranceVIFFactorToleranceVIF
Elevation0.5831.717SPI0.3123.205
Slope0.2204.545TWI0.3263.066
Aspect0.9921.008TPI0.7051.419
Plan curve0.6341.577Surface roughness0.3043.291
Profile curve0.6011.663TWI0.8381.193
Valley depth0.6071.646Slope height0.6031.658
Precipitation0.5971.675
Landslide Evaluation FactorsClassificationNumber of Landslide Points/ptsClassified Area/km IVCF
Elevation (m)190~94713,838139,307.130.767310.56164
947~1894586452,958.200.875910.61173
1894~2862249054,838.27−0.01553−0.01614
2862~364996470,326.52−1.21324−0.71253
3649~426129894,711.69−2.68492−0.93471
4261~71436197,758.19−4.30280−0.98709
Slope (°)0~7.566344133,091.410.031790.03280
7.56~16.337590122,597.760.293240.26646
16.33~25.095335117,720.17−0.01871−0.01941
25.09~34.77299792,670.65−0.35612−0.30963
34.77~77.10124543,090.92−0.46886−0.38542
Aspect−1~34.39199744,950.31−0.03860−0.03962
34.39~69.77230948,636.400.027760.02870
69.77~105.17273755,293.060.069530.07042
105.17~141.97282255,917.750.088880.08916
141.97~178.77259451,414.750.088590.08888
178.77~215.58218348,367.48−0.02281−0.02362
215.58~250.97209049,126.10−0.08191−0.08214
250.97~286.36225054,098.21−0.10456−0.10358
286.36~321.74220250,808.30−0.06338−0.06419
321.74~359.96232750,558.54−0.00324−0.00339
Plan curve−8.15~−0.27187046,660.81−0.14039−0.13645
−0.27~0.1415,420324,460.310.030080.03106
0.14~6.636225138,778.88−0.02773−0.02864
Profile curve−6.01~−0.24227770,903.92−0.36189−0.31372
−0.24~0.1914,813325,580.89−0.01353−0.01408
0.19~6.266425113,415.190.205720.19493
Valley depth−585~−276.43934104,816.84−1.64392−0.81403
−276.43~519.093355135,899.21−0.62489−0.47644
519.09~749.616953135,408.950.107450.10680
749.61~1004.40514981,580.200.313790.28235
1004.40~1344.12437039,889.340.865230.60704
1344.12~2508.88275412,305.461.579590.83232
Precipitation (mm)303.25~732.01155444,519.13−0.27837−0.25177
732.01~886.362685102,990.07−0.57023−0.44624
886.36~1017.855470129,631.35−0.08870−0.08862
1017.85~1149.346068104,914.830.226600.21256
1149.34~1309.40527086,028.960.284070.25924
1309.40~1761.03246641,845.090.245340.22808
SPI−13.81~−10.28166045,712.77−0.23991−0.22134
−10.28~−5.836010134,011.18−0.02884−0.02977
−5.83~−1.614883107,532.39−0.01638−0.01702
−1.61~0.446767155,616.51−0.05968−0.06057
0.44~3.52284354,786.870.117060.11582
3.52~15.26134911,705.660.914950.62847
TWI2.05~5.207430196,949.33−0.20177−0.18988
5.20~6.979709190,209.150.100580.10032
6.97~9.53375281,459.60−0.00216−0.00227
9.53~13.36170131,524.240.156110.15153
13.36~27.129209223.060.770570.56325
TPI−310.39~−14.75138625,236.750.174610.16796
−14.75~−6.766377111,620.640.214080.20203
−6.76~1.2310,301243,522.16−0.08648−0.08650
1.23~11.894663106,020.09−0.04748−0.04850
11.89~368.8079023,499.23−0.31620−0.28052
Surface roughness1~1.0614,575265,304.230.174130.16754
1.06~1.155766138,199.78−0.10102−0.10027
1.15~1.30247284,701.59−0.45841−0.37876
1.30~1.5662419,267.98−0.35435−0.30836
1.56~5.72741892.43−0.16586−0.15905
FVC0~0.5018018,681.15−1.56552−0.79872
0.50~0.73151657,141.71−0.55266−0.43616
0.73~0.889286189,037.800.063360.06436
0.88~112,528245,042.680.103330.10291
Slope height0~84.5518,467302,731.600.279710.25580
84.55~225.464140139,212.17−0.43873−0.36603
225.46~469.7079555,600.24−1.17102−0.69996
469.70~2395.4911312,355.98−1.61793−0.80909
ModelsSusceptibility LevelClassified Area/km Proportion of
Classified Area/%
Classified Area/km Proportion of the
Number of
Landslide Points/%
Density of Landslide
Points/(pts/km )
IVVery low58,864.9540.116520.0020.00088
Low88,686.6910.1741740.0070.00196
Moderate100,583.1070.1988030.0340.00798
High96,681.7580.19047030.2000.04864
Very high163,505.7220.32217,7600.7560.10862
IV-LRVery low97,641.3310.192610.0030.00062
Low94,576.9970.1862970.0130.00314
Moderate69,544.9010.1378960.0380.01288
High62,252.7410.12230320.1290.04870
Very high185,138.6110.36419,2250.8180.10384
IV-DTVery low73,635.9410.145550.0020.00075
Low76,111.4740.1501610.0070.00212
Moderate91,643.1120.1805320.0230.00581
High91,303.7040.18036840.1570.04035
Very high175,628.0020.34619,0600.8110.10852
CFVery low57,115.0220.1121290.0050.00226
Low113,965.8190.2246120.0260.00537
Moderate123,232.5360.24225850.1100.02098
High99,130.8100.19570030.2980.07064
Very high114,878.0450.22613,1630.5600.11458
CF-LRVery low97,641.3310.192610.0030.00069
Low94,576.9970.1862970.0130.00244
Moderate69,544.9010.1378960.0380.01512
High62,252.7410.12230320.1290.04765
Very high185,138.6110.36419,2250.8180.10299
CF-DTVery low73,635.9410.145550.0020.00166
Low76,111.4740.1501610.0070.00404
Moderate91,643.1120.1805320.0230.01515
High91,303.7040.18036840.1570.06848
Very high175,628.0020.34619,0600.8110.11182
The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

Zhang, J.; Qian, J.; Lu, Y.; Li, X.; Song, Z. Study on Landslide Susceptibility Based on Multi-Model Coupling: A Case Study of Sichuan Province, China. Sustainability 2024 , 16 , 6803. https://doi.org/10.3390/su16166803

Zhang J, Qian J, Lu Y, Li X, Song Z. Study on Landslide Susceptibility Based on Multi-Model Coupling: A Case Study of Sichuan Province, China. Sustainability . 2024; 16(16):6803. https://doi.org/10.3390/su16166803

Zhang, Jinming, Jianxi Qian, Yuefeng Lu, Xueyuan Li, and Zhenqi Song. 2024. "Study on Landslide Susceptibility Based on Multi-Model Coupling: A Case Study of Sichuan Province, China" Sustainability 16, no. 16: 6803. https://doi.org/10.3390/su16166803

Article Metrics

Article access statistics, further information, mdpi initiatives, follow mdpi.

MDPI

Subscribe to receive issue release notifications and newsletters from MDPI journals

COMMENTS

  1. A hill slope failure analysis: A case study of Malingoan ...

    This study focused on following objectives: (1) identify the typical hill slope parameters and valuranable soil properties; (2) application of Arc GIS for identification of vulnerability of hill slope failure; (3) suggest preventive measures to avoid landslide.

  2. Life and death of slow-moving landslides | Nature Reviews ...

    This Review discusses the recent advances in understanding slow-moving landslides and examines circumstances where slow-moving landslides can accelerate and fail catastrophically.

  3. Multi-source remote sensing-based landslide investigation ...

    This study presents a case study approach to fully leverage variety of multi-source remote sensing technologies for analyzing the characteristics of a landslide. The selected case is a landslide...

  4. Landslide case studies - British Geological Survey

    Landslide case studies. The landslides team at the BGS has studied numerous landslides. This work informs our geological maps, memoirs and sheet explanations and provides data for our National Landslide Database , which underpins much of our research.

  5. Comprehensive analysis of landslide stability and related ...

    In order to illustrate the presented method, the so-called Lanmuxi landslide in China is selected as a case study. Results demonstrated that the stress and strain majorly concentrated at the...

  6. Deadly Landslides in India Made Worse by Climate Change ...

    Extreme rainfall made 10 percent heavier by human-caused climate change triggered landslides that killed hundreds, according to a new study. Rescuers searching through mud and debris after ...

  7. Settlements and slides: a large landslide case study from the ...

    Abstract. The Central Cordillera of Luzon in the Philippines is home to some of the most complex geology, active plate margin tectonism, heaviest rainfall and steepest terrain in the world.

  8. Case Study: Oso, Washington, Landslide of March 22, 2014 ...

    This paper describes investigation, testing, analysis, and slope history used to determine the two-phase failure mechanism involved in the 2014 landslide near Oso, Washington. The first phase involves a slide mass located above the frequent landslides in the lower portion of the slope and extends to near the slope crest.

  9. Post-failure process and kinematic behavior of two landslides ...

    This paper presents a comprehensive and well-documented landslide case study involving detailed site investigation, field tests, seismic signal analyses, and advanced numerical analyses using the material point method to investigate the post-failure process and kinematic behavior of two landslides.

  10. Sustainability | Free Full-Text | Study on Landslide ...

    Landslides are among the most prevalent geological hazards and are characterized by their high frequency, significant destructive potential, and considerable incident rate. Annually, these events lead to substantial casualties and property losses. Thus, conducting landslide susceptibility assessments in the regions vulnerable to such hazards has become crucial. In recent years, the coupling of ...