Abstract. Quaternary deposits are complex and heterogeneous. They contain some of the most abundant and extensively used aquifers. In order to improve the knowledge of the spatial heterogeneity of such deposits, we acquired a large (1500 ha) and dense (20 m spacing) time domain electromagnetic (TDEM) data set in the upper Aare Valley, Switzerland (available at https://doi.org/10.5281/zenodo.4269887; Neven et al., 2020). TDEM is a fast and reliable method to measure the magnetic field directly related to the resistivity of the underground. In this paper, we present the inverted resistivity models derived from this acquisition. The depth of investigation ranges between 40 and 120 m, with an average data residual contained in the standard deviation of the data. These data can be used for many different purposes: from sedimentological interpretation of quaternary environments in alpine environments, geological and hydrogeological modeling, to benchmarking geophysical inversion techniques.
Abstract. Ground-penetrating radar (GPR) is widely used for determining mountain glacier thickness. However, this method provides thickness data only along the acquisition lines, and therefore interpolation has to be made between them. Depending on the interpolation strategy, calculated ice volumes can differ and can lack an accurate error estimation. Furthermore, glacial basal topography is often characterized by complex geomorphological features, which can be hard to reproduce using classical interpolation methods, especially when the field data are sparse or when the morphological features are too complex. This study investigates the applicability of multiple-point statistics (MPS) simulations to interpolate glacier bedrock topography using GPR measurements. In 2018, a dense GPR data set was acquired on the Tsanfleuron Glacier (Switzerland). These data were used as the source for a bedrock interpolation. The results obtained with the direct-sampling MPS method are compared against those obtained with kriging and sequential Gaussian simulations (SGSs) on both a synthetic data set – with known reference volume and bedrock topography – and the real data underlying the Tsanfleuron Glacier. Using the MPS modeled bedrock, the ice volume for the Scex Rouge and Tsanfleuron glaciers is estimated to be 113.9 ± 1.6 million cubic meters. The direct-sampling approach, unlike the SGS and kriging, allowed not only an accurate volume estimation but also the generation of a set of realistic bedrock simulations. The complex karstic geomorphological features are reproduced and can be used to significantly improve for example the precision of subglacial flow estimation.
In Quaternary deposits, the characterization of subsurface heterogeneity and its associated uncertainty is critical when dealing with groundwater resource management. The combination of different data types through joint inversion has proven to be an effective way to reduce final model uncertainty. Moreover, it allows the final model to be in agreement with a wider spectrum of data available on site. However, integrating them stochastically through an inversion is very time-consuming and resource expensive, due to the important number of physical simulations needed. The use of multi-fidelity models, by combining low-fidelity inexpensive and less accurate models with high-fidelity expensive and accurate models, allows one to reduce the time needed for inversion to converge. This multiscale logic can be applied for the generation of Quaternary models. Most Quaternary sedimentological models can be considered as geological units (large scale), populated with facies (medium scale), and finally completed by physical parameters (small scale). In this paper, both approaches are combined. A simple and fast time-domain EM 1D geophysical direct problem is used to first constrain a simplified stochastic geologically consistent model, where each stratigraphic unit is considered homogeneous in terms of facies and parameters. The ensemble smoother with multiple data assimilation (ES-MDA) algorithm allows generating an ensemble of plausible subsurface realizations. Fast identification of the large-scale structures is the main point of this step. Once plausible unit models are generated, high-fidelity transient groundwater flow models are incorporated. The low-fidelity models are populated stochastically with heterogeneous facies and their associated parameter distribution. ES-MDA is also used for this task by directly inferring the property values (hydraulic conductivity and resistivity) from the generated model. To preserve consistency, geophysical and hydrogeological data are inverted jointly. This workflow ensures that the models are geologically consistent and are therefore less subject to artifacts due to localized poor-quality data. It is able to robustly estimate the associated uncertainty with the final model. Finally, due to the simplification of both the direct problem and the geology during the low-fidelity part of the inversion, it greatly reduces the time required to converge to an ensemble of complex models while preserving consistency.
In most urbanized and agricultural areas of central Europe, the shallow underground is constituted of Quaternary deposits which are often the most extensively used layers (water pumping, shallow geothermic, material excavation). All these deposits are often complexly intertwined, leading to high spatial variability and high complexity. Geophysical data can be a fast and reliable source of information about the underground. Still, the integration of these data can be time-consuming, it lacks realistic interpolation in a full 3D space, and the final uncertainty is often not represented. In this study, we propose a new methodology to combine boreholes and geophysical data with uncertainty in an automatic framework. A spatially varying translator function that predicts the clay fraction from resistivity is inverted using boreholes description as control points. It is combined with a 3D stochastic interpolation framework based on a Multiple Points Statistics algorithm and Gaussian Random Function. This novel workflow allows incorporating robustly the data and their uncertainty and requires less user intervention than the already existing workflows. The methodology is illustrated for ground-based towed transient electromagnetic data (tTEM) and borehole data from the upper Aare valley, Switzerland. In this location, a 3D realistic high spatial resolution model of clay fraction was obtained over the whole valley. The very dense data set allowed to demonstrate the quality of the predicted values and their corresponding uncertainties using cross-validation.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.