An interactive Bayesian geostatistical inverse protocol for hydraulic tomography

Water Resources Research
By: , and 

Links

Abstract

Hydraulic tomography is a powerful technique for characterizing heterogeneous hydrogeologic parameters. An explicit trade-off between characterization based on measurement misfit and subjective characterization using prior information is presented. We apply a Bayesian geostatistical inverse approach that is well suited to accommodate a flexible model with the level of complexity driven by the data and explicitly considering uncertainty. Prior information is incorporated through the selection of a parameter covariance model characterizing continuity and providing stability. Often, discontinuities in the parameter field, typically caused by geologic contacts between contrasting lithologic units, necessitate subdivision into zones across which there is no correlation among hydraulic parameters. We propose an interactive protocol in which zonation candidates are implied from the data and are evaluated using cross validation and expert knowledge. Uncertainty introduced by limited knowledge of dynamic regional conditions is mitigated by using drawdown rather than native head values. An adjoint state formulation of MODFLOW-2000 is used to calculate sensitivities which are used both for the solution to the inverse problem and to guide protocol decisions. The protocol is tested using synthetic two-dimensional steady state examples in which the wells are located at the edge of the region of interest.

Additional publication details

Publication type Article
Publication Subtype Journal Article
Title An interactive Bayesian geostatistical inverse protocol for hydraulic tomography
Series title Water Resources Research
DOI 10.1029/2007WR006730
Volume 44
Issue 12
Year Published 2008
Language English
Publisher Wiley
Contributing office(s) Wisconsin Water Science Center
Description 19 p.
Larger Work Type Article
Larger Work Subtype Journal Article
Larger Work Title Water Resources Research