Probabilistic Exposure Risk Assessment with Advective-Dispersive Well Vulnerability Criteria

Abstract Time-related advection-based well-head protection zones are commonly used to manage the contamination risk of drinking water wells. According to current water safety plans advanced risk management schemes are needed to better control and monitor all possible hazards within catchments. The goal of this work is to cast the four advective–dispersive intrinsic well vulnerability criteria by Frind et al. [1] into a framework of probabilistic risk assessment framework. These criteria are: (i) arrival time, (ii) level of peak concentration, (iii) time until first arrival of critical concentrations and (iv) exposure time. Our probabilistic framework yields catchment-wide maps of probabilities to not comply with these criteria. This provides indispensable information for catchment managers to perform probabilistic exposure risk assessment and thus improves the basis for risk-informed well-head management. We resolve heterogeneity with high-resolution Monte Carlo simulations and use a new reverse formulation of temporal moment transport equations to keep computational costs low. Our method is independent of dimensionality and boundary conditions, and can account for arbitrary sources of uncertainty. It can be coupled with any method for conditioning on available data. For simplicity, we demonstrate the concept on a 2D example that includes conditioning on synthetic data.

[1]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

[2]  Ahmed E. Hassan,et al.  Evaluation of analytical solute discharge moments using numerical modeling in absolute and relative dispersion frameworks , 2002 .

[3]  P. Kitanidis Parameter Uncertainty in Estimation of Spatial Functions: Bayesian Analysis , 1986 .

[4]  Y. Rubin,et al.  Bayesian geostatistical design: Task‐driven optimal site investigation when the geostatistical model is uncertain , 2010 .

[5]  J. Shao,et al.  The jackknife and bootstrap , 1996 .

[6]  Roseanna M. Neupauer,et al.  Adjoint method for obtaining backward‐in‐time location and travel time probabilities of a conservative groundwater contaminant , 1999 .

[7]  R. Scheaffer,et al.  Mathematical Statistics with Applications. , 1992 .

[8]  Alberto Guadagnini,et al.  Time‐Related Capture Zones for Contaminants in Randomly Heterogeneous Formations , 1999 .

[9]  D. W. Pollock Semianalytical Computation of Path Lines for Finite‐Difference Models , 1988 .

[10]  Fabien Cornaton,et al.  Deterministic models of groundwater age, life expectancy and transit time distributions in advective-dispersive systems , 2003 .

[11]  David W. Pollock,et al.  Sources of Water to Wells for Transient Cyclic Systems , 1996 .

[12]  S. P. Neuman,et al.  Maximum likelihood Bayesian averaging of uncertain model predictions , 2003 .

[13]  W. Nowak,et al.  A modified Levenberg-Marquardt algorithm for quasi-linear geostatistical inversing , 2004 .

[14]  P. Kitanidis,et al.  Geostatistical inversing for large-contrast transmissivity fields , 2009 .

[15]  Matthijs van Leeuwen,et al.  Stochastic determination of well capture zones , 1998 .

[16]  A. Scheidegger General Theory of Dispersion in Porous Media , 1961 .

[17]  W. Nowak,et al.  Application of FFT-based Algorithms for Large-Scale Universal Kriging Problems , 2009 .

[18]  M. Stein,et al.  A Bayesian analysis of kriging , 1993 .

[19]  Peter J. Diggle,et al.  Bayesian methodology to stochastic capture zone determination: Conditioning on transmissivity measurements , 2002 .

[20]  Albert J. Valocchi,et al.  Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state , 2007 .

[21]  T. Ulrych,et al.  A full‐Bayesian approach to the groundwater inverse problem for steady state flow , 2000 .

[22]  A. Mohammad-Djafari A Matlab Program to Calculate the Maximum Entropy Distributions , 2001, physics/0111126.

[23]  David J. Cushman,et al.  Risk assessment for environmental contamination: an overview of the fundamentals and application of risk assessment at contaminated sites , 2001 .

[24]  E. Frind,et al.  Well Vulnerability: A Quantitative Approach for Source Water Protection , 2006, Ground water.

[25]  M. Willmann,et al.  Impact of log-transmissivity variogram structure on groundwater flow and transport predictions , 2009 .

[26]  Adrian P. Butler,et al.  Worth of head data in well-capture zone design: deterministic and stochastic analysis , 2004 .

[27]  R. Neupauer,et al.  Adjoint‐derived location and travel time probabilities for a multidimensional groundwater system , 2001 .

[28]  Vladimir Cvetkovic,et al.  Relative dispersion for solute flux in aquifers , 1998, Journal of Fluid Mechanics.

[29]  Felipe P. J. de Barros,et al.  A divide and conquer approach to cope with uncertainty, human health risk, and decision making in contaminant hydrology , 2011 .

[30]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[31]  Wolfgang Nowak,et al.  Probability density functions of hydraulic head and velocity in three‐dimensional heterogeneous porous media , 2008 .

[32]  P. Kitanidis,et al.  Characterization of mixing and dilution in heterogeneous aquifers by means of local temporal moments , 2000 .

