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.

Gaussian Process Modelling for Uncertainty Quantification in Convectively-Enhanced Dissolution Processes in Porous Media

A Surrogate Modelling Approach Based on Nonlinear Dimension Reduction for Uncertainty Quantification in Groundwater Flow Models

Residual Gaussian Process: A Tractable Nonparametric Bayesian Emulator for Multi-Fidelity Simulations

Uncertainty Quantification of Two-Phase Flow in Porous Media Via the Coupled-TgNN Surrogate Model

Uncertainty Quantification for Flow and Transport in Highly Heterogeneous Porous Media Based on Simultaneous Stochastic Model Dimensionality Reduction

Gaussian Process Emulators for Quantifying Uncertainty Inco2spreading Predictions in Heterogeneous Media

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

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

Uncertainty modeling and propagation for groundwater flow: a comparative study of surrogates

An Adaptive Gaussian Process‐based Method for Efficient Bayesian Experimental Design in Groundwater Contaminant Source Identification Problems

Propagation of Uncertainties in Density-Driven Flow

Deep Gaussian Process Emulation and Uncertainty Quantification for Large Computer Experiments

Monte Carlo Simulation for Variable-Density Groundwater Flow Through Reduced-Order Model Coupled with Gaussian Process

An Adaptive Gaussian Process-Based Iterative Ensemble Smoother for Data Assimilation

Uncertainty quantification for mineral precipitation and dissolution in fractured porous media

PDE-constrained Gaussian process surrogate modeling with uncertain data locations

Evaluating Two Sparse Grid Surrogates and Two Adaptation Criteria for Groundwater Bayesian Uncertainty Quantification

Propagation of Input Uncertainty in Presence of Model-Form Uncertainty: A Multifidelity Approach for Computational Fluid Dynamics Applications

Analytical uncertainty quantification approach based on adaptive generalized co‐Gaussian process model

An Efficient Approach for Quantifying Parameter Uncertainty in the SST Turbulence Model