The Regional Municipality of Waterloo in Ontario, Canada (population 250,000) depends on ground water for most of its water supply. The ground water is extracted from the Waterloo Moraine, an extensive and complex glacial aquifer system extending over a 400 km2 area. A methodology is being developed to inventory the ground water resource, to define its susceptibility to contamination, and to create the basis for optimal management and protection strategies. A key component of this methodology is a three-dimensional conceptual hydrogeologic model based on the geologic characteristics of the multiple-aquifer Moraine system. The steps in the development of the model include screening of the large database, interpretation and interpolation of the data to define the variable hydrostatigraphy and to generate consistent hydraulic conductivity functions, and model calibration. The numerical basis is a fully three-dimensional finite element model that provides flexibility and adaptability in representing the natural boundaries, the highly irregular stratigraphy, and the numerous wellfields. The model has the capability to automatically find the location of the water table consistent with given recharge and pumping conditions, and to direct recharge from low-permeability areas to higher-permeability areas. Capture zones generated by the model are found to be highly sensitive with respect to the geologic structure, in particular the presence or absence of windows in the aquitard units. Professional judgment is found to be an essential component of the modeling process.
[1]
Emil O. Frind,et al.
Thermal energy storage in an unconfined aquifer: 2. Model development, validation, and application
,
1992
.
[2]
Wolfgang Kinzelbach,et al.
Determination of groundwater catchment areas in two and three spatial dimensions
,
1992
.
[3]
Modeling Ground‐Water Remediation in the High Plains Aquifer
,
1994
.
[4]
E. Sudicky,et al.
A Noniterative Technique for the Direct Implementation of Well Bore Boundary Conditions in Three‐Dimensional Heterogeneous Formations
,
1995
.
[5]
Fred J. Molz,et al.
An Evaluation of Interpolation Methodologies for Generating Three-Dimensional Hydraulic Property Distributions from Measured Data
,
1995
.
[6]
Graham E. Fogg,et al.
Groundwater Flow and Sand Body Interconnectedness in a Thick, Multiple-Aquifer System
,
1986
.
[7]
E. Poeter,et al.
Inverse Models: A Necessary Next Step in Ground‐Water Modeling
,
1997
.
[8]
M. Anderson.
Hydrogeologic facies models to delineate large-scale spatial trends in glacial and glaciofluvial sediments
,
1989
.
[9]
D. W. Pollock.
Semianalytical Computation of Path Lines for Finite‐Difference Models
,
1988
.
[10]
David C. Schafer,et al.
Determining 3D Capture Zones in Homogeneous, Anisotropic Aquifers
,
1996
.
[11]
E. Sudicky.
A natural gradient experiment on solute transport in a sand aquifer: Spatial variability of hydraulic conductivity and its role in the dispersion process
,
1986
.