Modeling hypoxia in the Chesapeake Bay: Ensemble estimation using a Bayesian hierarchical model

Abstract Quantifying parameter and prediction uncertainty in a rigorous framework can be an important component of model skill assessment. Generally, models with lower uncertainty will be more useful for prediction and inference than models with higher uncertainty. Ensemble estimation, an idea with deep roots in the Bayesian literature, can be useful to reduce model uncertainty. It is based on the idea that simultaneously estimating common or similar parameters among models can result in more precise estimates. We demonstrate this approach using the Streeter–Phelps dissolved oxygen sag model fit to 29 years of data from Chesapeake Bay. Chesapeake Bay has a long history of bottom water hypoxia and several models are being used to assist management decision-making in this system. The Bayesian framework is particularly useful in a decision context because it can combine both expert-judgment and rigorous parameter estimation to yield model forecasts and a probabilistic estimate of the forecast uncertainty.

[1]  S. Cooper,et al.  Chesapeake Bay Watershed Historical Land Use: Impact on Water Quality and Diatom Communities , 1995 .

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

[3]  Dubravko Justic,et al.  Predicting the response of Gulf of Mexico hypoxia to variations in Mississippi River nitrogen load , 2003 .

[4]  B. Efron,et al.  Empirical Bayes on vector observations: An extension of Stein's method , 1972 .

[5]  James O. Berger,et al.  Statistical Analysis and the Illusion of Objectivity , 1988 .

[6]  G. Milliken Nonlinear Regression Analysis and Its Applications , 1990 .

[7]  Kenneth H. Reckhow,et al.  Importance of scientific uncertainty in decision making , 1994 .

[8]  B. Efron,et al.  Combining Possibly Related Estimation Problems , 1973 .

[9]  C. Stow,et al.  Predicting the frequency of water quality standard violations: a probabilistic approach for TMDL development. , 2002, Environmental science & technology.

[10]  R. Díaz Overview of hypoxia around the world. , 2001, Journal of environmental quality.

[11]  James D. Hagy,et al.  A simple model for forecasting the effects of nitrogen loads on Chesapeake Bay hypoxia , 2006 .

[12]  Michael R. Roman,et al.  Eutrophication of Chesapeake Bay: historical trends and ecological interactions , 2005 .

[13]  C. Morris Parametric Empirical Bayes Inference: Theory and Applications , 1983 .

[14]  Zehao Shen,et al.  Ecological applications of multilevel analysis of variance. , 2007, Ecology.

[15]  James D. Hagy,et al.  Hypoxia in Chesapeake Bay, 1950–2001: Long-term change in relation to nutrient loading and river flow , 2004 .

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

[17]  V. L. Johnson,et al.  Factors contributing to hypoxia in rivers, lakes, and streams , 2006 .

[18]  G. Brush,et al.  Long-Term History of Chesapeake Bay Anoxia , 1991, Science.

[19]  Stephen J. Ganocy,et al.  Bayesian Statistical Modelling , 2002, Technometrics.

[20]  L E Cronin,et al.  Chesapeake Bay Anoxia: Origin, Development, and Significance , 1984, Science.

[21]  B. Efron,et al.  Data Analysis Using Stein's Estimator and its Generalizations , 1975 .

[22]  Carl F. Cerco,et al.  Three‐Dimensional Eutrophication Model of Chesapeake Bay , 1993 .

[23]  Mark E. Borsuk,et al.  Approaches to Evaluate Water Quality Model Parameter Uncertainty for Adaptive TMDL Implementation1 , 2007 .

[24]  Craig A Stow,et al.  A Bayesian approach to retransformation bias in transformed regression. , 2006, Ecology.

[25]  Joseph Hilbe,et al.  Data Analysis Using Regression and Multilevel/Hierarchical Models , 2009 .

[26]  B. Efron,et al.  Stein's Estimation Rule and Its Competitors- An Empirical Bayes Approach , 1973 .

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

[28]  Dubravko Justic,et al.  Reducing hypoxia in the Gulf of Mexico: Advice from three models , 2004 .

[29]  Larry B. Crowder,et al.  Hypoxia-induced habitat shifts and energetic consequences in Atlantic croaker and brown shrimp on the Gulf of Mexico shelf , 2005 .

[30]  Larry B. Crowder,et al.  Habitat degradation from intermittent hypoxia: impacts on demersal fishes , 2005 .

[31]  C L Newcombe,et al.  OXYGEN-POOR WATERS OF THE CHESAPEAKE BAY. , 1938, Science.

[32]  George E. P. Box,et al.  The Bayesian estimation of common parameters from several responses , 1965 .

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

[34]  David Higdon,et al.  A Bayesian hierarchical model to predict benthic oxygen demand from organic matter loading in estuaries and coastal zones , 2001 .

[35]  Douglas M. Bates,et al.  Nonlinear Regression Analysis and Its Applications , 1988 .