Rain-on-snow events in cold regions frequently trigger slope failures. This study elucidates the instability mechanism of low-permeability silty clay slopes under combined rainfall and snowmelt conditions. A refined numerical model was established based on the sequential coupling of SEEP/W and SLOPE/W, utilizing the Morgenstern-Price method for stability analysis. A rigorous mesh sensitivity analysis confirmed that a locally refined mesh of 0.2 m with exponential time-stepping is essential to eliminate numerical dispersion at the wetting front. Simulation results indicate a significant time-lag effect in stability response; the critical failure time lags behind rainfall cessation (e.g., ~8 h for moderate rain) due to gravity-driven moisture redistribution. Spatially, the slope toe reaches saturation first, generating excess pore-water pressure and suggesting a tendency toward retrogressive instability. Furthermore, snowmelt superposition functions as a continuous hydraulic load, creating a base flow effect that advances the acceleration phase of failure by 1–2 h and further reduces the minimum safety factor. These findings highlight the critical role of the slope toe saturation and the necessity of considering snowmelt intensity in landslide early warning systems for cold regions.
Tang et al. (Sat,) studied this question.