A 3D seepage flow numerical simulation model was developed for seepage analysis of a landslide dam. A 3D seepage flow numerical simulation model coupled with a 2D surface flow and erosion/deposition model was developed for seepage analysis of a slope due to a rainfall event. The conventional water-phase (one-phase) seepage flow model assumes only water phase flow in seepage analysis, which is inadequate for unsaturated soil domains. Hence, a water-air two-phase seepage flow model that considers both water and air phase in the seepage flow process is also used for seepage analysis. Pore water pressure and moisture content data obtained by the seepage flow model were then used to analyze the stability of the slope. Janbu’s simplified method as well as the extended Spencer method was used for the stability analysis. Numerical simulation results and experimental measurements are satisfactorily in agreement.