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.

Use of stacked Gaussian processes regression method to improve prediction of groundwater solute transport model

Evaluation of Gaussian process regression kernel functions for improving groundwater prediction

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

Assessing Human Health Risk of Groundwater DNAPL Contamination by Quantifying the Model Structure Uncertainty

Improving prediction of groundwater quality in situations of limited monitoring data based on virtual sample generation and Gaussian process regression

Efficient Bayesian Inverse Modeling of Water Infiltration in Layered Soils

Development of a surrogate method of groundwater modeling using gated recurrent unit to improve the efficiency of parameter auto-calibration and global sensitivity analysis

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

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

Groundwater Contaminant Source Identification Via Bayesian Model Selection and Uncertainty Quantification

Comparison of Multiple Machine Learning Methods for Correcting Groundwater Levels Predicted by Physics-Based Models

An improved method of groundwater model structural uncertainty analysis

Efficient Simulation of Groundwater Solute Transport Using the Multipoint Flux Approximation Method with Arbitrary Polygon Grids

Assessment of Parametric Uncertainty for Groundwater Reactive Transport Modeling

A Spatiotemporal Regression-Kriging Model for Space-Time Interpolation: a Case Study of Chlorophyll-a Prediction in the Coastal Areas of Zhejiang, China.

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

Stochastic Modeling of Phosphorus Transport in the Three Gorges Reservoir by Incorporating Variability Associated with the Phosphorus Partition Coefficient.

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

Identifying the spatial pattern and driving factors of nitrate in groundwater using a novel framework of interpretable stacking ensemble learning

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

A New Efficient Approach to Model Stiff Biogeochemical Reactive Transport Scenarios Across Groundwater Systems