Bayesian calibration and sensitivity analysis for a karst aquifer model using active subspaces

Abstract

In this article, we perform a parameter study for a recently developed karst hydrological model. The study consists of a high-dimensional Bayesian inverse problem and a global sensitivity analysis. For the first time in karst hydrology, we use the active subspace method to find directions in the parameter space that dominate the Bayesian update from the prior to the posterior distribution in order to effectively reduce the dimension of the problem and for computational efficiency. Additionally, the calculated active subspace can be exploited to construct sensitivity metrics on each of the individual parameters and be used to construct a natural model surrogate. The model consists of 21 parameters to reproduce the hydrological behavior of spring discharge in a karst aquifer located in the Kerschbaum spring recharge area at Waidhofen a.d. Ybbs in Austria. The experimental spatial and time series data for the inference process were collected by the water works in Waidhofen. We show that this case study has implicit low-dimensionality, and we run an adjusted Markov chain Monte Carlo algorithm in a low-dimensional subspace to construct samples of the posterior distribution. The results are visualized and verified by plots of the posterior's push-forward distribution displaying the uncertainty in predicting discharge values due to the experimental noise in the data. Finally, a discussion provides hydrological interpretation of these results for the Kerschbaum area.

arXiv comment27 pages, 5 figures, 2 tables; 5 pages supplementary information

References in corpus (10)

• MCMC using Hamiltonian dynamics
• User-friendly tail bounds for sums of random matrices
• Active subspace methods in theory and practice: applications to kriging surfaces
• Dimension-independent likelihood-informed MCMC
• Exploiting Active Subspaces to Quantify Uncertainty in the Numerical Simulation of the HyShot II Scramjet
• Accelerating MCMC with active subspaces
• Solving Large-Scale PDE-constrained Bayesian Inverse Problems with Riemann Manifold Hamiltonian Monte Carlo
• Computing active subspaces with Monte Carlo
• Efficient parameter estimation for a methane hydrate model with active subspaces
• Goal-oriented adaptive surrogate construction for stochastic inversion

Cited by in corpus (3)

• Machine learning accelerates parameter optimization and uncertainty assessment of a land surface model
• Generalized bounds for active subspaces
• An efficient estimation of time-varying parameters of dynamic models by combining offline batch optimization and online data assimilation