Subsurface flow model for estimating 2D wetting pattern in drip irrigation

Abstract

Numerical model using 2D Richard’s equation in cylindrical polar coordinate system was developed to assess the soil wetting patterns in a drip irrigation system. A fully implicit finite‐difference scheme was used to solve equations with suitable initial and boundary conditions. The developed model was validated with an experimental result available in the literature, which shows a high level of concordance with low values of RMSE and R2 greater than 0.93 for all discharge rates. Furthermore, the versatility and applicability of the model were demonstrated for complex subsurface conditions considering homogeneous and heterogeneous soil profiles. Water ponding on the surface was observed in clay loam soil, contrasting with deep percolation in sandy loam soil. This underscores the significance of adjusting both discharge rate and irrigation frequency in accordance with prevailing soil conditions.