Bedrock-soil layer slopes (BSLSs) are widely distributed in nature. The existence of the interface between bedrock and soil layer (IBSL) affects the failure modes of the BSLSs, and the seismic action makes the failure modes more complex. In order to accurately evaluate the safety and its corresponding main failure modes of BSLSs under seismic action, a system reliability method combined with the upper bound limit analysis method and Monte Carlo simulation (MCS) is proposed. Four types of failure modes and their corresponding factors of safety (Fs) were calculated by MATLAB program coding and validated with case in existing literature. The results show that overburden layer soil's strength, the IBSL's strength and geometric characteristic, and seismic action have significant effects on BSLSs' system reliability, failure modes and failure ranges. In addition, as the cohesion of the inclination angle of the IBSL and the horizontal seismic action increase, the failure range of the BSLS gradually approaches the IBSL, which means that the damage range becomes larger. However, with the increase of overburden layer soil's friction angle, IBSL's depth and strength, and vertical seismic actions, the failure range gradually approaches the surface of the BSLS, which means that the failure range becomes smaller.

The frequent occurrence of earthquakes worldwide has rendered highway slope protection projects highly vulnerable to damage from seismic events and their secondary disasters. This severely hampers the smooth implementation of post-disaster rescue and recovery efforts. To address this challenge, this study proposes a comprehensive method for assessing seismic losses in slope protection projects, incorporating factors such as topography and elevation to enhance its universality. The method categorizes seismic losses into two main components: damage to protection structures and costs associated with landslide and rockfall clearance and transportation. This study estimates the cost range for common protection structures and clearance methods under general conditions based on widely recognized quota data in China. It establishes criteria for classifying the damage states of protection structures and provides loss ratio values based on real-world seismic examples and expert experience, constructing a model for assessing damage losses. Additionally, by summarizing the geometric characteristics of soil and rock accumulations on road surfaces, a method for estimating landslide volumes is proposed, considering the dynamic impact of slope gradients on clearance and transportation volumes, and a corresponding cost assessment model for clearance and transportation is developed. The feasibility and reliability of the proposed method are verified through two case studies. The results demonstrate that the method is easy to implement and provides a scientific basis for improving relevant standards and practices. It also offers an efficient and scientific tool for loss assessment to industry practitioners.

Mesh-free methods, such as the Smooth Particle Hydrodynamics (SPH) method, have recently been successfully developed to model the entire wetting-induced slope collapse process, such as rainfall-induced landslides, from the onset to complete failure. However, the latest SPH developments still lack an advanced unsaturated constitutive model capable of capturing complex soil behaviour responses to wetting. This limitation reduces their ability to provide detailed insights into the failure processes and to correctly capture the complex behaviours of unsaturated soils. This paper addresses this research gap by incorporating an advanced unsaturated constitutive model for clay and sand (CASM-X) into a recently proposed fully coupled seepage flow-deformation SPH framework to simulate a field-scale wetting-induced slope collapse test. The CASM-X model is based on the unified critical state constitutive model for clay and sand (CASM) and incorporates a void-dependent water retention curve and a modified suction-dependent compression index law, enabling the accurate prediction various unsaturated soil behaviours. The integration of the proposed CASM-X model in the fully coupled flow deformation SPH framework enables the successful prediction of a field-scale wetting-induced slope collapse test, providing insights into slope failure mechanisms from initiation to post-failure responses.

Seismic fragility denotes the probabilities of a system exceeding some prescribed damage levels under a range of seismic intensities. Classical seismic fragility studies in slope engineering usually construct fragility functions by making some assumptions for fragility curve shape, and always neglect spatial variability of soil materials. In this study, an assumption-free method on the basis of probability density evolution theory (PDET) is proposed for seismic fragility assessment of slopes. The random input earthquakes and spatially-variable soil parameters in slope are simultaneously quantified. By the proposed method, assumption-free fragility curves of a slope are established without limiting the fragility curve shape. The obtained fragility results are also compared with those from two classic parametric fragility methods (linear regression and maximum likelihood estimation) and Monte Carlo simulation. The results demonstrate that the proposed assumption-free method has potential to gives more rigorous and accurate fragility results than classical parametric fragility analysis methods. With the proposed method, more reliable fragility results can be obtained for slope seismic risk assessment.

