To　investigate　the　influence　of　macropore　parameters　on　the　non-uniform　migration　and　stability　of　slopes　under　rainfall,　a　solution　model　was　developed　based　on　the　two-domain　model　and　the　stability　coefficient　field　principle.　This　model　addressed　non-uniform　flow　and　slope　stability　under　rainfall　infiltration.　Using　the　COMSOL　Multiphysics　finite　element　platform,　a　corresponding　model　solving　program　was　created.　The　numerical　results　were　validated　through　indoor　rainfall　tests　on　macropore　soil　columns.　A　comparison　was　made　between　slope　volume　water　content　and　point　stability　coefficient　under　conditions　of　uniform　and　non-uniform　flow.　Subsequently,　the　impact　of　macropore　parameters　(namely,　the　proportion　of　macropore　domain　ωf,　the　ratio　of　water　conductivity　between　macropore　and　matrix　domain　μ,　and　the　macropore　empirical　parameter　rw)　on　slope　seepage　field　and　stability　coefficient　field　was　analyzed.　The　findings　revealed　that　compared　to　scenarios　without　macropores,　considering　macropores　led　to　a　7.7%　increase　in　volume　water　content　in　the　matrix　domain　and　a　5.1%　decrease　in　the　macropore　domain.　Additionally,　infiltration　depth　increased　by　83.3%　and　150.0%,　respectively,　and　the　shallow　instability　area　of　the　slope　expanded　by　3.9%.　Infiltration　depth　decreased　with　an　increase　in　ωf　for　both　the　matrix　and　macropore　domains.　Conversely,　with　an　increase　in　μ,　infiltration　depth　decreased　for　the　matrix　domain　and　increased　for　the　macropore　domain.　There　was　no　significant　relationship　observed　with　the　empirical　parameter　rw.　At　the　end　of　the　rainfall,　volume　water　content　in　the　matrix　domain　peaked,　while　the　macropore　domain　increased　with　higher　values　of　ωf　and　μ,　showing　minimal　impact　from　the　empirical　parameter　rw.　Water　exchange　was　categorized　into　negative　exchange　area,　positive　exchange　area,　and　no　exchange　area　along　the　profile.　The　equilibrium　depth　of　water　exchange　aligned　with　the　change　in　infiltration　depth　of　the　matrix　domain.　Both　the　negative　and　positive　exchange　areas　exhibited　peak　values　that　decreased　with　higher　ωf　and　increased　with　higher　μ　and　rw　values.　Under　varying　parameter　values,　the　slope　experienced　shallow　instability　failures.　Higher　values　of　ωf　and　μ　corresponded　to　deeper　instability　layers　and　lower　point　stability　coefficients,　indicating　that　macropores　were　detrimental　to　slope　stability.　©　2024　Sichuan　University.　All　rights　reserved.