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.

Bayesian Analysis of Uncertainties in Groundwater Numerical Simulation

Bayesian analysis for uncertainty and risk in a groundwater numerical model's predictions

Uncertainty Analysis in Statistical Modeling of Extreme Hydrological Events

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

Uncertainty analysis for decision-making in river basin management

Review of the Uncertainty Analysis of Groundwater Numerical Simulation

Uncertainty analysis of groundwater flow conceptual model

Uncertainty analysis for groundwater modeling

Numerical Approach for Water Distribution System Model Calibration Through Incorporation of Multiple Stochastic Prior Distributions.

Analysis on influence of aquifer heterogeneous simplification to groundwater numerical modeling

Uncertainty Assessment of Groundwater Prediction

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

Parameter Uncertainty Analysis on Hydrologic Frequency Models Based on Bayesian Theory

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

A Bayesian Approach for the Identification of Pollution Source in Water Quality Model Coupled with Hydrodynamics

The Characteristics of Probability Distribution of Groundwater Model Output Based on Sensitivity Analysis

Joint analysis of input and parametric uncertainties in watershed water quality modeling: A formal Bayesian approach

Quantification and Reduction of Groundwater Model Structural Uncertainty Based on Gaussian Process Regression

Risk analysis for reservoir flood control operation considering two-dimensional uncertainties based on Bayesian network

Uncertainty evaluation of groundwater modeling based on multimodel method

Reliability Analysis of the Groundwater Conceptual Model