TY - JOUR
T1 - Non‐Gaussian Parameter Inference for Hydrogeological Models Using Stein Variational Gradient Descent
AU - Ramgraber, Maximilian
AU - Weatherl, Robin
AU - Blumensaat, Frank
AU - Schirmer, Mario
PY - 2021/4
Y1 - 2021/4
N2 - The sustainable management of groundwater demands a faithful characterization of the subsurface. This, in turn, requires information which is generally not readily available. To bridge the gap between data need and availability, numerical models are often used to synthesize plausible scenarios not only from direct information but also from additional, indirect data. Unfortunately, the resulting system characterizations will rarely be unique. This poses a challenge for practical parameter inference: computational limitations often force modelers to resort to methods based on questionable assumptions of Gaussianity, which do not reproduce important facets of ambiguity such as Pareto fronts or multimodality. In search of a remedy, an alternative could be found in Stein Variational Gradient Descent (SVGD), a recent development in the field of statistics. This ensemble-based method iteratively transforms a set of arbitrary particles into samples of a potentially non-Gaussian posterior, provided the latter is sufficiently smooth. A prerequisite for this method is knowledge of the Jacobian, which is usually exceptionally expensive to evaluate. To address this issue, we propose an ensemble-based, localized approximation of the Jacobian. We demonstrate the performance of the resulting algorithm in two cases: a simple, bimodal synthetic scenario, and a complex numerical model based on a real world, prealpine catchment. Promising results in both cases—even when the ensemble size is smaller than the number of parameters—suggest that SVGD can be a valuable addition to hydrogeological parameter inference.
AB - The sustainable management of groundwater demands a faithful characterization of the subsurface. This, in turn, requires information which is generally not readily available. To bridge the gap between data need and availability, numerical models are often used to synthesize plausible scenarios not only from direct information but also from additional, indirect data. Unfortunately, the resulting system characterizations will rarely be unique. This poses a challenge for practical parameter inference: computational limitations often force modelers to resort to methods based on questionable assumptions of Gaussianity, which do not reproduce important facets of ambiguity such as Pareto fronts or multimodality. In search of a remedy, an alternative could be found in Stein Variational Gradient Descent (SVGD), a recent development in the field of statistics. This ensemble-based method iteratively transforms a set of arbitrary particles into samples of a potentially non-Gaussian posterior, provided the latter is sufficiently smooth. A prerequisite for this method is knowledge of the Jacobian, which is usually exceptionally expensive to evaluate. To address this issue, we propose an ensemble-based, localized approximation of the Jacobian. We demonstrate the performance of the resulting algorithm in two cases: a simple, bimodal synthetic scenario, and a complex numerical model based on a real world, prealpine catchment. Promising results in both cases—even when the ensemble size is smaller than the number of parameters—suggest that SVGD can be a valuable addition to hydrogeological parameter inference.
UR - http://www.scopus.com/inward/record.url?scp=85104852960&partnerID=8YFLogxK
U2 - 10.1029/2020WR029339
DO - 10.1029/2020WR029339
M3 - Article
SN - 0043-1397
VL - 57
JO - Water Resources Research
JF - Water Resources Research
IS - 4
M1 - e2020WR029339
ER -