Abstract

Quantitative risk assessment of landslides induced by slope failure is an important precondition for formulating effective disaster prevention, mitigation measures, and establishing a landslide risk warning system. In general, the location of the critical slip surface and the failure mode is unlikely to be predicted due to the spatial variability in soil. It remains a challenging task to effectively identify the critical slip surface and conduct the efficient risk assessment based on a three-dimensional (3-D) slope with spatial variability. Based on Monte Carlo simulation and the random field method, a quantitative risk evaluation method for slope failure considering the spatial variability of soil parameters is proposed in the study. Compared with a uniform soil slope, the landslide volume, the critical slip surface, and the factor of safety considering the spatial variability of soil are all uncertain; thus, the soil spatial variability has a significant effect on the failure mode and stability of the slope. By using the random finite element method, the critical slip surface of the slope is accurately identified, the corresponding landslide volume and slide distance are calculated, and the modified risk index for a landslide is further enriched, which can provide the reference basis for predicting the landslide deformation, quantitatively evaluating the landslide risk, and mitigating the landslide disaster.

1. Introduction

In mountain regions, landslides are one of the most dangerous natural hazards, which can result in devastating and fatal disasters to human society and urban agglomeration [14]. One of the most destructive landslide events is the 2020 Pettimudi landslide in India, which resulted in more than 66 fatalities and destroyed and buried severe village infrastructures. The occurrence dangers are affected or dominated by the concealment and irregularity of the critical slip surface of the slope, the spatial variability of geotechnical parameters, and the temporal and spatial variation of various boundary conditions [5, 6]; Qi and Li 2018 [7], where the spatial variability of geotechnical parameters is one of the most vital uncertainties affecting the slope stability [8]. Griffiths et al. [9] demonstrated that soils inherently exhibiting spatial variability can impact slope stability significantly on account of the complicated geological processes during their formation. Jiang et al. [10] illuminated the great necessity for quantitatively assessing the risk of landslide hazards induced by slope instability considering the soil spatial variability. Cheng et al. [11] expounded that many geotechnical designers should not only take into account the slope stability but also the failure consequence caused by slope instability from an engineering perspective.

