Appropriate vertical discretization of Richards' equation for two‐dimensional watershed‐scale modelling

A number of watershed-scale hydrological models include Richards' equation (RE) solutions, but the literature is sparse on information as to the appropriate application of RE at the watershed scale. In most published applications of RE in distributed watershed-scale hydrological modelling, coarse vertical resolutions are used to decrease the computational burden. Compared to point- or field-scale studies, application at the watershed scale is complicated by diverse runoff production mechanisms, groundwater effects on runoff production, runon phenomena and heterogeneous watershed characteristics. An essential element of the numerical solution of RE is that the solution converges as the spatial resolution increases. Spatial convergence studies can be used to identify the proper resolution that accurately describes the solution with maximum computational efficiency, when using physically realistic parameter values. In this study, spatial convergence studies are conducted using the two-dimensional, distributed-parameter, gridded surface subsurface hydrological analysis (GSSHA) model, which solves RE to simulate vadose zone fluxes. Tests to determine if the required discretization is strongly a function of dominant runoff production mechanism are conducted using data from two very different watersheds, the Hortonian Goodwin Creek Experimental Watershed and the non-Hortonian Muddy Brook watershed. Total infiltration, stream flow and evapotranspiration for the entire simulation period are used to compute comparison statistics. The influences of upper and lower boundary conditions on the solution accuracy are also explored. Results indicate that to simulate hydrological fluxes accurately at both watersheds small vertical cell sizes, of the order of 1 cm, are required near the soil surface, but not throughout the soil column. The appropriate choice of approximations for calculating the near soil-surface unsaturated hydraulic conductivity can yield modest increases in the required cell size. Results for both watersheds are quite similar, even though the soils and runoff production mechanisms differ greatly between the two catchments. Copyright © 2003 John Wiley & Sons, Ltd.

[1]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[2]  C. Downer.  Identification and modeling of important stream flow producing processes in watersheds , 2002 .

[3]  Eric F. Wood,et al.  NUMERICAL EVALUATION OF ITERATIVE AND NONITERATIVE METHODS FOR THE SOLUTION OF THE NONLINEAR RICHARDS EQUATION , 1991 .

[4]  S. T. Chu,et al.  INFILTRATION MODEL FOR SOIL PROFILES WITH A WATER TABLE , 1997 .

[5]  J. C. van Dam,et al.  Numerical simulation of infiltration, evaporation and shallow groundwater levels with the Richards equation , 2000 .

[6]  F. Melone,et al.  A unified model for infiltration and redistribution during complex rainfall patterns , 1997 .

[7]  D. L. Brakensiek,et al.  Estimation of Soil Water Properties , 1982 .

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

[9]  J. Crank.  Difference Methods for Initial-Value Problems 2nd edn , 1968 .

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

[11]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[12]  R. G. Hills,et al.  Algorithms for solving Richards' equation for variably saturated soils , 1992 .

[13]  P. E. O'connell,et al.  An introduction to the European Hydrological System — Systeme Hydrologique Europeen, “SHE”, 2: Structure of a physically-based, distributed modelling system , 1986 .

[14]  S. Sorooshian,et al.  Effective and efficient global optimization for conceptual rainfall‐runoff models , 1992 .

[15]  R. Horton.  The Rôle of infiltration in the hydrologic cycle , 1933 .

[16]  Florisa Melone,et al.  Modeling infiltration for multistorm runoff events , 1993 .

[17]  J. Monteith.  Evaporation and environment. , 1965, Symposia of the Society for Experimental Biology.

[18]  Glenn E. Moglen,et al.  Resolution Effects on Hydrologic Modeling Parameters and Peak Discharge , 2001 .

[19]  Stanley G. Benjamin,et al.  Performance of Different Soil Model Configurations in Simulating Ground Surface Temperature and Surface Fluxes , 1997 .

[20]  J. Clausen,et al.  Movement of Nitrogen through an Agricultural Riparian Zone: 1. Field Studies , 1993 .

[21]  P. E. O'connell,et al.  An introduction to the European Hydrological System — Systeme Hydrologique Europeen, “SHE”, 1: History and philosophy of a physically-based, distributed modelling system , 1986 .

[22]  G. Ling,et al.  The Courant and Peclet Number criteria for the numerical solution of the Richards Equation , 1993 .

[23]  Hatim O. Sharif,et al.  On the calibration and verification of two‐dimensional, distributed, Hortonian, continuous watershed models , 2000 .

[24]  A. Hjelmfelt,et al.  DEM Based Overland Flow Routing Model , 1998 .

[25]  D. Short,et al.  The practicability of using Richards' equation for general purpose soil-water dynamics models , 1995 .

[26]  J. Hutson,et al.  A retentivity function for use in soil–water simulation models , 1987 .

[27]  A. W. Benton,et al.  Soil survey of Tolland County, Connecticut. , 1966 .

[28]  J. Mullaney.  Hydrogeology And Ground-Water Flow At The Muddy Brook Riparian Zone, North-Central Connecticut , 1995 .

[29]  Fred L. Ogden,et al.  Prediction of runoff and soil moistures at the watershed scale: Effects of model complexity and parameter assignment , 2003 .

[30]  Fred L. Ogden,et al.  Green and Ampt Infiltration with Redistribution , 1997 .

[31]  Vijay P. Singh,et al.  CASC2D: a two-dimensional, physically-based, Hortonian hydrologic model. , 2002 .

[32]  David R. Miller,et al.  Deposition rates for sulfur and nitrogen to a hardwood forest in Northern Connecticut, U.S.A. , 1994 .

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

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

[35]  W. Sloan,et al.  UP Modelling System for large scale hydrology: deriving large-scale physically-based parameters for the Arkansas-Red River basin , 1999 .