Search Results
Geophysical investigation of a freshwater lens on the island of Langeoog, Germany – Insights from combined HEM, TEM and MRS data
2017, Costabel, Stephan, Siemon, Bernhard, Houben, Georg, Günther, Thomas
A multi-method geophysical survey, including helicopter-borne electromagnetics (HEM), transient electromagnetics (TEM), and magnetic resonance sounding (MRS), was conducted to investigate a freshwater lens on the North Sea island of Langeoog, Germany. The HEM survey covers the entire island and gives an overview of the extent of three freshwater lenses that reach depths of up to 45 m. Ground-based TEM and MRS were conducted particularly on the managed western lens to verify the HEM results and to complement the lithological information from existing boreholes. The results of HEM and TEM are in good agreement. Salt- and freshwater-bearing sediments can, as expected, clearly be distinguished due to their individual resistivity ranges. In the resistivity data, a large transition zone between fresh- and saltwater with a thickness of up to 20 m is identified, the existence of which is verified by borehole logging and sampling. Regarding lithological characterisation of the subsurface, the MRS method provides more accurate and reliable results than HEM and TEM. Using a lithological index derived from MRS water content and relaxation time, thin aquitard structures as well as fine and coarse sand aquifers can be distinguished. Complementing the existing borehole data with the lithology information estimated from MRS, we generate a map showing the occurrence of aquitard structures, which significantly improves the hydrogeological model of the island. Moreover, we demonstrate that the estimates of groundwater conductivity in the sand aquifers from geophysical data are in agreement with the fluid conductivity measured in the boreholes.
Electric resistivity and seismic refraction tomography: a challenging joint underwater survey at Äspö Hard Rock Laboratory
2017, Ronczka, Mathias, Hellman, Kristofer, Günther, Thomas, Wisén, Roger, Dahlin, Torleif
Tunnelling below water passages is a challenging task in terms of planning, pre-investigation and construction. Fracture zones in the underlying bedrock lead to low rock quality and thus reduced stability. For natural reasons, they tend to be more frequent at water passages. Ground investigations that provide information on the subsurface are necessary prior to the construction phase, but these can be logistically difficult. Geophysics can help close the gaps between local point information by producing subsurface images. An approach that combines seismic refraction tomography and electrical resistivity tomography has been tested at the Äspö Hard Rock Laboratory (HRL). The aim was to detect fracture zones in a well-known but logistically challenging area from a measuring perspective. The presented surveys cover a water passage along part of a tunnel that connects surface facilities with an underground test laboratory. The tunnel is approximately 100 m below and 20 m east of the survey line and gives evidence for one major and several minor fracture zones. The geological and general test site conditions, e.g. with strong power line noise from the nearby nuclear power plant, are challenging for geophysical measurements. Co-located positions for seismic and ERT sensors and source positions are used on the 450 m underwater section of the 700 m profile. Because of a large transition zone that appeared in the ERT result and the missing coverage of the seismic data, fracture zones at the southern and northern parts of the underwater passage cannot be detected by separated inversion. Synthetic studies show that significant three-dimensional (3-D) artefacts occur in the ERT model that even exceed the positioning errors of underwater electrodes. The model coverage is closely connected to the resolution and can be used to display the model uncertainty by introducing thresholds to fade-out regions of medium and low resolution. A structural coupling cooperative inversion approach is able to image the northern fracture zone successfully. In addition, previously unknown sedimentary deposits with a significantly large thickness are detected in the otherwise unusually well-documented geological environment. The results significantly improve the imaging of some geologic features, which would have been undetected or misinterpreted otherwise, and combines the images by means of cluster analysis into a conceptual subsurface model.
Impacts of a capillary barrier on infiltration and subsurface stormflow in layered slope deposits monitored with 3-D ERT and hydrometric measurements
2017, Hübner, Rico, Günther, Thomas, Heller, Katja, Noell, Ursula, Kleber, Arno
Identifying principles of water movement in the shallow subsurface is crucial for adequate process-based hydrological models. Hillslopes are the essential interface for water movement in catchments. The shallow subsurface on slopes typically consists of different layers with varying characteristics. The aim of this study was to draw conclusions about the infiltration behaviour, to identify water flow pathways and derive some general interpretations for the validity of the water movement on a hillslope with periglacial slope deposits (cover beds), where the layers differ in their sedimentological and hydrological properties. Especially the described varying influence of the basal layer (LB) as an impeding layer on the one hand and as a remarkable pathway for rapid subsurface stormflow on the other. We used a time lapse 3-D electrical resistivity tomography (ERT) approach combined with punctual hydrometric data to trace the spreading and the progression of an irrigation plume in layered slope deposits during two irrigation experiments. This multi-technical approach enables us to connect the high spatial resolution of the 3-D ERT with the high temporal resolution of the hydrometric devices. Infiltration through the uppermost layer was dominated by preferential flow, whereas the water flow in the deeper layers was mainly matrix flow. Subsurface stormflow due to impeding characteristic of the underlying layer occurs in form of "organic layer interflow" and at the interface to the first basal layer (LB1). However, the main driving factor for subsurface stormflow is the formation of a capillary barrier at the interface to the second basal layer (LB2). The capillary barrier prevents water from entering the deeper layer under unsaturated conditions and diverts the seepage water according to the slope inclination. With higher saturation, the capillary barrier breaks down and water reaches the highly conductive deeper layer. This highlights the importance of the capillary barrier effect for the prevention or activation of different flow pathways under variable hydrological conditions.
Noninvasive Estimation of Water Retention Parameters by Observing the Capillary Fringe with Magnetic Resonance Sounding
2014, Costabel, Stephan, Günther, Thomas
The magnetic resonance sounding (MRS) method is usually applied for delineation and characterization of aquifer system stratification. Its unique property, distinct from other hydrogeophysical methods, is the direct sensitivity to water content in the subsurface. The inversion of MRS data yields the subsurface water content distribution without need of a petrophysical model. Recent developments in instrumentation, i.e., decreased instrumental dead times and advanced noise cancellation strategies, enable the use of this method for investigating the vadose zone. A possible way to interpret MRS measurements with focus on water retention (WR) parameters is an inversion approach that directly provides WR parameters by modeling the capillary fringe (CF inversion). We have developed this kind of inversion further to account for different WR models and present a sensitivity study based on both synthetic and real field data. To assess the general applicability of the CF inversion, we analyzed the resolution properties for different measurement layouts and the parameter uncertainties for different realistic scenarios. Under moderate noise conditions and if the water table position is known, all WR parameters except the residual water content can be reliably estimated. The relative accuracy of the estimated pore distribution index estimation is better for larger CF. Small measurement loops of 5-m diameter achieve the best resolution for shallow investigation depths of <10 m.
pyGIMLi: An open-source library for modelling and inversion in geophysics
2017, Rücker, Carsten, Günther, Thomas, Wagner, Florian M.
Many tasks in applied geosciences cannot be solved by single measurements, but require the integration of geophysical, geotechnical and hydrological methods. Numerical simulation techniques are essential both for planning and interpretation, as well as for the process understanding of modern geophysical methods. These trends encourage open, simple, and modern software architectures aiming at a uniform interface for interdisciplinary and flexible modelling and inversion approaches. We present pyGIMLi (Python Library for Inversion and Modelling in Geophysics), an open-source framework that provides tools for modelling and inversion of various geophysical but also hydrological methods. The modelling component supplies discretization management and the numerical basis for finite-element and finite-volume solvers in 1D, 2D and 3D on arbitrarily structured meshes. The generalized inversion framework solves the minimization problem with a Gauss-Newton algorithm for any physical forward operator and provides opportunities for uncertainty and resolution analyses. More general requirements, such as flexible regularization strategies, time-lapse processing and different sorts of coupling individual methods are provided independently of the actual methods used. The usage of pyGIMLi is first demonstrated by solving the steady-state heat equation, followed by a demonstration of more complex capabilities for the combination of different geophysical data sets. A fully coupled hydrogeophysical inversion of electrical resistivity tomography (ERT) data of a simulated tracer experiment is presented that allows to directly reconstruct the underlying hydraulic conductivity distribution of the aquifer. Another example demonstrates the improvement of jointly inverting ERT and ultrasonic data with respect to saturation by a new approach that incorporates petrophysical relations in the inversion. Potential applications of the presented framework are manifold and include time-lapse, constrained, joint, and coupled inversions of various geophysical and hydrological data sets.