Human health risk assessment of CO2 leakage into overlying aquifers using a stochastic, geochemical reactive transport approach.

Increased human health risk associated with groundwater contamination from potential carbon dioxide (CO2) leakage into a potable aquifer is predicted by conducting a joint uncertainty and variability (JUV) risk assessment. The approach presented here explicitly incorporates heterogeneous flow and geochemical reactive transport in an efficient manner and is used to evaluate how differences in representation of subsurface physical heterogeneity and geochemical reactions change the calculated risk for the same hypothetical aquifer scenario where a CO2 leak induces increased lead (Pb(2+)) concentrations through dissolution of galena (PbS). A nested Monte Carlo approach was used to take Pb(2+) concentrations at a well from an ensemble of numerical reactive transport simulations (uncertainty) and sample within a population of potentially exposed individuals (variability) to calculate risk as a function of both uncertainty and variability. Pb(2+) concentrations at the well were determined with numerical reactive transport simulation ensembles using a streamline technique in a heterogeneous 3D aquifer. Three ensembles with variances of log hydraulic conductivity (σ(2)lnK) of 1, 3.61, and 16 were simulated. Under the conditions simulated, calculated risk is shown to be a function of the strength of subsurface heterogeneity, σ(2)lnK and the choice between calculating Pb(2+) concentrations in groundwater using equilibrium with galena and kinetic mineral reaction rates. Calculated risk increased with an increase in σ(2)lnK of 1 to 3.61, but decreased when σ(2)lnK was increased from 3.61 to 16 for all but the highest percentiles of uncertainty. Using a Pb(2+) concentration in equilibrium with galena under CO2 leakage conditions (PCO2 = 30 bar) resulted in lower estimated risk than the simulations where Pb(2+) concentrations were calculated using kinetic mass transfer reaction rates for galena dissolution and precipitation. This study highlights the importance of understanding both hydrologic and geochemical conditions when numerical simulations are used to perform quantitative risk calculations.

[1]  Martin J. Blunt,et al.  Streamline‐based simulation of carbon dioxide storage in a North Sea aquifer , 2006 .

[2]  Yoram Rubin,et al.  Impact of hydrogeological data on measures of uncertainty, site characterization and environmental performance metrics , 2012 .

[3]  Rainer Helmig,et al.  Numerical simulation of biodegradation controlled by transverse mixing , 1999 .

[4]  R. Wilkin,et al.  Geochemical impacts to groundwater from geologic carbon sequestration: controls on pH and inorganic carbon concentrations from reaction path and kinetic modeling. , 2010, Environmental science & technology.

[5]  Liange Zheng,et al.  Effect of dissolved CO2 on a shallow groundwater system: a controlled release field experiment. , 2013, Environmental science & technology.

[6]  Daniel M. Tartakovsky,et al.  Introduction to the special issue on uncertainty quantification and risk assessment , 2012 .

[7]  A. Lasaga.  Rate laws of chemical reactions , 1981 .

[8]  Carsten Vogt,et al.  Investigation of the geochemical impact of CO2 on shallow groundwater: design and implementation of a CO2 injection test in Northeast Germany , 2012, Environmental Earth Sciences.

[9]  Martin J. Blunt,et al.  Streamline‐based simulation of solute transport , 1999 .

[10]  Martin J. Blunt,et al.  Streamline-based simulation of advective–dispersive solute transport , 2004 .

[11]  Reed M. Maxwell,et al.  Evaluating effective reaction rates of kinetically driven solutes in large‐scale, statistically anisotropic media: Human health risk implications , 2012 .

[12]  D. McLaughlin,et al.  Macrodispersivity and Large-scale Hydrogeologic Variability , 2001 .

[13]  Rajesh J. Pawar,et al.  The impact of CO2 on shallow groundwater chemistry: observations at a natural analog site and implications for carbon sequestration , 2010 .

