Numerical Simulation on the Response Mechanism of Soil Water Migration to Mining Subsidence Cracks

Mining-induced subsidence has significantly altered the structure of the vadose zone in coal mining areas, where soil cracks act as preferential pathways controlling water infiltration and redistribution. In this study, a Hydrus-2D dual-domain seepage model incorporating geometric parameterization of cracks was developed to simulate water migration in the vadose zone of a typical subsidence area in the Ordos Basin. The model integrates field-measured crack geometry, soil texture, and rainfall characteristics to quantitatively analyze preferential flow formation under twelve combinations of crack width, soil type, and rainfall intensity. The results show that (i) crack width dominates preferential flow behavior, with wider cracks (≥5 cm) deepening the wetting front from approximately 107 cm to 144 cm within 120 h and sustaining high conductivity after rainfall; (ii) soil texture governs infiltration pathways, as sandy soils promote deeper wetting fronts (up to 99 cm, ~40% deeper than loam) and layered soils induce interface retention or “jump” infiltration; and (iii) rainfall intensity controls infiltration depth, with storm events producing wetting fronts more than four times deeper than those under light rain. Overall, this study demonstrates the feasibility and significance of integrating crack parameterization into vadose-zone hydrological modeling using Hydrus-2D, providing a quantitative basis for understanding rapid infiltration–migration–recharge processes and supporting ecological restoration and water resource management in arid and semi-arid mining regions.