Effects of soil heterogeneity on subsurface water movement in agricultural fields: A numerical study

This study aims to numerically examine the effects of soil heterogeneity on subsurface water movement in agricultural fields under wheat (Triticum aestivum L.) and rice (Oryza sativa L.) crop covers. The 3-D Richards equation that accounts for root water uptake was coupled with a random field generation method to produce realizations of soil hydraulic properties, and the effects of soil heterogeneity were analyzed under realistic atmospheric boundary conditions. Four soil types and three different soil heterogeneity levels were considered to investigate the response of soil heterogeneity on soil water dynamics. It was observed that the effect of soil heterogeneity varies across soil types and its effect was more pronounced in silty clay loam soils. Further, an increase in soil heterogeneity level can increase or decrease the mean soil moisture, hydraulic conductivity, and flow velocity based on land cover, upper boundary condition, and atmospheric demand, and its effect was higher during irrigation periods. A maximum difference of 30% in soil moisture values was observed near the surface soil layers due to an increase in soil heterogeneity on irrigation days. The mean soil moisture content under rice crop cover was always lower than that of wheat crop cover due to higher atmospheric demand in the former. Soil layering effect on soil moisture dynamics, which is related to capillary and hydraulic barrier formation, was found to depend on the layering configuration, available soil moisture, and land cover. The numerical experiments provided useful insights into the complex soil moisture dynamics in agricultural fields.