Simulation of groundwater age distributions

The objective of our work is to examine how to simulate the age of groundwater in such a way that it can be compared to actual measurements. We start by showing that computation of kinematic age, the one obtained by tracking water along streamlines, is ill posed in heterogeneous aquifers. This, together with its inability to account for mixing processes, makes it inadequate for comparison with age measurements, which are the result of some averaging of the age distribution in the water sample (the type of averaging depends on the measurement procedure). Therefore we go on to write the equations for the cumulative distribution function of residence time under transient flow conditions. This allows us to derive transient equations for the mean age, as well as for the higher-order moments of its distribution, which generalize previous results by others. These moments can be used for approximating age measurements, which need not be equal to the mean age of the water sample. Using both a synthetic and a real example, we show that mean age is an acceptable estimate of radiometric age measurements in many cases. Including a second-order correction (variance of residence time distribution) always improves results. On the other hand, higher-order approximations converge slowly for old (compared to half-life) waters, to the point that third-order approximations often worsen the results.

[1]  P. Cook,et al.  The importance of considering diffusion when using carbon-14 to estimate groundwater recharge to an unconfined aquifer , 1991 .

[2]  A. Bath,et al.  Palaeoclimatic trends deduced from the hydrochemistry of a triassic sandstone aquifer, United Kingdom , 1979 .

[3]  S. Gorelick,et al.  Multiple‐Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore‐Scale Heterogeneity , 1995 .

[4]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[5]  E. Simpson,et al.  Groundwater residence times and recharge rates using a discrete-state compartment model and 14C data , 1984 .

[6]  Charles F. Harvey,et al.  Temporal Moment‐Generating Equations: Modeling Transport and Mass Transfer in Heterogeneous Aquifers , 1995 .

[7]  P. Maloszewski,et al.  DETERMINING THE TURNOVER TIME OF GROUNDWATER SYSTEMS WITH THE AID OF ENVIRONMENTAL TRACERS 1. Models and Their Applicability , 1982 .

[8]  P. V. Danckwerts Continuous flow systems. Distribution of residence times , 1995 .

[9]  P. Maloszewski,et al.  Influence of matrix diffusion and exchange reactions on radiocarbon ages in fissured carbonate aquifers , 1991 .

[10]  J. Motyka,et al.  Matrix porosity as the most important parameter of fissured rocks for solute transport at large scales , 1994 .

[11]  X. Sanchez‐Vila,et al.  On matrix diffusion: formulations, solution methods and qualitative effects , 1998 .

[12]  E. Sudicky,et al.  Carbon 14 dating of groundwater in confined aquifers: Implications of aquitard diffusion , 1981 .

[13]  W. Kinzelbach,et al.  Continuous Groundwater Velocity Fields and Path Lines in Linear, Bilinear, and Trilinear Finite Elements , 1992 .

[14]  Daniel J. Goode,et al.  Direct Simulation of Groundwater Age , 1996 .

[15]  F. J. Pearson,et al.  The age of groundwater in the chalk of the London Basin , 1976 .

[16]  Kurt Roth,et al.  Transfer Functions and Solute Movement Through Soil: Theory and Applications , 1991 .

[17]  A. Zuber Chapter 1 – MATHEMATICAL MODELS FOR THE INTERPRETATION OF ENVIRONMENTAL RADIOISOTOPES IN GROUNDWATER SYSTEMS , 1986 .

[18]  P. Maloszewski,et al.  Principles and practice of calibration and validation of mathematical models for the interpretation of environmental tracer data in aquifers , 1993 .

[19]  M. E. Campana,et al.  Model-derived estimates of groundwater mean ages, recharge rates, effective porosities and storage in a limestone aquifer , 1985 .

[20]  E. Eriksson The Possible Use of Tritium' for Estimating Groundwater Storage , 1958 .

[21]  M. Llamas,et al.  Ground‐Water Age Distribution in Madrid Basin, Spain , 1982 .

[22]  W. Durner Transfer Functions and Solute Movement through Soil , 1992 .

[23]  S. Davis,et al.  Dating Groundwater: A Short Review , 1982 .

[24]  Charles F. Harvey,et al.  Mapping Hydraulic Conductivity: Sequential Conditioning with Measurements of Solute Arrival Time, Hydraulic Head, and Local Conductivity , 1995 .

[25]  Jesús Carrera,et al.  Coupled estimation of flow and solute transport parameters , 1996 .

[26]  P. V. Danckwerts Continuous flow systems , 1953 .