Stability analysis of loose accumulation slopes under rainfall:case study of a high‑speed railway in Southwest China

The high and steep slopes along a high-speed railway in the mountainous area of Southwest China are mostly composed of loose accumulations of debris with large internal pores and poor stability,which can easily induce adverse geological disasters under rainfall conditions.To ensure the smooth construction of the high-speed railway and the subsequent safe operation,it is necessary to master the stability evolution process of the loose accumulation slope under rainfall.This article simulates rainfall using the finite element analysis software’s hydromechanical coupling module.The slope stability under various rainfall situations is calculated and analysed based on the strength reduction method.To validate the simulation results,a field monitoring system is established to study the deformation characteristics of the slope under rainfall.The results show that rainfall duration is the key factor affecting slope stability.Given a constant amount of rainfall,the stability of the slope decreases with increasing duration of rainfall.Moreover,when the amount and duration of rainfall are constant,continuous rainfall has a greater impact on slope stability than intermittent rainfall.The setting of the field retaining structures has a significant role in improving slope stability.The field monitoring data show that the slope is in the initial deformation stage and has good stability,which verifies the rationality of the numerical simulation method.The research results can provide some references for understanding the influence of rainfall on the stability of loose accumulation slopes along high-speed railways and establishing a monitoring system.

Stability analysis of a slope containing water-sensitive mudstone considering different rainfall conditions at an open-pit mine

Mudstone, as a typical soft rock with wide distribution, has been endangering the slopes containing mudstone by its watersensitivity of swelling and weakening strength when encountering water. To comprehensively understand the water-sensitivity of mudstone and reveal its infuence on slope stability, we took the working slope containing water-sensitive mudstone of Shengli No.1 open-pit coal mine in Xilinhot, Inner Mongolia, China, as an example. Mudstone samples taken from the working slope were remodeled and saturated, and then triaxial tested to obtain the efective cohesion and efective internal friction angle. The flter paper method was used to obtain the soil–water characteristic curve of unsaturated mudstone. The pore structure of mudstone samples with diferent water contents were analyzed using the mercury intrusion porosimetry tests combined with the fractal dimension. The total pore content of the mudstone sample with lower water content is greater than that of the mudstone sample with higher water content. The mesopores are more in the mudstone sample with lower water content, while the small pores are more in the mudstone sample with higher water content. The variation of water content will change the complexity of mudstone pore structure. The higher the water content, the simpler the mudstone pore structure and the smoother the pore surface. Numerical calculations were conducted on the stability of the working slope under diferent rainfall conditions. The efective saturation on the mudstone layer surface changed and the plastic strain all occurred on the mudstone steps under diferent rainfall conditions. The key to preventing landslide of the slope containing water-sensitive mudstone in Shengli No.1 open-pit coal mine is to control the deformation and sliding of the mudstone layer.

Experimental study on instability mechanism and critical intensity of rainfall of high-steep rock slopes under unsaturated conditions

Two critical factors,namely intense precipitation and intricate excavation,can trigger rock mass disasters in mining operations.In this study,an indoor rainfall system was developed to precisely regulate the flow and intensity of precipitation.A large-scale model experiment was conducted on a self-designed physical simulation experiment platform to investigate the failure and instability of high-steep rock slopes under unsaturated conditions.The real-time reproduction of the progressive failure process in high-steep rock slopes enabled the determination of the critical rainfall intensity and revealed the mechanism underlying slope instability.Experiment results indicated that rainfall may be the primary factor contributing to rock mass instability,while continuous pillar mining exacerbates the extent of rock mass failure.The critical failure stage of high-steep rock slopes occurs at a rainfall intensity of 40 mm/h,whereas a rainfall exceeding 50 mm can induce critical instability and precipitation reaching up to 60 mm will result in slope failure.The improved region growing segmentation method(IRGSM)was subsequently employed for image recognition of rock mass deformation in underground mines.Herein an error comparison with the simple linear iterative cluster(SLIC)superpixel method and the original region growing segmentation method(ORGSM)showed that the average identification error in the X and Y directions by the method was reduced significantly(1.82%and 1.80%in IRGSM;4.70%and 6.26%in SLIC;9.45%and 12.40%in ORGSM).Ultimately,the relationship between rainfall intensity and failure probability was analyzed using the Monte Carlo method.Moreover,the stability assessment criteria of rock slope under unsaturated condition were quantitatively and accurately evaluated.