Study area: The Binggou and adjacent Yakou catchments in the northeastern Tibetan Plateau. Study focus: Hillslope flow paths were studied using hydrochemical data of various water types in the spring snowmelt and summer rainfall periods based on hydrochemical tracers and endmember mixing analysis. New hydrological insights for the study region: End-member mixing analysis confirmed the dominance of surface and near-surface runoff during the spring snowmelt. Specifically, the spring Binggou stream water had 61 % surface runoff, 22 % shallow groundwater, and 17 % near-surface runoff. The spring Yakou stream water had 64 % snowmelt, 25.5 % near-surface runoff, and 10.5 % riparian saturated soil water at a depth of 20 cm. The application of end-member mixing analysis failed in the summer rainfall period, and shallow subsurface flow contributed the most to the streamflow (similar to 100 %). The average acid-neutralizing capacity of the spring Yakou stream water was 611 mu eq/L, increasing to 841 mu eq/L in the summer, and for the Binggou stream water, the values were 747 mu eq/L and 1084 mu eq/L, respectively, indicating that the thawed soil layers had a significant buffering effect on stream water chemistry. This study revealed seasonal shifts in flow paths and stream sources, with a transition from surface to subsurface flow influenced by meteorological conditions and the active layer thickness. Future climate change may enhance subsurface flow recharge, leading to less diluted streamflow and stronger water-soil interactions.

In the loess tableland, gully slope instability induces severe soil erosion and land degradation, yet the synergistic effects of dominant vegetation under varying restoration modes combined with dynamic rainfall regimes and topographic variations on gully slope stabilization mechanisms remain inadequately quantified. Therefore, the dominant vegetation species under natural (NR) and artificial restoration (AR) was chosen as the object. Through field sampling, root-soil complex mechanical experiments, and numerical simulations, the protection effect of dominant vegetation under different restoration modes combination with rainfall and topographic variations was investigated. The result revealed significant differences in basic soil physical properties, root morphological characteristics, root and root-soil complex mechanical properties among five dominant vegetated plots under the different restoration modes (P < 0.05). The soil properties in the Scop plot under AR were slightly better than those in the other plots. The roots in the Spp plot developed better under NR. The shear strength of Lespedeza bicolor Turcz. was the highest under NR. The tensile strength of Digitaria sanguinalis (L.) Scop. was greatest under AR. The tensile force and tensile strength of single roots exhibited a significant positive linear correlation and a significant negative exponential correlation, with root diameter, respectively (P < 0.01). For the unstable gully slopes (F-s < 1.0), maximum displacement occurred at the slope foot, where tensile shear failure dominated, while the interior experienced compressive yielding. The grey relational analysis identified rainfall intensity as the primary destabilizing factor, followed by dominant vegetation species, slope height, and slope gradient. Notably, when rainfall intensity reaches or exceeds 0.06 m/h, or when slope height exceeds 20 m combined with long-duration rainfall, the regulatory impacts of dominant vegetation under different restoration modes on the gully slope stability are substantially diminished and become negligible. This study provides a theoretical basis for gully slope protection and ecological environmental construction in loess tableland.

