PurposeThis paper aims to develop an efficient algorithm combining straightforward response surface functions with Monte Carlo simulation to conduct seismic reliability analysis in a systematical way.Design/methodology/approachThe representative slip surfaces are identified and based on to calibrate multiple response surface functions with acceptable accuracy. The calibrated response surfaces are used to determine the yield acceleration in Newmark sliding displacement analysis. Then, the displacement-based limit state function is adopted to conduct seismic reliability analysis.FindingsThe calibrated response surface functions have fairly good accuracy in predicting the yield acceleration in Newmark sliding displacement analysis. The seismic reliability is influenced by such factors as PGA, spatial variability and threshold value. The proposed methodology serves as an effective tool for geotechnical practitioners.Originality/valueThe multiple sources of a seismic slope response can be effectively determined using the multiple response surface functions, which are easily implemented within geotechnical engineering.