PurposeFor gravity dams built on foundations with directional joint sets, the seepage in the foundation possesses anisotropic characteristics and may have adverse effects on the foundation stability. A methodology for system reliability analysis of gravity dam foundations considering anisotropic seepage and multiple sliding surfaces is proposed in this paper.Design/methodology/approachAnisotropic seepages in dam foundations are simulated using finite element method (FEM) with the equivalent continuum model (ECM), and their effect on dam foundation stability is involved by uplift pressures acting on the potential sliding surfaces. The system failure probability of the dam foundation is efficiently estimated using Monte Carlo method (MCM) combined with response surface method (RSM).FindingsThe case study shows that it is necessary to consider the possibly adverse effect of anisotropic seepage on foundation stability of gravity dams and the deterministic analysis of the foundation stability may be misleading. The system reliability analysis of the dam foundation is justified, as the uncertainties in shear strength parameters of the foundation rocks and joint sets as well as aperture, connectivity and spacing of the joint sets are quantified and the system effect of the multiple potential sliding surfaces on the foundation reliability is reasonably considered.Originality/value(1) A methodology is proposed for efficient system reliability analysis of foundation stability of gravity dams considering anisotropic seepage and multiple sliding surfaces (2) The influence of anisotropic seepage on the stability of gravity dam foundation is revealed (3) The influence of estimation errors of RSMs on the system reliability assessment of dam foundation is investigated.