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 for uncertainty and risk in a groundwater numerical model's predictions

Bayesian Analysis of Uncertainties in Groundwater Numerical Simulation

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

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

Uncertainty analysis of groundwater flow conceptual model

Uncertainty analysis for groundwater modeling

Analysis on influence of aquifer heterogeneous simplification to groundwater numerical modeling

State Estimation Based on Enhanced Bayesian Approach: Application in Water Distribution Systems

Uncertainty Assessment of Groundwater Prediction

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

Review of the Uncertainty Analysis of Groundwater Numerical Simulation

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

Parameter Uncertainty Analysis on Hydrologic Frequency Models Based on Bayesian Theory

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

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

Reliability Analysis of the Groundwater Conceptual Model

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

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

Uncertainty and Risk in Water Quality Modelling and Management

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