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.

Uncertainty Assessment of Groundwater Prediction

Uncertainty Analysis in Statistical Modeling of Extreme Hydrological Events

Decision-making in Water Management under Uncertainty

Uncertainty in High and Low Flows Due to Model Structure and Parameter Errors

Uncertainty analysis for decision-making in river basin management

Application of Uncertainty and Sensitivity Analysis in River Basin Management.

Uncertainty analysis of groundwater flow conceptual model

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

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

Review of the Uncertainty Analysis of Groundwater Numerical Simulation

Assessing Bayesian Model Averaging Uncertainty of Groundwater Modeling Based on Information Entropy Method

Uncertainty analysis for groundwater modeling

Bayesian Analysis of Uncertainties in Groundwater Numerical Simulation

Predictive Assessment of Groundwater Flow Uncertainty in Multiscale Porous Media by Using Truncated Power Variogram Model

Uncertainty evaluation of groundwater modeling based on multimodel method

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

Reliability Analysis of the Groundwater Conceptual Model

Sensitivity Analysis of the Probability Distribution of Groundwater Level Series Based on Information Entropy

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

Uncertainty of heterogeneous hydrogeological models in groundwater flow and land subsidence simulations – A case study in Huwei Town, Taiwan

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