Flume experiments to study fine-grain migration and its impact on slope stability

Fine grains migration is a primary cause of landslides and debris flows.This study investigates the effect of fine-grain migration on slope failure through flume experiments,focusing on the spatiotemporal characteristics and mechanisms of slope stability.A series of artificial rainfall flume experiments with varying rainfall intensities and slopes were conducted using soil samples collected from Wei Jia Gully.The experiments monitored pore-water pressure,grain migration,and failure sequences.Grain-size distribution parameters(μand Dc)were analyzed to understand the migration path and accumulation of fine grains.The experiments reveal that fine-grain migration significantly alters soil structure,leading to random blockage and interconnection of internal pore channels.These changes result in fluctuating pore-water pressure distributions and uneven fine-grain accumulation,critical factors in slope stability.Slope failures occur randomly and intermittently,influenced by fine-grain content in runoff and resulting pore-water pressure variations.This study highlights that fine-grain migration plays a vital role in slope stability,with significant implications for predicting and mitigating slope failures.The stochastic nature of fine-grain migration and its impact on soil properties should be incorporated into predictive models to enhance their accuracy and reliability.

Three-dimensional numerical analysis of plant-soil hydraulic interactions on pore water pressure of vegetated slope under different rainfall patterns

Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.However,the hydrological response of vegetated slopes,especially three-dimensional(3D)slopes covered with shrubs,under different rainfall patterns remains unclear and requires further investigation.To address this issue,this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil.Three series of numerical parametric studies are conducted to investigate the influences of slope inclination,rainfall pattern and rainfall duration.Four rainfall patterns(advanced,bimodal,delayed,and uniform)and two rainfall durations(4-h intense and 168-h mild rainfall)are considered to study the hydrological response of the slope.The computed results show that 17%higher transpiration-induced suction is found for a steeper slope,which remains even after a short,intense rainfall with a 100-year return period.The extreme rainfalls with advanced(PA),bimodal(PB)and uniform(PU)rainfall patterns need to be considered for the short rainfall duration(4 h),while the delayed(PD)and uniform(PU)rainfall patterns are highly recommended for long rainfall durations(168 h).The presence of plants can improve slope stability markedly under extreme rainfall with a short duration(4 h).For the long duration(168 h),the benefit of the plant in preserving pore-water pressure(PWP)and slope stability may not be sufficient.

Influence of joint spacing and rock characteristics on the toppling stability of cut rock slope through a simplified limit equilibrium method

Toppling failure of rock mass/soil slope is an important geological and environmental problem.Clarifying its failure mechanism under different conditions has great significance in engineering.The toppling failure of a cutting slope occurred in a hydropower station in Kyushu,Japan illustrates that the joint characteristic played a significant role in the occurrence of rock slope tipping failure.Thus,in order to consider the mechanical properties of jointed rock mass and the influence of geometric conditions,a simplified analytical approach based on the limit equilibrium method for modeling the flexural toppling of cut rock slopes is proposed to consider the influence of the mechanical properties and geometry condition of jointed rock mass.The theoretical solution is compared with the numerical solution taking Kyushu Hydropower Station in Japan as one case,and it is found that the theoretical solution obtained by the simplified analysis method is consistent with the numerical analytical solution,thus verifying the accuracy of the simplified method.Meanwhile,the Goodman-Bray approach conventionally used in engineering practice is improved according to the analytical results.The results show that the allowable slope angle may be obtained by the improved Goodman-Bray approach considering the joint spacing,the joint frictional angle and the tensile strength of rock mass together.

Infiltration,runoff,and slope stability behaviors of infinite slope with macropores based on an improved Green–Ampt model

