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.　©　2021　by　the　authors.　Licensee　MDPI,　Basel,　Switzerland.