DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » 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.« less

Authors:
ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1];  [2];  [1]
  1. University of Texas, Austin, TX (United States)
  2. 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
  • DOI: 10.1016/j.jcp.2015.10.008

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
  • DOI: 10.1029/2011wr010562

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
  • DOI: 10.5194/esd-5-15-2014

The global groundwater crisis
journal, October 2014


Land subsidence in China
journal, July 2005


Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 2. Uniqueness, Stability, and Solution Algorithms
journal, February 1986


Radar interferometry techniques for the study of ground subsidence phenomena: a review of practical issues through cases in Spain
journal, March 2013


Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 3. Application to Synthetic and Field Data
journal, February 1986


Recent progress on reservoir history matching: a review
journal, July 2010


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
  • DOI: 10.1137/12089586X

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
  • DOI: 10.1137/130934805

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
  • DOI: 10.1137/110845598

General Theory of Three‐Dimensional Consolidation
journal, February 1941

  • Biot, Maurice A.
  • Journal of Applied Physics, Vol. 12, Issue 2
  • DOI: 10.1063/1.1712886

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
  • DOI: 10.1137/090780717

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
  • DOI: 10.1016/j.jcp.2010.03.018

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
  • DOI: 10.1029/2010gl044571

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
  • DOI: 10.1002/2013jb010272

Three-dimensional deformation and strain induced by municipal pumping, Part 2: Numerical analysis
journal, November 2006


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
  • DOI: 10.1109/tgrs.2007.894440

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
  • DOI: 10.1145/3428447

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
  • DOI: 10.21105/joss.00940

On uncertainty quantification in hydrogeology and hydrogeophysics
journal, December 2017


A Reassessment of the Groundwater Inverse Problem
journal, May 1996

  • McLaughlin, Dennis; Townley, Lloyd R.
  • Water Resources Research, Vol. 32, Issue 5
  • DOI: 10.1029/96WR00160

Groundwater depletion: A global problem
journal, February 2005


Geometric MCMC for infinite-dimensional inverse problems
journal, April 2017


The Influence of Geologic Structures on Deformation due to Ground Water Withdrawal
journal, March 2008


Three-dimensional deformation and strain induced by municipal pumping, part 1: Analysis of field data
journal, March 2006


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
  • DOI: 10.1007/978-3-642-23099-8

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
  • DOI: 10.1029/2020WR027391

Inherent Limitations of Hydraulic Tomography
journal, September 2010


Groundwater-pumping optimization for land-subsidence control in Beijing plain, China
journal, January 2018


Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information
journal, February 1986


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
  • DOI: 10.1111/j.1467-9868.2011.00777.x

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
  • DOI: 10.1002/nag.1062

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
  • DOI: 10.1137/14098171x

On the importance of geological heterogeneity for flow simulation
journal, February 2006


Overcoming the problem of locking in linear elasticity and poroelasticity: an heuristic approach
journal, December 2008


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
  • DOI: 10.1214/15-BA969