Adjoint method for obtaining backward‐in‐time location and travel time probabilities of a conservative groundwater contaminant

Backward location and travel time probabilities can be used to determine the prior location of contamination in an aquifer. For a contaminant particle that was detected in an aquifer, the backward location probability is the probability of where the particle was located at some prior time. Backward travel time probability is the probability of when the particle was located at some position upgradient of the detection. These probabilities can be used to improve characterization of known sources of groundwater contamination, to identify previously unknown contamination sources, and to delineate capture zones. For simple model domains, backward probabilities can be obtained heuristically from a forward model of contaminant transport. For multidimensional problems and complex domain geometries, the heuristic approach is difficult to implement and verify. The adjoint method provides a formal approach for obtaining backward probabilities for all model domains and geometries. We formally show that the backward model probabilities are adjoint states of resident concentration. We provide a methodology for obtaining the governing equations and boundary and final conditions for these probabilities. The approach is illustrated using a one-dimensional, semi-infinite domain that mimics flow to a production well, and these results are compared to equivalent probabilities derived heuristically.

[1]  T. L. Saaty Modern nonlinear equations , 1981 .

[2]  R. W. Andrews,et al.  Sensitivity Analysis for Steady State Groundwater Flow Using Adjoint Operators , 1985 .

[3]  A. Zuber,et al.  On the physical meaning of the dispersion equation and its solutions for different initial and boundary conditions , 1978 .

[4]  W. Yeh,et al.  Sequential estimation of aquifer parameters , 1988 .

[5]  William W.-G. Yeh,et al.  Coupled inverse problems in groundwater modeling - 1. Sensitivity analysis and parameter identification. , 1990 .

[6]  William W.-G. Yeh,et al.  Variational sensitivity analysis, data requirements, and parameter identification in a leaky aquifer system , 1990 .

[7]  David A. Chin,et al.  Risk Management in Wellhead Protection , 1994 .

[8]  W. Yeh,et al.  Identification of Parameter Structure in Groundwater Inverse Problem , 1985 .

[9]  John M. Mulvey,et al.  Contaminated groundwater remediation design using simulation, optimization, and sensitivity theory: 1. Model development , 1988 .

[10]  John L. Wilson,et al.  Illustration and Verification of Adjoint Sensitivity Theory for Steady State Groundwater Flow , 1985 .

[11]  William A. Jury,et al.  Simulation of solute transport using a transfer function model , 1982 .

[12]  L. Townley,et al.  Computationally Efficient Algorithms for Parameter Estimation and Uncertainty Propagation in Numerical Models of Groundwater Flow , 1985 .

[13]  G. Dagan Flow and transport in porous formations , 1989 .

[14]  S. P. Neuman A statistical approach to the inverse problem of aquifer hydrology: 3. Improved solution method and added perspective , 1980 .

[15]  Chia-Shyun Chen,et al.  Analytical solution for aquifer decontamination by pumping , 1988 .

[16]  M. V. Genuchten,et al.  Flux-Averaged and Volume-Averaged Concentrations in Continuum Approaches to Solute Transport , 1984 .

[17]  A. Bagtzoglou,et al.  Application of particle methods to reliable identification of groundwater pollution sources , 1992 .

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

[19]  T. Skaggs,et al.  Recovering the release history of a groundwater contaminant , 1994 .