#### Abstract

Taking the high and steep slope of a ski resort as the background, in order to ensure the structural safety of the slope in the subsequent complex rainfall environment, the calculation principle of unsaturated soil slope stability is summarized, and the parameter interaction analysis is proposed. The numerical simulation method is used to study the influence of various rainfall types on slope stability; combined with a multiparameter combination, the slope stability is analyzed. The results show that different types of rainfall have a significant impact on the safety factor of the slope, the safety factor of the slope under the front rainfall type is the smallest, the minimum value is 1.38, and the front rainfall mostly occurs in the late summer, so it is necessary to avoid construction in the late summer or do a good job of safety monitoring in the late summer. Through setting monitoring points, it is found that the pore water pressure of the three positions from small to large is in the order of slope toe<slope middle<slope top, and the pore water pressure at the slope toe is basically 0 kPa. At this time, the matric suction at the slope toe is the minimum, so the protection should be strengthened. In the combination analysis of the four parameters, the direct effect and coupling effect of C15 are at a small level, which indicates that when the four parameters fluctuate together, the mutual coupling is not the effect of a single parameter, which provides a theoretical reference for the follow-up related research.

#### 1. Introduction

With the development of the national economy, a large number of railways, highways, water conservancies, mines, towns, and other facilities are built, especially in the construction of hills and mountainous areas. In human engineering activities, more and more slopes are excavated and filled, and the slope height is increasing [1–4]. Due to the complex geological conditions in China and with the increasing scale of human transformation of nature, engineering disasters caused by improper design and construction are frequent. In addition to the inherent factors such as geology and geomorphology, the complex natural environment such as rainfall, freeze-thaw cycle, and earthquake will have an inestimable impact on the slope stability.

Domestic and foreign scholars have studied slope stability in complex environments by means of statistical analysis, field tests, theoretical analysis, and numerical simulation and achieved a series of research results [5–8].

Tang et al. [9] studied parameters of the soil-water characteristic curve with a homogeneous slope as an example and on slope stability. Lin et al. [10] obtained the conclusion that the low rainfall intensity and long duration rainfall will increase the pore water pressure of deep soil, resulting in sliding failure and large-scale landslide. Yang et al. [11, 12] studied the influence order of heavy rain, rainstorm, and torrential rain on the pore water pressure change of the slope with a weak interlayer and the variation law of the safety factor and plastic zone under different rainfall conditions. Zhou et al. [13] showed that the safety factor of the slope decreased continuously in the process of rainfall infiltration, but after the rain stopped, the safety factor of the slope decreased gradually. Based on the saturated-unsaturated seepage theory and strength reduction method, Xu et al. [14] analyzed the soil slope, which showed that with the increase of rainfall intensity, the change range of the safety factor was larger. Wang et al. [15] combined with the saturated-unsaturated seepage theory studied the change law of the slope seepage field under different rainfall types. In the above studies, the mechanical parameters of soil material are set as fixed values, and the spatial variability of soil material is not considered.

Luo et al. [16] took 1022 landslides in the National Geopark induced by the Jiuzhaigou earthquake in 2017 as sample data and cross-validated the spatial prediction accuracy of the initial model by using the Bayesian probability method and generalized additive statistical model. Ping et al. [17] took a natural loess slope in Jingyang as the research object, carried out artificially simulated rainfall experiments of the large loess slope under different rainfall intensities, and summarized the water infiltration law and deformation failure mode of the loess slope. Liang et al. [18, 19] used 3D printing technology to make plant rhizomes, studied the influence of slope height on vegetation to improve slope stability under earthquake action, and found that vegetation is a very effective method to reduce earthquake sliding for slopes with smaller height (such as low height embankment along traffic infrastructure). Chen et al. [20] used the integral limit equilibrium variational method to give a detailed theoretical derivation and solution process, obtained the expressions of critical slip surface function and stress function, and gave the solution method of the safety factor. Chen [21] vigorously developed the limit analysis method and successfully applied it to the stability analysis of saturated homogeneous soil. Based on the limit equilibrium method, Han et al. [22] proposed a three-dimensional loading stability analysis method for the anisotropic heterogeneous slope, which can be regarded as the improvement and expansion of the solution, and studied the influence of the friction angle, anisotropic factor, heterogeneous factor, slope angle, width depth ratio, top load, and other factors on slope stability. Cai et al. [23], based on the principle of saturated-unsaturated seepage, taking the low highway slope in Yunnan Province as the research object, calculated and analyzed the slope stability and reliability under different rainfall intensities by using GeoStudio finite element software.

