Two models for horizontally averaged unsaturated flow have been developed from two different approaches in the first (Chen et al., this issue) of these companion papers. In this paper the results from both the spatially horizontally averaged Richards equation (SHARE) model and the averaged Green-Ampt model are compared with the results from a three-dimensional finite difference model of unsaturated flow which is perceived as the reference solution. The results of the averaged Green-Ampt model show very good agreement with the averaged results from the three-dimensional model, while SHARE model results are applicable only when fluctuations in soil parameters are small with respect to their mean values. It is also shown that methods of simple parameter averaging (arithmetic or geometric averages) with the local Richards equation does not yield meaningful results in heterogeneous soils. This study suggests that spatially horizontally averaged simplified models (such as the averaged Green-Ampt model) are attractive alternatives to perturbation models (such as the SHARE model) in heterogeneous fields. Due to their simplicity in formulation, accuracy in predicting average behaviors, and minimal requirement of computer effort, the spatially horizontally averaged simplified models can be easily implemented in large-scale models, such as atmospheric mesoscale models. boundary and initial conditions. The average saturation at each depth was obtained by areal integration of the local scale water saturation over the specified area under the assumption of independent vertical soil columns at field scales. Relative occurrence frequency curves of the spatially varying parameters were required for the solution of the average quantities. The second model developed was the spatially horizontally averaged Richards equation (SHARE) model. This model was developed by using spatial averaging and regular perturbation techniques under the assumption of no source or sink in the study area. The SHARE model was expressed as a system of two coupled one-dimensional partial differential equations in terms of the mean saturation and the cross-covariance of saturation and saturated hydrau-
[1]
G. Destouni.
Prediction uncertainty in solute flux through heterogeneous soil
,
1992
.
[2]
George F. Pinder,et al.
Application of the Digital Computer for Aquifer Evaluation
,
1968
.
[3]
Gedeon Dagan,et al.
Solute Dispersion in Unsaturated Heterogeneous Soil at Field Scale: I. Theory
,
1979
.
[4]
Gedeon Dagan,et al.
Unsaturated flow in spatially variable fields: 1. Derivation of models of infiltration and redistribution
,
1983
.
[5]
S. Planton,et al.
A Simple Parameterization of Land Surface Processes for Meteorological Models
,
1989
.
[6]
William H. Press,et al.
Numerical recipes
,
1990
.
[7]
David Russo,et al.
Determining soil hydraulic properties by parameter estimation: On the selection of a model for the hydraulic properties
,
1988
.
[8]
David Russo,et al.
Stochastic analysis of simulated vadose zone solute transport in a vertical cross section of heterogeneous soil during nonsteady water flow
,
1991
.
[9]
Rao S. Govindaraju,et al.
Spatial averaging of unsaturated flow equations under infiltration conditions over areally heterogeneous fields: 1. Development of models
,
1994
.