Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.256 IF 4.256
  • IF 5-year value: 4.819 IF 5-year 4.819
  • CiteScore value: 4.10 CiteScore 4.10
  • SNIP value: 1.412 SNIP 1.412
  • SJR value: 2.023 SJR 2.023
  • IPP value: 3.97 IPP 3.97
  • h5-index value: 58 h5-index 58
  • Scimago H index value: 99 Scimago H index 99
Volume 22, issue 10 | Copyright
Hydrol. Earth Syst. Sci., 22, 5485-5508, 2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 24 Oct 2018

Research article | 24 Oct 2018

Contributions to uncertainty related to hydrostratigraphic modeling using multiple-point statistics

Adrian A. S. Barfod1,2, Troels N. Vilhelmsen2, Flemming Jørgensen1, Anders V. Christiansen2, Anne-Sophie Høyer1, Julien Straubhaar3, and Ingelise Møller1 Adrian A. S. Barfod et al.
  • 1Department of Groundwater and Quaternary Geological Mapping, Geological Survey of Denmark & Greenland (GEUS), C.F. Møllers Allé 8, 8000 Aarhus C, Denmark
  • 2Hydrogeophysics Group, Department of Geoscience, Aarhus University, C.F. Møllers Allé 4, 8000 Aarhus C, Denmark
  • 3Centre d'Hydrogéologie et de Géothermie (CHYN), Université de Neuchâtel, Neuchâtel, Switzerland

Abstract. Forecasting the flow of groundwater requires a hydrostratigraphic model, which describes the architecture of the subsurface. State-of-the-art multiple-point statistical (MPS) tools are readily available for creating models depicting subsurface geology. We present a study of the impact of key parameters related to stochastic MPS simulation of a real-world hydrogeophysical dataset from Kasted, Denmark, using the snesim algorithm. The goal is to study how changes to the underlying datasets propagate into the hydrostratigraphic realizations when using MPS for stochastic modeling. This study focuses on the sensitivity of the MPS realizations to the geophysical soft data, borehole lithology logs, and the training image (TI). The modeling approach used in this paper utilizes a cognitive geological model as a TI to simulate ensemble hydrostratigraphic models. The target model contains three overall hydrostratigraphic categories, and the MPS realizations are compared visually as well as quantitatively using mathematical measures of similarity. The quantitative similarity analysis is carried out exhaustively, and realizations are compared with each other as well as with the cognitive geological model.

The results underline the importance of geophysical data for constraining MPS simulations. Relying only on borehole data and the conceptual geology, or TI, results in a significant increase in realization uncertainty. The airborne transient electromagnetic SkyTEM data used in this study cover a large portion of the Kasted model area and are essential to the hydrostratigraphic architecture. On the other hand, the borehole lithology logs are sparser, and 410 boreholes were present in this study. The borehole lithology logs infer local changes in the immediate vicinity of the boreholes, thus, in areas with a high degree of geological heterogeneity, boreholes only provide limited large-scale structural information. Lithological information is, however, important for the interpretation of the geophysical responses. The importance of the TI was also studied. An example was presented where an alternative geological model from a neighboring area was used to simulate hydrostratigraphic models. It was shown that as long as the geological settings are similar in nature, the realizations, although different, still reflect the hydrostratigraphic architecture. If a TI containing a biased geological conceptualization is used, the resulting realizations will resemble the TI and contain less structure in particular areas, where the soft data show almost even probability to two or all three of the hydrostratigraphic units.

Publications Copernicus
Short summary
The focus of this study is on the uncertainty related to using multiple-point statistics (MPS) for stochastic modeling of the upper 200 m of the subsurface. The main research goal is to showcase how MPS methods can be used on real-world hydrogeophysical data and show how the uncertainty related to changing the underlying MPS setup propagates into the finalized 3-D subsurface models.
The focus of this study is on the uncertainty related to using multiple-point statistics (MPS)...