15 results
Search Results
Now showing 1 - 10 of 15
- ItemGeophysical investigation of a freshwater lens on the island of Langeoog, Germany – Insights from combined HEM, TEM and MRS data(Amsterdam [u.a.] : Elsevier Science, 2017) Costabel, Stephan; Siemon, Bernhard; Houben, Georg; Günther, ThomasA 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.
- ItemEvaluation of spectral induced polarization field measurements in time and frequency domain(Amsterdam [u.a.] : Elsevier Science, 2020) Martin, Tina; Günther, Thomas; Orozco, Adrian Flores; Dahlin, TorleifSpectral induced polarization (SIP) measurements have been demonstrated to correlate with important parameters in hydrogeological and environmental investigations. Although SIP measurements were often collected in the frequency domain (FDIP), recent developments have demonstrated the capabilities to solve for the frequency-dependence of the complex conductivity through measurements collected in the time domain (TDIP). Therefore, the aim of our field investigations is a comparison of the measured frequency-dependence at a broad frequency range resolved through FDIP and TDIP. In contrast to previous studies, we conducted measurements with different instruments and measuring technologies for both FDIP and TDIP. This allows for investigating the robustness of different measurements and assessing various sources of errors, for the assessment of the advantages and drawbacks from different measuring techniques. Our results demonstrate that data collected through different instruments are consistent. Apparent resistivity measurements as well as the inversion results revealed quantitatively the same values for all instruments. The measurements of the IP effect are also comparable, particularly FDIP readings in the low frequencies (< 10 Hz) revealed to be quantitatively the same for different instruments. TDIP measurements are consistent for data collected with both devices. As expected, the spatial distribution of the values is also consistent for low frequency data (in FDIP) and late times measurements in TDIP (> 0.1 s). However, data quality for higher frequencies in FDIP (i.e., early times in TDIP) show larger variations, which reflects the differences between the instruments to deal with the electromagnetic contamination of the IP data. Concluded in general, the different instruments and measuring techniques can provide consistent responses for varying signal-to-noise ratio and measuring configurations. © 2020 The Authors
- ItemComparison of novel semi-airborne electromagnetic data with multi-scale geophysical, petrophysical and geological data from Schleiz, Germany(Amsterdam [u.a.] : Elsevier Science, 2020) Steuer, Annika; Smirnova, Maria; Becken, Michael; Schiffler, Markus; Günther, Thomas; Rochlitz, Raphael; Yogeshwar, Pritam; Mörbe, Wiebke; Siemon, Bernhard; Costabel, Stephan; Preugschat, Benedikt; Ibs-von Seht, Malte; Zampa, Luigi Sante; Müller, FranzIn the framework of the Deep Electromagnetic Sounding for Mineral EXploration (DESMEX) project, we carried out multiple geophysical surveys from regional to local scales in a former mining area in the state of Thuringia, Germany. We prove the applicability of newly developed semi-airborne electromagnetic (EM) systems for mineral exploration by cross-validating inversion results with those of established airborne and ground-based investigation techniques. In addition, supporting petrophysical and geological information to our geophysical measurements allowed the synthesis of all datasets over multiple scales. An initial regional-scale reconnaissance survey was performed with BGR's standard helicopter-borne geophysical system deployed with frequency-domain electromagnetic (HEM), magnetic and radiometric sensors. In addition to geological considerations, the HEM results served as base-line information for the selection of an optimal location for the intermediate-scale semi-airborne EM experiments. The semi-airborne surveys utilized long grounded transmitters and two independent airborne receiver instruments: induction coil magnetometers and SQUID sensors. Due to the limited investigation depth of the HEM method, local-scale electrical resistivity tomography (ERT) and long-offset transient electromagnetic (LOTEM) measurements were carried out on a reference profile, enabling the validation of inversion results at greater depths. The comparison of all inversion results provided a consistent overall resistivity distribution. It further confirmed that both semi-airborne receiver instruments achieve the bandwidth and sensitivity required for the investigation of the resistivity structure down to 1 km depth and therewith the detection of deeply seated earth resources. A 3D geological model, lithological and geophysical borehole logs as well as petrophysical investigations were integrated to interpret of the geophysical results. Distinct highly-conductive anomalies with resistivities of less than 10 Om were identified as alum shales over all scales. Apart from that, the petrophysical investigations exhibited that correlating geophysical and geological information using only one single parameter, such as the electrical resistivity, is hardly possible. Therefore, we developed a first approach based on clustering methods and self-organizing maps (SOMs) that allowed us to assign geological units at the surface to a given combination of geophysical and petrophysical parameters, obtained on different scales. © 2020 The Authors
- ItemClassification of slag material by spectral induced polarization laboratory and field measurements(Amsterdam [u.a.] : Elsevier Science, 2021) Martin, Tina; Günther, Thomas; Weller, Andreas; Kuhn, KerstinHistorical slag dumps are of increasing interest due to economic, environmental or archaeological reasons. Geophysical investigations can help accessing the potential reuse of slag material to recover metallic raw material or for the estimation of the hazard potential of the buried slag material due to dissolution occurrence. In our study, we have investigated various slag material in the laboratory with the spectral induced polarization (SIP) method, obtained from different historical slag dumps, located in the Harz Mountains, Germany. We also present SIP results from field measurements at a historical slag dump where most of the slag samples reveal high amounts of iron, zinc, silica, and barium. Our results reveal a discrimination between three different slag grades (low, medium, high) by using the imaginary conductivity σ″ at a medium frequency (1–10 Hz) in both laboratory and field. Furthermore, additional information is obtained by a classification based on the spectral polarization behaviour and considering the field frequency range (0.1 Hz – 100 Hz). Five different types of spectra (ascending, descending, constant, maximum and minimum type) can be discriminated and recognized in the laboratory and in distinct areas of the slag dump. Even though a direct comparison between the laboratory and field results still needs to be proven, the buried slag material can be differentiated from the surrounding material by the polarization magnitude.
- ItemMonitoring freshwater–saltwater interfaces with SAMOS – installation effects on data and inversion(Oxford : Wiley, 2020) Ronczka, Mathias; Günther, Thomas; Grinat, Michael; Wiederhold, HelgaA major problem for the freshwater supply of coastal regions is the intrusion of saltwater into aquifers. Due to extensive extraction of freshwater to suffice increasing drinking water demands and/or in periods of reduced groundwater recharge, the equilibrium state may be disturbed. The result is an upconing or movement of the fresh–saline groundwater interface, which reduces the local drinking water resources at coastal regions or islands. The saltwater monitoring system (SAMOS) is a vertical electrode chain installed in a backfilled borehole. It provides a solution to observe the transition zone in detail, both temporally and spatially. We present monitoring data of the first year from three locations - with different geological conditions that show disturbances in the resistivity distribution that result from the drilling processes. A clayey backfilling, for example, can lead to beam-like artefacts, and a mixed fluid within the backfilling changes its bulk resistivity, both leading to misinterpretations. We performed data inversion under cylindrically symmetrical conditions in full-space in order to separate these resistivity artefacts from the undisturbed background. Data inversion reveals that it is possible to separate drilling effects on the resistivity distribution from the undisturbed background. Thus, an interpretation of the natural transition zones can be made immediately after the installation. © 2020 The Authors. Near Surface Geophysics published by John Wiley & Sons Ltd on behalf of European Association of Geoscientists and Engineers.
- ItemStructural style and neotectonic activity along the Harz Boundary Fault, northern German: a multimethod approach integrating geophysics, outcrop data and numerical simulations(Berlin ; Heidelberg : Springer, 2020) Müller, Katharina; Polom, Ulrich; Winsemann, Jutta; Steffen, Holger; Tsukamoto, Sumiko; Günther, Thomas; Igel, Jan; Spies, Thomas; Lege, Thomas; Frechen, Manfred; Franzke, Hans‑Joachim; Brandes, ChristianWe present new evidence for neotectonic activity along the Harz Boundary Fault, a Cretaceous reverse fault that represents a key structure in northern Germany. For the fault analysis, we use a multimethod approach, integrating outcrop data, luminescene dating, shear wave seismics, electrical resistivity tomography (ERT) and numerical simulations. A recent sinkhole at the SSW-ward dipping and WNW–ESE striking Harz Boundary Fault exposes a NNE-ward dipping and WNW–ESE striking planar fault surface that cuts through unconsolidated debris-flow deposits thus pointing to young Lateglacial tectonic activity. The fault shows a polyphase evolution with initial normal fault movement and a later reactivation as an oblique fault with reverse and strike-slip components. A shear wave seismic profile was acquired to analyse the geometry of the fault and show that the Harz Boundary Fault is steeply dipping and likely has branches. Partly, these branches propagate into overlying alluvial-fan deposits that are probably Pleniglacial to Lateglacial in age. The outcrop data in combination with the seismic data give evidence for a splay fault system with steep back-thrusts. One of these back-thrusts is most likely the NNE-ward dipping fault that is exposed in the sinkhole. The lateral extent of the fault was mapped with electrical resistivity tomography (ERT) profiles. The timing of fault movement was estimated based on optically stimulated luminescence dating of the faulted debris-flow deposits using both quartz and feldspar minerals. Consistent feldspar and quartz ages indicate a good bleaching of the sediment prior to deposition. The results imply fault movements post-dating ~ 15 ka. Numerical simulations of glacio isostatic adjustment (GIA)-related changes in the Coulomb failure stress regime at the Harz Boundary Fault underpin the assumption that the fault was reactivated during the Lateglacial due to stress changes induced by the decay of the Late Pleistocene (Weichselian) Fennoscandian ice sheet. © 2020, The Author(s).
- ItemA remote-control datalogger for large-scale resistivity surveys and robust processing of its signals using a software lock-in approach(Göttingen : Copernicus Publ., 2018) Oppermann, Frank; Günther, ThomasWe present a new versatile datalogger that can be used for a wide range of possible applications in geosciences. It is adjustable in signal strength and sampling frequency, battery saving and can remotely be controlled over a Global System for Mobile Communication (GSM) connection so that it saves running costs, particularly in monitoring experiments. The internet connection allows for checking functionality, controlling schedules and optimizing pre-amplification. We mainly use it for large-scale electrical resistivity tomography (ERT), where it independently registers voltage time series on three channels, while a square-wave current is injected. For the analysis of this time series we present a new approach that is based on the lock-in (LI) method, mainly known from electronic circuits. The method searches the working point (phase) using three different functions based on a mask signal, and determines the amplitude using a direct current (DC) correlation function. We use synthetic data with different types of noise to compare the new method with existing approaches, i.e. selective stacking and a modified fast Fourier transformation (FFT)-based approach that assumes a 1∕f noise characteristics. All methods give comparable results, but the LI is better than the well-established stacking method. The FFT approach can be even better but only if the noise strictly follows the assumed characteristics. If overshoots are present in the data, which is typical in the field, FFT performs worse even with good data, which is why we conclude that the new LI approach is the most robust solution. This is also proved by a field data set from a long 2-D ERT profile.
- ItemcustEM: Customizable finite-element simulation of complex controlled-source electromagnetic data(Tulsa, Okla. : SEG, 2019) Rochlitz, Raphael; Skibbe, Nico; Günther, ThomasWe have developed the open-source toolbox custEM (customizable electromagnetic modeling) for the simulation of complex 3D controlled-source electromagnetic (CSEM) problems. It is based on the open-source finite-element library FEniCS, which supports tetrahedral meshes, multiprocessing, higher order polynomials, and anisotropy. We use multiple finite-element approaches to solve the time-harmonic Maxwell equations, which are based on total or secondary electric field and gauged potential formulations. In addition, we develop a secondary magnetic field formulation, showing superior performance if only magnetic fields are required. Using Nédélec basis functions, we robustly incorporate the current density on the edges of the mesh for the total field formulations. The latter enable modeling of CSEM problems taking topography into account. We evaluate semianalytical 1D layered-earth solutions with the pyhed library, supporting arbitrary configurations of dipole or loop sources for secondary field calculations. All system matrices have been modified to be symmetric and solved in parallel with the direct solver MUMPS. Aside from the finite-element kernel, mesh generation, interpolation, and visualization modules have been implemented to simplify and automate the modeling workflow. We prove the capability of custEM, including validation against analytic-solutions, crossvalidation of all implemented approaches, and results for a model with 3D topography with four examples. The object-oriented implementation allows for customizable modifications and additions or to use only submodules designed for special tasks, such as mesh generation or matrix assembly. Therefore, the toolbox is suitable for crossvalidation with other codes and as the basis for developing 3D inversion routines.
- ItemImpacts of a capillary barrier on infiltration and subsurface stormflow in layered slope deposits monitored with 3-D ERT and hydrometric measurements(Munich : EGU, 2017) Hübner, Rico; Günther, Thomas; Heller, Katja; Noell, Ursula; Kleber, ArnoIdentifying 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.
- ItempyGIMLi: An open-source library for modelling and inversion in geophysics(Amsterdam [u.a.] : Elsevier Science, 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.