WEC-C: a distributed, deterministic catchment model -- theory, formulation and testing

WEC-C is a distributed, deterministic, catchment-scale water flow and solute transport model containing a number of innovations not present in comparable models. For example, it allows for the imposition of catchment topological changes resulting from land uses such as surface mining. In addition. it features preferential vertical flow in the vadose zone modelled using a dual continuum approach. Numerically, it differs from many existing models in that it employs a direct linkage between the vertical and lateral solvers of its split- solver scheme and, due to its use of explicit solvers, is stable regardless of the form of the soil moisture characteristic curves. The WEC-C model framework is a rectangular grid of uniform cell size in the lateral plane combined with a system of soil layers, of variable thickness, in the vertical direction. Within this structure, the governing equations for flow and transport are discretised and solved. All parameters are defined locally in each computational cell so that all available data on catchment variability can be incorporated directly into the model. This paper describes the formulation of WEC-C, which is based on operator splitting with first- order accurate solvers for both the vertical and lateral flow and transport models. WEC-C was subjected to four stringent, synthetic tests designed to evaluate its accuracy by comparison with available analytical and numerical solutions. These showed that there were scale issues associated with the model, and that induced numeric dispersion of solutes could be significant. Nonetheless, it is suggested that WEC-C is still useful as a distributed catchment model.

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

[2]  Philip Broadbridge,et al.  Constant rate rainfall infiltration: A versatile nonlinear model, I. Analytic solution , 1988 .

[3]  Keith Beven,et al.  A Discussion of Distributed Hydrological Modelling , 1990 .

[4]  D. A. Barry,et al.  Exact Solutions for Water Infiltration With an Arbitrary Surface Flux or Nonlinear Solute Adsorption , 1991 .

[5]  J. Parlange,et al.  Free-surface flow in porous media and periodic solution of the shallow-flow approximation , 1984 .

[6]  R. Noyé,et al.  Numerical Solutions of Partial Differential Equations , 1983 .

[7]  R. R. van der Ploeg,et al.  A Numerical Study of the Effects of Noncapillary-Sized Pores Upon Infiltration 1 , 1979 .

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

[9]  Philip Broadbridge,et al.  Reply [to Comment on Constant rate rainfall infiltration: A versatile nonlinear model: 1. Analytic , 1988 .

[10]  Y. Tsang,et al.  On thermohydrologic conditions near high‐level nuclear wastes emplaced in partially saturated fractured tuff: 1. Simulation studies with explicit consideration of fracture effects , 1990 .

[11]  M. Sharma,et al.  Areal distribution of infiltration parameters and some soil physical properties in lateritic catchments , 1987 .

[12]  D. A. Barry,et al.  Capillarity correction to periodic solutions of the shallow flow approximation , 2013 .

[13]  Philip Broadbridge,et al.  Reply [to “Comment on ‘Constant rate rainfall infiltration: A versatile nonlinear model: 2. Applications of solutions’ by I. White and P. Broadbridge”] , 1989 .

[14]  D. A. Barry,et al.  Note on common mixing cell models , 1994 .

[15]  Malcolm R. Davidson Numerical calculation of saturated-unsaturated infiltration in a cracked soil , 1985 .

[16]  J.-Y. Parlange,et al.  Comment on “Constant rate rainfall infiltration: A versatile nonlinear model: 1. Analytic solution” by P. Broadbridge and I. White , 1988 .

[17]  Arnold Verruijt,et al.  Modeling Groundwater Pollution , 1987 .

[18]  James O. Duguid,et al.  Flow in fractured porous media , 1977 .

[19]  E. O'Loughlin Prediction of Surface Saturation Zones in Natural Catchments by Topographic Analysis , 1986 .

[20]  V. Singh,et al.  Computer Models of Watershed Hydrology , 1995 .

[21]  Philip Broadbridge,et al.  Constant rate rainfall infiltration: A versatile nonlinear model, II. Applications of solutions , 1988 .

[22]  Gedeon Dagan,et al.  SECOND-ORDER THEORY OF SHALLOW FREE-SURFACE FLOW IN POROUS MEDIA , 1967 .

[23]  G. Strang On the Construction and Comparison of Difference Schemes , 1968 .

[24]  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 .

[25]  J.-Y. Parlange,et al.  Exact nonlinear solution for constant flux infiltration , 1988 .

[26]  F. D. Whisler,et al.  Analysis of radial, steady-state solution, and solute flow. , 1970 .

[27]  L. Lapidus,et al.  Mathematics of Adsorption in Beds. VI. The Effect of Longitudinal Diffusion in Ion Exchange and Chromatographic Columns , 1952 .

[28]  M. V. Genuchten,et al.  A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media , 1993 .

[29]  Keith Beven,et al.  Distributed Hydrological Modelling , 1998 .

[30]  Cass T. Miller,et al.  Temporal discretisation errors in non-iterative split-operator approaches to solving chemical reaction/groundwater transport models , 1996 .

[31]  Thomas A. McMahon,et al.  Physically based hydrologic modeling: 2. Is the concept realistic? , 1992 .

[32]  Cass T. Miller,et al.  Development of split‐operator, Petrov‐Galerkin Methods to simulate transport and diffusion problems , 1993 .

[33]  D. A. Barry,et al.  A class of exact solutions for Richards' equation , 1993 .

[34]  Roger Moore Computer simulation techniques in hydrology: By George Fleming , 1976 .

[35]  W. G. Gray,et al.  Finite Element Simulation in Surface and Subsurface Hydrology , 1977 .

[36]  K. Terzaghi,et al.  Soil mechanics in engineering practice , 1948 .

[37]  Jens Christian Refsgaard,et al.  Application of the SHE to catchments in India Part 2. Field experiments and simulation studies with the SHE on the Kolar subcatchment of the Narmada River , 1992 .