Next Article in Journal
The Identification and Quantification of Hidden Hazards in Small Scale Reservoir Engineering Based on Deep Learning: Intelligent Perception for Safety of Small Reservoir Projects in Jiangxi Province
Previous Article in Journal
Dual Method for Comprehensive Evaluation of Sustainable Water Resources’ Utilization Capacity in Huangshui River in Yellow River Basin, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Modified Method for Evaluating the Stability of the Finite Slope during Intense Rainfall

1
Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Water 2024, 16(20), 2877; https://doi.org/10.3390/w16202877
Submission received: 6 September 2024 / Revised: 7 October 2024 / Accepted: 8 October 2024 / Published: 10 October 2024
(This article belongs to the Section Hydrology)

Abstract

:
The Green–Ampt (GA) model is a widely used analytical method to calculate the depth of the wetting front during intense rainfall. However, it neglects the existence of the transition layer and the seepage parallel to the slope surface. Therefore, a modified stratified Green–Ampt (MSGA) model is proposed. A process to assess the stability of the finite slope during a rainfall event is demonstrated by combining the MSGA model and the limit equilibrium method. In the case of the Liangshuijing landslide, the factor of safety presents a negative correlation with the depth of the wetting front. The factor of safety obtained by the stratified Green–Ampt (SGA) model is smaller than that calculated by the MSGA model, and the gap between the factor of safety based on the two methods widens with time. The moving speed of the wetting front accelerates with the increase in the length of the slope surface, and the size effect becomes apparent when the length is short. In the initial stage of infiltration, the effect of the seepage parallel to the slope surface is small. The effect of the seepage cannot be neglected at the latter stage. The result calculated by the MSGA model agrees well with the measured result in the test.

1. Introduction

Landslides, one of the most serious and widespread geological disasters, pose a severe threat to the safety of the public. The main induced factor of the landslide is heavy rainfall [1,2,3]. During the rainfall infiltration process, the gravity of the soils and the seepage force will increase continuously with the development of the wetting front. In addition, the internal friction angle and cohesion of the soils will decrease. Numerous studies have shown that the factor of safety of the slope reduces gradually with rainfall infiltration [4,5,6]. Rainwater infiltration can greatly affect the stability of the slope. Therefore, it is essential to simulate the rainfall infiltration process precisely.
Richards’ equation model and GA are the two most popular methods used to simulate the process of rainfall infiltration [7,8]. Richards’ equation is a numerical model derived from the mass conservation law and Darcy’s law, which has been used widely [9,10,11]. However, using Richards’ equation to simulate the infiltration process requires many implicit parameters, which are sometimes unavailable [12]. Moreover, the numerical solution requires constant iteration and a large amount of calculation [13]. The GA model is an analytic method, another popular method used to simulate rainfall infiltration [14,15]. Since the concept of the GA model is clear and concise, the GA model is widely used [16,17].
The GA model was initially proposed to simulate the rainfall infiltration process in the plane. Then, it was extended to the slope by Chen and Young [18]. Wang [19] improved the GA model by considering the seepage parallel to the slope surface. The original GA model assumes that the soil above the wetting front is saturated without considering the existence of the transition layer. Yao [20] modified the GA model by using an elliptic curve to describe the water content of the transition layer. Jiang and Liu [21] proposed an improved Green–Ampt model considering the heterogeneity of the soil. However, previous research did not consider both the transition layer and seepage parallel to the slope surface at the same time. Therefore, a modified stratified Green–Ampt (MSGA) model is proposed that takes the transition layer and the seepage parallel to the slope surface into account at the same time. A stability evaluation method considering the seepage force for a finite slope during a rainfall event is given.
In this paper, the MSGA model is proposed to simulate rainfall infiltration in a finite slope considering both the transition layer and the seepage parallel to the slope surface. A process to evaluate the time-dependent stability of the finite slope is demonstrated based on the MSGA model and limit equilibrium method. This paper is organized as follows: The background of the GA model and one-dimensional Richards’ equation is introduced in Section 2. The proposed MSGA model and time-dependent method for evaluating the finite slope stability are demonstrated in Section 3. Section 4 gives an example of using the modified method for evaluating the stability of a finite slope during intense rainfall. Several discussions and the validation of the MSGA model are included in Section 5. The main conclusions are summarized in Section 6: Conclusions.

2. Theoretical Background

2.1. GA Model

Green and Ampt [7] proposed the GA model based on Darcy’s law to describe the infiltration process. The GA model assumes that the soil profile is homogeneous and that the initial water content is uniformly distributed. It also assumes that there is a wetting front in the profile. The soil above the wetting front is saturated, while the soil below the wetting front is unsaturated. The water content of the unsaturated soil is equal to the initial water content value. Figure 1a shows the original GA model. However, the original GA model did not consider the influence of the slope dip angle. Chen and Yong [18] modified the GA model by extending it from the plane to the slope. The modified model is shown in Figure 1b. When β is 0 degrees, the modified model degenerates into the original model. The original model is a special case of the modified model. In this paper, the infiltration equation of the GA model has considered the effect of the dip angle of the slope surface, which can be expressed by Equation (1).
i = K s Z f cos 2 β + H cos β + S f Z f cos β
where i represents the infiltration rate for ponded infiltration (m/h), Ks represents the saturated hydraulic conductivity (m/h), Zf represents the depth of the wetting front in the vertical direction (m), H represents the depth of the ponding water in the vertical direction (m), Sf represents the suction head of the wetting front (m), and β represents the dip angle of the slope surface (°).
The GA model is by far one of the best models available to describe infiltration during rainfall events [22]. Since the GA model is concise and simple, it has been widely modified and applied to slopes for evaluating the stability of landslide masses [23].

