摘要:AbstractThis study aims to extend the multivariate adaptive regression splines (MARS)‒Monte Carlo simulation (MCS) method for reliability analysis of slopes in spatially variable soils. This approach is used to explore the influences of the multiscale spatial variability of soil properties on the probability of failure (Pf) of the slopes. In the proposed approach, the relationship between the factor of safety and the soil strength parameters characterized with spatial variability is approximated by the MARS, with the aid of Karhunen–Loève expansion. MCS is subsequently performed on the established MARS model to evaluatePf. Finally, a nominally homogeneous cohesive‒frictional slope and a heterogeneous cohesive slope, which are both characterized with different spatial variabilities, are utilized to illustrate the proposed approach. Results showed that the proposed approach can estimate thePfof the slopes efficiently in spatially variable soils with sufficient accuracy. Moreover, the approach is relatively robust to the influence of different statistics of soil properties, thereby making it an effective and practical tool for addressing slope reliability problems concerning time-consuming deterministic stability models with low levels ofPf. Furthermore, disregarding the multiscale spatial variability of soil properties can overestimate or underestimate thePf. Although the difference is small in general, the multiscale spatial variability of the soil properties must still be considered in the reliability analysis of heterogeneous slopes, especially for those highly related to cost effective and accurate designs.Graphical abstractDisplay OmittedHighlights•Efficient slope reliability analysis using MARS and K-L expansion is proposed.•Influence of soil statistics and multi-scale spatial variability is fully studied.•Neglecting the soil multi-scale spatial variability would give unreasonablePf.