[14]  Timothy R. Ginn,et al.  Stochastic‐Convective Transport with Nonlinear Reaction: Mathematical Framework , 1995 .

[15]  Reed M. Maxwell,et al.  Streamline-based simulation of virus transport resulting from long term artificial recharge in a heterogeneous aquifer , 2003 .

[16]  Reed M. Maxwell,et al.  Using streamlines to simulate stochastic reactive transport in heterogeneous aquifers: Kinetic metal release and transport in CO2 impacted drinking water aquifers , 2013 .

[17]  Daniel M. Tartakovsky,et al.  Assessment and management of risk in subsurface hydrology: A review and perspective , 2013 .

[18]  J. Daniels,et al.  Analysis of Uncertainty and Variability in Exposure to Characterize Risk: Case Study Involving Trichloroethylene Groundwater Contamination at Beale Air Force Base in California , 2000 .

[19]  Carl I. Steefel,et al.  Multidimensional, multicomponent, subsurface reactive transport in nonuniform velocity fields: code verification using an advective reactive streamtube approach , 1998 .

[20]  Martin J. Blunt,et al.  A 3D Field-Scale Streamline-Based Reservoir Simulator , 1997 .

[21]  T.-C. Jim Yeh,et al.  Stochastic modelling of groundwater flow and solute transport in aquifers , 1992 .

[22]  S. Carroll,et al.  Speciation and fate of trace metals in estuarine sediments under reduced and oxidized conditions, Seaplane Lagoon, Alameda Naval Air Station (USA) , 2002, Geochemical transactions.

[23]  Liange Zheng,et al.  Changes in the chemistry of shallow groundwater related to the 2008 injection of CO2 at the ZERT field site, Bozeman, Montana , 2010 .

[24]  Little,et al.  Comment on “ Potential Impacts of Leakage from Deep CO 2 Geosequestration on Overlying Freshwater Aquifers ” , 2011 .

[25]  Sookyun Wang,et al.  Dissolution of a mineral phase in potable aquifers due to CO2 releases from deep formations; effect of dissolution kinetics , 2004 .

[26]  C. S. Simmons,et al.  Stochastic-Convective Transport with Nonlinear Reaction: Biodegradation With Microbial Growth , 1995 .

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

[28]  E J Calabrese,et al.  Distinguishing outdoor soil ingestion from indoor dust ingestion in a soil pica child. , 1992, Regulatory toxicology and pharmacology : RTP.

[29]  T. Ginn.  Stochastic-convective transport with nonlinear reactions and mixing: finite streamtube ensemble formulation for multicomponent reaction systems with intra-streamtube dispersion. , 2001, Journal of contaminant hydrology.

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

[31]  Steven F. Carle,et al.  Contamination, risk, and heterogeneity: on the effectiveness of aquifer remediation , 2008 .

[32]  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 .

[33]  Albert J. Valocchi,et al.  A multidimensional streamline-based method to simulate reactive solute transport in heterogeneous porous media , 2010 .

[34]  H. Christopher Frey,et al.  Characterization and Simulation of Uncertain Frequency Distributions: Effects of Distribution Choice, Variability, Uncertainty, and Parameter Dependence , 1998 .

[35]  Y. Rubin.  Transport in heterogeneous porous media: Prediction and uncertainty , 1991 .

[36]  Hari S Viswanathan,et al.  A system model for geologic sequestration of carbon dioxide. , 2009, Environmental science & technology.

[37]  R. Spear,et al.  Integrating uncertainty and interindividual variability in environmental risk assessment. , 1987, Risk analysis : an official publication of the Society for Risk Analysis.

[38]  T E McKone,et al.  Uncertainties in health-risk assessment: an integrated case study based on tetrachloroethylene in California groundwater. , 1992, Regulatory toxicology and pharmacology : RTP.