Quantitative risk assessment of slope instability plays a significant role in slope design and risk mitigation, and most of the previous studies have paid close attention to quantitative risk assessment methods in geotechnical engineering (e.g., [1214]. Huang et al. [15] developed an efficient framework for quantitative risk assessment of slope instability where the failure consequence for each failure mode of the slope is assessed respectively; however, the influence of the correlation between slip surfaces has not been considered effective by the method. Li and Chu [16] proposed a quantitative risk assessment approach for two-dimensional (2-D) slope lope failure that can identify a series of representative slip surfaces by limit equilibrium method, where this method only considers the slope failure risk with the horizontal spatial variability of soil parameters. The typical investigations on risk assessment of landslide hazards due to slope failure are summarized in Table 1.

Table 1 illustrates that the landslide hazard risk assessment mainly includes two elements. One is the possibility of slope instability, namely failure probability. The other is the impact of landslides on people, property, society, and the environment, which can be summarized as failure consequences. Obviously, the failure probability and the failure consequences are directly related to the failure mode of the slope. The effective identification of the failure mode of the slope is a key step in the quantitative assessment of landslide risk. Due to the spatial variability of soil parameters, multiple failure modes may exist in the slope and the effective identification of the critical slip surface is a complicated problem. Overall, how to quantify the risk assessment of landslide hazards induced by slope failure considering the soil spatial variability remains an intricate challenge.

Previous researches on slope stability analysis and risk assessment considering soil spatial variability were dominated by 2-D analysis (e.g., [24]; Li et al. 2016; [25]. However, the 2-D model cannot reflect the failure mode of slope more realistically in terms of the shape, location, and length of the critical slip surface. Liu et al. [26] illuminated that slope failure will occur in areas with low soil shear strength due to the spatial variability of soil parameters, which directly affects the shape and position of the slip surfaces of the slope. Meanwhile, how to effectively identify the potential and critical failure slip surface of a 3-D slope and quantitatively calculate the landslide volume is still a complicated problem. Thus, it is necessary to investigate the stability analysis and risk assessment of 3-D slope with the spatial variability of soil properties.

Several studies (e.g., [22, 27, 28] have made great efforts to estimate 3-D slope reliability. Xiao et al. [5] proposed an auxiliary random finite element method to quantify the risk of 3-D slope considering spatial variability of soil properties. Vanmarcke et al. [29] pioneered the quantitative risk assessment of 3-D slope that is regarded as an extension of 2-D slope reliability analysis based on the first-order second-moment method. However, previous studies on the quantitative risk assessment of slope failure require a presumed critical failure surface of the slope, which becomes impractical when the soil spatial variability is taken into account.

In this study, the problem that effectively identifies the exact location of slip surface and estimate the probability of failure in a 3-D slope consisting of spatially random soils can be solved by the finite element analysis (denoted by FEA), where the irregular slip surface caused by the spatial variation can be captured automatically in a 3-D slope without a prescribed slip surface and the volume of slide mass can therefore be calculated exactly. Based on Monte Carlo simulation and the random field method, a quantitative risk evaluation method for slope failure considering the spatial variability of soil parameters is proposed. By calculating the corresponding landslide volume and travel distance, a modified risk index is proposed to evaluate the landslide hazard, where the risk index can be used for a regional landslide hazard assessment. This research aims to provide a practical method for quantitative risk assessment of slope failure, formulate the risk prevention, and control plan of landslide disasters.

2. Finite Element Model

The quantitative risk assessment of the 3-D slope in spatial variability soil is conducted by the finite element method. Figure 1(a) illustrates the 3-D slope model that consists of 43,200 eight-node brick elements and the soil undrained shear strength (denoted by cu) is characterized by a random field. The corresponding 2-D plane strain model is shown in Figure 2(a), where the plane strain boundary conditions are also in accord with the 3-D finite element model. Since the properties of materials are strictly positive, a lognormal random field is applied to embody the random variables so, that guarantees all soil properties are nonnegative. Therefore, a 3-D lognormal random field is adopted to characterize the undrained shear strength, where the lognormal random field is obtained from a Gaussian random field that is generated by the modified linear estimation method (see [30]).

Figure 1(b) illustrates the contour of the maximum principal plastic strain of a typical sample. A small element size is used to capture the potential slip surface by the finite-element program in each simulation, where the sensitivity of results to mesh refinement in a 3-D slope is verified by Griffiths and Marquez [2]. The spatial correlation length along the vertical directions (denoted by Θy) and the spatial correlation length along the horizontal direction (denoted as Θx) refer to the values in Li et al. [31]. Xiao et al. [8] demonstrated that the failure mode and stability of 3-D slopes can be significantly influenced by the horizontal spatial variability. Thus, the Θx is considered to explore its effect on the overall performance of landslides. The soil is regarded as an elastic-perfectly-plastic material based on the Mohr-Coulomb yield criterion. The parameters displayed in Table 2 are assumed to be deterministic except for the undrained shear strength, which is referred to as those reported by Liu and Shields [32].

Due to the spatial randomness of the soil, the potential slip surfaces and the failure modes varied from one simulation (or realization) to another. For this reason, 600 realizations of the lognormal random field for the 3-D full-model were generated in this study; that is, Monte-Carlo simulations were conducted to examine the performance of those models [33]. For each simulation, the factor of safety (FS) is calculated by the strength reduction method, and the volume of the slide mass is calculated by summing the volume of elements whose displacement is greater than a specific value that is set at 0.2 m in this study. The displacement of element and volume of slide mass can be calculated by equations (1)–(3). Thus, 600 FSs and the corresponding volume of the slide mass can be obtained for each scenario (i.e., each implementation of random field), where the volume of the slide mass can be considered as an “equivalent quantity” to characterize the failure consequence of slope in landslide risk assessment [18]. For the models shown in Figure 3, the elements with nongrey color formed the slide region of the landslide, which suggests the failure mode of the slope is dominated by the spatial variability of soil parameters. The calculation procedure of the quantitative risk assessment of a 3-D slope instability considering the spatially variable soils parameters is schematically shown in Figure 4.where uk denotes the deformation displacement of the kth element greater than a specific value of displacement; uki is the displacement of the ith node of the kth element; xki,yki, and zki are the coordinate of ith node of the kth element, respectively; ΔVk is the volume of the kth element; V is the volume of slide mass.

3. Results

The irregular slip surface caused by the spatial variation of soil parameters may fail at any location of the slope, which can lead to additional difficulty in identifying the exact location of the critical slip surface and calculating the volume of slide mass. However, the critical slip surface can be automatically identified by the finite element method in each simulation prior to a prescribed slip surface. Thus, it is reasonable to accurately evaluate the failure probability of the slope and the corresponding failure consequence that depends on the area (2-D model) or volume (3-D model) of the slide mass. To evaluate the risk of landslide hazard induced by slope failure with the spatial variables in soil, 600 realizations of different random strength fields are generated so that can obtain 600 values of the FSs and the corresponding consequence of slope failure. In the study, the landslide volume and the corresponding travel distance are regarded as the evaluation indicators of the failure consequences of slope instability, and probabilistic parametric studies on landslide hazard risk are conducted for the cases where the Θx value is presented in Table 2.

3.1. Comparison between 2-D and 3-D Slope Stability Analyses

To preferably assess the landslide hazard induced by the slope failure with the spatial variability of soil, we herein introduce a concept of the standardized landslide volume, which is defined as the ratio of slide volume of a slope with spatially random strength over the slide volume of a slope with uniform strength. Thus, the standardized landslide volume can be treated as a simple indicator to characterize the consequence of the landslide hazard. Figure 5 displays the histogram of the 600 values of the standardized landslide volume on the basis of 2-D and 3-D analyses with various Θx values. The histograms suggest that the spatial variability of soil has a prominent influence on the perspective of the standardized landslide volume compared to the deterministic model (a blue arrow) whether it is a 3-D or a 2-D analysis. The reason is that the material parameter characteristic values of nonuniform slopes are different in different regions, the region with stronger or weaker soil strength is uncertain from one simulation to the next, however, slope failure will occur in the region with a low shear strength of soil, which can directly affect the location and shape of the critical slip surface and result in the irregular failure modes of slope and these uncertainties will inevitably affect the volume of the slide mass (see Figure 3). This result also emphasizes the necessity of considering the spatial variation in soils in slope stability analysis.

For cases presented in Table 2, the standardized landslide volume range of the 2-D model is larger, and the range of standardized landslide volume of 3-D is more concentrated. The findings can be explained that the 3-D model can more realistically reflect the location and shape of the critical slip surface and failure consequences of the slope compared to the 2-D model. The range of the standardized landslide volume in Figure 5 is overtly wider when Θx is 1000 m compared with those of the models that Θx is 2 m for both 2-D and 3-D models. These findings indicate that the location of the critical slip surface of the slope with a larger value of Θx will be tough to be identified which results in the difficulty of quantitative assessment of the landslide volume. This phenomenon can be reflected by the COV of the calculated results summarized in Table 3. For instance, when Θx is 1000 m in the 3-D model, the COV of the standardized landslide volume is as high as 0.19. Such a high COV value indicates that the volume of the slide mass is uneasy to be evaluated, which can make the landslide risk assessment sophisticated in reality. The reason is that a larger value of Θx indicates a stronger correlation between neighboring zones and the possibility of a local failure becomes greater on account of the inherent existence of randomness in soil (Li et al. 2021).

The difference between the 2-D model and the 3-D model (e.g., the standardized landslide volume) is inversely proportional to the Θx. The failure mechanism of landslides can be factually reflected by the 3-D slope in terms of the shape, position, and length of the critical slip surface compared to the 2-D model. For a conservative estimation, the infinitely large horizontal spatial correlation length is recommended to quantitatively evaluate the landslide hazard taking into account the inherent spatial variability in soils. The tail of the histograms in Figure 5 lies on the left-hand side of the distribution; that is, the data are left-skewed, which can be verified by the skewness of all cases with random soils. The skewness of the calculated results summarized in Table 3 is generally negative. This finding suggests that some smaller portion of the slide (i.e., local failure) may come up in reality. From the perspective of standardized landslide volume, the study reveals that the spatial variability of soil parameters has played a significant role in the quantitative assessment of landslide risk.

The empirical cumulative distributions of FS of all the cases are plotted in Figure 6. Figure 6 illustrates that the stability of the slope can be dominated by the spatial variability of soil parameters when compared to the deterministic model (black dashed line). However, the traditional slope stability analysis method demands a presumed critical failure surface of the slope, which becomes impractical when the soil spatial variability is taken into account. The spatial variability of soil parameters can result in the critical failure surface being irregular and unknown, which makes the probability of a slope local failure become greater. The slip surface of local failure is small, and the energy required for failure is low, which can result in a small FS for slope. Moreover, Figure 6 clearly illuminates the conservativeness of the 2-D analysis from the perspective of FS, which is in accord with the results of previous researchers [2, 23].

3.2. Stability Analyse and Failure Consequence for 3-D Slope

For a 3-D slope with spatially random strength, the slide volume varies from one simulation (or realization) to another. The variables plotted in Figure 7 can examine the relationship between the slide volume fraction (defined as slide volume over the total volume) and the stability (FS) of the 3-D model. As Figure 7 shows that the deterministic result is located near the most likely occurrence point in terms of the slide volume fraction, however, the FS calculated by the deterministic FEA is significantly greater than the results of random analysis, which can be verified by the mode of the calculated results in Table 3. Especially, for the case shown in Figure 7(a), where the spatial correlation length in both directions is set at 2 m, the FS obtained by random FEA are all lower than the deterministic result. The FS obtained from random FEA increases to a certain extent with the increases of Θx but still less than the deterministic result in general, which can manifest the assumption that a uniform soil is likely to overestimate the stability of a slope in terms of FS.

The spatial variability of soil properties has a remarkable effect on the slide volume fraction, especially as shown in Figure 7(d). The COV of the slide volume fraction can reach up to 0.19 and the slide volume fraction change from 15.2% to 75.8% when Θx is 1000 m. The larger COV indicates that the evaluation of slide volume of the slope is of complexity, which increases the difficulty of quantitative assessment of the risk of landslide hazard. The slide volume fraction of random finite element analysis is still different from the deterministic result with the decrease of Θx. Overall, the spatial variability of soil parameters can impact observably the slope stability and failure mode. Figure 7 indicates that the slide volume fraction can be treated as being independent of the FS, which can make the landslide hazard assessment be complicated and changeable. In addition, when the horizontal spatial correlation length is greater than the model size, the results of cases show display similar characteristics, which can be demonstrated by the histogram in Figures 5 and 6.

4. Risk Assessment

Quantitative risk assessment of slope instability is an important basis for landslide disaster risk management and the formulation of various disaster prevention and mitigation measures. The main purpose of carrying out a special landslide disaster risk assessment is to conduct risk management and control, mitigate or eliminate landslide risks, and avoid the hazards caused by the landslide disasters to the safety of personnel and property in complex geological conditions. The effect of spatial variability of soil properties on slope risk may be profound on account of the failure consequences associated with different slope failure modes generally differ. Li et al. (2021) illustrate that the traditional equation may not be directly applicable to a 3-D slope risk assessment with a large number of potential slip surfaces due to the spatial variability of soil properties [11]. When the spatial variability exists in the undrained shear strength, the potential slip surface is unknown and makes the consequence of slope failure also uncertain. However, the failure probability and the failure consequences are directly connected with the failure mode of the slope; thus, the effective identification of the critical slip surface is a crucial step in the quantitative risk assessment of landslide hazard. The potential and irregular slip surface can be automatically captured by the finite element software so that the quantitative risk assessment of slope failure with the spatially variable soils can be implemented [3]. Since the 2-D model cannot reflect the location and shape of the critical slip surface, the failure mode of the slope is a large-scale 3-D problem. Thus, a quantitative risk assessment method for 3-D slope failure considering the spatially variable soils is proposed based on Monte Carlo simulation and the random field method, which can provide valuable risk assessment information for project designers.

Huang et al. [15] elaborated that the failure consequence of slope depends on the landslide volume. Meanwhile, Scheidegger [34] demonstrated that the landslide volume is one of the predominant factors on landslide travel distance. The prediction of landslide travel distance is a significant prerequisite for formulating rational strategies to improve disaster prevention and relocation of residents and infrastructures facilities nearby. Thus, the travel distance can be regarded as the landslide hazard intensity factor, and the modified risk index R for landslide hazard is proposed in this study, where R combines the probability of slope instability and the corresponding failure consequence.

The empirical formula for predicting the value of travel distance is as follows:where H is the maximum fall elevation of the slope; L denotes the travel distance; V is the slide volume.where Ci denotes the consequences of landslide hazard; Li denotes the travel distance of ith; Lt is the travel distance induced by the total failure of the slope; Ri and FSi are the risk index and the ith simulation of FS, respectively; FSmin is the minimum FS calculated out of the Monte-Carlo simulations. In this manner, we can readily evaluate the average risk and maximum risk out of Monte-Carlo simulations, which can be employed to categorize the landslide hazard. Thus, it is viable to effectively evaluate the risk induced by slope failure, which can result in the landslide risk assessment is not complicated. By doing so, the risk index can make the serious consequences of landslide hazards be easy-to-understand, which enhances the risk perception and communication for geotechnical practitioners.

Figure 8 demonstrates the histogram of R of the cases listed in Table 4, where the risk index R is generally less than 0.95. The acceptable value of R should be determined by the practitioner based on the specific site in a particular geological environment. The results in the histograms can provide powerful evidence that the horizontal spatial variability of soil parameters has a vital role in determining the slope risk, which can be proved by the COV of the calculated results summarized in Figure 8. Figure 9 shows the empirical cumulative distribution functions of the 3-D slope model under various Θx values and can obtain a risk assessment table (see Table 4). According to the table, different risk levels can be categorized. By this means, more quantitative information for risk assessment can be provided rather than merely giving an FS or the consequence of slope failure.

5. Conclusions

This study proposed a simple but effective approach for quantitative risk assessment of slope instability considering the spatial variability of soil properties based on Monte Carlo simulation and the random field method. The critical slip surface of the slope is accurately identified the corresponding landslide volume and travel distance are calculated by the random finite element method. By comparing the 2-D model and 3-D model, the effect of horizontal spatial variability on slope stability and risk assessment of landslides is explored. Three main findings can be summarized:(1)Compared to a uniform soil slope, the landslide volume, the critical slip surface, and FS of slope considering the spatial variability of soil are all uncertain, which cast lights the role that the soil with inherently spatial variability has a vital effect on the stability of the slope and landslide hazard assessment(2)A longer horizontal spatial length can lead to a larger range of slope failure consequences. For a conservative estimation, it is suggested that the infinite horizontal spatial correlation length should be considered in the quantitative risk assessment of landslide hazards(3)A modified risk index proposed in this study combines the probability of slope instability and corresponding failure consequence, which can be adapted to the risk assessment of landslide hazard under various working conditions, and provide a reference for the emergency plan of landslide hazard and formulating the risk prevention and control plan of landslide disaster

However, the current study merely focuses on the effect of the spatial variability of soil properties on the landslide risk assessment, and the landslide hazard induced by rainfall or earthquakes is not conducted. These issues are equally important for the quantitative assessment of landslide disaster risk and are worthy of further study. In addition, a more complete and specific assessment of the damage consequences of landslide disasters is a future research objective.

Data Availability

The data are available on request to the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Key Research and Development Program of China (Grant No. 2018YFC1505005) and the Natural Science Foundation of Hunan Province, China (Grant No. 2021JJ40201), and Science and Technology Progress and Innovation Project of Transport Department of Hunan Province (202009).