Infiltration–runoff–slope instability mechanism of macropore slope under heavy rainfall is unclear.This paper studied its instability mechanism with an improved Green–Ampt(GA)model considering the dual-porosity(i.e.,matrix and macropore)and ponding condition,and proposed the infiltration equations,infiltration–runoff coupled model,and safety factor calculation method.Results show that the infiltration processes of macropore slope can be divided into three stages,and the proposed model is rational by a comparative analysis.The wetting front depth of the traditional unsaturated slope is 17.2%larger than that of the macropore slope in the early rainfall stage and 27%smaller than that of the macropore slope in the late rainfall stage.Then,macropores benefit the slope stability in the early rainfall but not in the latter.Macropore flow does not occur initially but becomes pronounced with increasing rainfall duration.The equal depth of the wetting front in the two domains is regarded as the onset criteria of macropore flow.Parameter analysis shows that macropore flow is delayed by increasing proportion of macropore domain(ω_(f)),whereas promoted by increasing ratio of saturated permeability coefficients between the two domains(μ).The increasing trend of ponding depth is sharp at first and then grows slowly.Finally,when rainfall duration is less than 3 h,ωf andμhave no significant effect on the safety factor,whereas it decreases with increasingωf and increases with increasingμunder longer duration(≥3 h).With the increase ofω_(f),the slope maximum instability time advances by 10.5 h,and with the increase ofμ,the slope maximum instability time delays by 3.1 h.

Three-dimensional stability calculation method for high and large composite slopes formed by mining stope and inner dump in adjacent open pits

The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.

Classification and rating of disintegrated dolomite strata for slope stability analysis

Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gaining increasing prominence.This is primarily due to their unique properties,including low strength and loose structure.Current methods for evaluating slope stability,such as basic quality(BQ)and slope stability probability classification(SSPC),do not adequately account for the poor integrity and structural fragmentation characteristic of disintegrated dolomite.To address this challenge,an analysis of the applicability of the limit equilibrium method(LEM),BQ,and SSPC methods was conducted on eight disintegrated dolomite slopes located in Baoshan,Southwest China.However,conflicting results were obtained.Therefore,this paper introduces a novel method,SMRDDS,to provide rapid and accurate assessment of disintegrated dolomite slope stability.This method incorporates parameters such as disintegrated grade,joint state,groundwater conditions,and excavation methods.The findings reveal that six slopes exhibit stability,while two are considered partially unstable.Notably,the proposed method demonstrates a closer match with the actual conditions and is more time-efficient compared with the BQ and SSPC methods.However,due to the limited research on disintegrated dolomite slopes,the results of the SMRDDS method tend to be conservative as a safety precaution.In conclusion,the SMRDDS method can quickly evaluate the current situation of disintegrated dolomite slopes in the field.This contributes significantly to disaster risk reduction for disintegrated dolomite slopes.

Evaluation of slope stability through rock mass classification and kinematic analysis of some major slopes along NH-1A from Ramban to Banihal, North Western Himalayas

The network of Himalayan roadways and highways connects some remote regions of valleys or hill slopes,which is vital for India’s socio-economic growth.Due to natural and artificial factors,frequency of slope instabilities along the networks has been increasing over last few decades.Assessment of stability of natural and artificial slopes due to construction of these connecting road networks is significant in safely executing these roads throughout the year.Several rock mass classification methods are generally used to assess the strength and deformability of rock mass.This study assesses slope stability along the NH-1A of Ramban district of North Western Himalayas.Various structurally and non-structurally controlled rock mass classification systems have been applied to assess the stability conditions of 14 slopes.For evaluating the stability of these slopes,kinematic analysis was performed along with geological strength index(GSI),rock mass rating(RMR),continuous slope mass rating(CoSMR),slope mass rating(SMR),and Q-slope in the present study.The SMR gives three slopes as completely unstable while CoSMR suggests four slopes as completely unstable.The stability of all slopes was also analyzed using a design chart under dynamic and static conditions by slope stability rating(SSR)for the factor of safety(FoS)of 1.2 and 1 respectively.Q-slope with probability of failure(PoF)1%gives two slopes as stable slopes.Stable slope angle has been determined based on the Q-slope safe angle equation and SSR design chart based on the FoS.The value ranges given by different empirical classifications were RMR(37-74),GSI(27.3-58.5),SMR(11-59),and CoSMR(3.39-74.56).Good relationship was found among RMR&SSR and RMR&GSI with correlation coefficient(R 2)value of 0.815 and 0.6866,respectively.Lastly,a comparative stability of all these slopes based on the above classification has been performed to identify the most critical slope along this road.