[33]  U. Epa,et al.  Guiding Principles for Monte Carlo Analysis , 1997 .

[34]  J. M. Shafer,et al.  Assessment of Uncertainty in Time‐Related Capture Zones Using Conditional Simulation of Hydraulic Conductivity , 1991 .

[35]  Tomas Öberg,et al.  A Review of Probabilistic Risk Assessment of Contaminated Land (12 pp) , 2005 .

[36]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[37]  Sabine Attinger,et al.  Temporal behavior of a solute cloud in a heterogeneous porous medium: 1. Point‐like injection , 2000 .

[38]  E. Jaynes Information Theory and Statistical Mechanics , 1957 .

[39]  Keith Beven,et al.  Stochastic capture zone delineation within the generalized likelihood uncertainty estimation methodology: Conditioning on head observations , 2001 .

[40]  P. Kitanidis Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[41]  J. Wishart,et al.  Methods of Statistical Analysis , 1954 .

[42]  Peter K. Kitanidis,et al.  The concept of the Dilution Index , 1994 .

[43]  Vassilios A. Tsihrintzis,et al.  Delineation of groundwater protection zones by the backward particle tracking method: theoretical background and GIS-based stochastic analysis , 2008 .

[44]  Lorna Fewtrell,et al.  Water safety plans: managing drinking-water quality from catchment to consumer , 2005 .

[45]  Stephen E. Silliman,et al.  Utility of Simple Models for Capture Zone Delineation in Heterogeneous Unconfined Aquifers , 2000 .

[46]  Peter J. Diggle,et al.  Bayesian Geostatistical Design , 2006 .

[47]  W. Kinzelbach,et al.  Determination of a well head protection zone by stochastic inverse modelling , 1998 .

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

[49]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[50]  Harrie-Jan Hendricks Franssen,et al.  Joint estimation of transmissivities and recharges—application: stochastic characterization of well capture zones , 2004 .

[51]  Y. Rubin Applied Stochastic Hydrogeology , 2003 .

[52]  Keith Beven,et al.  Bayesian methodology for stochastic capture zone delineation incorporating transmissivity measurements and hydraulic head observations , 2003 .

[53]  T. Ulrych,et al.  Minimum relative entropy and probabilistic inversion in groundwater hydrology , 1998 .

[54]  Alberto Guadagnini,et al.  Probabilistic estimation of well catchments in heterogeneous aquifers , 1996 .

[55]  David N. Lerner,et al.  How Uncertain Is Our Estimate of a Wellhead Protection Zone? , 1998 .

[56]  Andres Alcolea,et al.  Pilot points method incorporating prior information for solving the groundwater flow inverse problem , 2006 .

[57]  Terje Aven,et al.  Some reflections on uncertainty analysis and management , 2010, Reliab. Eng. Syst. Saf..

[58]  P. Bickel,et al.  Obstacles to High-Dimensional Particle Filtering , 2008 .

[59]  G. Dagan Solute transport in heterogeneous porous formations , 1984, Journal of Fluid Mechanics.

[60]  G. Meek Mathematical statistics with applications , 1973 .

[61]  John Skilling,et al.  Maximum Entropy and Bayesian Methods , 1989 .

[62]  Eulogio Pardo-Igúzquiza,et al.  Bayesian Inference of Spatial Covariance Parameters , 1999 .

[63]  Yoram Rubin,et al.  A risk‐driven approach for subsurface site characterization , 2008 .

[64]  Rajandrea Sethi,et al.  An automatic, stagnation point based algorithm for the delineation of Wellhead Protection Areas , 2008 .

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

[66]  Ramana V. Grandhi,et al.  A Bayesian approach for quantification of model uncertainty , 2010, Reliab. Eng. Syst. Saf..

[67]  Jennifer A Hoeting,et al.  Model selection for geostatistical models. , 2006, Ecological applications : a publication of the Ecological Society of America.

[68]  Sabine Attinger,et al.  Temporal behavior of a solute cloud in a heterogeneous porous medium: 2. Spatially extended injection , 2000 .

[69]  R. Baierlein Probability Theory: The Logic of Science , 2004 .

[70]  H. H. Franssen,et al.  A comparison of seven methods for the inverse modelling of groundwater flow. Application to the characterisation of well catchments , 2009 .

[71]  Roko Andricevic,et al.  Probabilistic capture zone delineation based on an analytic solution. , 2002, Ground water.

[72]  Yoram Rubin,et al.  The concept of comparative information yield curves and its application to risk‐based site characterization , 2009 .

[73]  M. Riva,et al.  Probabilistic study of well capture zones distribution at the Lauswiesen field site. , 2006, Journal of contaminant hydrology.

[74]  Alberto Guadagnini,et al.  Delineation of Source Protection Zones Using Statistical Methods , 2005 .

[75]  Adrian E. Raftery,et al.  Bayesian model averaging: a tutorial (with comments by M. Clyde, David Draper and E. I. George, and a rejoinder by the authors , 1999 .

[76]  Wolfgang Kinzelbach,et al.  Semianalytical uncertainty estimation of well catchments: Conditioning by head and transmissivity data , 2004 .