Modeling multidomain hydraulic properties of shrink‐swell soils

Shrink‐swell soils crack and become compacted as they dry, changing properties such as bulk density and hydraulic conductivity. Multidomain models divide soil into independent realms that allow soil cracks to be incorporated into classical flow and transport models. Incongruously, most applications of multidomain models assume that the porosity distributions, bulk density, and effective saturated hydraulic conductivity of the soil are constant. This study builds on a recently derived soil shrinkage model to develop a new multidomain, dual‐permeability model that can accurately predict variations in soil hydraulic properties due to dynamic changes in crack size and connectivity. The model only requires estimates of soil gravimetric water content and a minimal set of parameters, all of which can be determined using laboratory and/or field measurements. We apply the model to eight clayey soils, and demonstrate its ability to quantify variations in volumetric water content (as can be determined during measurement of a soil water characteristic curve) and transient saturated hydraulic conductivity, K s (as can be measured using infiltration tests). The proposed model is able to capture observed variations in K s of one to more than two orders of magnitude. In contrast, other dual‐permeability models assume that K s is constant, resulting in the potential for large error when predicting water movement through shrink‐swell soils. Overall, the multidomain model presented here successfully quantifies fluctuations in the hydraulic properties of shrink‐swell soil matrices, and are suitable for use in physical flow and transport models based on Darcy's Law, the Richards Equation, and the advection‐dispersion equation.