Dimensionality and scales of preferential flow in soils of Shale Hills hillslope simulated using HYDRUS

Preferential flow (PF) processes are governed by subsurface soil structures at various scales. Still, model validation and mechanistic understanding of PF are very lacking. We hypothesize that PF at hillslope and larger scales cannot be described and quantified when neglecting small‐scaled spatially variable processes and simplifying the model dimensionality. The objective was to learn from comparing simulation results of multidimensional (1D, 2D, and 3D) and multiscale (pedon, catena, and catchment) modeling approaches with comprehensive datasets, and so as to evaluate PF simulations based on the Richards’ equation (solved by the HYDRUS software). Results showed limited alignment between 1D simulations and soil moisture data, mainly affected by vertical changes in porosity, permeability, and precipitation features. 2D and 3D simulations outperformed 1D models. 3D simulations provided satisfactory description of PF dynamics at the pedon scale, considering accurate representations of soil and bedrock structures for three dimensions (vertical, horizontal, and surrounding area). In 2D simulations at the pedon scale, models incorporating dual‐porosity and anisotropy of soils yielded more accurate predictions of water dynamics than single‐porosity and isotropic models. Furthermore, the application of 2D simulation at the catena scale identify PF pathways owing to the enhanced representation of the hydraulic connectivity between different locations along the slope. The results confirmed the significance of multidimensional and multiscale modeling approaches for PF simulations in hillslope hydrology. Considering the complexity and parameterization of 2D and 3D “bottom‐up” physically based models in representing spatial variability within and between soil profiles and/or underlying bedrock geology, the results contribute to creating a modeling framework applicable to identify the PF processes and thus their implications in managing water resources. Core Ideas The 3D model was better than 2D model to detect the preferential flow. A dual‐porosity or anisotropy model produced more accurate predictions than a single‐porosity or isotropy model. A model domain that considered fractured bedrock performed better than without it.