Approaches to Evaluate Water Quality Model Parameter Uncertainty for Adaptive TMDL Implementation1

Abstract:  The National Research Council recommended Adaptive Total Maximum Daily Load implementation with the recognition that the predictive uncertainty of water quality models can be high. Quantifying predictive uncertainty provides important information for model selection and decision-making. We review five methods that have been used with water quality models to evaluate model parameter and predictive uncertainty. These methods (1) Regionalized Sensitivity Analysis, (2) Generalized Likelihood Uncertainty Estimation, (3) Bayesian Monte Carlo, (4) Importance Sampling, and (5) Markov Chain Monte Carlo (MCMC) are based on similar concepts; their development over time was facilitated by the increasing availability of fast, cheap computers. Using a Streeter-Phelps model as an example we show that, applied consistently, these methods give compatible results. Thus, all of these methods can, in principle, provide useful sets of parameter values that can be used to evaluate model predictive uncertainty, though, in practice, some are quickly limited by the “curse of dimensionality” or may have difficulty evaluating irregularly shaped parameter spaces. Adaptive implementation invites model updating, as new data become available reflecting water-body responses to pollutant load reductions, and a Bayesian approach using MCMC is particularly handy for that task.

[1]  Kenneth H. Reckhow,et al.  Engineering approaches for lake management , 1983 .

[2]  Robert C. Spear,et al.  Large simulation models: calibration, uniqueness and goodness of fit , 1997 .

[3]  Adrian F. M. Smith,et al.  Bayesian computation via the gibbs sampler and related markov chain monte carlo methods (with discus , 1993 .

[4]  U. Aswathanarayana,et al.  Assessing the TMDL Approach to Water Quality Management , 2001 .

[5]  K. Beven,et al.  On constraining the predictions of a distributed model: The incorporation of fuzzy estimates of saturated areas into the calibration process , 1998 .

[6]  J. Geweke,et al.  Bayesian Inference in Econometric Models Using Monte Carlo Integration , 1989 .

[7]  S. Chib,et al.  Understanding the Metropolis-Hastings Algorithm , 1995 .

[8]  Mark E. Borsuk,et al.  On Monte Carlo methods for Bayesian inference , 2003 .

[9]  S. Carpenter,et al.  Ecosystem experiments : Frontiers in biology : ecology , 1995 .

[10]  Keith Beven,et al.  Equifinality and the problem of robust calibration in nitrogen budget simulations , 1999 .

[11]  C. S. Holling Adaptive Environmental Assessment and Management , 2005 .

[12]  David W. Dilks,et al.  Development of Bayesian Monte Carlo techniques for water quality model uncertainty , 1992 .

[13]  M. B. Beck,et al.  Water quality modeling: A review of the analysis of uncertainty , 1987 .

[14]  G. Hornberger,et al.  Approach to the preliminary analysis of environmental systems , 1981 .

[15]  G. Casella,et al.  Explaining the Gibbs Sampler , 1992 .

[16]  Keith Beven,et al.  The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[17]  T. Page,et al.  Predictive Capability in Estimating Changes in Water Quality: Long-Term Responses to Atmospheric Deposition , 2004 .

[18]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[19]  Keith Beven,et al.  Equifinality, sensitivity and predictive uncertainty in the estimation of critical loads , 1999 .

[20]  K. Beven,et al.  TOPOGRAPHIC CONTROLS ON SUBSURFACE STORM FLOWAT THE HILLSLOPE SCALE FOR TWO HYDROLOGICALLY DISTINCT SMALL CATCHMENTS , 1997 .

[21]  Keith Beven,et al.  On the sensitivity of soil-vegetation-atmosphere transfer (SVAT) schemes: equifinality and the problem of robust calibration , 1997 .

[22]  F. Pappenberger,et al.  Ignorance is bliss: Or seven reasons not to use uncertainty analysis , 2006 .

[23]  Keith Beven,et al.  Hydrological processes—Letters. Topographic controls on subsurface storm flow at the hillslope scale for two hydrologically distinct small catchmetns , 1997 .

[24]  Walter R. Gilks,et al.  A Language and Program for Complex Bayesian Modelling , 1994 .

[25]  D. Rubin Using the SIR algorithm to simulate posterior distributions , 1988 .

[26]  Keith Beven,et al.  UNCERTAINTY IN THE ESTIMATION OF CRITICAL LOADS: A PRACTICAL METHODOLOGY , 1997 .