The unsaturated seepage field coupled with heavy rainfall-induced surface flow mainly accounts for the slope instability. If the slope contains macropores, the coupled model and solution process significantly differ from the traditional one (without macropores). Most of the studies on the variation of the water field under the coupled effect of runoff and seepage on the slope did not consider the macropore structure. In this paper, two coupled Richards equations were used to describe the MF (Macropore Flow), and along with the kinematic wave equation, they were applied to establish a coupled model of SR (Slope Runoff) and MF. The numerical solving of the coupled model was realized by the COMSOL PDE finite element method, and an innovative laboratory test was conducted to verify the numerical results. The effects of different factors (i.e., rainfall intensity, rainfall duration, saturated conductivity, and slope roughness coefficient) on water content and ponding depth with and without macropores were compared and analyzed. The results show that infiltration is more likely to happen in MF than UF (Unsaturated Flow, without macropore). The depths of the saturation zone and the wetting front of MF are obviously greater than those of UF. When SR occurs, rainfall duration has the most significant influence on infiltration. When macropores are considered, the ponding depth is smaller at the beginning of rainfall, while the effects are not obvious in the later period. Rain intensity and roughness coefficient have significant influences on the ponding depth. Therefore, macropores should not be ignored in the analysis of the slope seepage field.