It is not difficult to find that at present, the research on slope stability at home and abroad is mostly on common slopes [24–27], but the study on the stability of high and steep slopes is still less, especially in the complex rainfall environment. In view of this, this paper takes a ski resort as the engineering background and uses GeoStudio finite element analysis software to study the influence of various rainfall types (forward type, center front type, back front type, and average type) on the stability of high and steep slopes. Meanwhile, the best coupling combination is proposed by using the method of parameter combination analysis, and the internal mechanism of the influence of rainfall infiltration on slope stability is analyzed, providing guidance and guarantee for subsequent construction and operation.

#### 2. Calculation Principle of Slope Stability

##### 2.1. Unsaturated Fluid-Structure Interaction and Variance Theory

Rainfall infiltration is a typical unsaturated fluid-solid coupling process. The incompleteness of soil makes the infiltration water flow along the pores and produces seepage water pressure, which forms the strain field of soil structure in the form of seepage volume force and then forms the changes of displacement and volume strain. Rainfall infiltration affects the permeability coefficient and seepage field of soil. It can be seen from the basic principle of hydraulics [28] that the seepage volume force of the continuous porous rock and soil medium is in a positive proportion to the hydraulic gradient, which can be expressed as where is a constant determined by soil structure and fluid properties; is the total rainfall head; is the seepage volume force of rainfall infiltration; , , and are the components of in , , and directions, respectively; , , and are the angle between the seepage volume force and the component value; and , , and are the hydraulic gradients of volume elements in , , and directions, respectively.

In the numerical analysis of geotechnical engineering, when the finite element method is used to solve the stress field, the seepage volume force of the element can be transformed into the equivalent nodal load through formulas (2) and (3) below: where is the infiltration volume; is the shape function matrix of the slope element; and and are the equivalent nodal forces produced by the seepage volume force and increment value, respectively.

In the seepage stress coupling analysis, the empirical relationship between the porosity and the permeability coefficient is adopted. where and are the initial state of slope soil and the porosity after rainfall infiltration, respectively; and and are the permeability coefficient before and after infiltration.

Seepage volume force and other external loads cause the change of the soil stress field and consolidation deformation and the change of the void ratio and porosity, then affect the permeability coefficient, and finally cause the change of the seepage field [29].

The two-dimensional steady seepage field under the action of the stress field can be characterized as where is the element seepage coefficient matrix formed by nodes; is the total head; is a known function of ; and are the normal directions of and , respectively; is the water head boundary; is the flow boundary; is the boundary of infiltration basin; and is the unit seepage.

The relationship between the permeability coefficient, saturation, and matric suction can be expressed as where is the current permeability coefficient; is the saturated permeability coefficient of soil; and are the air pressure and water pressure in the slope soil, where and ; and , , and are the material coefficients. where and are the saturation and residual saturation, respectively, where ; and is the maximum saturation, where .

Under the action of the seepage field, the function characterization model of the stress field is as follows: where is the surface force distribution on the boundary; is the slope stiffness matrix; is the node load increment caused by the external load; is the node load increment caused by seepage volume force; is the displacement change; and are the node displacement and initial node displacement; and and are the displacement and stress boundary, respectively.

It can be seen from equations (7) and (8) that the matrix function expression of seepage stress coupling is as follows: where is the water distribution of the rainfall infiltration field.

##### 2.2. Stability Theory of the Unsaturated Soil Slope

Landslide refers to the phenomenon that the slope loses its original stable state and moves down the slope along a sliding surface as a whole. Landslide is mainly caused by the shear stress on one side of the soil (landslide surface) exceeding its shear strength. The slope soil under the complex rainfall environment is mostly unsaturated soil. For unsaturated soil, the main reason for slope instability is the decrease or loss of soil matric suction caused by rainfall. Therefore, the research on the mechanism between unsaturated soil shear strength and matric suction is of great significance.

The shear strength formula of saturated soil was first proposed by Coulomb: where is the shear strength; is the cohesive force; is the normal stress; and is the internal friction angle.

On the basis of the Mohr-Coulomb shear strength formula, the concept of matric suction is introduced and different shear strength formulas of unsaturated soil are established.