A vector sum analysis method for stability evolution of expansive soil slope considering shear zone damage softening

Slope stability analysis is a classical mechanical problem in geotechnical engineering and engineering geology.It is of great significance to study the stability evolution of expansive soil slopes for engineering construction in expansive soil areas.Most of the existing studies evaluate the slope stability by analyzing the limit equilibrium state of the slope,and the analysis method for the stability evolution considering the damage softening of the shear zone is lacking.In this study,the large deformation shear mechanical behavior of expansive soil was investigated by ring shear test.The damage softening characteristic of expansive soil in the shear zone was analyzed,and a shear damage model reflecting the damage softening behavior of expansive soil was derived based on the damage theory.Finally,by skillfully combining the vector sum method and the shear damage model,an analysis method for the stability evolution of the expansive soil slope considering the shear zone damage softening was proposed.The results show that the shear zone subjected to large displacement shear deformation exhibits an obvious damage softening phenomenon.The damage variable equation based on the logistic function can be well used to describe the shear damage characteristics of expansive soil,and the proposed shear damage model is in good agreement with the ring shear test results.The vector sum method considering the damage softening behavior of the shear zone can be well applied to analyze the stability evolution characteristics of the expansive soil slope.The stability factor of the expansive soil slope decreases with the increase of shear displacement,showing an obvious progressive failure behavior.

Early warning system for shallow landslides using rainfall threshold and slope stability analysis

A combined cluster and regression analysis were performed for the first time to identify rainfall threshold that triggers landslide events in Amboori, Kerala, India. Amboori is a tropical area that is highly vulnerable to landslides. The 2, 3, and 5-day antecedent rainfall data versus daily rainfall was clustered to identify a cluster of critical events that could potentially trigger landslides. Further, the cluster of critical events was utilized for regression analysis to develop the threshold equations. The 5-day antecedent(xvariable) vs. daily rainfall(y-variable) provided the best fit to the data with a threshold equation of y = 80.7-0.1981 x. The intercept of the equation indicates that if the 5-day antecedent rainfall is zero, the minimum daily rainfall needed to trigger the landslide in the Amboori region would be 80.7 mm. The negative coefficient of the antecedent rainfall indicates that when the cumulative antecedent rainfall increases, the amount of daily rainfall required to trigger monsoon landslide decreases. The coefficient value indicates that the contribution of the 5-day antecedent rainfall is～20% to the landslide trigger threshold. The slope stability analysis carried out for the area, using Probabilistic Infinite Slope Analysis Model(PISA-m), was utilized to identify the areas vulnerable to landslide in the region. The locations in the area where past landslides have occurred demonstrate lower Factors of Safety(FS) in the slope stability analysis. Thus, rainfall threshold analysis together with the FS values from slope stability can be suitable for developing a simple, cost-effective, and comprehensive early-warning system for shallow landslides in Amboori and similar regions.

Stability analysis of unsaturated soil slope during rainfall infiltration using coupled liquid-gas-solid three-phase model

