Dimension reduction for integrating data series in Bayesian inversion of geostatistical models
- This study explores methods with which multidimensional data, e.g. time series, can be effectively incorporated into a Bayesian framework for inferring geostatistical parameters. Such series are difficult to use directly in the likelihood estimation procedure due to their high dimensionality; thus, a dimension reduction approach is taken to utilize these measurements in the inference. Two synthetic scenarios from hydrology are explored in which pumping drawdown and concentration breakthrough curves are used to infer the global mean of a log-normally distributed hydraulic conductivity field. Both cases pursue the use of a parametric model to represent the shape of the observed time series with physically-interpretable parameters (e.g. the time and magnitude of a concentration peak), which is compared to subsets of the observations with similar dimensionality. The results from both scenarios highlight the effectiveness for the shape-matching models to reduce dimensionality from 100+ dimensions down to less than five. The modelsThis study explores methods with which multidimensional data, e.g. time series, can be effectively incorporated into a Bayesian framework for inferring geostatistical parameters. Such series are difficult to use directly in the likelihood estimation procedure due to their high dimensionality; thus, a dimension reduction approach is taken to utilize these measurements in the inference. Two synthetic scenarios from hydrology are explored in which pumping drawdown and concentration breakthrough curves are used to infer the global mean of a log-normally distributed hydraulic conductivity field. Both cases pursue the use of a parametric model to represent the shape of the observed time series with physically-interpretable parameters (e.g. the time and magnitude of a concentration peak), which is compared to subsets of the observations with similar dimensionality. The results from both scenarios highlight the effectiveness for the shape-matching models to reduce dimensionality from 100+ dimensions down to less than five. The models outperform the alternative subset method, especially when the observations are noisy. This approach to incorporating time series observations in the Bayesian framework for inferring geostatistical parameters allows for high-dimensional observations to be faithfully represented in lower-dimensional space for the non-parametric likelihood estimation procedure, which increases the applicability of the framework to more observation types. Although the scenarios are both from hydrogeology, the methodology is general in that no assumptions are made about the subject domain. Any application that requires the inference of geostatistical parameters using series in either time of space can use the approach described in this paper.…
Author details: | Heather SavoyORCiD, Falk HeßeORCiDGND |
---|---|
DOI: | https://doi.org/10.1007/s00477-019-01697-9 |
ISSN: | 1436-3240 |
ISSN: | 1436-3259 |
Title of parent work (English): | Stochastic environmental research and risk assessment |
Publisher: | Springer |
Place of publishing: | New York |
Publication type: | Article |
Language: | English |
Date of first publication: | 2019/07/03 |
Publication year: | 2019 |
Release date: | 2021/01/14 |
Tag: | Bayesian inference; Dimension reduction; Geostatistics; Stochastic hydrogeology |
Volume: | 33 |
Issue: | 7 |
Number of pages: | 18 |
First page: | 1327 |
Last Page: | 1344 |
Funding institution: | National Science FoundationNational Science Foundation (NSF) [EAR-1011336]; Office of Science of the U.S. Department of EnergyUnited States Department of Energy (DOE) [DE-AC02-05CH11231] |
Organizational units: | Mathematisch-Naturwissenschaftliche Fakultät / Institut für Umweltwissenschaften und Geographie |
DDC classification: | 5 Naturwissenschaften und Mathematik / 55 Geowissenschaften, Geologie / 550 Geowissenschaften |
Peer review: | Referiert |