Frequency-domain Electromagnetic Measurements Research Articles

A comparison between Kalman ensemble generator and geostatistical frequency-domain electromagnetic inversion: The impacts on near-surface characterization

The spatial distribution of the physical properties of the first meters beneath the earth’s surface is often complex due to its highly dynamic nature and small-scale heterogeneities resulting from natural and anthropogenic processes. Therefore, obtaining numerical 3D models that accurately describe the spatial distribution of these properties is often challenging yet essential for different fields such as environmental assessment and remediation, geoarchaeological conservation, and precision agriculture. Frequency-domain electromagnetic (FDEM) induction methods have proven their potential to image these properties in high (spatial) detail because FDEM measurements are sensitive to two key soil properties: electrical conductivity and magnetic susceptibility. Predicting subsurface properties from FDEM data requires solving an ill-posed and nonlinear inverse problem with multiple solutions. Recently, there has been a rapid growth of FDEM inversion methods, which may be broadly divided into probabilistic and deterministic methods. We compare two stochastic FDEM inversion approaches: the Kalman ensemble generator (KEG) and another one formulated as an iterative geostatistical FDEM inversion. Both methods are applied to a synthetic data set with spatially heterogeneous physical properties of interest, mimicking a real landfill mining site. The predicted models are compared with the reference models in terms of histogram and variogram models reproduction and in their ability to quantify spatial uncertainty. The results indicate the ability of both methods to predict the reference values. Although the KEG is computationally efficient, it struggles to reproduce the extreme values. In contrast, the geostatistical inversion approach ensures the reproduction of the imposed histograms and variogram models in the predicted models. As the prior information is included in both inversion methods in different ways, the pointwise variance models computed from all of the posterior models have different information. The synthetic data set is available to the community, so it can be used as a benchmark for other FDEM inversion methods.

Understanding Ancient Landscapes in the Venetian Plain through an Integrated Geoarchaeological and Geophysical Approach

This paper reports the results of the multidisciplinary study carried out in the SE area of Ceggia, in the eastern part of the Venetian Plain. The area has been characterized, since ancient times, by numerous morphological transformation, due to the presence of lagoon and marshes, and interested by repeated reclamation. Aerial and satellite images have identified many natural and anthropogenic traces. From a geophysical point of view, electrical resistivity tomography (ERT) combined with frequency-domain electromagnetic measurements (FDEM) can help to discriminate the spatial distribution of different buried structures in conductive systems. The electrical conductivity is, in fact, directly related to the soil moisture content. The multidisciplinary approach adopted in this context, with the results obtained thanks to the contribution of aerial and satellite images, historical cartography, archaeological survey, geophysical measurements, geomorphological characterization, and 14C dating, allow us to suggest a possible interpretation of the different traces highlighted in the studied area. This approach suggests a potentially useful and replicable methodology to study similar evidence, such as along the North Adriatic coast and in broad sectors of the Po Valley. The key issue, in this kind of system, lies, in fact, in the possibility to date and compare traces visible on the surface by remote sensing, establishing their interest from an archaeological and geomorphological point of view using an integration of field measurements. At the end of this research, the classification of the different anomalies found in this hydraulic variable context, thanks to the multidisciplinary approach here adopted, suggest new hypotheses for reading the complex history of this understudied area.

Investigating the Earth Fill Embankment of the Lotsane Dam for Internal Defects Using Time-lapse Resistivity Imaging and Frequency Domain Electromagnetics