Generally, most soil slope failures are induced by rainfall infiltration, a process that involves interactions between the liquid phase, gas phase,and solid skeleton in an unsaturated soil slope. In this study, a loosely coupled liquid-gas-solid three-phase model, linking two numerical codes,TOUGH2/EOS3, which is used for water-air two-phase flow analysis, and FLAC^(3D), which is used for mechanical analysis, was established. The model was validated through a documented water drainage experiment over a sandy column and a comparison of the results with measured data and simulated results from other researchers. The proposed model was used to investigate the features of water-air two-phase flow and stress fields in an unsaturated soil slope during rainfall infiltration. The slope stability analysis was then performed based on the simulated water-air two-phase seepage and stress fields on a given slip surface. The results show that the safety factor for the given slip surface decreases first, then increases, and later decreases until the rainfall stops. Subsequently, a sudden rise occurs. After that, the safety factor decreases continually and reaches its lowest value, and then increases slowly to a steady value. The lowest value does not occur when the rainfall stops, indicating a delayed effect of the safety factor. The variations of the safety factor for the given slip surface are therefore caused by a combination of pore-air pressure, matric suction, normal stress, and net normal stress.

Influences of rainfall infiltration on stability of accumulation slope by in-situ monitoring test

In order to improve the understanding of the fundamental mechanism of rainfall infiltration induced landslides in accumulation slope and to clarify some important characteristics of slope performance,artificial rainfall simulation tests and field synthetic monitoring were carried out on a typical accumulation slope of Shangrui Freeway in Guizhou Province,China.The monitoring results show that the most accumulation landslides caused by rainfall infiltration are shallow relaxation failure,whose deformation zone lies within the top 0-4 m soil layer.The deformation of slope gradually reduces from the surface,where the greatest deformation lies in,to the deep part of slope.The average percentage of infiltration during the first 2 h is 86%,and then it reduces gradually with time because of the increase of the surface runoff.The average percentage of infiltration drop to a relatively stable value(50%)after 6 h.Rainfall infiltration causes obvious increase of pore-water pressure,which may result in a reduction of shear strength due to a decrease in effective stress and wetting-induced softening.The double-effect of rainfall infiltration is the main reason of rainfall infiltration induced landslides in accumulation slope.

Infiltration regulation and stability analysis of soil slope under sustained and small intensity rainfall

Rainfall infiltration depth and mode can severely influence slope stability.With the sustained rainfall,the influenced region of slope gradually expands.By using the Green-Ampt model to the soil slope,infiltration regulation was discussed under sustained and small intensity rainfall.And the infiltration rate of unsaturated soil was proposed according to the saturated infiltration theory.Because of the changing of initial moisture content in depth of slope,the saturated or unsaturated infiltration rate and depth could also be changeable with the sustained rainfall infiltration.Based on the principle of strength reduction,the calculation model of slope safety factor was established under different initial moisture contents and infiltration modes.Then,the slope stability was quantitatively analyzed through software FLAC3D.The calculation results of soil slope engineering show that there is a shorter period for slope stability under different initial moisture contents and unsaturated infiltration ways at the slope wetting front.The stability period of slope is 33.3%according to different initial moisture contents of wetting front less than that of the same initial moisture content of wetting front.And the slope is easier to fail under the unsaturated infiltration.The results agree well with the actual situation under sustained and small intensity rainfall.

Surficial Stability Analysis of Unsaturated Loess Slopes Subjected to Rainfall Infiltration Effects

According to theory of unsaturated soil strength and Green-Ampt model, an analysis method based on limit equilibrium theory is introduced to consider rainfall infiltration effects in loess slope stability analysis. The relationships between wetting band depth and surficial stability of slopes are analyzed. It is found that the infiltration adds to the weight of the soil and at the same time reduces the shear strength provided by matric suction of the soil. The wetting band depth plays a key role in the stability of slopes. The minimum rainfall intensity and the minimum rainfall duration needed to infiltrate to the wetting front depth are calculated based on the Green-Ampt model. The method in this paper will contribute to the predication of slope stability considering rainfall characteristics.

Stability analysis of an unsaturated soil slope considering rainfall infiltration based on the Green-Ampt model

