skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Soil moisture estimation using tomographic ground penetrating radar in a MCMC–Bayesian framework

Abstract

In this paper, we focus on a hydrogeological inverse problem specifically targeting monitoring soil moisture variations using tomographic ground penetrating radar (GPR) travel time data. Technical challenges exist in the inversion of GPR tomographic data for handling non-uniqueness, nonlinearity and high-dimensionality of unknowns. We have developed a new method for estimating soil moisture fields from crosshole GPR data. It uses a pilot-point method to provide a low-dimensional representation of the relative dielectric permittivity field of the soil, which is the primary object of inference: the field can be converted to soil moisture using a petrophysical model. We integrate a multi-chain Markov chain Monte Carlo (MCMC)–Bayesian inversion framework with the pilot point concept, a curved-ray GPR travel time model, and a sequential Gaussian simulation algorithm, for estimating the dielectric permittivity at pilot point locations distributed within the tomogram, as well as the corresponding geostatistical parameters (i.e., spatial correlation range). We infer the dielectric permittivity as a probability density function, thus capturing the uncertainty in the inference. The multi-chain MCMC enables addressing high-dimensional inverse problems as required in the inversion setup. The method is scalable in terms of number of chains and processors, and is useful for computationally demanding Bayesian model calibrationmore » in scientific and engineering problems. The proposed inversion approach can successfully approximate the posterior density distributions of the pilot points, and capture the true values. The computational efficiency, accuracy, and convergence behaviors of the inversion approach were also systematically evaluated, by comparing the inversion results obtained with different levels of noises in the observations, increased observational data, as well as increased number of pilot points.« less

Authors:
 [1];  [2];  [3];  [2];  [4];  [2]
  1. Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Washington State Univ., Richland, WA (United States)
  2. Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
  3. Sandia National Lab. (SNL-CA), Livermore, CA (United States)
  4. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Publication Date:
Research Org.:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Sandia National Lab. (SNL-CA), Livermore, CA (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21)
OSTI Identifier:
1459653
Alternate Identifier(s):
OSTI ID: 1496976
Report Number(s):
PNNL-SA-116707; SAND-2018-7074J
Journal ID: ISSN 1436-3240; PII: 1571
Grant/Contract Number:  
AC05-76RL01830; AC04-94AL85000
Resource Type:
Accepted Manuscript
Journal Name:
Stochastic environmental research and risk assessment
Additional Journal Information:
Journal Volume: 32; Journal Issue: 8; Journal ID: ISSN 1436-3240
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
54 ENVIRONMENTAL SCIENCES; Tomographic ground penetrating radar; Soil moisture; Multi-chain Markov chain Monte Carlo; Bayesian

Citation Formats