We investigated the internal structure of the Lotsane Dam for zones that may be prone to seepage and internal erosion using the electrical resistivity imaging (ERI) and the frequency domain electromagnetic (FDEM) methods. Time-lapse ERI measurements were also made for a period of 8 months in order to monitor the temporal evolution of defective zones. Results from both the FDEM and ERI measurements show two main layers. The first is an upper conductive layer varying in thickness from 10 to 25 m which is related to the clay core embankment. Situated beneath this upper conductive layer is a highly resistive crystalline basement on which the dam was founded. Furthermore, the ERI and FDEM measurements revealed the presence of fractures and possible zones of weakness within the dam foundation. Time-lapse ERI measurements revealed resistivity increases in the observed possible defective zones, including proximal to the spillway and at the embankment-foundation interface. The long-term resistivity variation may be indicating change in material properties in those portions of the dam, and may evolve to destabilize the structural integrity of the dam and or develop into preferential seepage pathways with time. The identified anomalous zones are good indicators that the embankment integrity is at risk and we suggest continuous geophysical monitoring of Lotsane Dam structure in order to ensure dam safety and integrity on the long term.

Remote Sensing, Archaeological, and Geophysical Data to Study the Terramare Settlements: The Case Study of Fondo Paviani (Northern Italy)

During the Middle and Recent Bronze Age, the Po Plain and, more broadly Northern Italy were populated by the so-called “Terramare”, embanked settlements, surrounded by a moat. The buried remains of these archaeological settlements are characterized by the presence of a system of palaeo-environments and a consequent natural gradient in soil moisture content. These differences in the soil are often firstly detectable on the surface during the seasonal variations, with aerial, satellite, and Laser Imaging Detection and Ranging (LIDAR) images, without any information on the lateral and in-depth extension of the related buried structures. The variation in the moisture content of soils is directly related to their differences in electrical conductivity. Electrical resistivity tomography (ERT) and frequency domain electromagnetic (FDEM), also known as electromagnetic induction (EMI) measurements, provide non-direct measurements of electrical conductivity in the soils, helping in the reconstruction of the geometry of different buried structures. This study presents the results of the multidisciplinary approach adopted to the study of the Terramare settlement of Fondo Paviani in Northern Italy. Remote sensing and archaeological data, collected over about 10 years, combined with more recent ERT and FDEM measurements, contributed to the analysis of this particular, not yet wholly investigated, archaeological site. The results obtained by the integrated multidisciplinary study here adopted, provide new useful, interesting information for the archaeologists also suggesting future strategies for new studies still to be conducted around this important settlement.

Assessment of spatial representativity of X-ray tomography to study Vertical Flow Treatment wetlands

French Vertical Flow (VF) treatment wetlands receive raw wastewater and provide simultaneous sludge and wastewater treatment. For proper sludge handling, the treatment wetland must be designed adequately and specific operational conditions must be maintained. When these conditions are not met, accumulation of biosolids may lead to clogging. Filtration in French VF Treatment wetlands is governed by mechanisms at the pore-scale. They must be better understood to predict reliably biosolid accumulation. X-ray Computed Tomography (Xray-CT) is a promising technique to characterize in detail the morphology of the filtering media in treatment wetlands. In order to set a solid basis for the use of Xray-CT, the spatial representativity of measurements must be assessed. This issue is addressed in this study by successively analyzing spatial properties at the filter scale using Frequency Domain Electromagnetic Measurements (FDEMs), and at the pore scale using Xray-CT. A map of the electric conductivity at the surface of a French VF Treatment wetland is obtained by FDEM that indicates a homogeneous distribution of biosolids to which electrical conductivity is highly correlated. Different morphological properties were computed from Xray-CT after phase segmentation: phase volume fraction profiles, Specific Surface Area profiles and pore size distributions. Samples show several similarities of pore scale properties obtained by Xray-CT independently of the sampling region and especially the same vertical gradients. FDEM measurements and Xray-CT analysis are in agreement to indicate a good influent distribution at the surface of a full-scale mature French VF Treatment wetland. A criterion to define the limits of the deposit layer and gravel layer is introduced. This division allows to compare layers independently. Finally, a 2D-REV analysis suggests that the selected sample diameter of 5 cm is large enough to be representative of the heterogeneous distribution of phases at the pore-scale as long as no Phragmites are present.

Integrated application of 2D resistivity and electromagnetic methods to investigate a metallic-sulfide deposit in Soap Gulch, Montana. A case study

