An efficient Bayesian framework incorporating the Stochastic Response Surface Method for the uncertainty analysis of subsurface and multimedia transport models
S. Balakrishnan (2,3), A. Roy (3), M. Ierapetritou (2), G.P. Flach (1), P.G. Georgopoulos (2,3)
1. Savannah River Technology Center; 2. Department of Chemical and Biochemical Engineering, Rutgers University; 3. Environmental and Occupational Health Sciences Institute, UMDNJ - R.W. Johnson Medical School and Rutgers University
Comprehensive uncertainty analyses of complex models of environmental and biological systems are essential but often not feasible due to the computational resources they require. Traditional Monte Carlo methods, e.g. involving standard or Latin Hypercube sampling, for propagating uncertainty and developing probability densities of model outputs, may require performing a large number of model simulations that can be prohibitive in the case of 3D dynamic finite element or finite difference models. An alternative approach is provided by the Stochastic Response Surface Method (SRSM), which is a computationally efficient framework that facilitates uncertainty analysis through the determination of statistically equivalent reduced models; SRSM is applied here to subsurface transport modeling. Past applications of SRSM included the uncertainty analysis of atmospheric chemistry and transport models, groundwater models, and physiologically based pharmacokinetic models [Isukapalli et al., 1998, Risk Analysis, 18(3), 351-363;Isukapalli et al., 2000, Risk Analysis, 20(5), 591-602]. The present work further demonstrates how the SRSM framework can be used in a Bayesian manner not only to characterize but also to reduce uncertainties, by incorporating observational information in estimates of model parameters. This is achieved by the combined application of the SRSM and Markov Chain Monte Carlo (MCMC) methods. Current and ongoing work includes the application of this novel framework to the groundwater and unsaturated zone modeling of the General Separations Area (GSA) of the US DOE Savannah River Site. Subsurface processes are simulated with the Flow And Contaminant Transport (FACT) code Flach, G.P. and M.K. Harris, 2000, WSRC-TR-96-00399). FACT is transient three-dimensional, finite element code designed to simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media. Furthermore, FACT is designed specifically to handle complex multi-layer and/or heterogeneous aquifer systems in an efficient manner while accommodating a wide range of boundary conditions (Hamm, L.L and S.E. Aleman, 2000, WSRC-TR-99-00282). The present work demonstrates the feasibility of performing systematic uncertainty analysis and parameter optimization with a complex model such as FACT, by employing the combined SRSM/MCMC approach.
Acknowledgments: This work has been funded in part by the US Environmental Protection Agency under Cooperative Agreement # EPAR-827033 to the Environmental and Occupational Health Sciences Institute; and by a grant to the Institute for Responsible Management, Consortium for Risk Evaluation with Stakeholder Participation from the US Department of Energy, Instrument DE-FG2600NT 40938. The viewpoints expressed in this work are solely the responsibility of the authors and do not necessarily reflect the views of the US Department of Energy, the US Environmental Protection Agency, or their contractors.