Based on the principle of saturated infiltration and the Green-Ampt model,an unsaturated infiltration model for a soil slope surface was established for either constant moisture content,or depth-varying moisture content and the slope.Infiltration parameters in the partially saturated slope were revealed under sustained rainfall.Through analysis of the variation of initial moisture content in the slope,the ponding time,infiltration depth,and infiltration rate were deduced for an unsaturated soil slope subject to rainfall infiltration.There is no ponded water on the surface of the slope under sustained low-intensity rainfall.The results show that the infiltration parameters of an unsaturated slope are influenced by the initial moisture content and the wetting front saturation,the soil cohesion and rainfall intensity under sustained rainfall.More short-term slope failures can occur with the decrease of cohesion of the soil of the slope.The ponding time and infiltration depth differ considering constant or different initial moisture content respectively in the soil slope.Then,best-fit curves of the infiltration rate,ponding time,and infiltration depth to the wetting front saturation were obtained with constant or different initial moisture contents.And the slope failure time is roughly uniform when subject to a rainfall intensity I>5 mm/h.

The mechanism of slope instability due to rainfall-induced structural decay of earthquake-damaged loess

Natural loess slopes are characterized by a strong geological structure,which is an important factor in maintaining slope stability.The magnitude and duration of the earthquake may disturb the soil structure at different levels degrees,locally changing the arrangement between soil particles.The process of rainfall humidification weakens the cementation between soil particles,and the disturbance and humidification change the structural state of the soil,which in turn causes sliding of the slope along with the decay of soil mechanical properties.As slope instability is often the result of a series of post-earthquake ripple effects,it is of great scientific significance to study the mechanism of slope instability due to the structural decay of earthquake-damaged loess exacerbated by rainfall.In this paper,the impact of structural decay of loess on slope stability is simulated by GEOSTUDIO software under three conditions:pre-earthquake rainfall,post-earthquake rainfall and earthquake,taking the landslide in Buzi Village,Min County,Gansu Province as an example.The comparative analysis of the calculation results shows that the structural properties of the slope without earthquake disturbance are influenced by infil-tration amount.When it is fully saturated,the structural properties are similar to those of saturated soil,and the safety factor is reduced by 12.9%.In addition,the earthquake intensity and duration have different degrees of structural damage to the soil.When the structure is fully damaged,it is similar to that of remodelled soil,and the safety factor is reduced by 45.84%.Notably,the process of the earthquake and the following humidification generates the most serious damage to the loess structure,with a reduction in the safety factor of up to 56.15%.The quantitative analysis above obviously illustrates that the post-earthquake rainfall causes the most severe damage to structural loess slopes,and the resulting landslide hazard should not be underestimated.

Upper Bound Solution of Soil Slope Stability under Coupling Effect of Rainfall and Earthquake

Based on the limit analysis upper bound method,a new model of soil slope collapse has been proposed which consists of two rigid block zones and a plastic shear zone.Soil slope was induced failure by coupling effect of rainfall and earthquake,and these blocks were also incorporated horizontal earthquake force and vertical gravitate.The velocities and forces were analyzed in three blocks,and the expression of velocity discontinuities was obtained by the principle of incompressibility.The external force work for the blocks,the internal energy of the plastic shear zone and the velocity discontinuous were solved.The present stability ratios are compared to the prevenient research,which shows the superiority of the mechanism and rationality of the analysis.The critical height of the soil slope can provide theoretical basis for slope support and design.

Effects of intense rainfall on stability of infinite terraced slope

For fully understanding the hydrological dynamics of an infinite terraced slope, the infiltration process was studied by employing the Green and Ampt infiltration model. The limit equilibrium method and the Mohr-Coulomb failure criterion were adopted to derive a stability model for the infinite terraced slope subjected to an intense rainfall. Numerical simulation was performed for verifying its applicability. The results of numerical simulation indicate that a set of stepped wetting fronts are found during infiltration, and the infiltration of terraced slope covered by coarse-textured soils can be approximated as one-dimensional infiltration. The potential sliding surfaces from the numerical method are all parallel to the slope line, and the proposed model and framework can provide an approximate method of estimating how the infiltration affects the stability of an infinite terraced slope.