Predicting cumulative surface slope displacements induced by rainfall infiltration is crucial for accurately assessing the risks to potentially affected infrastructure. In this paper the numerical modelling of the case history of Miscano slope is presented. Plaxis 2D code has been used adopting two constitutive laws: the linear elastoplastic model (Mohr-Coulomb, MC) and the Hardening Soil with small strain stiffness (HSsmall). The aim is to test the suitability of these constitutive laws in predicting the hydro-mechanical behaviour of clayey soil slope. Based on long-term field measurements, the parameters of MC and HSsmall have been determined by back analysing the first-year field measurements in terms of cumulative surficial horizontal displacements and pore water pressure. Subsequently, the numerical models have been validated against the analogous field measurements collected from the second year. The numerical models predict with a good agreement the field measurements for both years. In terms of cumulative surficial horizontal displacements, the HSsmall underestimates the field measurements by 21.2% at the end of the first year, while that based on MC exhibits a 32.8% overestimation. Moreover, the initialization procedure clearly affects the cumulative surficial horizontal displacements results obtained with both the HSsmall and MC models for the second year. In fact, the best results have been achieved when the second-year net rainfall have been applied starting from the initial phase used to generate the lithostatic stress state.

Flow instability impacts negatively on hydraulic structures. Changes in water pressure or the periodic impact of water flows cause structural damage to channels. The rapid increase in water depth leads to overflows or sprays, which erode soil adjacent to channels. In this study, flow instability was examined through the basis of theories and experiments. The theoretical discriminants for flow instability were inferred by Vedernikov number and the effect of slopes on the Froude number was considered. A rectangular cross- channel was selected for the experiments. The experimental results were compared with theories, it was shown that when the flow conditions were on the margin of instability, the discriminant established by this study is able to accurately determine the occurrence of instability. Through this new discriminant, the discrepancy which appears in traditional method can be avoided. The presented results are ideal for channel design and offer new approaches for flow instability prevention.

This study investigated how soil properties affect levee erosion and foundation scouring by evaluating the behavior of loose and cohesive (mixed) soils beneath a rigid crest under overflow conditions and analyzing flow dynamics within the scoured hole to understand the scouring mechanism. Four cases were examined with varying overtopping depths (Od): LS-FS, LS-FM, and LM-FS, at Od = 2 cm, and LS-FM at Od = 3 cm, where 'L' stands for levee, 'F' for foundation, 'S' for sand (#8), and 'M' for mixed soil (20% silt + 80% sand #8). The results revealed distinct differences among the cases. Notably, erosion of the back slope in the LM-FS case was delayed fourfold compared to LS-FS. In the LS-FM case, breaching of the levee body was delayed by 1.6 times compared to the LS-FS case with a 2 cm overtopping depth. Moreover, different scour hole geometries with complex flow patterns occurred in different timespans. Particle image velocimetry (PIV) was utilized on two physical scoured hole models to analyze the flow behavior within these scoured holes. The PIV analysis revealed the formation of twin eddies, moving in opposite directions and shaped by the nappe flow jet, which was instrumental in the development of the scour holes. This study found that foundation cohesion is more essential than the levee body in delaying levee breaches under rigid crest. Additionally, it revealed the role of twin eddies, especially the levee-side eddy, in increasing the size of the scoured hole upstream and causing levee breaches.

In the dynamic response analysis of slopes, the displacement of sliding surfaces is an important indicator for assessing stability. However, due to the uniform dynamic parameters of the Newmark slide block method, it is difficult to accurately analyze the displacements of large-scale slopes. To address this issue, the spatial distribution characteristics of dynamic parameters need to be considered to accurately analyze the stability of slopes. Under the combined action of rainfall and reservoir water level change, the Shiliushubao old landslide in the Three Gorges Reservoir area remains stable. To investigate the seismic stability of slopes, simulated seismic waves were employed. Firstly, the dynamic triaxial test, designed with cyclic loading, was employed to investigate the variation rules of the dynamic parameters of slope soil, and to establish a functional relationship. Then, the stochastic seismic motion model was used to generate artificially seismic waves in the Three Gorges Reservoir Area. Finally, to assess the stability of the old landslide, finite element software, GeoStudio 2018 was used to obtain the spatial distribution characteristics of the dynamic parameters and to calculate the permanent displacements of the reservoir bank slope by inputting random ground motion loads and dynamic characteristic functions. It is demonstrated that under the most unfavorable working conditions of heavy rainfall and the earthquake in the specific region, the permanent displacement of the Shiliushubao old landslide will be less than the critical permanent displacement, the old landslide is not to experience instability again.