Numerical simulation of infiltration, evaporation and shallow groundwater levels with the Richards equation

Abstract Analysis of water and solute movement in unsaturated–saturated soil systems would greatly benefit from an accurate and efficient numerical solution of the Richards equation. Recently the mass balance problem has been solved by proper evaluation of the water capacity term. However, the Darcy fluxes as calculated by various numerical schemes still deviate significantly due to differences in nodal spacing and spatial averaging of the hydraulic conductivity K. This paper discusses a versatile, implicit, backward finite difference scheme which is relatively easy to implement. Special attention is given to the selection of a head or flux controlled top boundary condition during the iterative solution of the Richards equation. The stability of the scheme is shown for extreme events of infiltration, evaporation and rapidly fluctuating, shallow groundwater levels in case of two strongly non-linear soils. For nodal distances of 5 cm, arithmetic means of K overestimate the soil water fluxes, while geometric means of K underestimate these fluxes. At smaller nodal distances, arithmetic means of K converge faster to the theoretical solution than geometric means. In case of nodal distances of 1 cm and arithmetic averages of K, errors due to numerical discretization are small compared to errors due to hysteresis and horizontal spatial variability of the soil hydraulic functions.

[1]  David Russo,et al.  Estimation of finite difference interblock conductivities for simulation of infiltration into initially dry soils , 1992 .

[2]  Donald L. Suarez,et al.  Two-dimensional transport model for variably saturated porous media with major ion chemistry , 1994 .

[3]  R. Feddes,et al.  Simulation of field water use and crop yield , 1978 .

[4]  D. L. Baker.  Darcian Weighted Interblock Conductivity Means for Vertical Unsaturated Flow , 1995 .

[5]  Bruno Brunone,et al.  Numerical analysis of one-dimensional unsaturated flow in layered soils , 1998 .

[6]  J.J.B. Bronswijk,et al.  Modelling soil water dynamics in the unsaturated zone — State of the art , 1988 .

[7]  M. Th. van Genuchten,et al.  A comparison of numerical solutions of the one-dimensional unsaturated—saturated flow and mass transport equations , 1982 .

[8]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[9]  U. Hornung,et al.  Truncation errors in the numerical solution of horizontal diffusion in saturated/unsaturated media , 1983 .

[10]  Reinder A. Feddes,et al.  Simulation model of the water balance of a cropped soil: SWATRE , 1983 .

[11]  A. W. Warrick,et al.  Numerical approximations of darcian flow through unsaturated soil , 1991 .

[12]  Cass T. Miller,et al.  Robust solution of Richards' equation for nonuniform porous media , 1998 .

[13]  Reinder A. Feddes,et al.  SWACROP: A water management and crop production simulation model. , 1992 .

[14]  M. Th. van Genuchten,et al.  HYSWASOR — Simulation Model of Hysteretic Water and Solute Transport in the Root Zone , 1993 .

[15]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[16]  P. J. Ross,et al.  Efficient numerical methods for infiltration using Richards' equation , 1990 .

[17]  Nicholas Kouwen,et al.  Hysteretic effects on net infiltration , 1983 .

[18]  E. E. Miller,et al.  Physical Theory for Capillary Flow Phenomena , 1956 .

[19]  Lehua Pan,et al.  A Transformed Pressure Head‐Based Approach to Solve Richards' Equation for Variably Saturated Soils , 1995 .

[20]  Klaas P. Groen.  Pesticide leaching in polders : field and model studies on cracked clays and loamy sand , 1997 .

[21]  J.H.M. Wösten,et al.  The use of pedotransfer in soil hydrology research in Europe , 1997 .

[22]  Joop G Kroes,et al.  User's guide of SWAP version 2.0 : Simulation of water flow, solute transport and plant growth in the Soil-Water-Atmosphere-Plant environment , 1997 .

[23]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[24]  R. Feddes,et al.  Salinization and crop transpiration of irrigated fields in Pakistan's Punjab , 1997 .

[25]  Attila Nemes,et al.  Using existing soil data to derive hydraulic parameters for simulation models in environmental studies and in land use planning; final report on the European Union funded project, 1998 , 1998 .

[26]  P. Berg.  Long-term simulation of water movement in soils using mass-conserving procedures , 1999 .

[27]  Jan W. Hopmans,et al.  Stochastic analysis of soil water regime in a watershed , 1989 .

[28]  P. Milly,et al.  A mass-conservative procedure for time-stepping in models of unsaturated flow , 1985 .

[29]  B. Mohanty,et al.  A new convergence criterion for the modified Picard iteration method to solve the variably saturated flow equation , 1996 .

[30]  Robert S. Mansell,et al.  Advances in Infiltration , 1984 .

[31]  Jack C. Parker,et al.  Development and evaluation of closed-form expressions for hysteretic soil hydraulic properties , 1987 .

[32]  Michel Vauclin,et al.  A note on estimating finite difference interblock hydraulic conductivity values for transient unsaturated flow problems , 1979 .

[33]  Lehua Pan,et al.  Finite element methods for modeling water flow in variably saturated porous media: Numerical oscillation and mass‐distributed schemes , 1996 .

[34]  Gedeon Dagan,et al.  Unsaturated flow in spatially variable fields: 2. Application of water flow models to various fields , 1983 .

[35]  Randel Haverkamp,et al.  A Comparison of Numerical Simulation Models For One-Dimensional Infiltration1 , 1977 .

[36]  A. J. Desbarats.  An Interblock Conductivity Scheme for Finite Difference Models of Steady Unsaturated Flow in Heterogeneous Media , 1995 .