Non‐Gaussian Parameter Inference for Hydrogeological Models Using Stein Variational Gradient Descent

• Published in: Water resources research 2021-04, Vol.57 (4), p.n/a
• Main Authors: Ramgraber, Maximilian, Weatherl, Robin, Blumensaat, Frank, Schirmer, Mario
• Format: Article
• Language: English
• Subjects: Algorithms; Approximation; Bayesian statistics; Catchment area; Computer applications; Fronts; Geology; Groundwater; Hydrogeological models; Hydrogeology; Inference; Mathematical models; modeling; non‐Gaussian; Numerical models; Parameters; Statistical methods; stein variational gradient descent; Sustainability management
• ISSN: 0043-1397; 1944-7973
• DOI: 10.1029/2020WR029339

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. Key Points This manuscript explores Stein Variational Gradient Descent for parameter inference in hydrogeological models We introduce an ensemble‐based Jacobian approximation and test the algorithm in a synthetic and a real test case The algorithm performs well in high dimensional and multimodal, non‐Gaussian settings