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.

Debris flows are a type of natural disaster induced by vegetation-water-soil coupling under external dynamic conditions. Research on the mechanism by which underground plant roots affect the initiation of gulley debris flows is currently limited. To explore this mechanism, we designed 14 groups of controlled field-based simulation experiments. Through monitoring, analysis, calculation, and simulation of the changes in physical parameters, such as volumetric water content, pore-water pressure, and matric suction, during the debris flow initiation process, we revealed that underground plant roots change the pore structure of soil masses. This affects the response time of pore-water pressure to volumetric water content, as well as hydrological processes within soil masses before the initiation of gully debris flows. Underground plant roots increase the peak volumetric water content of rock and soil masses, reduce the rates of increase of volumetric water content and pore-water pressure, and increase the dissipation rate of pore-water pressure. Our results clarify the influence of underground roots on the initiation of gulley debris flows, and also provide support for the initiation warning of gully debris flow. When the peak value of stable volumetric water content is taken as the early warning value, the early warning time of soil with underground plant roots is delayed by 534 to 1253 s. When the stable peak value of pore-water pressure is taken as the early warning value, the early warning time of soil with underground plant roots is delayed by 193 to 1082 s. This study provides a basis for disaster prevention and early warning of gully debris flows in GLP, and also provides ideas and theoretical basis under different vegetation-cover conditions area similar to GLP.

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.

The overconsolidation ratio (OCR) is a critical factor in determining the mechanical behaviour of overconsolidated clays. On the basis of the three requirements for the peak strength line, a continuous and smooth peak strength line is constructed from the perspective of the peak stress ratio, and then a new yield function for overconsolidated clays is developed. The developed yield function in the stress space is characterized by an elliptical curve. The evolution of the developed yield function in the stress space is captured by a new hardening parameter, which is constructed by integrating the proposed peak strength surface with the subloading surface concept. By combining the developed yield function with the non-orthogonal plastic flow rule, a non-orthogonal elastoplastic constitutive model of overconsolidated clays is established to consider the influence of the OCR on strength and deformation. The proposed model requires seven material parameters, all of which have a clear physical meaning and can be easily determined via conventional laboratory tests. Three typical stress paths are employed to demonstrate the essential features of the proposed model. The effectiveness of the proposed model is confirmed by comparing the experimental data with corresponding model predictions.

Loose sandy soil layers are prone to liquefaction under strong earthquakes, causing damage to civil engineering structures inside or upon the liquefied ground. According to the present Japanese design guideline on liquefaction countermeasures for river levees, the entire depth of the liquefiable subsoil below river embankments should be improved. However, this approach is not economical against deep liquefiable subsoil. To rationalize the design approach, this contribution investigated the performance of a floating-type cement treatment method, in which only the shallower part of the liquefiable subsoil is reinforced. A series of centrifuge shaking table model tests was conducted under a 50g environment. The depth of improvement (cement treatment) was varied systematically, and the effect of the sloping ground was examined. The experimental results revealed that the settlements of river embankments can be reduced linearly by increasing the depth of improvement. Moreover, the acceleration of embankments can be reduced drastically by the vibration-isolation effect between the cement-treated soil and the liquefiable soil. These effects contribute to the safe retention of the embankment shape even when the liquefied sloping ground causes lateral flows. Towards practical implementation, discussions on the effect of permeability on cement-treated soil were expanded. Furthermore, the stress acting on cement-treated soil during shaking was measured using an acrylic block to explain the occurrence of cracks in the soil. (c) 2025 Japanese Geotechnical Society. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http:// creativecommons.org/licenses/by-nc-nd/4.0/).

Study region: Indus Basin Study focus: Meteorological droughts can result in hydrological and soil moisture droughts with severe consequences for food production. In the Indus basin there are strong upstream-downstream linkages and upstream droughts may have strong downstream impacts. This study identifies periods of meteorological, hydrological and soil moisture drought in the Indus Basin for the period 1981-2010, analyses drought propagation and evaluates the role of meltwater in mitigating drought. We used outputs from a cryosphere-hydrology model (SPHY) and a crop-hydrology model (LPJmL), analysed the Standardized Precipitation Evapotranspiration Index (SPEI), the Standardized Streamflow Index (SSI), Soil Moisture Anomaly Index (SMAI) and crop yield, which are used as drought indicators to identify periods of drought, analyse drought propagation and its impacts. New hydrological insights for the region: Propagation of meteorological drought to hydrological drought and hydrological drought to soil moisture drought shows varied patterns and lag times. There were slightly more periods of soil moisture drought when meltwater was not available than when meltwater was available for irrigation. Our results show that identifying the link between soil moisture drought and yield anomaly remains challenging due to differences in temporal resolution of the data. Nevertheless, the results highlight the critical role of meltwater in mitigating yield variability, especially in the more downstream areas. This provides insight into the potential consequences of future cryosphere degradation for food production in the future.

