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 analysis of groundwater flow conceptual model

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

Reliability Analysis of the Groundwater Conceptual Model

Uncertainty analysis for groundwater modeling

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

Review of the Uncertainty Analysis of Groundwater Numerical Simulation

Bayesian Analysis of Uncertainties in Groundwater Numerical Simulation

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

Uncertainty Assessment of Groundwater Prediction

Uncertainty evaluation of groundwater modeling based on multimodel method

Analysis on influence of aquifer heterogeneous simplification to groundwater numerical modeling

Uncertainty Evaluation of a Groundwater Conceptual Model by Using a Multimodel Averaging Method

Assessing Groundwater Modeling Uncertainty By Model Averaging Method

Evaluating Groundwater Interbasin Flow Using Multiple Models and Multiple Types of Data

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

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

Incorporating interpretation uncertainties from deterministic 3D hydrostratigraphic models in groundwater models

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

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

Application of Multiscale Finite Element Method in the Uncertainty Qualification of Large-Scale Groundwater Flow