Practical Postcalibration Uncertainty Analysis: Yucca Mountain, Nevada

The values of parameters in a groundwater flow model govern the precision of predictions of future system behavior. Predictive precision, thus, typically depends on an ability to infer values of system properties from historical measurements through calibration. When such data are scarce, or when their information content with respect to parameters that are most relevant to predictions of interest is weak, predictive uncertainty may be high, even if the model is "calibrated." Recent advances help recognize this condition, quantitatively evaluate predictive uncertainty, and suggest a path toward improved predictive accuracy by identifying sources of predictive uncertainty and by determining what observations will most effectively reduce this uncertainty. We demonstrate linear and nonlinear predictive error/uncertainty analyses as applied to a groundwater flow model of Yucca Mountain, Nevada, the United States' proposed site for disposal of high-level radioactive waste. Linear and nonlinear uncertainty analyses are readily implemented as an adjunct to model calibration with medium to high parameterization density. Linear analysis yields contributions made by each parameter to a prediction's uncertainty and the worth of different observations, both existing and yet-to-be-gathered, toward reducing this uncertainty. Nonlinear analysis provides more accurate characterization of the uncertainty of model predictions while yielding their (approximate) probability distribution functions. This article applies the above methods to a prediction of specific discharge and confirms the uncertainty bounds on specific discharge supplied in the Yucca Mountain Project License Application.

[1]  R. Hunt,et al.  Are Models Too Simple? Arguments for Increased Parameterization , 2007, Ground water.

[2]  D. M. Ely,et al.  A method for evaluating the importance of system state observations to model predictions, with application to the Death Valley regional groundwater flow system , 2004 .

[3]  John Doherty,et al.  Predictive error analysis for a water resource management model , 2007 .

[4]  John Doherty,et al.  Two statistics for evaluating parameter identifiability and error reduction , 2009 .

[5]  B. W. Arnold SATURATED ZONE FLOW AND TRANSPORT MODEL ABSTRACTION , 2004 .

[6]  J. Doherty,et al.  Calibration‐constrained Monte Carlo analysis of highly parameterized models using subspace techniques , 2009 .

[7]  John Doherty,et al.  Predictive error dependencies when using pilot points and singular value decomposition in groundwater model calibration , 2008 .

[8]  John Doherty,et al.  Efficient nonlinear predictive error variance for highly parameterized models , 2006 .

[9]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

[10]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[11]  G. Zyvoloski,et al.  A numerical model for thermo-hydro-mechanical coupling in fractured rock , 1997 .

[12]  J. Doherty,et al.  The cost of uniqueness in groundwater model calibration , 2006 .

[13]  George A. Zyvoloski,et al.  Finite element methods for geothermal reservoir simulation , 1983 .

[14]  J. Doherty,et al.  Role of the calibration process in reducing model predictive error , 2005 .

[15]  D. M. Ely,et al.  Preliminary evaluation of the importance of existing hydraulic-head observation locations to advective-transport predictions, Death Valley regional flow system, California and Nevada , 2001 .

[16]  C. Welty,et al.  Stochastic analysis of transverse dispersion in density‐coupled transport in aquifers , 2003 .

[17]  D. Sweetkind,et al.  Death Valley regional groundwater flow system, Nevada and California : hydrogeologic framework and transient groundwater flow model , 2010 .

[18]  P. Kitanidis On the geostatistical approach to the inverse problem , 1996 .

[19]  Richard L. Cooley,et al.  Simultaneous confidence and prediction intervals for nonlinear regression models with application to a groundwater flow model , 1987 .

[20]  J. Doherty,et al.  A hybrid regularized inversion methodology for highly parameterized environmental models , 2005 .

[21]  C. Tiedeman,et al.  Methods for using groundwater model predictions to guide hydrogeologic data collection, with application to the Death Valley regional groundwater flow system , 2003 .

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

[23]  W. Menke Geophysical data analysis : discrete inverse theory , 1984 .

[24]  SATURATED ZONE SITE-SCALE FLOW MODEL , 2004 .