We 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.

No abstract

We have developed a novel semiairborne frequency-domain electromagnetic (EM) system and successfully tested it within the DESMEX project. The semiairborne approach relies on the fact that part of the system is positioned on the ground and the rest is airborne. This allows us to take advantage of ground and airborne techniques. In particular, a high-moment transmitter can be installed on the earth’s surface, which enables us to inject and induce strong EM fields in the subsurface. Moreover, galvanic coupling is possible, which is an advantage if additional ground stations are deployed. The airborne receivers allow easier, significantly faster, and more uniform spatial coverage of the study area than the ground receivers. In our implementation, transmitters and electric field receivers are installed on the ground. Magnetic field sensors, such as commercially available fluxgate, total field magnetometers, and newly developed induction coils, are installed on a helicopter-towed bird. First, we describe the results of a semiairborne survey performed in a selected area with ancient mining located in the Saxothuringian zone near Schleiz, Germany. A 3D semiairborne inversion model represents several conductive anomalies, which agree well with the outcrop of alum shale formations at the surface. In addition, the shallow parts of the semiairborne model are compared with the result of an independent helicopter-borne survey, which consists of stepwise 1D models.

Transient electromagnetics (TEM) is a well-established method for mineral, groundwater, and geothermal exploration. Superconducting quantum interference device (SQUID)-based magnetic-field receivers used for TEM have quantitative advantages and higher sensitivity compared with commonly used induction coils. Special applications are deep soundings with target depths [Formula: see text] and settings with conductive overburden. However, SQUIDs have rarely been applied for TEM measurements in environments with significant anthropogenic noise. We compared a low-temperature SQUID with a commercially available induction coil in an area affected by anthropogenic noise. We acquired four fixed-loop data sets with totally 61 receiver stations close to Bad Frankenhausen, Germany. The high sensitivity of the SQUID enables low noise levels, which lead to longer high-quality transient data compared with the induction coil. The effect of anthropogenic and natural noise sources is more critical for the coil than for the SQUID data. In the vicinity of the transmitter loop, systematic distortion of the coil signals occurs at early times, most probably caused by sferic interferences. We have developed 1D inversion results of both receivers that matched well in general. However, the SQUID-based models were more consistent and showed greater depths of investigation. This led to a superior resolution of deeper layers and even enabled a potential detection of thin conducting targets at up to a 500 m depth. Moreover, we find that the SQUID data inversion revealed multidimensional effects within the conductive overburden. In this regard, we applied forward modeling to analyze systematic differences between inversion results of SQUID and coil data. We determine that low-temperature SQUIDs have the potential to significantly improve the reliability of subsurface models in suburban environments. Nevertheless, we recommend combined application of both types of receivers.