The shear strength formula of unsaturated soil proposed by Bishop [30] in 1960 can be expressed as follows: where is the effective cohesion; is the pore gas pressure; is the pore water pressure; is a parameter related to the type and saturation of soil; and is the effective internal friction angle. Although the value of is related to the type and saturation of soil, its physical meaning is not clear and cannot be obtained simply by experiment or theory, so it is difficult to popularize it in practical engineering.

In 1978, Fredlund et al. [31] put forward the shear strength formula of unsaturated soil with normal stress and suction as variables. where is the rate at which shear strength increases with suction .

Comparing equation (11) with equation (12), it can be found that equation (11) simply regards matric suction as the increase of effective stress, while equation (12) considers the influence of matric suction on the shear strength of unsaturated soil separately, which is more accurate than equation (12). However, and cannot be measured by a simple test, and they are difficult to be applied to practical engineering.

Vanapalli et al. [32, 33] put forward the empirical model of shear strength of unsaturated soil, innovatively connecting with the change of water content in soil, in which the relationship between and water content in soil can be expressed as follows: where is the volume moisture content; is the residual volume moisture content; and is the saturated volumetric moisture content, and the further shear strength of unsaturated soil can be expressed as

The shear strength formula of unsaturated soil is established by formula (14), and can be determined by the characteristic curve of soil and water, which can be widely used in practical engineering. However, the study of slope stability involves not only the design and construction of the existing subgrade slope but also an important research direction of the stability of the natural slope. How to prevent the landslide and how to warn early the quantity of landslide earthwork are inevitable problems in actual engineering. At present, the early warning of the slope is mainly through two methods, namely, field detection and finite element model prediction.

Based on the study of the soil-water characteristic curve, van Genuchten obtained the relationship between unsaturated soil-water content and matric suction in the form of power function (V-G model): where is the matric suction; and and are the fitting parameters. By substituting the V-G model into equation (14), the following results are obtained:

The GeoStudio finite element analysis software used in this paper is based on the V-G model described in equation (16).

##### 2.3. Analysis Principle of Influencing Factors

The soil saturated permeability coefficient, rainfall intensity, Mohr-Coulomb parameters of material properties, and friction angle are the key internal and external factors to determine the stability of the soil slope. In order to analyze the influence of the above four parameters on slope stability, it is necessary to analyze their interaction. Finally, the quantization degree of the output parameters under the input parameters and the interaction between the parameters are realized by using the effect of the variance of the input parameters on the variance of the output results [34]. Taking the pore pressure, saturation, safety factor, and effective stress of the slope model under rainfall infiltration as output parameters and taking the saturated permeability coefficient of soil, rainfall intensity, Mohr-Coulomb parameters of material properties, and friction angle as input parameters, the global quantitative analysis is carried out. The parameter is defined as the input parameter of the model , and the relationship between the input parameter and the output parameter can be defined as where is the influence degree of the change of parameter on the output value , that is, the sensitivity of the output value to the change of the parameter, where is the output parameter (saturated permeability coefficient of soil, rainfall intensity, Mohr-Coulomb parameter of material properties, and friction angle); is the output parameter (pore pressure, saturation, safety factor, and effective stress); and is the difference before and after the change of the same input parameter.

The model of input parameters is established, , and the input parameters are valued at levels. Multiple sampling is used to replace repeated sampling, and the mean value of output parameters determined by multiple sampling calculation is used to represent the direct effect of input parameters on output, and the standard deviation of statistics is used to represent the interaction between input parameters.

#### 3. Engineering Examples

##### 3.1. Project Overview

This paper takes the high and steep slope of a ski resort as an example. The slope height is 130 m, and the longitudinal length is 350 m. The slope is divided into two steps, the first step is the starting area, and the slope is 0.5. The slope of the second step is the auxiliary slide with a slope of 0.26. According to the slope classification, the slope in the starting area is high and steep. The water level elevation on the left side of the slope is 50 m, and that on the right side is 40 m. The schematic diagram of the slope is shown in Figure 1. Taking the high and steep slope in the starting area as the research object, this paper studies the influence of various rainfall conditions on the stability of high and steep slopes.

##### 3.2. Numerical Simulation

###### 3.2.1. Model Parameters

According to the geological exploration report, the soil layer of the first step slope is mainly completely weathered siltstone, and its corresponding physical and mechanical parameters are shown in Table 1. The soil-water characteristic curve adopts the V-G model: , , , , and , and the soil-water characteristic curve and function permeability curve are shown in Figure 2.

###### 3.2.2. Rainfall Parameters