Geoelectrical methods involving electrical resistivity tomography (ERT), self-potential (SP), frequency domain electromagnetic (FDEM), and very low frequency (VLF) methods have been used to provide valuable information in locating a known sulfide ore body in Soap Gulch, Montana. The study develops basis of comparison for the geophysical techniques employed. Ranges of resistivity along the area have been established using interpreted ERT which can help to understand the subsurface distribution of sulfides in the area. A sulfide body was delineated from the survey area corresponding to anomalously low resistivity values on the ERT section, negative SP, and high apparent current density zone in VLF. Depth to the localized ore zone ranges approximately from 10 to 20 m. FDEM data reflect the conductivity distribution of the shallow subsurface (less than 6 m deep); hence, the delineated sulfide zone had minimal contribution to FDEM measurements. The results of the study show that SP, VLF, and ERT methods provide significant information in localizing ore bodies. The survey revealed that the resistivity values obtained from ERT profile corroborate the FDEM, SP, and VLF from the area.

Toward real-time three-dimensional mapping of surficial aquifers using a hybrid modeling approach

A hybrid modeling approach is proposed for near real-time three-dimensional (3D) mapping of surficial aquifers. First, airborne frequency-domain electromagnetic (FDEM) measurements are numerically inverted to obtain subsurface resistivities. Second, a machine-learning (ML) algorithm is trained using the FDEM measurements and inverted resistivity profiles, and borehole geophysical and hydrogeologic data. Third, the trained ML algorithm is used together with independent FDEM measurements to map the spatial distribution of the aquifer system. Efficacy of the hybrid approach is demonstrated for mapping a heterogeneous surficial aquifer and confining unit in northwestern Nebraska, USA. For this case, independent performance testing reveals that aquifer mapping is unbiased with a strong correlation (0.94) among numerically inverted and ML-estimated binary (clay-silt or sand-gravel) layer resistivities (5–20 ohm-m or 21–5,000 ohm-m), and an intermediate correlation (0.74) for heterogeneous (clay, silt, sand, gravel) layer resistivities (5–5,000 ohm-m). Reduced correlation for the heterogeneous model is attributed to over-estimating the under-sampled high-resistivity gravels (about 0.5 % of the training data), and when removed the correlation increases (0.87). Independent analysis of the numerically inverted and ML-estimated resistivities finds that the hybrid procedure preserves both univariate and spatial statistics for each layer. Following training, the algorithms can map 3D surficial aquifers as fast as leveled FDEM measurements are presented to the ML network.

Frequency domain electromagnetic and ground penetrating radar investigation of ephemeral streams: case study near the Southern High Plains, Texas

A frequency domain electromagnetic (FDEM) and ground penetrating radar (GPR) study was conducted on an ephemeral stream in north-central Texas to determine if FDEM and GPR measurements can be combined to determine the electrical characteristics of current and ancient stream channels. GPR data were collected at several frequencies to image sedimentary structures of different scale lengths, and to determine the formation porosity and water content of stream sediments. FDEM measurements were collected using Geonics EM31 and EM34 loop–loop instruments on a profile along the current stream channel and five profiles perpendicular to the channel. The results indicate that the greater spatial resolution of the EM31 mapped the current and possible ancient channels better than the EM34, however, the EM34 provided depth information on the formations underlying the channel sediments that the EM31 could not image. GPR measurements taken along a point bar deposit with 200, 100 and 25 MHz antennae indicated that the higher frequency antenna better resolved channel structures including laminar bedding, trough scours and cross-bedding, however, lower frequency antenna (25 MHz) imaged sedimentary structures within the underlying channel sediments. Common midpoint GPR measurements collected along the point bar deposit were used to estimate the sediment formation porosity (26%) and the water content, during a dry period, of the unsaturated (12%) and saturated (26 or 100% of the pore space) sediments. The combined results indicated that the FDEM data should be collected first because of the speed and ease of measurements. The FDEM data (especially the EM31) pointed to the locations of possible sedimentary structures, which can then be resolved by using different frequency GPR measurements.