Abstract:The validity of using Gaussian assumptions for model residuals in uncertainty quantification of a groundwater reactive transport model was evaluated in this study. Least squares regression methods explicitly assume Gaussian residuals, and the assumption leads to Gaussian likelihood functions, model parameters, and model predictions. While the Bayesian methods do not explicitly require the Gaussian assumption, Gaussian residuals are widely used. This paper shows that the residuals of the reactive transport model are non-Gaussian, heteroscedastic, and correlated in time; characterizing them requires using a generalized likelihood function such as the formal generalized likelihood function developed by Schoups and Vrugt (2010). For the surface complexation model considered in this study for simulating uranium reactive transport in groundwater, parametric uncertainty is quantified using the least squares regression methods and Bayesian methods with both Gaussian and formal generalized likelihood functions. While the least squares methods and Bayesian methods with Gaussian likelihood function produce similar Gaussian parameter distributions, the parameter distributions of Bayesian uncertainty quantification using the formal generalized likelihood function are non-Gaussian. In addition, predictive performance of formal generalized likelihood function is superior to that of least squares regression and Bayesian methods with Gaussian likelihood function. The Bayesian uncertainty quantification is conducted using the differential evolution adaptive metropolis (DREAM((ZS))) algorithm; as a Markov chain Monte Carlo (MCMC) method, it is a robust tool for quantifying uncertainty in groundwater reactive transport models. For the surface complexation model, the regression-based local sensitivity analysis and Morris-and DREAM((ZS))-based global sensitivity analysis yield almost identical ranking of parameter importance. The uncertainty analysis may help select appropriate likelihood functions, improve model calibration, and reduce predictive uncertainty in other groundwater reactive transport and environmental modeling.

Global sensitivity analysis in Monte Carlo radiation transport

Polynomial chaos expansion for global sensitivity analysis applied to a model of radionuclide migration in a randomly heterogeneous aquifer

Preliminary Study of Generalized Sensitivity Analysis with Continuous-Energy Monte Carlo Code RMC

Global sensitivity analysis of natural convection in porous enclosure: Effect of thermal dispersion, anisotropic permeability and heterogeneity

Application of Stochastic and Deterministic Techniques for Uncertainty Quantification and Sensitivity Analysis of Energy Systems

Comparative Assessment of Two Global Sensitivity Approaches Considering Model and Parameter Uncertainty

Global sensitivity analysis of computer models with functional inputs

A geostatistics‐informed hierarchical sensitivity analysis method for complex groundwater flow and transport modeling

Uncertainty quantification of contaminant transport and risk assessment with conditional stochastic collocation method

Improvement of Sensitivity and Uncertainty Analysis Capabilities of Generalized Response in Monte Carlo Code RMC

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

A Universal Adjoint-Weighted Algorithm for Geometric Sensitivity Analysis of K-Eigenvalue Based on Continuous-Energy Monte Carlo Method

Global variance reduction method for global Monte Carlo particle transport simulations of CFETR

Generalized Sensitivity Analysis of Responses to Nuclear Data in 2-D/1-D Transport Code KYADJ

Development of Sensitivity Analysis Capabilities of Generalized Responses to Nuclear Data in Monte Carlo Code RMC

Quasi-Global Sensitivity Analysis of Supersonic Computational Fluid Dynamics Simulations

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

Geometric Sensitivity Analysis of Generalized Response Function with RMC Code

Efficient Evaluation of Small Failure Probability in High-Dimensional Groundwater Contaminant Transport Modeling Via A Two-Stage Monte Carlo Method

Mixed-Variable Global Sensitivity Analysis For Knowledge Discovery And Efficient Combinatorial Materials Design

Advanced methodology for uncertainty propagation in computer experiments with large number of inputs