The project site belongs to the warm temperate semihumid monsoon climate zone, with significant continental monsoon climate characteristics and four distinct seasons. It is cold and dry in winter, with northwest wind prevailing. In summer, the rainfall is concentrated and southeast wind prevails. It is dry and windy in spring and cool and rainless in autumn. In addition, due to the influence of topography, precipitation is unevenly distributed not only in four seasons but also in different places. The annual difference between the rainy area and the less rainy area is 100-200 mm. Spring precipitation accounts for 5%-15% of annual precipitation. Summer precipitation accounts for 70%-80% of annual precipitation. The precipitation in autumn accounts for 11%-15% of annual precipitation, slightly more than that in spring. Winter precipitation is very little, accounting for only 1%-4% of annual precipitation. The average annual precipitation of the project site from 1956 to 2008 is 559 mm. According to the rainfall grade specified by the meteorological department (Table 2), combined with the local rainfall, the rainfall parameters are defined as four types of rainfall [35], namely, forward type, center type, rear type, and average type. The duration of rainfall is 7 days. Considering the influence of rainfall infiltration, the working conditions of 7 days after rainfall are also included in the time range. The rainfall process is shown in Figure 3.

**(a) Forward type (referred to as F)**

**(b) Center type (referred to as C)**

**(c) Back forward type (referred to as B)**

**(d) Average type (referred to as A)**

#### 4. Result Analysis

##### 4.1. Natural Slope Stability Analysis

In order to better analyze the influence of different rainfall types on slope stability, the stability of the natural slope is analyzed firstly. The cloud diagram of pore water pressure under the natural state is shown in Figure 4. It can be seen that the saturated soil is below the initial water level, the pore water pressure is positive, and the unsaturated soil is above the initial water level. Under the action of capillary water, the capillary water rises, and the pore water pressure is negative. It can be seen from the equation that when the pore water pressure is negative, the matric suction is greater than zero, and the shear strength of the soil is enhanced. At this time, the initial water level is at the foot of the slope, and because the safety factor of the slope is 1.5, the slope can be considered to be stable.

##### 4.2. Slope Stability Analysis under Various Rainfall Conditions

###### 4.2.1. Safety Factor

The influence of rainfall infiltration on slope stability has been a key issue in geotechnical engineering. The safety factor is a parameter which can reflect the slope stability directly. Figure 5 shows the slope safety coefficient under different rainfall types. It can be seen from Figure 5 that in the rainfall stage, the safety coefficient of the slope decreases with the continuous rainfall and the change range is more intense. After the rain stops on the 7th day, the rainwater will continue to infiltrate, the water level will gradually decrease, and the change range of the safety coefficient is slow and gradually stable. In addition, the law of the change of the safety coefficient under different rainfall types can be seen from the comparison of the rainfall stage: forward type<average type<center front type<back front type, which is because the rainfall of the forward type is much larger than those of the other three types in the initial stage of rainfall, and the soil mass of the slope can reach the saturation state in a short time, which leads to the rapid decrease of matric suction and slope stability, The macroperformance is that the safety coefficient is reduced rapidly. The theory infers that the variation law of the safety coefficient of the average type and center front type is also verified. The safety coefficient of the first three days of the front-type rainfall is always greater than that of the average rainfall condition. However, with the increasing rainfall of the center front type, the gap between the two is smaller and smaller. On the fourth day, the rainfall intensity of the center front type and the average type is 70 mm. At this time, the safety coefficient of the center front type is almost equal to that of the average type. In the next three days, with the decrease of the rainfall intensity of the center front type, the safety coefficient of the front-type rainfall condition gradually exceeds that of the average one. The whole change of the safety factor can be divided into two stages, the first stage is the sharp decline stage, and the four rainfall modes show the same law. With the rainfall no longer changing, the safety coefficient is gradually stable. Although the total rainfall is consistent, the final safety factor decreases in the same trend, but the change of the safety factor caused by different rainfall types in the process of rainfall is different, which is also the significance of this study.

In general, the safety factors of the four rainfall types from small to large are F<A<B<C, and the minimum safety factor is 1.38. In engineering practice, the front-type rainfall mostly occurs at the end of summer, so the construction unit should avoid the construction at the end of summer, and the slope engineering in use should be well monitored at the end of summer.

###### 4.2.2. Pore Water Pressure

