Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty
Abstract
Abstract Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. Here, we develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)‐based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one‐dimensional problem and demonstrate it numerically for a three‐dimensional aquifer model. We design a generalized preconditioned Crank‐Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low‐rank‐based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstratedmore »
- Authors:
-
- University of Texas, Austin, TX (United States)
- Washington University, St. Louis, MO (United States)
- Publication Date:
- Research Org.:
- Univ. of Texas, Austin, TX (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC); National Science Foundation (NSF); Ministry of Education Saudi Arabia
- OSTI Identifier:
- 1978568
- Alternate Identifier(s):
- OSTI ID: 1832644
- Grant/Contract Number:
- SC0019303; CBET-1508713; ACI-1550593; DMS-2012453; DE‐SC0019303
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Water Resources Research
- Additional Journal Information:
- Journal Volume: 57; Journal Issue: 11; Journal ID: ISSN 0043-1397
- Publisher:
- American Geophysical Union (AGU)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 54 ENVIRONMENTAL SCIENCES
Citation Formats
Alghamdi, Amal, Hesse, Marc A., Chen, Jingyi, Villa, Umberto, and Ghattas, Omar. Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty. United States: N. p., 2021.
Web. doi:10.1029/2021wr029775.
Alghamdi, Amal, Hesse, Marc A., Chen, Jingyi, Villa, Umberto, & Ghattas, Omar. Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty. United States. https://doi.org/10.1029/2021wr029775
Alghamdi, Amal, Hesse, Marc A., Chen, Jingyi, Villa, Umberto, and Ghattas, Omar. Fri .
"Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty". United States. https://doi.org/10.1029/2021wr029775. https://www.osti.gov/servlets/purl/1978568.
@article{osti_1978568,
title = {Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty},
author = {Alghamdi, Amal and Hesse, Marc A. and Chen, Jingyi and Villa, Umberto and Ghattas, Omar},
abstractNote = {Abstract Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. Here, we develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)‐based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one‐dimensional problem and demonstrate it numerically for a three‐dimensional aquifer model. We design a generalized preconditioned Crank‐Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low‐rank‐based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two‐dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.},
doi = {10.1029/2021wr029775},
journal = {Water Resources Research},
number = 11,
volume = 57,
place = {United States},
year = {Fri Aug 06 00:00:00 EDT 2021},
month = {Fri Aug 06 00:00:00 EDT 2021}
}
Works referenced in this record:
Dimension-independent likelihood-informed MCMC
journal, January 2016
- Cui, Tiangang; Law, Kody J. H.; Marzouk, Youssef M.
- Journal of Computational Physics, Vol. 304
Nonsustainable groundwater sustaining irrigation: A global assessment: NONSUSTAINABLE GROUNDWATER SUSTAINING IRRIGATION
journal, January 2012
- Wada, Yoshihide; van Beek, L. P. H.; Bierkens, Marc F. P.
- Water Resources Research, Vol. 48, Issue 6
Global modeling of withdrawal, allocation and consumptive use of surface water and groundwater resources
journal, January 2014
- Wada, Y.; Wisser, D.; Bierkens, M. F. P.
- Earth System Dynamics, Vol. 5, Issue 1
The global groundwater crisis
journal, October 2014
- Famiglietti, J. S.
- Nature Climate Change, Vol. 4, Issue 11
Land subsidence in China
journal, July 2005
- Xue, Yu-Qun; Zhang, Yun; Ye, Shu-Jun
- Environmental Geology, Vol. 48, Issue 6
Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 2. Uniqueness, Stability, and Solution Algorithms
journal, February 1986
- Carrera, Jesus; Neuman, Shlomo P.
- Water Resources Research, Vol. 22, Issue 2
Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain
journal, March 2013
- Tomás, R.; Romero, R.; Mulas, J.
- Environmental Earth Sciences, Vol. 71, Issue 1
Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 3. Application to Synthetic and Field Data
journal, February 1986
- Carrera, Jesus; Neuman, Shlomo P.
- Water Resources Research, Vol. 22, Issue 2
Recent progress on reservoir history matching: a review
journal, July 2010
- Oliver, Dean S.; Chen, Yan
- Computational Geosciences, Vol. 15, Issue 1
A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems Part I: The Linearized Case, with Application to Global Seismic Inversion
journal, January 2013
- Bui-Thanh, Tan; Ghattas, Omar; Martin, James
- SIAM Journal on Scientific Computing, Vol. 35, Issue 6
A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems, Part II: Stochastic Newton MCMC with Application to Ice Sheet Flow Inverse Problems
journal, January 2014
- Petra, Noemi; Martin, James; Stadler, Georg
- SIAM Journal on Scientific Computing, Vol. 36, Issue 4
A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion
journal, January 2012
- Martin, James; Wilcox, Lucas C.; Burstedde, Carsten
- SIAM Journal on Scientific Computing, Vol. 34, Issue 3
General Theory of Three‐Dimensional Consolidation
journal, February 1941
- Biot, Maurice A.
- Journal of Applied Physics, Vol. 12, Issue 2
Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations
journal, January 2011
- Flath, H. P.; Wilcox, L. C.; Akçelik, V.
- SIAM Journal on Scientific Computing, Vol. 33, Issue 1
A fully coupled 3-D mixed finite element model of Biot consolidation
journal, June 2010
- Ferronato, Massimiliano; Castelletto, Nicola; Gambolati, Giuseppe
- Journal of Computational Physics, Vol. 229, Issue 12
Global depletion of groundwater resources: GLOBAL GROUNDWATER DEPLETION
journal, October 2010
- Wada, Yoshihide; van Beek, Ludovicus P. H.; van Kempen, Cheryl M.
- Geophysical Research Letters, Vol. 37, Issue 20
Joint inversion in coupled quasi-static poroelasticity: COUPLED POROELASTIC INVERSION
journal, February 2014
- Hesse, Marc A.; Stadler, Georg
- Journal of Geophysical Research: Solid Earth, Vol. 119, Issue 2
Three-dimensional deformation and strain induced by municipal pumping, Part 2: Numerical analysis
journal, November 2006
- Burbey, Thomas J.
- Journal of Hydrology, Vol. 330, Issue 3-4
Submillimeter Accuracy of InSAR Time Series: Experimental Validation
journal, May 2007
- Ferretti, Alessandro; Savio, Giuliano; Barzaghi, Riccardo
- IEEE Transactions on Geoscience and Remote Sensing, Vol. 45, Issue 5
hIPPYlib: An Extensible Software Framework for Large-Scale Inverse Problems Governed by PDEs: Part I: Deterministic Inversion and Linearized Bayesian Inference
journal, April 2021
- Villa, Umberto; Petra, Noemi; Ghattas, Omar
- ACM Transactions on Mathematical Software, Vol. 47, Issue 2
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet
journal, September 2015
- Isaac, Tobin; Petra, Noemi; Stadler, Georg
- Journal of Computational Physics, Vol. 296
hIPPYlib: An Extensible Software Framework for Large-Scale Inverse Problems
journal, October 2018
- Villa, Umberto; Petra, Noemi; Ghattas, Omar
- Journal of Open Source Software, Vol. 3, Issue 30
On uncertainty quantification in hydrogeology and hydrogeophysics
journal, December 2017
- Linde, Niklas; Ginsbourger, David; Irving, James
- Advances in Water Resources, Vol. 110
A Reassessment of the Groundwater Inverse Problem
journal, May 1996
- McLaughlin, Dennis; Townley, Lloyd R.
- Water Resources Research, Vol. 32, Issue 5
Groundwater depletion: A global problem
journal, February 2005
- Konikow, Leonard F.; Kendy, Eloise
- Hydrogeology Journal, Vol. 13, Issue 1
Geometric MCMC for infinite-dimensional inverse problems
journal, April 2017
- Beskos, Alexandros; Girolami, Mark; Lan, Shiwei
- Journal of Computational Physics, Vol. 335
The Influence of Geologic Structures on Deformation due to Ground Water Withdrawal
journal, March 2008
- Burbey, Thomas J.
- Ground Water, Vol. 46, Issue 2
Three-dimensional deformation and strain induced by municipal pumping, part 1: Analysis of field data
journal, March 2006
- Burbey, Thomas J.; Warner, Sandy M.; Blewitt, Geoff
- Journal of Hydrology, Vol. 319, Issue 1-4
Automated Solution of Differential Equations by the Finite Element Method
book, January 2012
- Logg, Anders; Mardal, Kent-Andre; Wells, Garth
- Lecture Notes in Computational Science and Engineering
Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate
journal, October 2020
- Alghamdi, Amal; Hesse, Marc A.; Chen, Jingyi
- Water Resources Research, Vol. 56, Issue 10
Inherent Limitations of Hydraulic Tomography
journal, September 2010
- Bohling, Geoffrey C.; Butler, James J.
- Ground Water, Vol. 48, Issue 6
Groundwater-pumping optimization for land-subsidence control in Beijing plain, China
journal, January 2018
- Qin, Huanhuan; Andrews, Charles B.; Tian, Fang
- Hydrogeology Journal, Vol. 26, Issue 4
Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information
journal, February 1986
- Carrera, Jesus; Neuman, Shlomo P.
- Water Resources Research, Vol. 22, Issue 2
An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach: Link between Gaussian Fields and Gaussian Markov Random Fields
journal, August 2011
- Lindgren, Finn; Rue, Håvard; Lindström, Johan
- Journal of the Royal Statistical Society: Series B (Statistical Methodology), Vol. 73, Issue 4
On the causes of pressure oscillations in low-permeable and low-compressible porous media: PRESSURE OSCILLATIONS IN LOW-PERMEABLE POROUS MEDIA
journal, July 2011
- Haga, Joachim Berdal; Osnes, Harald; Langtangen, Hans Petter
- International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 36, Issue 12
Algorithms for Kullback--Leibler Approximation of Probability Measures in Infinite Dimensions
journal, January 2015
- Pinski, F. J.; Simpson, G.; Stuart, A. M.
- SIAM Journal on Scientific Computing, Vol. 37, Issue 6
On the importance of geological heterogeneity for flow simulation
journal, February 2006
- Eaton, Timothy T.
- Sedimentary Geology, Vol. 184, Issue 3-4
Overcoming the problem of locking in linear elasticity and poroelasticity: an heuristic approach
journal, December 2008
- Phillips, Phillip Joseph; Wheeler, Mary F.
- Computational Geosciences, Vol. 13, Issue 1
On Bayesian A- and D-Optimal Experimental Designs in Infinite Dimensions
journal, September 2016
- Alexanderian, Alen; Gloor, Philip J.; Ghattas, Omar
- Bayesian Analysis, Vol. 11, Issue 3