[39]  F. O. Hoffman,et al.  Propagation of uncertainty in risk assessments: the need to distinguish between uncertainty due to lack of knowledge and uncertainty due to variability. , 1994, Risk analysis : an official publication of the Society for Risk Analysis.

[40]  Glenn E. Hammond,et al.  Elucidating geochemical response of shallow heterogeneous aquifers to CO2 leakage using high-performance computing: Implications for monitoring of CO2 sequestration , 2013 .

[41]  M. Blunt,et al.  Simulating Flow in Heterogeneous Systems Using Streamtubes and Streamlines , 1996 .

[42]  Kate Maher,et al.  The dependence of chemical weathering rates on fluid residence time , 2009 .

[43]  R. Maxwell,et al.  Stochastic environmental risk analysis: an integrated methodology for predicting cancer risk from contaminated groundwater , 1999 .

[44]  V. Cvetkovic,et al.  Evaluation of Risk from Contaminants Migrating by Groundwater , 1996 .

[45]  Dmitri Kavetski,et al.  Development of a hybrid process and system model for the assessment of wellbore leakage at a geologic CO2 sequestration site. , 2008, Environmental science & technology.

[46]  John L. Wilson,et al.  Efficient and accurate front tracking for two‐dimensional groundwater flow models , 1991 .

[47]  Carl I. Steefel,et al.  Fluid-rock interaction: A reactive transport approach , 2009 .

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

[49]  A. Lasaga.  Chemical kinetics of water‐rock interactions , 1984 .

[50]  Susan D. Pelmulder,et al.  On the development of a new methodology for groundwater‐Driven health risk assessment , 1998 .

[51]  Jery R. Stedinger,et al.  Probabilistic risk and uncertainty analysis for bioremediation of four chlorinated ethenes in groundwater , 2007 .

[52]  Steven F. Carle,et al.  Analysis of groundwater migration from artificial recharge in a large urban aquifer: A simulation perspective , 1999 .

[53]  John E. McCray,et al.  A quantitative methodology to assess the risks to human health from CO2 leakage into groundwater , 2010 .

[54]  S. Carroll,et al.  Rock−Water Interactions Controlling Zinc, Cadmium, and Lead Concentrations in Surface Waters and Sediments, U.S. Tri-State Mining District. 2. Geochemical Interpretation , 1998 .

[55]  H. Viswanathan,et al.  Comparison of streamtube and three-dimensional models of reactive transport in heterogeneous media , 2004 .

[56]  Jens Birkholzer,et al.  On mobilization of lead and arsenic in groundwater in response to CO2 leakage from deep geological storage , 2009 .

[57]  D. K. Smith,et al.  On the evaluation of groundwater contamination from underground nuclear tests , 2002 .

[58]  Liange Zheng,et al.  Evaluation of Potential Changes in Groundwater Quality in Response to CO2 Leakage from Deep Geologic Storage , 2010 .

[59]  Yue Hao,et al.  Geochemical detection of carbon dioxide in dilute aquifers , 2009, Geochemical transactions.

[60]  C. Steefel,et al.  A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution rea , 1994 .

[61]  E. R. Siirila,et al.  A new perspective on human health risk assessment: development of a time dependent methodology and the effect of varying exposure durations. , 2012, The Science of the total environment.

[62]  Giehyeon Lee,et al.  Geochemical implications of gas leakage associated with geologic CO2 storage--a qualitative review. , 2013, Environmental science & technology.

[63]  Y. Rubin,et al.  A methodology to integrate site characterization information into groundwater‐driven health risk assessment , 1999 .

[64]  A. Navarre‐Sitchler,et al.  Kinetic metal release from competing processes in aquifers. , 2012, Environmental science & technology.

[65]  Critical chemical reaction rates for multicomponent groundwater contamination models , 1987 .

[66]  Philippe Gouze,et al.  Effective non-local reaction kinetics for transport in physically and chemically heterogeneous media. , 2011, Journal of contaminant hydrology.