The safety factor is a quantitative description of slope stability. According to the calculation principle of unsaturated soil slope stability described in Section 2, matric suction has an important influence on the shear strength of unsaturated soil. Under the influence of rainfall, the matric suction is mainly related to the pore water pressure. In order to better study the influence of rainfall on slope stability, monitoring points are arranged at the top, middle, and foot of the slope to explore the changes of pore water pressure in the slope under different rainfall types. The results are shown in Figure 6.

**(a) Pore pressure at the top of the slope**

**(b) Pore pressure in the middle of the slope**

**(c) Pore pressure at the toe of the slope**

According to Figure 6, the pore water pressure of each monitoring point in the slope increases with the increase of rainfall. After the rain, the pore water pressure decreases with the rain infiltration in the soil and finally tends to be stable. The pore water pressure at the top of the slope reaches a stable value on the fourth day, and the pressure of the slope all reaches a stable value on the sixth day. It can also reflect the process of rainwater seepage and the reaction time of pore water pressure. In addition, it can be seen from Figure 6(c) that under three rainfall conditions, the initial pore water pressure shows a trend of decreasing first and then rising. This is caused by the rainfall intensity being less than the permeability intensity. The permeability coefficient is assumed to be 20 mm/h (480 mm/d), but the total rainfall reduction in the initial two days of C, B, and A is less than the daily infiltration rainfall of the slope, while the pore water pressure starts to rise on the second day. It can be seen that the rainfall intensity and permeability coefficient are not simple linear addition and subtraction to slope stability but a response aging problem. Figure 7 shows the change chart of pore water pressure at each point of the slope under the forward rainfall condition. The change law is consistent with Figure 6. The pore water pressure increases first and then decreases, and the increase and decrease tend to be stable. The pore water pressure in three positions is arranged from small to large: the toe of the slope<the middle of the slope<the top of the slope, but it is worth noting that the pore water pressure at the toe of the slope is basically maintained at about 0 kPa. At this time, the matric suction value is the smallest, and the protection measures for the toe should be strengthened. Figure 8 is the pore water pressure cloud chart of the slope under the forward rainfall condition. The representative 2, 5, 8, and 14 days are selected for analysis in the cloud chart of 0-14 days. It can be seen from Figure 8 that the pore water pressure near the slope surface rises rapidly when rainfall occurs, and rainwater penetrates down the slope, and the growth rate of pore water pressure from top to bottom gradually slows down. With the continuous slope saturated with rainfall time, on the 14th day, pore water pressure in the main area of the slope was 0 kPa, the matric suction value was the smallest, and the overall stability of the slope was poor.

**(a) Lasted two days**

**(b) Lasted five days**

**(c) Lasted eight days**

**(d) Lasted fourteen days**

##### 4.3. Stability Analysis of the Slope with Parameter Interaction

According to the initial values of the saturated permeability coefficient, rainfall intensity, Mohr-Coulomb parameter of material properties, and friction angle, the unsaturated parameter is 0.018, rainfall intensity is 20 mm/h, cohesion is 15 kPa, and internal friction angle is 30°. The setting ranges of the four parameters are ±10%, ±10%, ±10%, and ±5%. Based on the Morris quantitative analysis method and statistical principle, there are 15 combinations of the four parameters, as shown in Table 3. Using different parameter combinations and distribution interval values, the slope stability model can be calculated, and the action values of the four input parameters on the slope stability characteristics can be determined.

Under the unsaturated seepage condition, the shear strength parameters of soil which determine the external precipitation of fluid-solid coupling are not independent, and there are interactions among parameters, which change pore pressure, saturation, displacement, and plastic strain which characterize the stability of slopes. The results of finite element simulation based on the Morris input parameter combination can determine the change characteristics of four stability parameters under the action of 15 groups of parameters.

It can be seen from Figure 9(a) that the first four groups (1, 2, 3, and 4) of the 15 groups of parameters are the independent effects of the four parameters, and the abscissa is the standard deviation of the parameters, indicating the coupling interaction between the parameters. The ordinate is the mean value of parameter combination, which indicates the direct effect of parameters on pore pressure characteristics of slope stability. The mean value and standard deviation of C2>C1>C3>C4; that is, the rainfall intensity has the most obvious effect on the pore pressure characteristics of slope stability. The standard deviation of C6, C7, C8, and C9 is greater than the mean value and standard deviation of C5; that is, the former parameter combination has a positive effect on the pore pressure characteristics of slope stability, C4 has a negative effect on the pore pressure characteristics, and the interaction coupling between them reduces the impact on the pore pressure characteristics of slope stability.

