Probabilistic Analysis of Rainfall-Induced Slope Stability Using KL Expansion and Polynomial Chaos Kriging Surrogate Model

Rainfall infiltration is one of the main factors inducing slope instability, while the spatial heterogeneity and uncertainty of soil parameters have profound impacts on slope response characteristics and stability evolution. Traditional deterministic analysis methods struggle to reveal the dynamic risk evolution process of the system under heavy rainfall. Therefore, this paper proposes an uncertainty analysis framework combining Karhunen–Loève Expansion (KLE) random field theory, Polynomial Chaos Kriging (PCK) surrogate modeling, and Monte Carlo simulation to efficiently quantify the probabilistic characteristics and spatial risks of rainfall-induced slope instability. First, for key strength parameters such as cohesion and internal friction angle, a two-dimensional random field with spatial correlation is constructed to realistically depict the regional variability of soil mechanical properties. Second, a PCK surrogate model optimized by the LARS algorithm is developed to achieve high-precision replacement of finite element calculation results. Then, large-scale Monte Carlo simulations are conducted based on the surrogate model to obtain the probability distribution characteristics of slope safety factors and potential instability areas at different times. The research results show that the slope enters the most unstable stage during the middle of rainfall (36–54 h), with severe system response fluctuations and highly concentrated instability risks. Deterministic analysis generally overestimates slope safety and ignores extreme responses in tail samples. The proposed method can effectively identify the multi-source uncertainty effects of slope systems, providing theoretical support and technical pathways for risk early warning, zoning design, and protection optimization of slope engineering during rainfall periods.