Preferential flow and soil strength degradation induced by desiccation cracks are important causes for expansive clay slope instability. The cyclic opening and closing of desiccation cracks during drying-wetting processes incessantly alters preferential flow paths and soil strength. Quantify the impact of desiccation crack dynamics on slope hydrology and stability remains a major unresolved challenge. To bridge this gap, we developed the first slope-scale hydro-mechanical model that couples weather-driven crack evolution with preferential flow while incorporating the deterioration effect on soil strength. This unified approach is a major contribution to our capacity to model the integration of hydrological processes and mechanical degradation of soil strength induced by dynamic cracks. The hydrological part adopted a dynamic dual-permeability model (dynamic DPM) and was validated by a physical slope model test. The dynamic DPM was then integrated into a set of numerical slope stability analyses under one-year atmospheric conditions. The groundwater level, water balance, pore water distribution, crack evolution and slope stability were investigated in the case of dynamic cracks and fixed cracks. The hydrological results showed that the slope model with dynamic cracks retained more water and higher groundwater level than that with fixed cracks. The narrowing of desiccation cracks slows down slope drainage process, resulting in a rapid build-up of pore water pressure due to preferential flow, which emerges as an often overlooked and significant factor contributing to slope instability. Conversely, fixed and well-connected cracks in soils enhance water drainage and thus benefit slope stability. The mechanical results revealed that the irreversible deterioration effect induced by crack dynamics on soil strength persistently degrades long-term slope stability. These findings provide new insights into failure mechanisms in cracked soil slopes, and show the importance of the integration of dynamic crack properties into climate-resilient slope design. Also, our results underscore the importance of understanding and quantifying the physical behavior of soil structures for soil hydrological response and slope stability assessment.

Preferential flow induced by desiccation cracks (PF-DC) has been proven to be an important hydrological effect that could cause various geotechnical engineering and ecological environment problems. Investigation on the PF-DC remains a great challenge due to the soil shrinking–swelling behavior. This work presents an experimental and numerical study of the PF-DC considering the dynamic changes of desiccation cracks. A soil column test was conducted under wetting–drying cycles to investigate the dynamic changes of desiccation cracks and their hydrological response. The ratios between the crack area and soil matrix area (crack ratio), crack aperture and depth were measured. The soil water content, matrix suction and water drainage were monitored. A new dynamic dual-permeability preferential flow model (DPMDy) was developed, which includes physically consistent functions in describing the variation of both porosity and hydraulic conductivity in crack and matrix domains. Its performance was compared to the single-domain model (SDM) and rigid dual-permeability model (DPM) with fixed crack ratio and hydraulic conductivity. The experimental results showed that the maximum crack ratio and aperture decreased when the evaporation intensity was excessively raised. The self-closure phenomenon of cracks and increased surficial water content was observed during low-evaporation periods. The simulation results showed that the matrix evaporation modeled by the DPMDy is lower than that of the SDM and DPM, but its crack evaporation is the highest. Compared to the DPM, the DPMDy simulated a faster pressure head building-up process in the crack domain and higher water exchange rates from the crack to the matrix domain during rainfall. Using a fixed crack ratio in the DPM, whether it is the maximum or the average value from the experiment data, will overestimate the infiltration fluxes of PF-DC but underestimate its contribution to the matrix domain. In conclusion, the DPMDy better described the underlying physics involving crack evolution and hydrological response with respect to the SDM and DPM. Further improvement of the DPMDy should focus on the hysteresis effect of the soil water retention curve and soil deformation during wetting–drying cycles.

Quantitative investigation on the preferential flow induced by desiccation cracks (PF-DC) remains a great challenge due to the soil shrinking-swelling behavior. This work presents a series of comparative numerical studies to investigate the accuracy and substitutability of different models in simulating the water flux, hydrological response and crack evolution induced by PF-DC. As a comparative study, an effective dynamic dual-permeability model (DDPM) we recently developed and validated was regarded as a benchmark model. Three numerical experiments were conducted to (i) compare the difference among the single-domain model (SDM), rigid dual-permeability model (RDPM) and DDPM; (ii) test the sensitivity of the DDPM to the shrinking-swelling parameters; (iii) test the rationality of a “lighter” dynamic DPM (LDPM) only considering the proportion changes of each domain while neglecting the variation of hydraulic properties. The results showed that compared to the DDPM, the SDM overestimated the water content under low-rainfall intensity while underestimating the water content under high-intensity rainfall and failed to capture the early increase of water content in deep soils induced by PF-DC. The RDPM greatly overestimated the total water content and water storage capacity of the crack domain, which was not suggested to be used in the surface runoff or flood forecast. The DDPM is overall not sensitive to the shrinking-swelling parameters, indicative of relatively loose accuracy requirements in measuring the soil shrinking-swelling parameters. The LDPM can be a tentative alternative option for the DDPM, but it is better not to use it to evaluate the surface runoff or use it under long-term extreme drought. In conclusion, the prediction errors without considering crack evolution and variation of hydraulic properties of each domain (RDPM) are the highest, then followed by the only considering crack evolution (LDPM) and uncertainties of shrinking-swelling parameters.