**(a) Pore pressure**

**(b) Saturation**

**(c) Displacement**

**(d) Equivalent plastic strain**

It can be seen from Figure 9(b) that the mean value and standard deviation of C2>C1>C3>C4; that is, the rainfall intensity has the most obvious influence on the pore pressure characteristics of slope stability, and the coupling interaction between rainfall intensity and other three parameters is weak.

The direct effect of the combination of two parameters is shown. The order of parameter combinations from high to low is C6>C7>C5>C10>C8>C9, and the smallest two groups of C8 and C9 have the same rainfall intensity. Combined with the direct effect of single-parameter analysis, the influence of rainfall intensity on saturation is negative. According to the standard deviation, the coupling effect of parameters from high to low is C8>C9>C5>C10>C3, and the slope stability analysis after reinforcement is C6>C7.

It can be seen from Figure 9(c) that among the direct effects represented by the mean value, C2>C4>C3>C1, and the standard deviation from large to small is C4>C3>C2>C1; that is, the coupling effect between the saturation coefficient and the other three terms is the smallest, and the coupling effect between the internal friction angle and the other three terms has the greatest impact on the displacement characteristics of slope stability.

According to the mean distribution of two-parameter combination, C10>C6>C5>C8>C9>C7. From the mean value of a single parameter, it can be seen that there is a negative coupling effect between the saturated permeability coefficient and the internal friction angle, which leads to the decrease of the direct and coupling effects of C7 on the displacement variation. Under the multiparameter action, the rainfall intensity, cohesion, and friction angle are positively coupled, so the mean value and standard deviation of C13 are large. When the three factors work together with the saturated permeability coefficient, the direct effect and coupling effect of C15 on the displacement change are significantly reduced; that is, the saturated permeability coefficient has a negative effect on the displacement change of characteristic points of slope stability.

It can be seen from Figure 9(d) that in terms of the direct effect represented by the mean value, C3>C2>C4>C1, and the standard deviation is C3>C4>C1>C2 in turn. It can be seen from the comprehensive mean value and standard deviation coefficient that the saturated permeability coefficient has the weakest influence on the equivalent plastic strain of slope stability.

From the mean value of the two-parameter combination, C6>C10>C5>C8>C9>C9, and the standard deviation from high to low is C10>C6>C8>C9>C7>C5. From the comprehensive mean value and standard deviation, it can be seen that the combination of the saturated permeability coefficient, internal friction angle, and cohesion has a more significant effect on the plastic strain characteristics of slope stability than other two-parameter combinations, and the cohesion has a positive coupling effect with them. There is a negative relationship between rainfall intensity and other three parameters.

The results show that the direct effect and coupling effect of C15 are at a small level, which indicates that the interaction coupling effect is not the effect of a single parameter when the four parameters fluctuate together.

#### 5. Conclusion

In view of the stability of high and steep slopes under complex conditions, this paper analyzes the stability of slopes in a ski resort under four different rainfall conditions by using GeoStudio finite element analysis software. The coupling analysis is carried out by using the method of parameter combination analysis, and the following conclusions are drawn: (1)For the same slope, the higher the initial rainfall, the faster the soil saturation speed and the faster the safety coefficient decreases. Different rainfall types have a significant impact on the slope safety coefficient. Qianfeng-type rainfall (the minimum safety factor) occurs at the end of summer, so it is necessary to avoid construction at the end of summer or do safety monitoring at the end of summer(2)The pore water pressure of each monitoring point increases with the increase of rainfall. At the slope toe, the initial pore water pressure shows a trend of decreasing first and then rising because the rainfall intensity is less than the permeability strength. In addition, the influence of the rainfall intensity and permeability coefficient on slope stability is not simple linear addition and subtraction(3)For different types of rainfall, the trend and value of the safety factor decrease are basically the same. However, considering the timeliness, the impact of different rainfall types on the safety coefficient of the slope is different during rainfall. For the large rainfall, such as the forward type, the sudden drop occurs in the early stage of rainfall, and the latter with small rainfall is slowly decreased in the early stage. This also provides the theoretical basis of time effectiveness for slope support and reinforcement(4)The direct and coupling effects of the saturated permeability coefficient, cohesion, internal friction angle, and external rainfall intensity on the stability of slopes are at a low level. With the joint changes of the four parameters, the interaction coupling is stacked by algebra with the effect not equal to the single parameter

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.