In this paper, a numerical procedure for probabilistic slope stability analysis is presented. This procedure extends the traditional limit equilibrium method of slices to a probabilistic approach that accounts for the uncertainties and spatial variation of the soil
strength parameters. In this study, two-dimensional random fields were generated based on a Karhunen-Loeve expansion in a fashion consistent with a specified marginal distribution function and an autocorrelation function. A Monte Carlo simulation was then used to determine the statistical response based on the generated random fields. This approach makes no assumption about the critical failure surface. Rather, the critical failure surface corresponding to the input random fields of soil properties is searched during the process of analysis. A series of analyses was performed to verify the application potential of the proposed method and to study the effects of uncertainty due to the spatial heterogeneity on the stability of slope. The results show that the proposed method can efficiently consider the various failure mechanisms caused by the spatial variability of soil property in the probabilistic slope stability assessment.