2.2. Distribution of Water Content Based on the SGA Model

The GA model assumes that there are only saturated and natural layers in the water content profile. However, this is not consistent with observations. According to the observation from Colman and Bodman [24], there is a transition layer between the saturated layer and the natural layer. The water content of the soil in the transition layer gradually changes from saturated to the initial water content. According to the experiments and theoretical analysis, it was found that the thickness of the saturated layer approximates to half of the depth of the wetting front [25,26]. Additionally, Wang [25] proposed to use the ellipse function to describe the transition layer, which has been used and validated by other researchers [20,27].
Considering the existence of the transition layer, Yao [20] proposed the SGA model to describe the infiltration of rainwater. The SGA model assumes that there are three layers in the soil water content profile, and the thickness of the saturated layer equals half of the depth of the wetting front. In addition, it assumes that the soil water content varies with depth as an ellipse function. The distribution of the water content in the soil profile can be defined by Equations (2) and (3):
θ ( Z ) = { θ s , 0 Z Z s θ i + ( θ s θ i ) 1 ( Z Z s Z s ) 2 , Z s Z Z f θ i , Z Z f
Z s = Z f 2
where θ denotes water content, θi denotes initial water content, θs denotes saturated water content, Zs denotes the thickness of the saturated layer (m), and Z denotes the depth (m). According to the SGA model, the water content distribution is shown in the Figure 2:

2.3. One-Dimensional Richards’ Equation

The one-dimensional Richards’ equation is widely used to simulate the infiltration of rainwater in unsaturated soil, and the equation is solved by the numerical method. In the case of one-dimensional vertical infiltration on the plane, Richards’ equation can be expressed as follows:
θ t = z ( D θ z ) d K d θ θ z
where z denotes the vertical spatial coordinate (m), D denotes soil water diffusivity (m2/h), K denotes the permeability coefficient (m/h), and t denotes time (h). According to Chen and Young [18], when the one-dimensional Richards’ equation is used on the slope, it can be expressed as follows:
θ t = z * ( D θ z * ) d K d θ θ z * cos β
where z* denotes the coordinate in the normal direction to the slope surface (m), and β denotes the dip angle of the slope (°). The diagram of the coordinates is shown in Figure 3.

3. Methodology

3.1. MSGA Model for a Finite Slope

The water exchange among the slope consists of the recharge from the rainfall infiltration and the discharge of the seepage parallel to the slope surface. For the infinite slope, since the volume of the recharge from the infiltration of rainwater is much greater than the discharge of the seepage parallel to the slope surface, it is reasonable to neglect the influence of seepage parallel to the slope surface. However, the seepage parallel to the slope surface should be considered for the finite slope. Few expansions of the GA model have considered both the existence of the transition layer and the seepage parallel to the slope surface for the finite slope. Therefore, the MSGA model is proposed. While soil parameters are widely recognized to exhibit spatial variability, most of the modified GA models assume that soil is homogeneous for simplicity [19,28,29]. The MSGA model assumes that soil is homogeneous. The MSGA model considers the seepage parallel to the slope based on the SGA model. The schematic diagram of the MSGA model for a finite slope during a rainfall event is shown in Figure 4. Usually, the ponding water head is small compared to the suction head of the wetting front, so the value of the ponding water head is assumed to be equal to 0 m.
The infiltration process during intense rainfall can be divided into two stages depending on the relation between rainfall intensity and infiltration capacity. At stage 1, the rainfall intensity is less than the infiltration capacity. At stage 2, the rainfall intensity is greater than the infiltration capacity. The infiltration ability of the slope can be expressed as follows:
i a = K s Z f cos 2 β + S f Z f cos β
where ia denotes infiltration ability (m/h).
Since the permeability coefficient has an exponential relation to the water content, the permeability coefficient in the saturated state is much greater than that in the unsaturated state. Therefore, the MSGA model only considers the seepage in the saturated layer at stage 2.
At stage 1, which represents the initial stage of the rainfall infiltration, the shallow soil of the slope is unsaturated. Rainwater infiltrates into the soil entirely because the infiltration capacity is greater than the rainfall intensity. Therefore, the infiltration rate is dependent on the rainfall intensity, which can be calculated by:
i = q cos β
in which q denotes the rainfall intensity (m/h).
With the infiltration of rainwater, the state of shallow soil changes from unsaturated to saturated. When the wetting front forms, this is called the critical time. The infiltration ability (ia) is equal to the projection of the rainfall intensity on the slope surface (qcosβ) at the critical time. At this time, the depth of the wetting front can be calculated by combing Equations (6) and (7), which is expressed as follows:
Z f c = S f ( q / K s 1 ) cos 2 β
where Zfc is the depth of the wetting front at the critical time (cm). Based on Equations (2) and (3), the cumulative infiltration depth is defined by:
I = ( θ s θ i ) Z s cos β + Z s Z f ( θ θ i ) cos β d z = ( θ s θ i ) Z f 2 cos β + Z f / 2 Z f ( θ s θ i ) cos β 1 ( z Z f / 2 Z f / 2 ) 2 d z = ( 4 + π ) ( θ s θ i ) 8 cos β Z f
where I is the cumulative infiltration depth (m). The cumulative infiltration depth at the critical time and the critical time can be calculated as follows:
I c = ( 4 + π ) ( θ s θ i ) 8 cos β Z f c = ( 4 + π ) ( θ s θ i ) S f 8 ( q / K s 1 ) cos β
t c = I c q cos β = ( 4 + π ) ( θ s θ i ) S f 8 q ( q / K s 1 ) cos 2 β
where Ic is the cumulative infiltration depth at the critical time (m), and tc is the critical time (h).
Stage 1 changes to stage 2 when the critical time is reached. At stage 2, the shallow soil is saturated. The moving speed of the wetting front is controlled by the rainfall infiltration and the seepage parallel to the slope surface. The effect of the rainfall infiltration on the moving speed of the wetting front can be expressed as follows:
( d Z f d t ) 1 = 8 ( θ s θ i ) ( 4 + π ) cos β ( d I d t ) 1 = 8 K s ( Z f cos 2 β + S f ) ( θ s θ i ) ( 4 + π ) Z f cos 2 β
According to the law of conservation of mass and Darcy’s law, the following equation can be obtained:
K s sin β Z s cos β d t = ( 4 + π ) ( θ s θ i ) 8 L d Z f cos β
where L is the length of the surface for the finite slope (m). Therefore, the effect of the seepage on the moving speed of the wetting front can be expressed as follows:
( d Z f d t ) 2 = 4 K s Z f sin β ( θ s θ i ) ( 4 + π ) L
The actual infiltration speed of the wetting front is calculated by combining Equations (12) and (14), which is expressed as follows:
d Z f d t = 8 K s ( θ s θ i ) ( 4 + π ) ( Z f cos 2 β + S f Z f cos 2 β Z f sin β 2 L )
Since the expression for Zf obtained by integrating Equation (15) is complicated, it is suggested to use the numerical method to calculate Zf. The depth of the wetting front at stage 2 is calculated by the Runge-Kutta 4th-order method, with the initial condition that Zf = Zfc and t = tc.

3.2. Stability Analysis for Various Sliding Surfaces

Slope failure induced by rainfall frequently occurs, and the shallow failure is always parallel to the slope surface [11,30]. The shallow failure may occur along the wetting front or the interface between the saturated layer and the transition layer. Besides shallow failure, failure also occurs along the surface of the bedrock.
For the slip parallel to the slope surface along the wetting front, the factor of safety can be calculated by:
F s = ( G s + G t ) cos β tan φ + ( u a u w ) tan φ b L + c L ( G s + G t ) sin β + P
in which Fs denotes the factor of safety, Gs denotes the gravity of the soil in the saturated layer (N/m), Gt denotes the gravity of the soil in the transition layer (N/m), φ′ denotes the effective friction angle of the natural soil (°), φb is the angle that denotes the contribution of matrix suction against shear strength (°), c′ denotes the effective cohesion of the natural soil (Pa), ua denotes pore air pressure (Pa), uw denotes pore water pressure (Pa), and P denotes the seepage force (N/m).
According to Vanapalli [31], tanφb can be expressed as follows:
tan φ b = Θ tan φ
in which Θ = (θθr)/(θsθr) denotes relative water content. For the slip parallel to the slope surface along the wetting front, combing Equations (16) and (17), the factor of safety can be calculated by:
F s = tan φ [ ( G s + G t ) cos β + ( u a u w ) Θ i L ] + c L ( G s + G t ) sin β + P
in which Θi = (θiθr)/(θsθr).
For the slip parallel to the slope surface along the interface between the transition layer and saturated layer, the factor of safety can be calculated by:
F s = G s cos β tan φ s + c s L G s sin β + P
in which φs′ represents the effective friction angle of the saturated soil (°), and cs′ represents the effective cohesion of the saturated soil (Pa).
For the slip along the surface of the bedrock, the factor of safety can be calculated by:
F s = ( G s + G t + G n ) cos β tan φ + c L ( G s + G t + G n ) sin β + P
in which Gs denotes the gravity of the soil in the natural layer (N/m).

3.3. Variation in Fs for a Natural Slope during a Rainfall Event

Natural slopes usually have irregularly curved surfaces and slip surfaces, which makes it difficult to calculate the depth of the wetting front for natural slopes during intense rainfall. Without the depth of the wetting front, the stability of the natural slope cannot be evaluated by the limit equilibrium method.
The slice method is adopted to calculate the depth of the wetting front. The sides of the slices are vertical. The curved slope surface and slip surface of each slice are replaced by straight lines. Each slice can be deemed as a finite slope. The simplified model of the natural slope is shown in Figure 5. Then, the MSGA model is used to calculate the depth of the wetting front of each slice.
According to Equation (15), the depth of the wetting front during intense rainfall is influenced by the slice’s dip angle and the slice’s length. Therefore, different slices have different depths of the wetting front. The depth of the wetting front should be continuous. This kind of discontinuity can be reduced by increasing the number of slices.
A slice is divided into three layers based on the MSGA model, which is shown in Figure 6a. The gravity of the ith slice can be expressed as follows:
G i = G i s + G i t + G i n = γ s A i s + G i t + γ i A i n
where Gi represents the gravity of the ith slice. Gis represents the gravity of the saturated layer of the ith slice (N/m), Git represents the gravity of the transitional layer of the ith slice (N/m), and Gin represents the gravity of the natural layer of the ith slice (N/m). γs represents the unit weight of the saturated soils (N/m3), and γi represents the unit weight of the natural soils (N/m3). Ais represents the area of the saturated layer of the ith slice (m2), and Ain represents the area of the natural layer of the ith slice (m2).
According to Yao [20], Git can be expressed as follows:
G i t = { [ L 1 Z f 3 Z f 2 sin β 1 4 3 Z f 2 cos β 1 4 tan ( α 1 β 1 ) ] cos β 1 γ d [ 4 ( 1 + θ i ) + π ( θ s θ i ) 8 ] , i = 1 L i cos β i γ d Z f [ 4 ( 1 + θ i ) + π ( θ s θ i ) 8 ] , 1 < i < n [ L n Z f 3 Z f 2 sin β n 4 3 Z f 2 cos β n 4 tan ( β n α n ) ] cos β n γ d [ 4 ( 1 + θ i ) + π ( θ s θ i ) 8 ] , i = n
where Li denotes the length of the slope surface of the ith slice (m), βi denotes the dip angle of the slope surface of the ith slice (°), and αi denotes the dip angle of the slip surface of the ith slice (°).
In the saturated layer, the seepage parallel to the slope surface would cause the seepage force parallel to the slope surface, which is expressed as follows:
P i = γ w A i s sin β i
where Pi is the seepage force of the ith slice (N/m), and γw is the unit weight of water (N/m3). The force diagram of the ith soil slice is shown in Figure 6b.
The factor of safety can be calculated by:
F s = [ G i cos α i + P i sin ( β i α i ) + E i 1 sin ( α i 1 α i ) ] tan φ + E i + c l i G i sin α i + P i cos ( β i α i ) + E i 1 cos ( α i 1 α i )
where Fs represents the factor of safety, Ei represents the thrust between the ith slice and the (i + 1)th slice (N/m), and li represents the length of the slip surface (m). Equation (24) can be written as follows:
E i = F s G i sin α i ( G i cos a i tan φ + c l i ) + P i [ cos ( β i α i ) sin ( β i α i ) tan φ ] , i = 1 E i = F s G i sin α i + E i 1 [ F s cos ( α i 1 α i ) sin ( α i 1 α i ) tan φ ] ( G i cos a i tan φ + c l i ) + P i [ cos ( β i α i ) sin ( β i α i ) tan φ ] , i 2
To be consistent with reality, it is assumed that there is no tensile stress between adjacent slices. Therefore, the value of Ei is replaced by 0 when the calculation result of Ei is smaller than 0, which can be expressed as follows:
E i = { 0 , E i < 0 E i , E i 0
The process for calculating the stability of the natural slope is illustrated as follows: Firstly, an arbitrarily small value of Fs is input. Then, for i = 1 to i = n, the values of Ei can be calculated in sequence. If En is not equal to 0, then Fs will be reset with a greater value, and the value of Ei will also be recalculated. Fs needs to be iterated until En equals 0. When En equals 0, the corresponding Fs is the actual factor of safety of the slope. The process of the proposed method is shown in Figure 7.

4. Case Study and Results

4.1. Background

The Liangshuijing landslide is selected as an example to assess the stability of a natural slope during a rainfall event. The Liangshuijing landslide is located in the Three Gorges reservoir area of Chongqing, adjacent to Guling Town. The location of the landslide is shown in Figure 8b. The Liangshuijing landslide was first discovered in 2003. In early April 2009, it accelerated [32]. The landslide has been monitored for many years, which provides a rich dataset [33]. The sliding body consists of a Quaternary soil layer (Q4del). The bedrock is comprised of interbedded sandstone and mudstones (J2s).
The height of the landslide front is about 100 m, the height of the rear edge is about 319.5 m, the relative height difference is about 221.5 m, the longitudinal length of the plane is about 434 m, and the transverse width is 358 m. The total area of the Liangshuijing landslide is 11.82 × 104 m2, and the total volume is 407.79 × 104 m3. The average thickness of the sliding body is about 34.5 m. The topographic map of the Liangshuijing landslide is shown in Figure 8a.
The Liangshuijing landslide is located in a subtropical monsoon humid climate area. Rainfall in the Liangshuijing landslide area is mainly distributed from March to August each year, and the total rainfall in these six months can account for more than 65% of the total annual rainfall. The average annual rainfall is 1436.5 mm, with the maximum daily rainfall reaching 189.6 mm.

4.2. Application of the MSGA Model

The 1-1′ profile is chosen as an example to study the process of calculating the safety factor of a natural slope during a rainfall event. Figure 9 shows the 1-1′ profile. According to the boundary and geometry of the slide body, the slide body is divided into eight pieces. From top to bottom, the numbers of the slices are from 1st to 8th.
The physical and mechanical parameters of the soil of the sliding body in the Liangshuijing landslide are shown in Table 1. The geometric elements and infiltration parameters of each slice in the 1-1′ profile are listed in Table 2. The rainfall intensity is set to 0.04 m/h. The MSGA model and SGA model are used to calculate the factor of safety of the slope. The variation in the factor of safety of the slope with time is shown in Figure 10.
The 5th slice has the smallest dip angle of the slope surface among the slices, and it reaches the critical state earliest. The difference between the MSGA model and the SGA model appears at tc5th = 2.39 h.
The slope of the curve of the average depth of the wetting front represents the average moving speed of the wetting front. The infiltration rate of the wetting front calculated by the MSGA model is fast at first, but it becomes slower and slower with time. At last, the depth of the wetting calculated by the MSGA model tends to be a constant. The explanation for the tendency might be that it reaches a state of balance. In the state of balance, the recharge from the rainwater infiltration just compensates for the discharge of the seepage parallel to the slope surface. The advancing speed of the wetting front calculated by the SGA model also decreases constantly. The gap between the speed calculated by the two methods widens with time. The factor of safety shows a negative correlation with the average depth of infiltration. The factors of safety calculated by the two models change from 1.195 at the initial time to 1.141 and 1.037 at 40 h, respectively.

5. Discussion

5.1. The Applicability of the MSGA Model

In order to verify the accuracy of the MSGA model, it is used to calculate the depth of the wetting front under the M-3 test condition proposed by Lin [34]. The results calculated by the MSGA model are compared with the measured data. In addition, the GA model and SGA model are also used to compare with the MSGA model.
According to Ma [35], the results of the Richards equation solved by HYDRUS-1D are not very consistent with the measured results of the experiments on the flat surface. The HYDRUS-1D code severely underestimates the advancing depth of the wetting front. However, Richards’ equation is widely used to simulate the infiltration of the rainwater for the slope. Therefore, it is necessary to study the applicability of Richards’ equation to the slope. The Runge-Kutta 4th-order method is used to solve the Richards equation under the M-3 test condition. The results derived from Richards’ equation are compared with the MSGA model, GA model, and SGA model.
The soil water retention curve of the silt is described by the Van Genuchten model, which is expressed as follows [36]:
{ θ θ r θ s θ r = ( 1 + | α h | n ) m , h > 0 θ = θ s , h 0
in which θr represents the residual water content, and h represents the soil water pressure head (m). a, m, and n represent empirical parameters (m = 1 − 1/n). The relation between the unsaturated hydraulic conductivity (K) and soil water pressure head (h) is expressed as follows [37]:
K ( h ) = K s { 1 ( α h ) m n [ 1 + ( α h ) n ] m } 2 [ 1 + ( α h ) n ] m l
in which l is an empirical parameter equal to 0.5 for most soils.
The M-3 test model conducted by Lin is shown in Figure 11. In the M-3 test, three tensiometers are installed in the model. Tensiometers 1, 2, and 3 are at a depth of 0.4 m, 0.1 m, and 0.167 m, respectively. The slope angle of the model is 33.7°. The rainfall intensity equals 0.04 m/h. Table 3 shows the soil parameters.
The depths of the wetting front calculated by different methods and measured by the tensiometers are shown in Figure 12. The times of the wetting front arriving at a depth of 0.1 m, 0.167 m, and 0.4 m are recorded in Table 4.
The depths of the wetting front calculated by the MSGA model and SGA model are consistent with the measured results. The accuracy of the GA model is worse than that of the MSGA model and SGA model. The difference between the results calculated by the MSGA model and the SGA model widens over time.
At the first 0.7 h, the depth calculated by Richards’ equation is close to the results calculated by different methods and measured by the tensiometer. Then, the depth calculated by Richards’ equation diverges greatly from the measured results. The Richards equation solved by the Runge-Kutta 4th-order method severely underestimates the moving speed of the wetting front, which is in accordance with the results calculated by HYDRUS-1D in the plane [35].

5.2. Effects of the Factors for Infiltration

Besides the properties of the soil, the rainwater infiltration process based on the MSGA model is affected by the length of the slope surface and the dip of the slope surface. A series of parametric studies are performed to assess the effect and relative contribution of the controlling factors on infiltration. The parameters of the soil and the rainfall intensity in the parametric studies are set as follows: θs = 0.45, θi = 0.15, Ks = 0.016 m/h, Sf = 0.09 m, and q = 0.04 m/h.
Four lengths of the slope surface, L (1, 3, 5, and 10 m with an inclination of 45° from the horizon), are used for comparison. The result is shown in Figure 13. The depth of infiltration is positively correlated with the length of the slope surface. The size effect becomes apparent when L is short. When L is large enough, the MSGA model degenerates into the SGA model.
Three dip angles of the slope surface, β (15, 30, and 45° with a length of the slope surface of 2 m), are used for comparison. The result is shown in Figure 14. As mentioned above, for the MSGA model, the moving speed of the wetting front is controlled by the recharge of the rainfall infiltration and the discharge of the seepage parallel to the slope surface.
At the initial stage, with the same dip of the slope surface, the depths of the wetting front calculated by MSGA and SGA are similar. The explanation for this might be that the recharge of the rainfall infiltration dominates the infiltration rate of the wetting front. The effect of the seepage parallel to the slope surface, which is neglected in the SGA model, is small at the initial stage. Therefore, the difference between the two methods is very small. With the same dip angle of the slope surface, the depth of the wetting front calculated by the MSGA model is less than the result calculated by the SGA model. It indicates that seepage parallel to the slope inhibits the development of rainfall infiltration depth. The infiltration depth increases with the increase in slope inclination, which is consistent with the conclusion of Chen and Young [18].
At the latter stage, the depth of the wetting front calculated by the MSGA model decreases with the increase in slope inclination, because the seepage parallel to the slope surface dominates the moving speed of the wetting front at the latter stage. The volume of the seepage parallel to the slope surface increases with the increasing depth of the wetting front, according to Equation (15). Therefore, from the initial stage to the latter stage, the effect of the seepage parallel to the slope surface increases constantly.

5.3. Stability for Different Sliding Surfaces

To compare the factor of safety for different sliding surfaces, a series of studies are performed on a finite slope, which is shown in Figure 15. The parameters of the soil are set as yd = 13.45 kN/m3, θs = 0.45, θi = 0.15, θr = 0.10, c′ = 19 kPa, cs′ = 12 kPa, φ′ = 19.8°, φs′ = 13.1°, ua-uw = 5 kPa. Two dip angles of the slope surface, α (15° and 25°), are used for comparison. The length of the slope surface is 5 meters. The method to calculate the factor of safety for various sliding surfaces is illustrated in Section 3.2: Stability analysis for various sliding surfaces.
The variation in factors of safety based on different sliding surfaces is shown in Figure 15. Compared with the factors of safety for sliding along the interface and the wetting front, the factor of safety for sliding along the bedrock surface decreases slowly with the depth of the wetting front. The factors of safety for sliding along the bedrock surface and the wetting front are equal when the wetting front reaches the bedrock surface. At the initial stage, the bedrock surface is the most dangerous sliding surface. With the advancing of the wetting front, the interface becomes the most dangerous sliding surface. When α = 15°, the factor of safety for sliding along the interface drops below 1 firstly. When α = 25°, the factor of safety for sliding along the bedrock surface drops below 1 firstly. The sliding surface of a slope is largely determined by the dip angle of α. The results are shown in Figure 16.

6. Conclusions

This paper proposes the MSGA model to simulate rainfall infiltration. The MSGA model considers the existence of the transition layer and the seepage parallel to the slope surface. This paper uses the Liangshuijing landslide as an example to demonstrate how to calculate the factor of safety of a natural slope. Also, this paper validates the accuracy of the MSGA model and discusses the applicability of Richards’ equation. Finally, this paper discusses the effects of the factors for infiltration. According to the analyses conducted, the conclusions are obtained as follows:
(1) The factor of safety of the Liangshuijing landslide decreases constantly with time during an intense rainfall event. The factor of safety presents a negative correlation with the depth of the wetting front. The factor of safety calculated by the SGA model is smaller than that calculated by the MSGA model, and the gap between the results widens with time.
(2) The MSGA model can more accurately simulate the process of rainfall infiltration compared with the SGA model and Richards’ equation. Richards’ equation solved by the Runge-Kutta 4th-order method evidently underestimates the advancing speed of the wetting front.
(3) The length of the slope surface presents a positive correlation with the advancing speed of the wetting front, and the size effect becomes apparent when the length is short. The recharge of the rainwater infiltration and seepage parallel to the slope surface controls the moving speed of the wetting front. At the initial stage, the effect of the seepage parallel to the slope surface is small. As the infiltration depth increases, the volume of the seepage parallel to the slope surface increases. At the latter stage, the seepage parallel to the slope surface dominates the moving speed of the wetting front.
In this paper, it is assumed that the soil is homogeneous. The spatial variability and the heterogeneity of the soil parameters are not considered. In addition, the 2D plane strain model is used in this article and most of the previous research. However, few studies have focused on the application of the GA model to 3D slope stability analysis due to the complexity of the calculation, while actual slopes are usually 3D. This may affect the reliability of the result. Therefore, further studies are required to be undertaken in the future to solve the problem mentioned above.

Author Contributions

Methodology, X.W. and W.X.; Software, X.W. and W.X.; Validation, X.W.; Writing—original draft, X.W.; Writing—review & editing, W.R., S.C. and L.L.; Funding acquisition, W.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (Grant No. 2021YFC1523400).

Data Availability Statement

All data analyzed during this study are included in this article.

Conflicts of Interest

The authors have no relevant financial or non-financial interests to disclose.

Code Availability

The custom code used during the current study is available from the corresponding author.

References

  1. Chinkulkijniwat, A.; Yubonchit, S.; Horpibulsuk, S.; Jothityangkoon, C.; Jeeptaku, C.; Arulrajah, A. Hydrological responses and stability analysis of shallow slopes with cohesionless soil subjected to continuous rainfall. Can. Geotech. J. 2016, 53, 2001–2013. [Google Scholar] [CrossRef]
  2. Dou, H.Q.; Han, T.C.; Gong, X.N.; Zhang, J. Probabilistic slope stability analysis considering the variability of hydraulic conductivity under rainfall infiltration-redistribution conditions. Eng. Geol. 2014, 183, 1–13. [Google Scholar] [CrossRef]
  3. Wang, G.H.; Sassa, K. Pore-pressure generation and movement of rainfall-induced landslides: Effects of grain size and fine-particle content. Eng. Geol. 2003, 69, 109–125. [Google Scholar] [CrossRef]
  4. Rahardjo, H.; Ong, T.H.; Rezaur, R.B.; Leong, E.C. Factors controlling instability of homogeneous soil slopes under rainfall. J. Geotech. Geoenviron. 2007, 133, 1532–1543. [Google Scholar] [CrossRef]
  5. Tozato, K.; Dolojan, N.L.J.; Touge, Y.; Kure, S.; Moriguchi, S.; Kawagoe, S.; Kazama, S.; Terada, K. Limit equilibrium method-based 3D slope stability analysis for wide area considering influence of rainfall. Eng. Geol. 2022, 308, 106808. [Google Scholar] [CrossRef]
  6. Xu, J.S.; Du, X.L.; Zhao, X.; Li, L.Y. Analytical Stability Analysis of Rainfall-Infiltrated Slopes Based on the Green-Ampt Model. Int. J. Geomech. 2023, 23, 04022277. [Google Scholar] [CrossRef]
  7. Green, W.H.; Ampt, G.A. Studies on soil physics Part I—The flow of air and water through soils. J. Agr. Sci. 1911, 4, 1–24. [Google Scholar] [CrossRef]
  8. Richards, L.A. Capillary conduction of liquids through porous mediums. Phys.-J. Gen. Appl. P 1931, 1, 318–333. [Google Scholar] [CrossRef]
  9. Barontini, S.; Ranzi, R.; Bacchi, B. Water dynamics in a gradually nonhomogeneous soil described by the linearized Richards equation. Water Resour. Res. 2007, 43, W08411. [Google Scholar] [CrossRef]
  10. Elmaloglou, S.; Diamantopoulos, E. The effect of intermittent water application by surface point sources on the soil moisture dynamics and on deep percolation under the root zone. Comput. Electron. Agr. 2008, 62, 266–275. [Google Scholar] [CrossRef]
  11. Rahimi, A.; Rahardjo, H.; Leong, E.C. Effect of hydraulic properties of soil on rainfall-induced slope failure. Eng. Geol. 2010, 114, 135–143. [Google Scholar] [CrossRef]
  12. Montrasio, L.; Valentino, R. A model for triggering mechanisms of shallow landslides. Nat. Hazard. Earth Sys. 2008, 8, 1149–1159. [Google Scholar] [CrossRef]
  13. Ogden, F.L.; Saghafian, B. Green and Ampt infiltration with redistribution. J. Irrig. Drain. Eng. 1997, 123, 386–393. [Google Scholar] [CrossRef]
  14. Dolojan, N.L.J.; Moriguchi, S.; Hashimoto, M.; Terada, K. Mapping method of rainfall-induced landslide hazards by infiltration and slope stability analysis A case study in Marumori, Miyagi, Japan, during the October 2019 Typhoon Hagibis. Landslides 2021, 18, 2039–2057. [Google Scholar] [CrossRef]
  15. Li, C.J.; Guo, C.X.; Yang, X.G.; Li, H.B.; Zhou, J.W. A GIS-based probabilistic analysis model for rainfall-induced shallow landslides in mountainous areas. Environ. Earth Sci. 2022, 81, 432. [Google Scholar] [CrossRef]
  16. Loáiciga, H.A.; Huang, A. Ponding analysis with Green-and-Ampt infiltration. J. Hydrol. Eng. 2007, 12, 109–112. [Google Scholar] [CrossRef]
  17. Wei, X.Y.; Liu, Y.B.; Ren, W.Z.; Cai, S.M.; Xu, W.H. Fully analytical method for probabilistic stability analysis of the rainfall-induced landslide considering variability of saturated hydraulic conductivity and rainfall uncertainty. Bull. Eng. Geol. Environ. 2024, 83, 374. [Google Scholar] [CrossRef]
  18. Chen, L.; Young, M.H. Green-Ampt infiltration model for sloping surfaces. Water Resour. Res. 2006, 42, W07420. [Google Scholar] [CrossRef]
  19. Wang, D.J.; Tang, H.M.; Zhang, Y.H.; Li, C.D.; Huang, L. An improved approach for evaluating the time-dependent stability of colluvial landslides during intense rainfall. Environ. Earth Sci. 2017, 76, 321. [Google Scholar] [CrossRef]
  20. Yao, W.M.; Li, C.D.; Zhan, H.B.; Zeng, J.B. Time-dependent slope stability during intense rainfall with stratified soil water content. Bull. Eng. Geol. Env. 2019, 78, 4805–4819. [Google Scholar] [CrossRef]
  21. Jiang, S.H.; Liu, X.; Ma, G.T.; Rezania, M. Stability analysis of heterogeneous infinite slopes under rainfall-infiltration by means of an improved Green-Ampt model. Can. Geotech. J. 2024, 61, 1560–1573. [Google Scholar] [CrossRef]
  22. Chu, S.T. Infiltration during an Unsteady Rain. Water Resour. Res. 1978, 14, 461–466. [Google Scholar] [CrossRef]
  23. Xiong, S.; Yao, W.M.; Li, C.D. Stability evaluation of multilayer slopes considering runoff in the saturated zone under rainfall. Eur. J. Environ. Civ. Eng. 2019, 25, 1718–1732. [Google Scholar] [CrossRef]
  24. Colman, E.A.; Bodman, G.B. Moisture and energy conditions during downward entry of water into moist and layered Soils. Soil. Sci. Soc. Am. J. 1945, 9, 3–11. [Google Scholar] [CrossRef]
  25. Wang, W.Y.; Wang, Z.R.; Wang, Q.J.; Zhang, J.F. Improvement and evaluation of the Green-Ampt model in loess soil. J. Hydraul. Eng. 2003, 34, 30–34. (In Chinese) [Google Scholar]
  26. Zhang, J.; Han, T.; Dou, H.; Ma, S. Stability of loess slope considering infiltration zonation. J. Cent. South Univ. (Sci. Technol.) 2014, 45, 4355–4361. (In Chinese) [Google Scholar]
  27. Peng, Z.; Huang, J.; Wu, J.; Guo, H. Modification of Green-Ampt model based on the stratification hypothesis. Adv. Water Sci. 2012, 23, 59–66. (In Chinese) [Google Scholar]
  28. Cheng, P.; Liu, Y.; Li, Y.P.; Yi, J.T. A large deformation finite element analysis of uplift behaviour for helical anchor in spatially variable clay. Comput. Geotech. 2022, 141, 104542. [Google Scholar] [CrossRef]
  29. Chen, X.J.; Fang, P.P.; Chen, Q.N.; Hu, J.; Yao, K.; Liu, Y. Influence of cutterhead opening ratio on soil arching effect and face stability during tunnelling through non-uniform soils. Undergr. Space 2024, 17, 45–59. [Google Scholar] [CrossRef]
  30. Cho, S.E. Infiltration analysis to evaluate the surficial stability of two-layered slopes considering rainfall characteristics. Eng. Geol. 2009, 105, 32–43. [Google Scholar] [CrossRef]
  31. Vanapalli, S.K.; Fredlund, D.G.; Pufahl, D.E.; Clifton, A.W. Model for the prediction of shear strength with respect to soil suction. Can. Geotech. J. 1996, 33, 379–392. [Google Scholar] [CrossRef]
  32. Wang, H.L.; Xu, W.Y. Stability of Liangshuijing landslide under variation water levels of Three Gorges Reservoir. Eur. J. Environ. Civ. Eng. 2013, 17, S158–S173. [Google Scholar] [CrossRef]
  33. Jiang, P.; Chen, J.J. Displacement prediction of landslide based on generalized regression neural networks with K-fold cross-validation. Neurocomputing 2016, 198, 40–47. [Google Scholar] [CrossRef]
  34. Lin, H.-C. The Study on Mechanism and Numerical Analysis of Rainfall-Induced Soil Slope Failure; Tsinghua University: Beijing, China, 2007. (In Chinese) [Google Scholar]
  35. Ma, Y.; Feng, S.Y.; Su, D.Y.; Gao, G.Y.; Huo, Z.L. Modeling water infiltration in a large layered soil column with a modified Green-Ampt model and HYDRUS-1D. Comput. Electron. Agr. 2010, 71, S40–S47. [Google Scholar] [CrossRef]
  36. Vangenuchten, M.T. A Closed-Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils. Soil. Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef]
  37. Vangenuchten, M.T.; Leij, F.J.; Yates, S.R. RETC Code for Quantifying the Hydraulic Functions of Unsaturated Soils; U.S Salinity Laboratory, U.S. Department of Agriculture, Agricultural Research Service: Riverside, CA, USA, 1991.
Figure 1. Diagram of the GA model. (a) Original model. (b) Modified model.
Figure 1. Diagram of the GA model. (a) Original model. (b) Modified model.
Water 16 02877 g001
Figure 2. Distribution of the water content based on the SGA model.
Figure 2. Distribution of the water content based on the SGA model.
Water 16 02877 g002
Figure 3. Definition of coordinates.
Figure 3. Definition of coordinates.
Water 16 02877 g003
Figure 4. Diagram of the MSGA model.
Figure 4. Diagram of the MSGA model.
Water 16 02877 g004
Figure 5. Diagram of the sliced division of the natural slope.
Figure 5. Diagram of the sliced division of the natural slope.
Water 16 02877 g005
Figure 6. The ith slice. (a) Geometric profile. (b) Forces diagram.
Figure 6. The ith slice. (a) Geometric profile. (b) Forces diagram.
Water 16 02877 g006
Figure 7. Flow chart for the proposed method.
Figure 7. Flow chart for the proposed method.
Water 16 02877 g007
Figure 8. (a) Topographic map of the Liangshuijing landslide. (b) Location of the study area.
Figure 8. (a) Topographic map of the Liangshuijing landslide. (b) Location of the study area.
Water 16 02877 g008
Figure 9. Cross section along 1-1′ profile of Liangshuijing landslide.
Figure 9. Cross section along 1-1′ profile of Liangshuijing landslide.
Water 16 02877 g009
Figure 10. Variation in the factor of safety and the average depth of the wetting front during a rainfall event.
Figure 10. Variation in the factor of safety and the average depth of the wetting front during a rainfall event.
Water 16 02877 g010
Figure 11. Sketch of the M-3 test model.
Figure 11. Sketch of the M-3 test model.
Water 16 02877 g011
Figure 12. The depths of the wetting front were calculated by different methods and measured by tensiometer.
Figure 12. The depths of the wetting front were calculated by different methods and measured by tensiometer.
Water 16 02877 g012
Figure 13. Effect of the length of the slope surface on infiltration.
Figure 13. Effect of the length of the slope surface on infiltration.
Water 16 02877 g013
Figure 14. Effect of the dip angle of the slope surface on infiltration.
Figure 14. Effect of the dip angle of the slope surface on infiltration.
Water 16 02877 g014
Figure 15. The model for comparison of safety for different sliding surfaces.
Figure 15. The model for comparison of safety for different sliding surfaces.
Water 16 02877 g015
Figure 16. Variation in factors of safety based on different sliding surfaces (a) α = 15°, (b) α = 25°.
Figure 16. Variation in factors of safety based on different sliding surfaces (a) α = 15°, (b) α = 25°.
Water 16 02877 g016
Table 1. Soil parameters of the Liangshuijing landslide.
Table 1. Soil parameters of the Liangshuijing landslide.
ParameterValue
Saturated hydraulic conductivity Ks0.011 m/h
Saturated water content θs0.330
Natural water content θi0.285
Suction at the wetting front Sf6 m
Natural unit weight of the soil γi23 kN/m3
Saturated unit weight of the soil γs23.8 kN/m3
Effective cohesion c19.17 kPa
Effective friction angle φ25.03°
Table 2. Geometric elements and infiltration parameters of each slice in 1-1′ profile.
Table 2. Geometric elements and infiltration parameters of each slice in 1-1′ profile.
Number of the Sliceα (°)β (°)L (m)Zfc (m)tc (h)
1st3520612.632.64
2nd4032433.233.24
3rd3431773.163.17
4th2615272.492.50
5th239312.382.39
6th1822562.712.71
7th1426812.872.89
8th1038503.743.76
Table 3. Parameters of silt.
Table 3. Parameters of silt.
ks (m/h)θsθrθiSf (m)a (1/m)nm
0.0160.4050.10.10.092.0201.5870.3699
Table 4. Comparison of the arrival time.
Table 4. Comparison of the arrival time.
MethodArriving Time at a Certain Depth (h)
0.1 m0.167 m0.4 m
Tensiometer0.681.264.06
GA0.771.404.41
SGA0.681.263.94
MSGA0.681.264.04
Richards’ equation0.801.957.02
Disclaimer/Publisher’s Note: 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

MDPI and ACS Style

Wei, X.; Ren, W.; Xu, W.; Cai, S.; Li, L. A Modified Method for Evaluating the Stability of the Finite Slope during Intense Rainfall. Water 2024, 16, 2877. https://doi.org/10.3390/w16202877

AMA Style

Wei X, Ren W, Xu W, Cai S, Li L. A Modified Method for Evaluating the Stability of the Finite Slope during Intense Rainfall. Water. 2024; 16(20):2877. https://doi.org/10.3390/w16202877

Chicago/Turabian Style

Wei, Xiaoyang, Weizhong Ren, Wenhui Xu, Simin Cai, and Longwei Li. 2024. "A Modified Method for Evaluating the Stability of the Finite Slope during Intense Rainfall" Water 16, no. 20: 2877. https://doi.org/10.3390/w16202877

APA Style

Wei, X., Ren, W., Xu, W., Cai, S., & Li, L. (2024). A Modified Method for Evaluating the Stability of the Finite Slope during Intense Rainfall. Water, 16(20), 2877. https://doi.org/10.3390/w16202877

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop