Abstract. Airborne electromagnetic (AEM) methods supply data over large areas in a cost-effective way. We used Artificial Neural Networks (ANN) to classify the geophysical signal into a meaningful geological parameter. By using examples of known relations between ground-based geophysical data (in this case electrical conductivity, EC, from electrical cone penetration tests) and geological parameters (presence of glacial till), we extracted learning rules that could be applied to map the presence of a glacial till using the EC profiles from the airborne EM data. The saline groundwater in the area was obscuring the EC signal from the till but by using ANN we were able to extract subtle and often non-linear, relations in EC that were representative of the presence of the till. The ANN results were interpreted as the probability of having till and showed a good agreement with drilling data. The glacial till is acting as a layer that inhibits groundwater flow, due to its high clay-content, and is therefore an important layer in hydrogeological modelling and for predicting the effects of climate change on groundwater quantity and quality.
Abstract. A direct current (DC) resistivity and time domain induced polarization (TDIP) survey was undertaken at a decommissioned landfill site situated in Hørløkke, Denmark, for the purpose of mapping the waste deposits and to discriminate important geological units that control the hydrology of the surrounding area. It is known that both waste deposits and clay have clear signatures in TDIP data, making it possible to enhance the resolution of geological structures compared to DC surveys alone. Four DC/TDIP profiles were carried out crossing the landfill, and another seven profiles in the surroundings provide a sufficiently dense coverage of the entire area. The whole dataset was inverted using a 1-D laterally constrained inversion scheme, recently implemented for TDIP data, in order to use the entire decay curves for reconstructing the electrical parameters of the soil in terms of the Cole-Cole polarization model. Results show that it is possible to resolve both the geometry of the buried waste body and key geological structures. In particular, it was possible to find a silt/clay lens at depth that correlates with the flow direction of the pollution plume spreading out from the landfill and to map a shallow sandy layer rich in clay that likely has a strong influence on the hydrology of the site. This interpretation of the geophysical findings was constrained by borehole data, in terms of geology and gamma ray logging. The results of this study are important for the impact of the resolved geological units on the hydrology of the area, making it possible to construct more realistic scenarios of the variation of the pollution plume as a function of the climate change.
ABSTRACT This study uses time‐domain induced polarization data for the delineation and characterization of the former landfill site at Eskelund, Denmark. With optimized acquisition parameters combined with a new inversion algorithm, we use the full content of the decay curve and retrieve spectral information from time‐domain IP data. Thirteen IP/DC profiles were collected in the area, supplemented by el‐log drilling for accurate correlation between the geophysics and the lithology. The data were inverted using a laterally constrained 1D inversion considering the full decay curves to retrieve the four Cole‐Cole parameters. For all profiles, the results reveal a highly chargeable unit that shows a very good agreement to the findings from 15 boreholes covering the area, where the extent of the waste deposits was measured. The thickness and depth of surface measurements were furthermore validated by el‐log measurements giving in situ values, for which the Cole‐Cole parameters were computed. The 3D shape of the waste body was pinpointed and well‐defined. The inversion of the IP data also shows a strong correlation with the initial stage of the waste dump and its composition combining an aerial map with acquired results.
The principle of equivalence is known to cause nonuniqueness in interpretations of direct current (DC) resistivity data. Low- or high-resistivity equivalences arise when a thin geologic layer with a low/high resistivity is embedded in a relative high-/low-resistivity background formation causing strong resistivity-thickness correlations. The equivalences often make it impossible to resolve embedded layers. We found that the equivalence problem could be significantly reduced by combining the DC data with full-decay time-domain induced polarization (IP) measurements. We applied a 1D Markov chain Monte Carlo algorithm to invert synthetic DC data of models with low- and high-resistivity equivalences. By applying this inversion method, it is possible to study the space of equivalent models that have an acceptable fit to the observed data, and to make a full sensitivity analysis of the model parameters. Then, we include a contrast in chargeability into the model, modeled in terms of spectral Cole-Cole IP parameters, and invert the DC and IP data in combination. The results show that the addition of IP data largely resolves the DC equivalences. Furthermore, we present a field example in which DC and IP data were measured on a sand formation with an embedded clay layer known from a borehole drilling. Inversion results show that the DC data alone do not resolve the clay layer due to equivalence problems, but by adding the IP data to the inversion, the layer is resolved.
Key to the effective management of natural ecosystems that characterise the floodplains of the Murray Basin in south eastern Australia, and maintenance of river health (a critical water resource) in a setting where severe salinisation is an ever-present threat, requires a sound understanding of surface water-groundwater processes. This paper presents results from an examination of hydrogeophysics, specifically airborne electromagnetics (AEM) data acquired by the SkyTEM time domain helicopter EM system, as a means for improving our knowledge of spatial patterns associated with inter-aquifer mixing where groundwater flow is complex. In the south-eastern part of the Murray Basin, AEM data shows considerable promise as a means for understanding of groundwater quality and its lateral variability. In the Bookpurnong and Loxton irrigation areas the high moment capability of SkyTEM permits us to investigate variations in the quality of groundwater at depth (>100m), which in turn allows us to visualise how groundwater may be moving across aquitards and within particular aquifer systems.
Adiabatic half-passage (AHP) pulses show great promise for significantly enhancing the signal-to-noise ratio of the surface nuclear magnetic resonance (NMR) free-induction decay measurement. Performing an AHP requires that the frequency sweep terminates when the transmit frequency is equal to the Larmor frequency, a condition that demands accurate knowledge of the true Larmor frequency. If the frequency sweep is terminated at an incorrect frequency, i.e., with an unknown offset between the transmit and Larmor frequency at the end of the pulse, the net excitation is affected and it can differ from that predicted by modeling that assumes a 0 Hz offset at the end of the sweep. Surface NMR surveys using a traditional single-frequency pulse have previously been shown to display degraded performance in the presence of an uncertain Larmor frequency estimate; the AHP pulse is also likely susceptible to such degraded performance. To ensure that reliable results can be produced by AHP pulses in the presence of an uncertain Larmor frequency estimate, we have developed an approach that adapts the frequency-cycling scheme for use with AHP pulses. We hypothesize that data collected using two similar AHP pulses, each with the exact same frequency sweep but where one sweeps toward the Larmor frequency from higher frequencies and the other from lower frequencies, can be stacked in such a manner that the impact of an unknown frequency offset is significantly reduced. We present synthetic and field results to demonstrate that frequency-cycling AHP pulse surface NMR data can ensure reliable performance even in the presence of an uncertain Larmor frequency estimate.
Approximately 50,000 ground based TEM soundings have been carried out in Denmark during<br>the last decade for groundwater investigations. This number will increase in the future and therefore we<br>have developed a new helicopter time-domain electromagnetic (TEM) system, SkyTEM.<br>In the development and design of the SkyTEM system it has been an unchangeable demand that<br>the data quality of the SkyTEM system should be the same or better than the data quality from ground<br>based systems obtained by e.g. the Protem 47 system (40 x 40 m central loop configuration). Because<br>the amounts of data produced from the system are very high, new concepts for the processing and<br>inversion of TEM data have been developed.<br>In this paper we discuss the processing of the data produced by the SkyTEM system – GPS data,<br>transmitter frame angle and altitude data, transmitter status parameters and the transient decays. We<br>conclude the paper by presenting the results from a 50 km2 large groundwater survey west of Århus<br>which demonstrates the high resolution capabilities of the SkyTEM system.
Abstract. A DC resistivity (DC) and Time Domain Induced Polarization (TDIP) survey was undertaken at a decommissioned landfill site situated in Hørløkke, Denmark, for the purpose of mapping the waste deposits and to discriminate important geological units that control the hydrology of the surrounding area. It is known that both waste deposits and clay have clear signatures in TDIP data, making possible to enhance the resolution of geological structures, when compared to DC surveys alone. Four DC/TDIP profiles were carried out crossing the landfill and another seven profiles in the surroundings, giving a dense coverage over the entire area. The whole dataset was inverted using a 1-D Laterally Constrained Inversion scheme, recently implemented for IP data, in order to use the entire decay curves for reconstructing the electrical parameters of the soil in terms of the Cole-Cole polarization model. Results show that it is possible to both resolve the geometry of the buried waste body and key geological structures. In particular, it was possible to find a silt/clay lens at depth, which correlates with the flow direction of the pollution plume spreading out from the landfill, and to map a shallow sandy layer rich in clay that likely has a strong influence on the hydrology of the site. This interpretation of the geophysical findings was constrained by boreholes data, in terms of geology and gamma ray logging. The results of this study are important for the impact that the resolved geological units have in the hydrology of the area, making it possible to construct more realistic scenarios of the variation of the pollution plume as a function of the climate change.
ABSTRACT Accurate modelling of the conductivity structure of mineralisations can often be difficult. In order to remedy this, a parametric approach is often used. We have developed a parametric thin‐sheet code, with a variable overburden. The code is capable of performing inversions of time‐domain airborne electromagnetic data, and it has been tested successfully on both synthetic data and field data. The code implements an integral solution containing one or more conductive sheets, buried in a half‐space with a laterally varying conductive overburden. This implementation increases the area of applicability compared to, for example, codes operating in free space, but it comes with a significant increase in computational cost. To minimise the cost, the code is parallelised using OpenMP and heavily optimised, which means that inversions of field data can be performed in hours on multiprocessor desktop computers. The code models the full system transfer function of the electromagnetic system, including variable flight height. The code is demonstrated with a synthetic example imitating a mineralisation buried underneath a conductive meadow. As a field example, the Valen mineral deposit, which is a graphite mineral deposit located in a variable overburden, is successfully inverted. Our results match well with previous models of the deposit; however, our predicted sheet remains inconclusive. These examples collectively demonstrate the effectiveness of our thin‐sheet code.