This study presents a method for remediating soils contaminated by organic pollutants through the selective blocking of pores. This technique is based on the use of yield stress fluids, specifically concentrated biopolymer solutions, which, due to their distinctive rheological properties, preferentially flow through high-conductance flow paths. Following the injection of yield stress fluid, its presence redirects subsequent water flow towards the pores that are typically unswept during standard waterflooding. Laboratory experiments at the pore scale were conducted to validate this method and confirm previous findings from core-flooding experiments. Aqueous xanthan gum solutions were used as microscopic blocking agents in well-characterized micromodels exhibiting microscopic heterogeneities in pore size. The impact of polymer concentration, soil wettability and operating conditions (injection pressure and flow rate) on the residual pollutant saturation following treatment was analyzed, enabling the optimization of the remediation strategy. The use of xanthan gum as a blocking agent led to a significant improvement in pollutant removal compared to conventional waterflooding, delivering consistently better results across all cases studied. The method demonstrated strong performance in water-wet medium, with the average polymer concentration yielding the highest efficiency in pollutant removal.

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.

Seepage problems in half-space domains are crucial in hydrology, environmental, and civil engineering, involving groundwater flow, pollutant transport, and structural stability. Typical examples include seepage through dam foundations, coastal aquifers, and levees under seepage forces, requiring accurate numerical modeling. However, existing methods face challenges in handling complex geometries, heterogeneous media, and anisotropic properties, particularly in multi-domain half-spaces. This study addresses these challenges by extending the modified scaled boundary finite element method (SBFEM) and using this method to explore steady seepage problems in complex half-space domain. In the modified SBFEM framework, segmented straight lines or curves, parallel to the far-field infinite boundary, are introduced as scaling lines, with a one-dimensional discretization applied to them, thereby reducing computational costs.Then the weighted residual method is applied to obtain the modified SBFEM governing equations and boundary conditions of steady-state seepage problem according to the Laplace diffusion equation and Darcy's law. Furthermore, the steady seepage matrix at infinity is obtained by solving the eigenvalue problem of Schur decomposition and then the 4th-order Runge-Kutta algorithm is used to iteratively solve until the seepage matrix at the boundary lines is reached. Comparisons between the present numerical results and solutions available in the published work have been conducted to demonstrate the efficiency and accuracy of this method. At the same time, the influences of the geometric parameters and complex half-space domain on the seepage flow characteristics in complex half-space domain are investigated in detail.

Cohesion provided by pore ice is a critical component influencing the mechanical behavior of frozen soil, as it not only cements soil particles together but also shares the external loads with them. In view the crucial role of cohesion in developing an elastoplastic model for frozen soil, this paper employs triaxial tensile strength (TTS) to characterize cohesion and proposes a TTS degradation expression driven by plastic shear strain. By directly incorporating TTS into the yield function, a framework for a Non-Orthogonal Elastoplastic (NOEP) constitutive model that accounts for cohesion degradation in frozen soil is developed. Furthermore, a hardening parameter incorporating TTS is introduced and used in conjunction with the modified yield function to determine the magnitude of the plastic strain increment. The non-orthogonal plastic flow rule is used to determine the direction of the plastic strain increment based on the modified yield function. Ultimately, by combining the elastic strain increment determined by Hooke's rule, a NOEP constitutive model incorporating cohesion degradation for frozen soil is established. The validity and rationality of the proposed NOEP model in representing the stress-strain relationship of frozen soil are confirmed through comparisons with test results of frozen soil under the triaxial compression conditions. The proposed constitutive model provides a more comprehensive and precise representation of frozen soil's response to external loading, enhancing the understanding of its shear deformation behavior and providing a robust theoretical foundation for engineering design and construction in cold regions.