Model predicts the impact of root system architecture on soil water infiltration.
There is strong experimental evidence that root systems substantially change the saturated hydraulic conductivity of soil. However, the mechanisms by which roots affect soil hydraulic properties remain largely unknown. In this work, we made the hypothesis that preferential soil moisture transport occurs along the axes of roots, and that this is what changes a soil's saturated hydraulic conductivity. We modified Richards' equation to incorporate the preferential flow of soil moisture along the axes of roots. Using the finite element method and Bayesian optimisation, we developed a pipeline to calibrate our model with respect to a given root system. When applied to simulated root systems, the pipeline successfully predicted the pore-water pressure profiles corresponding to saturated hydraulic conductivity values, observed by Leung et al. (2018), for soils vegetated by willow and grass. Prediction accuracy improved for root systems with more realistic architectures, therefore suggesting that changes in saturated hydraulic conductivity are a result of roots enabling preferential soil moisture transport along their axes. The model proposed in this work improves our ability to predict moisture transport through vegetated soil and could help optimise irrigation, forecast flood events and plan landslide prevention strategies.