The geotechnical properties of municipal solid waste (MSW) in landfills vary considerably depending on the composition, time, and the rate of waste density. This variability in geotechnical properties leads to many uncertainties in the analysis of landfill slope stability. This study, using probabilistic methods investigates the slope stability and the probability of failure of a Kahrizak landfill under the conditions of spatial variability of physical and geotechnical properties of MSW. To achieve this goal, the random field theory has been used with a finite difference numerical method in the framework of the Monte Carlo simulation. MSW’s shear strength parameters such as cohesion and friction angle as well as the unit weight of layers are considered as random variables. An extensive literature review was conducted to address the probable variation of above-mentioned parameters in Kahrizak landfill. The results of several different laboratory researches on MSW samples collected from Kahrizak landfill were also employed in the modelings. Output results are presented in the form of probability distributions of the factor of safety as well as the probability of failure corresponding to these distributions. Moreover, the effect of various parameters such as coefficient of variation and correlation distance of input parameters on these output results has been investigated. The results show that considering the spatial variability in probabilistic methods would help to determine the various mechanisms affecting the performance and the probability of failure of the landfill slope. In addition, the output of such probabilistic analyses can be used as guidance for engineers to design safe and reliable landfill.