Bao, Jie, Hou, Zhangshuan, Ray, Jaideep, Huang, Maoyi, Swiler, Laura, and Ren, Huiying. Soil moisture estimation using tomographic ground penetrating radar in a MCMC–Bayesian framework. United States: N. p., 2018. Web. doi:10.1007/s00477-018-1571-8.
Bao, Jie, Hou, Zhangshuan, Ray, Jaideep, Huang, Maoyi, Swiler, Laura, & Ren, Huiying. Soil moisture estimation using tomographic ground penetrating radar in a MCMC–Bayesian framework. United States. doi:10.1007/s00477-018-1571-8.
Bao, Jie, Hou, Zhangshuan, Ray, Jaideep, Huang, Maoyi, Swiler, Laura, and Ren, Huiying. Tue . "Soil moisture estimation using tomographic ground penetrating radar in a MCMC–Bayesian framework". United States. doi:10.1007/s00477-018-1571-8. https://www.osti.gov/servlets/purl/1459653.
@article{osti_1459653,
title = {Soil moisture estimation using tomographic ground penetrating radar in a MCMC–Bayesian framework},
author = {Bao, Jie and Hou, Zhangshuan and Ray, Jaideep and Huang, Maoyi and Swiler, Laura and Ren, Huiying},
abstractNote = {In this paper, we focus on a hydrogeological inverse problem specifically targeting monitoring soil moisture variations using tomographic ground penetrating radar (GPR) travel time data. Technical challenges exist in the inversion of GPR tomographic data for handling non-uniqueness, nonlinearity and high-dimensionality of unknowns. We have developed a new method for estimating soil moisture fields from crosshole GPR data. It uses a pilot-point method to provide a low-dimensional representation of the relative dielectric permittivity field of the soil, which is the primary object of inference: the field can be converted to soil moisture using a petrophysical model. We integrate a multi-chain Markov chain Monte Carlo (MCMC)–Bayesian inversion framework with the pilot point concept, a curved-ray GPR travel time model, and a sequential Gaussian simulation algorithm, for estimating the dielectric permittivity at pilot point locations distributed within the tomogram, as well as the corresponding geostatistical parameters (i.e., spatial correlation range). We infer the dielectric permittivity as a probability density function, thus capturing the uncertainty in the inference. The multi-chain MCMC enables addressing high-dimensional inverse problems as required in the inversion setup. The method is scalable in terms of number of chains and processors, and is useful for computationally demanding Bayesian model calibration in scientific and engineering problems. The proposed inversion approach can successfully approximate the posterior density distributions of the pilot points, and capture the true values. The computational efficiency, accuracy, and convergence behaviors of the inversion approach were also systematically evaluated, by comparing the inversion results obtained with different levels of noises in the observations, increased observational data, as well as increased number of pilot points.},
doi = {10.1007/s00477-018-1571-8},
journal = {Stochastic environmental research and risk assessment},
number = 8,
volume = 32,
place = {United States},
year = {2018},
month = {6}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:

Works referenced in this record:

Estimating flow parameter distributions using ground-penetrating radar and hydrological measurements during transient flow in the vadose zone
journal, June 2004


Geostatistical inversing for large-contrast transmissivity fields
journal, April 2008

  • Zanini, Andrea; Kitanidis, Peter K.
  • Stochastic Environmental Research and Risk Assessment, Vol. 23, Issue 5
  • DOI: 10.1007/s00477-008-0241-7

Measurement of Soil Moisture Content at Microwave Frequencies
journal, January 2015


Reservoir-parameter identification using minimum relative entropy-based Bayesian inversion of seismic AVA and marine CSEM data
journal, November 2006

  • Hou, Zhangshuan; Rubin, Yoram; Hoversten, G. Michael
  • GEOPHYSICS, Vol. 71, Issue 6
  • DOI: 10.1190/1.2348770

Hydrogeological characterization of the south oyster bacterial transport site using geophysical data
journal, October 2001

  • Hubbard, Susan S.; Chen, Jinsong; Peterson, John
  • Water Resources Research, Vol. 37, Issue 10
  • DOI: 10.1029/2001WR000279

Eikonal solver in the celerity domain
journal, July 2005


Finite‐difference calculation of traveltimes in three dimensions
journal, May 1990


Observations and stochastic modeling of soil moisture control on evapotranspiration in a Californian oak savanna: SOIL MOISTURE CONTROL ON ET
journal, August 2008

  • Chen, Xingyuan; Rubin, Yoram; Ma, Siyan
  • Water Resources Research, Vol. 44, Issue 8
  • DOI: 10.1029/2007WR006646

A fast sweeping method for Eikonal equations
journal, May 2004


A Bayesian peaks-over-threshold analysis of floods in the Itajaí-açu River under stationarity and nonstationarity
journal, November 2015

  • Silva, Artur Tiago; Portela, Maria Manuela; Naghettini, Mauro
  • Stochastic Environmental Research and Risk Assessment, Vol. 31, Issue 1
  • DOI: 10.1007/s00477-015-1184-4

Fast Sweeping Algorithms for a Class of Hamilton--Jacobi Equations
journal, January 2003

  • Tsai, Yen-Hsi Richard; Cheng, Li-Tien; Osher, Stanley
  • SIAM Journal on Numerical Analysis, Vol. 41, Issue 2
  • DOI: 10.1137/S0036142901396533

Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging
journal, October 2002


Ground-Penetrating Radar for High-Resolution Mapping of soil and rock Stratigraphy1
journal, July 1989


Stochastic estimation of facies using ground penetrating radar data
journal, November 2003

  • Moysey, S.; Caers, J.; Knight, R.
  • Stochastic Environmental Research and Risk Assessment (SERRA), Vol. 17, Issue 5
  • DOI: 10.1007/s00477-003-0152-6

Ground‐penetrating radar velocity tomography in heterogeneous and anisotropic media
journal, November 1997

  • Vasco, Don W.; Peterson, John E.; Lee, Ki Ha
  • GEOPHYSICS, Vol. 62, Issue 6
  • DOI: 10.1190/1.1444276

Distributed Soil Moisture from Crosshole Ground-Penetrating Radar Travel Times using Stochastic Inversion
journal, January 2013


Migration of common‐shot gathers
journal, February 1986


Soil moisture initialization for climate prediction: Characterization of model and observation errors
journal, January 2005


The assessment of high recharge areas using DO indicators and recharge potential analysis: a case study of Taiwan’s Pingtung plain
journal, September 2014

  • Tsai, Jui-Pin; Chen, Yu-Wen; Chang, Liang-Cheng
  • Stochastic Environmental Research and Risk Assessment, Vol. 29, Issue 3
  • DOI: 10.1007/s00477-014-0941-0

Applications of algebraic reconstruction techniques to crosshole seismic data
journal, October 1985

  • Peterson, John E.; Paulsson, Bjorn N. P.; McEvilly, Thomas V.
  • GEOPHYSICS, Vol. 50, Issue 10
  • DOI: 10.1190/1.1441847

A Bayesian approach for inverse modeling, data assimilation, and conditional simulation of spatial random fields
journal, October 2010

  • Rubin, Yoram; Chen, Xingyuan; Murakami, Haruko
  • Water Resources Research, Vol. 46, Issue 10
  • DOI: 10.1029/2009WR008799

Inference of multi-Gaussian relative permittivity fields by probabilistic inversion of crosshole ground-penetrating radar data
journal, July 2017


DRAM: Efficient adaptive MCMC
journal, December 2006


Dimensionality reduction for efficient Bayesian estimation of groundwater flow in strongly heterogeneous aquifers
journal, November 2016

  • Mara, Thierry A.; Fajraoui, Noura; Guadagnini, Alberto
  • Stochastic Environmental Research and Risk Assessment, Vol. 31, Issue 9
  • DOI: 10.1007/s00477-016-1344-1

Heterogeneous aquifer characterization from ground-penetrating radar tomography and borehole hydrogeophysical data using nonlinear Bayesian simulations
journal, July 2011

  • Dubreuil-Boisclair, Camille; Gloaguen, Erwan; Marcotte, Denis
  • GEOPHYSICS, Vol. 76, Issue 4
  • DOI: 10.1190/1.3571273

Ray‐based synthesis of bistatic ground‐penetrating radar profiles
journal, January 1995


Dynamic inversion for hydrological process monitoring with electrical resistance tomography under model uncertainties: A DYNAMIC INVERSION APPROACH
journal, April 2010

  • Lehikoinen, A.; Huttunen, J. M. J.; Finsterle, S.
  • Water Resources Research, Vol. 46, Issue 4
  • DOI: 10.1029/2009WR008470

Soil moisture content estimation using ground-penetrating radar reflection data
journal, June 2005


Simulation of ground‐penetrating radar waves in a 2-D soil model
journal, July 1996

  • Casper, David A.; Kung, K. ‐J. Samuel
  • GEOPHYSICS, Vol. 61, Issue 4
  • DOI: 10.1190/1.1444025

Computerized geophysical tomography
journal, January 1979


Efficient MCMC for Climate Model Parameter Estimation: Parallel Adaptive Chains and Early Rejection
journal, September 2012

  • Solonen, Antti; Ollinaho, Pirkka; Laine, Marko
  • Bayesian Analysis, Vol. 7, Issue 3
  • DOI: 10.1214/12-BA724

Crosshole Radar Tomography in a Fluvial Aquifer near Boise, Idaho
journal, September 2006

  • Clement, W. P.; Barrash, W.
  • Journal of Environmental & Engineering Geophysics, Vol. 11, Issue 3
  • DOI: 10.2113/JEEG11.3.171

Geophysical-hydrological identification of field permeabilities through Bayesian updating
journal, August 1993

  • Copty, Nadim; Rubin, Yoram; Mavko, Gary
  • Water Resources Research, Vol. 29, Issue 8
  • DOI: 10.1029/93WR00745

Ground Water Model Calibration Using Pilot Points and Regularization
journal, March 2003


Pre-inversion Corrections and Analysis of Radar Tomographic Data
journal, March 2001

  • Peterson, Jr., John E.
  • Journal of Environmental and Engineering Geophysics, Vol. 6, Issue 1
  • DOI: 10.4133/JEEG6.1.1

Smart pilot points using reversible-jump Markov-chain Monte Carlo: SMART PILOT POINTS
journal, May 2016

  • Jiménez, S.; Mariethoz, G.; Brauchler, R.
  • Water Resources Research, Vol. 52, Issue 5
  • DOI: 10.1002/2015WR017922

Laboratory Rainfall-Induced Slope Failure with Moisture Content Measurement
journal, May 2007


Application of the pilot point method to the identification of aquifer transmissivities
journal, October 1991


Forward modeling of ground-penetrating radar data using digitized outcrop images and multiple scenarios of water saturation
journal, June 2001

  • Kowalsky, M. B.; Dietrich, P.; Teutsch, G.
  • Water Resources Research, Vol. 37, Issue 6
  • DOI: 10.1029/2001WR900015

Bayesian Calibration of the Community Land Model Using Surrogates
journal, January 2015

  • Ray, J.; Hou, Z.; Huang, M.
  • SIAM/ASA Journal on Uncertainty Quantification, Vol. 3, Issue 1
  • DOI: 10.1137/140957998

Estimation of permeable pathways and water content using tomographic radar data
journal, November 1997

  • Hubbard, S. S.; Peterson, J. E.; Majer, E. L.
  • The Leading Edge, Vol. 16, Issue 11
  • DOI: 10.1190/1.1437539