Magnetotelluric Measurements Research Articles

The structure of a continental intraplate volcanic system and controls from shear zones: Insights into the central Hoggar Cenozoic volcanic province, Northwest Africa, from electrical resistivity images

Continental intraplate volcanic systems, with their locations far from plate tectonic boundaries, are not well understood: the crustal and lithospheric mantle structure of these systems remain enigmatic and there is no consensus on the mechanisms that cause melt generation and ascent. The Cenozoic saw the development of numerous volcanic provinces on the African plate, including within the Central Hoggar, located in Northwest Africa, part of the Tuareg shield. The magmatic activity began at approximately 34 Ma and continued throughout the Quaternary. In order to understand the origins and potential mechanisms that generated the intraplate volcanic activity in the Central Hoggar we aim to image the subsurface architecture, in terms of electrical resistivity, from the surface to the lithospheric mantle. To do so we use magnetotelluric measurements from 40 locations to generate a 3-D electrical resistivity model, over an area of about 100 km by 160 km. Low-resistivity features (i.e., conductors) are observed in the crust that are narrow, linear structures congruent with the boundaries of terranes and prominent fault zones (e.g., Azrou N’Fad). They likely reflect the Pan-African mega-shear zones, which were reactivated throughout the tectonic evolution of the region. The model reveals that these faults are lithospheric-scale. The low-resistivity features likely represent the signatures of past fluid pathways and mineralization. A deeper low-resistivity feature is observed in the upper lithospheric mantle directly beneath the Manzaz and Atakor volcanic districts. It may represent local, small-scale metasomatism of the sub-continental lithospheric mantle, and low-percent partial melting, that sits above a regional, large-scale asthenospheric upwelling associated with the Hoggar swell. It is likely the origin point of the fluids responsible for the overlying anomalies. The results highlight the control of the lithospheric-scale, mega-shear zones on the spatial distribution of the recent Cenozoic volcanic activity, which was influenced by the location of pre-existing structural weaknesses.

Crustal conductivity footprint of the Miocene porphyry copper polymetallic deposits in the Gangdese metallogenic belt, Tibetan Plateau

New evidence worldwide has linked the surface locations of mineral deposits and their crustal-scale electrical conductivity footprint. We examine the relationship between the Gangdese Miocene porphyry copper deposits, Tibetan Plateau, and the electrical conductivity signature from a three-dimensional model generated from 311 magnetotelluric measurements. The distribution of electrical resistivity throughout the crust and the conductance within the mid-lower crust (depth range of 25–70 km) is analyzed. The results clearly show that the large and ultra-large Miocene porphyry copper deposits coincide spatially with conductive zones and areas of very-high conductance (>10,000S) in the mid-lower crust. Computations are undertaken to determine the influence of water-bearing silicate melts and alkali-bearing (Na+ and K+) fluids on conductivity. Based on this, the bulk conductivity is interpretated to be caused by a system of alkali-rich volatile-rich partial melt. The alkali-rich volatile-rich magmatic-hydrothermal fluids facilitate the migration and concentration of metal ions originating in deep areas. The volumes necessary are much less than partial melt alone and can thus help to reconcile large conductivity variations with small seismic velocity variations. The electrical structure indicates the magma source area of anatexis in the lower crust, a multi-stage magmatic system with large mid-crustal and small upper-crustal magma reservoirs, and complex pathways related to rift zones. We determine that the conductive zones in the mid-lower crust have an influence on the development of the mineralization and the location of the mineral belt.

The Pannon LitH2Oscope magnetotelluric array in the Pannonian Basin

The Pannonian Basin is one of the best natural laboratories in the world to study the lithospheric response to continental extension and subsequent tectonic inversion. Here we address the topic of lithospheric structure by a combined geochemical and magnetotelluric analysis, which has been carried out in the framework of the Pannon LitH2Oscope project. The main objective was to detect the resistivity distribution over the entire lithosphere by magnetotelluric measurements, considering the lithological resistivity properties and relate the results to the structure and evolution of the Pannonian Basin. The Pannon LitH2Oscope MT array was used to estimate the depth of the Lithosphere-Asthenosphere Boundary (LAB), considering the legacy MT data and compared to previous estimates for the region. Using the MT and geomagnetic response functions, major structural zones of the Pannonian basin, such as the Mid-Hungarian Shear Zone or fault systems like the Makó Trough and the Békés Basin, were also imaged. In addition, we used the apparent resistivity soundings to compare 1D resistivity models computed from geochemistry and obtained from field MT measurements. This comparison provided new constrains for the composition, fluid and melt content variations at the local lithosphere-asthenosphere boundary. The Pannon LitH2Oscope MT dataset and the results presented in this paper provide input for more complex 3D inversions and further investigations of the lithospheric structure in the Carpathian-Pannonian region.

Validating a UK Geomagnetically Induced Current Model Using Differential Magnetometer Measurements

AbstractExtreme space weather can damage ground‐based infrastructure such as power lines, railways and gas pipelines through geomagnetically induced currents (GICs). Modeling GICs requires knowledge about the source magnetic field and the electrical conductivity structure of the Earth to calculate ground electric fields during enhanced geomagnetic activity. The electric field, in combination with detailed information about the power grid topology, enable the modeling of GICs in high‐voltage (HV) power lines. Directly monitoring GICs in substations is possible with a Hall probe, but scarcely realized in the UK. Therefore we deployed the differential magnetometer method (DMM) to measure GICs at 12 sites in the UK power grid. The DMM includes the installation of two fluxgate magnetometers, one directly under a power line affected by GICs, and one as a remote site. The difference in recordings of the magnetic field at each instrument yields an estimate of the GICs in the respective power line segment via the Biot‐Savart law. We collected data across the UK in 2018–2022, monitoring HV line segments where previous research indicated high GIC risk. We recorded magnetometer data during several smaller storms that allow detailed analysis of our GIC model. For the ground electric field computations we used recent magnetotelluric (MT) measurements recorded close to the DMM sites. Our results show that there is strong agreement in both amplitude and signal shape between measured and modeled line and substation GICs when using our HV model and the realistic electric field estimates derived from MT data.

Physics-Informed Deep Learning Inversion with Application to Noisy Magnetotelluric Measurements

Despite demonstrating exceptional inversion production for synthetic data, the application of deep learning (DL) inversion methods to invert realistic magnetotelluric (MT) measurements, which are inevitably contaminated by noise in acquisition, poses a significant challenge. Hence, to facilitate DL inversion for realistic MT measurements, this work explores developing a noise-robust MT DL inversion method by generating targeted noisy training datasets and constructing a physics-informed neural network. Different from most previous works that only considered the noise of one fixed distribution and level, we propose three noise injection strategies and compare their combinations to mitigate the adverse effect of measurement noise on MT DL inversion results: (1) add synthetic relative noise obeying Gaussian distribution; (2) propose a multiwindow Savitzky–Golay (MWSG) filtering scheme to extract potential and possible noise from the target field data and then introduce them into training data; (3) create an augmented training dataset based on the former two strategies. Moreover, we employ the powerful Swin Transformer as the backbone network to construct a U-shaped DL model (SwinTUNet), based on which a physics-informed SwinTUNet (PISwinTUNet) is implemented to further enhance its generalization ability. In synthetic examples, the proposed noise injection strategies demonstrate impressive inversion effects, regardless of whether they are contaminated by familiar or unfamiliar noise. In a field example, the combination of three strategies drives PISwinTUNet to produce considerably faithful reconstructions for subsurface resistivity structures and outperform the classical deterministic Occam inversions. The experimental results show that the proposed noise-robust DL inversion method based on the noise injection strategies and physics-informed DL architecture holds great promise in processing MT field data.

Caprock integrity at Çanakkale-Tuzla hydrothermal system inferred from magnetotelluric modeling using particle swarm optimization

The impermeable caprock within a geothermal system serves the purpose of effectively sealing the reservoir, resulting in the elevation of both pressure and temperature. This sealing mechanism plays a crucial role in the long-term preservation of the system, while also contributing to its overall sustainability. Caprock failure subsequent to seismic activity near a geothermal site can lead to the permeation of the caprock structure, resulting in diminished sealing capabilities and a decline in the reservoir temperature. Additionally, this process alters the geochemical composition of the water by creating a hydrothermal mixture zone that disrupts the resistivity structure of the caprock, which is typically characterized by low resistivity values due to its substantial clay content and mineral alteration. This study focuses on investigating the integrity of the caprock at ơnakkale-Tuzla geothermal field in Turkey, where the water temperature and conductivity were reported to have decreased after a moderate magnitude earthquake and subsequent aftershocks. For this purpose, we performed magnetotelluric (MT) measurements, a method known for its sensitivity to geochemical reactions. These measurements were conducted along two parallel profiles that encompassed a total of 32 stations. Particle swarm optimization (PSO) technique was employed to overcome subtle difficulties associated with conventional inversion methods when modeling the MT data of complex formations. This is the first study that overcomes the difficulties emanating from the caprock failure by modeling MT data using PSO. Our modeling approach produced resistivity images that we interpreted as the signature of the failed caprock following the earthquake at the study site. Our results appear to confirm the documented geochemical changes or hydrothermal mixture zone about caprock structure.

Influence of Tectonic and Geological Structure on GIC in Southern South Island, New Zealand

AbstractAs part of a 5‐year project to assess the risk posed by geomagnetically induced currents (GIC) to the New Zealand electrical transmission network, long‐period magnetotelluric (MT) measurements have been made at 62 sites in southern South Island of New Zealand, a region where there was an absence of previous MT data. The data are largely 3‐dimensional in character, but show distinct features that can be related to the known tectonic and geological structure. In this work we focus on how the measured MT impedance tensors, and a simple interpretation of conductivity structure, can be used to assess the influence of tectonic and geological structure on GIC. We use the impedance tensors to calculate the magnitudes and orientations of induced electric fields in response to various orientations of inducing magnetic field. The electric fields so calculated are then used in a simplified model of the transmission network to calculate GIC at grounded substations. Our results confirm that tectonic/geological structure in the lower South Island and the resulting electrical conductivity variations have important impacts on the GIC magnitude. In the south‐west, smaller induced electric fields, associated with the higher conductivity in that region, lead to much reduced GIC at a substation in that area. In contrast, higher electric fields occurring in a NW‐SE band across the center of the region, contribute to much larger GIC in Dunedin city. Our results thus help explain the observed GIC reported at transformers in the region.

Geothermal Resource Mapping in Northern Botswana Inferred from Three-Dimensional Magnetotelluric Inversion

A set of magnetotelluric (MT) data collected from 28 stations at the Kasane Hot Spring in northern Botswana was used to derive a 3-D model of the electrical resistivity distribution around deep geologic structures that we associated with geothermal resources. The dimensionality analysis results revealed that 3-D modeling should be used to determine the electrical structure. The resistivity models revealed a conductive layer with an average thickness of 200 m, representing the overlying sediments of the Proterozoic volcanic rocks. The thick high-resistivity zone (>100 Ωm) below the conductive layer can be associated with Mesoproterozoic bedrock. The MT measurements in this area show a tube-shaped conductive anomaly that could serve as a fluid pathway feeding the hot spring. A fracture-controlled meteoric fluid circulation presumably determines the existence of the Kasane Hot Spring system.

Application of Electromagnetic Methods for Reservoir Monitoring with Emphasis on Carbon Capture, Utilization, and Storage

The Controlled-Source ElectroMagnetic (CSEM) method provides crucial information about reservoir fluids and their spatial distribution. Carbon dioxide (CO2) storage, enhanced oil recovery (EOR), geothermal exploration, and lithium exploration are ideal applications for the CSEM method. The versatility of CSEM permits its customization to specific reservoir objectives by selecting the appropriate components of a multi-component system. To effectively tailor the CSEM approach, it is essential to determine whether the primary target reservoir is resistive or conductive. This task is relatively straightforward in CO2 monitoring, where the injected fluid is resistive. However, for scenarios involving brine-saturated (water-wet) or oil-wet (carbon capture, utilization, and storage—CCUS) reservoirs, consideration must also be given to conductive reservoir components. The optimization of data acquisition before the survey involves analyzing target parameters and the sensitivity of multi-component CSEM. This optimization process typically includes on-site noise measurements and 3D anisotropic modeling. Based on our experience, subsequent surveys tend to proceed smoothly, yielding robust measurements that align with scientific objectives. Other critical aspects to be considered are using magnetotelluric (MT) measurements to define the overall background resistivities and integrating real-time quality assurance during data acquisition with 3D modeling. This integration allows the fine tuning of acquisition parameters such as acquisition time and necessary repeats. As a result, data can be examined in real-time to assess subsurface information content while the acquisition is ongoing. Consequently, high-quality data sets are usually obtained for subsequent processing and initial interpretation with minimal user intervention. The implementation of sensitivity analysis during the inversion process plays a pivotal role in ensuring that the acquired data accurately respond to the target reservoirs’ expected depth range. To elucidate these concepts, we present an illustrative example from a CO2 storage site in North Dakota, USA, wherein the long-offset transient electromagnetic method (LOTEM), a variation of the CSEM method, and the MT method were utilized. This example showcases how surface measurements attain appropriately upscaled log-scale sensitivity. Furthermore, the sensitivity of the CSEM and MT methods was examined in other case histories, where the target reservoirs exhibited conductive properties, such as those encountered in enhanced oil recovery (EOR), geothermal, and lithium exploration applications. The same equipment specifications were utilized for CSEM and MT surveys across all case studies.

Fast parallel IGA-ADS solver for time-dependent Maxwell's equations

We propose a simulator for time-dependent Maxwell's equations with linear computational cost. We employ B-spline basis functions as considered in the isogeometric analysis (IGA). We focus on non-stationary Maxwell's equations defined on a regular patch of elements. We employ the idea of alternating-directions splitting (ADS) and employ a second-order accurate time-integration scheme for the time-dependent Maxwell's equations in a weak form. After discretization, the resulting stiffness matrix exhibits a Kronecker product structure. Thus, it enables linear computational cost LU factorization. Additionally, we derive a formulation for absorbing boundary conditions (ABCs) suitable for direction splitting. We perform numerical simulations of the scattering problem (traveling pulse wave) to verify the ABC. We simulate the radiation of electromagnetic (EM) waves from the dipole antenna. We verify the order of the time integration scheme using a manufactured solution problem. We then simulate magnetotelluric measurements. Our simulator is implemented in a shared memory parallel machine, with the GALOIS library supporting the parallelization. We illustrate the parallel efficiency with strong and weak scalability tests corresponding to non-stationary Maxwell simulations.

Modelling of subsurface geological structure related to anomalous magnetotelluric phase data exceeding 90° in Central India

SUMMARY We note anomalous phase (>90°) in magnetotelluric (MT) measurements at several stations along the Narmada river course normal to the Tan Shear Zone (TSZ) in central India. We made efforts to derive a possible cause of the anomalous phase by comparing the results of numerical modelling based on the tectonics of the TSZ with the measured impedance tensor data. A zone of multiple reactivations of the TSZ leads to the formation of a damage zone and conjugate Riedel shear faults parallel and normal to the TSZ, respectively. The conjugate Riedel shear faults act as a pathway for the Narmada River course. The multiple tectonic activities associated with the TSZ induce surficial heterogeneity that results in distortion in the MT data. The Mohr circle and phase tensor analysis establish the effect of distortion in the data. To understand the cause of distortion as well as anomalous phase in the present data, we carried out synthetic 3-D modelling. We modelled the damage zone and conjugate Riedel shear faults by a cross shape. From the results, we noted that the near-surface 3-D heterogeneity across the TSZ critically distorted the YX component of the MT impedance tensor due to the strong current channelling. When the current flows across the conjugate Riedel shear zone (Narmada river course) charges are accumulated along the boundaries of near-surface heterogeneity. The charges associated with the current channelling brought out the reversion of the electric field direction that is reflected in the form of an anomalous phase in the YX component of the impedance tensor in measured data.

Simulation and Observations of Audio Magnetotelluric Measurements over Water-Covered Areas

Electromagnetic (EM) surveys play a significant role in mineral exploration. However, the EM method often faces limitations when investigating minerals in areas covered by rivers, lakes, or other water bodies. This paper introduces audio magnetotelluric (AMT) observation technology that utilizes separated electric and magnetic channels to deal with this challenge over water-covered areas. The study analyzes and discusses the characteristics of the relative error of the magnetic field through forward simulation. The observation and profile experiments were conducted at the estuary of a river in Liaoning Province, China, and high-quality data in the river and the pseudo-geoelectric section of the underwater space were successfully obtained. The results demonstrate the feasibility and effectiveness of the AMT observation technology over water-covered areas, emphasizing the importance of locating the magnetic channel in a quiet zone at a certain distance from the shore. This configuration helps reduce the influence of resistivity differences between water and shore, ultimately improving data quality and accuracy. The research suggests that the AMT observation technology, utilizing separated electric and magnetic channels, has the potential for further improvement and can serve as a valuable guide for mineral exploration over water-covered areas.

Magnetotelluric Power Line Noise Removal Using Temporally Varying Sinusoidal Subtraction of the Grid Utility Frequency

The magnetotelluric method relies on variations of natural electromagnetic fields, which in the vicinity of human settlements are persistently distorted by anthropogenic electromagnetic noise. A large source of noise to the magnetotelluric response is caused by the harmonic oscillations of the power network utility frequency centered on 50/60 Hz along with the associated higher harmonics. Removing this type of noise is essential for high frequency magnetotelluric measurements used for shallow surveys. There are a large number of approaches for how to treat power line noise in magnetotelluric signals, however, commonly used methods do not take into account time variations/instabilities of the utility frequency. That is not serious problem in vicinity of well balanced grid networks, but can cause issues in regions with larger utility frequency variations. Under such conditions, commonly used methods loose more of the natural signal, which is undesirable especially in case of very noisy datasets. Hence, we adopted approach for removing of power line noise with respect to time variations of the utility frequency and applied it to magnetotelluric signals to preserve more of natural signal. The method is based on modelling of the grid network harmonic oscillations by the optimum utility frequency and its integer multiples. The resulting sum of sinusoidal signals is subsequently subtracted from recorded data and only particular noise frequencies are removed from the original signal with high precision, while frequency ranges around power line harmonics are cleaned.

A decade of successful magnetotelluric surveys for delineating shallow geothermal reservoirs beneath nonvolcanic hot springs in Thailand

Magnetotelluric (MT) applications to probe the shallow reservoirs of nonvolcanic hot springs in Thailand have become ever more convincing during the past few decades. In 2013, the first modern MT survey was deployed over the Mae Chan hot spring. A shallow hot water reservoir within a network of fractured granite was imaged as a low-resistivity zone beneath the hot spring surrounded by a resistive granite heat source. With the success of the first MT deployment, in 2014 and 2015, MT measurements were expanded to five hot springs with the objective of 1 km deep drilling where the Muang Rae hot spring was chosen out of these five hot springs to assess the volume of hot fluid for a future geothermal power plant. The conductive zone was again detected at a depth of 200–600 m beneath the Muang Rae hot spring, which our team interpreted as a shallow hot fluid reservoir. However, others interpreted it as just the clay cap with no hot fluid inside, which occurs in most volcanic systems. Therefore, a 1 km well was drilled next to the conductor, instead of through it as we suggested. The well did not find any hot water flowing. In 2018, an MT survey with densely spaced stations was used around the Mae Chan hot spring to validate the shallow conductive zone beneath the hot spring through which five wells (<200 m) were drilled. The depth positions where the hot fluid was found match well with the conductor location at most wells. This is a good indicator that in the Thailand nonvolcanic system, the conductive zone beneath the hot spring is genuinely the shallow reservoir, and not the clay cap zone as in the volcanic systems.

The geothermal magmatic system at the northern Rahat volcanic field, Saudi Arabia, revealed from 3D magnetotelluric inversion

Harrat Rahat (HR) volcanic field is one of the most prominent volcanic areas in western Saudi Arabia, as it is a promising site for geothermal energy. HR has been intensively studied from a geologic point of view. Detailed inspections distinguished 12 stages of volcanic eruption with characteristic volcanics and age intervals. The area is seismologically active and suffered a swarm of earthquake events in 1999. Magnetotelluric measurements were conducted across the northern HR for deep electrical imaging to provide a comprehensive image of the geothermal source and its impact on volcanic eruptions and seismicity. Acquisition sites (65 in total) provided a square grid with ∼5 km intervals. The recording at each site continued for 20 h on average and covered a broad band of frequencies from 1 mHz to 300 Hz. Phase tensor representations and the calculated skew angles indicated three-dimensional subsurface structures. 3D inversion was then carried out using a start model of half space with 100Ωm. The inverted model shows a main magma chamber with very low resistivities in its interior (3 Ωm). Protrusions branch upwards, or first laterally and then upwards, from the main magma reservoir. They represent magma pathways to the surface and account for various eruption locations and stages. The magma chamber is overlain by a continuous zone of higher resistivity (>500 Ωm) that depicts the cooled-down lava sheets. The tabulated and spatially mapped historical eruptions are projected on the results of the 3D MT inversion in the form of slices cut at various depths to obtain clues about the details observed in the model. The epicenters of the 1999 swarm showed no direct spatial matching with features of the inverted model. Similarly, the foci positions are concentrated in the subsurface away from the geothermal source. This confirms the model, as the epicenters usually exist at zones of faults and fractures in the crust, even if the source of stresses that induce seismicity originates away from them. The rocks near the thermal source are ductile and cannot be a source of an elastic rupture. We conclude that the thermal magma reservoir is shallow (10–25 Km), of a large spatial extent (∼1000 km2), impacts surface features and induces earthquake activities.

Prospecting for ophiolite-type chromite deposit in Sartohay, West Junggar (NW China): Constraints from geological and geophysical data

Chromite is a critical mineral that is of great significance for the development of metallurgical, chemical, and other industries. In order to study the deep ore prospecting methods of ophiolite-type chromitite, we have carried out plane-profile geological and mineral surveys, as well as high-precision gravity, magnetic, and controlled source audio-frequency magnetotelluric (CSAMT) measurements in the Sartohay, West Junggar. Based on the analysis, interpretation, and inference of multi-source geological information, a three-dimensional (3D) geological model of the typical mining area in Sartohay District was established. The results show that Sartohay ophiolite was occurred in the form of rootless tectonic blocks and generally displayed deformation and metamorphism (e.g., oblique thrust, transverse extrusion, and lateral strike slip). Sartohay ophiolite is in tectonic contact with surrounding strata. Several northwest dipping high-angle thrust faults are developed at the contact sites and inside the ophiolite rock mass, showing a thrust imbricated structure. Meanwhile, ophiolite in the modeling area is a NW-trending monoclinic rock mass. The base surface of the ophiolite is obliquely distributed in an “N-like” shape, and the chromitite bearing lithofacies zones also have a NW-trending distribution. Geophysical data reveal that the studied mining area is characterized by a layered electrical structure in the vertical direction. The upper layer displays a high-resistivity agglomerate and a low-resistivity band distributed alternately, while the lower layer is a complete high-resistivity block. Combined with the deposit characteristics in Sartohay, it is considered that the morphological occurrence, geological structure, lithology, and lithofacies of ophiolite control the spatial distribution of related chromitites. The dunite-harzburgite lithofacies belts, highly strain zones, and strong tectonic alteration belts of the ophiolite rock mass are significant sites for the chromite ore bodies/groups occurrence. This study provides approaches for researching on the production laws of ophiolite-type chromitite and supplys a prospecting pattern for deep concealed orebodies.

Relationship of the Crustal Structure, Rheology, and Tectonic Dynamics Beneath the Lhasa‐Gangdese Terrane (Southern Tibet) Based on a 3‐D Electrical Model

AbstractThe tectonic dynamics of the Lhasa‐Gangdese terrane, southern Tibet, are still unclear and open questions persist regarding the structure, physical properties, and rheology of the lithosphere. A three‐dimensional electrical resistivity model was generated beneath the Lhasa terrane, 79°−94°E and 28°–33°N, an area of ∼1,500 by ∼600 km, using data from 537 magnetotelluric measurements. To overcome computational challenges, we present an approach in which multiple high‐resolution models from overlapping subregions are combined by means of an error‐weighted averaging scheme to construct a continuous electrical resistivity model. Based on the electrical structure, the rheology of the mid‐lower crust was investigated by constraining the melt fraction, in order to explore controls on dynamic mechanisms from a whole‐system perspective. The models shed light on the vertical and lateral migration of crustal materials, magma transportation, and continental deformation in the Lhasa terrane. The model shows resistive features in the south that may represent the Indian plate, and conductive features above this that may represent the migration of materials beneath the Lhasa terrane, both of which vary substantially along the Indus‐Yarlung Zangbo suture. In the mid‐lower crust, conductive features (with a likely corresponding decrease of the effective viscosity) may be related to underplating, magma reservoirs, mineralization sources, and vertical motion of buoyant materials. A weakened mid‐lower crust has consequences for surface deformation and may have contributed to the formation of (north‐south) rift zones. Furthermore, the inhomogeneous distribution of weakened areas in the mid‐lower crust has implications for the local lateral migration of materials.

Numerical modeling of the co-seismic electromagnetic signals observed during the 2014 Mw 4.9 Hammam Melouane earthquake, Northern Algeria

The seismic activity in northern Algeria provides the opportunity to observe earthquake-related electromagnetic signals through magnetotelluric (MT) measurements. Through an MT station, we measured the two telluric and three orthogonal magnetic field components. On December 23rd, 2014, a moderate earthquake (Mw = 4.9) occurred in the village of Hammam Melouane (Algeria), located 36 km northeast of the MT site. This event allowed the observation of a co-seismic electromagnetic signal in the MT time series. The analysis of the electromagnetic time series using Fourier transform revealed the spectral characteristics of these co-seismic electromagnetic signals with a high spectral amplitude that started at the frequency of 0.3 Hz. Then, we simulate the electric and magnetic responses of the Hammam Melouane earthquake based on the electrokinetic mechanism. We conducted the simulations using a point source model with a realistic layered earth model and compared the synthetic signals to the observed data in the frequency band 0.3–6 Hz. The result shows that the synthetic electric signals match the measured data in both amplitude and waveform for the parts of the recordings after the arrival of the S wave with a phase shift. Similarly, the synthetic magnetic field agrees with the observed data in the waveforms with a phase shift, but their amplitudes are smaller than the observed data. In the latter part of the records (>18 s), a poor disagreement can be seen between the synthetic and observed electromagnetic field signals in their shape and amplitudes. We investigated the rotation-induced magnetic field effect as a possible mechanism for the generation of additional co-seismic magnetic signals. The synthetic magnetic fields resulting from the magnetometer rotation modeling are in agreement with the observed magnetic data in terms of amplitudes after the arrival of seismic waves. The main conclusion of this study reveals that the electrokinetic effect is responsible for the generation of co-seismic electric fields while the rotation of magnetic sensors is likely the main cause of co-seismic magnetic fields.

Inferring the magmatic roots of volcano-geothermal systems in the Rotorua Caldera and Okataina Volcanic Centre from magnetotelluric models

To investigate the underlying magmatic systems that drive high-temperature geothermal fields and volcanism in the north-central Taupo Volcanic Zone (TVZ), we use a combination of land and lake-bottom magnetotelluric (MT) measurements. In total, MT data from 399 sites (including 34 lake-bottom locations) across 2800 km2 covering the Rotorua Caldera and Okataina Volcanic Centre are inverted to image the 3-D resistivity structure of the crust to 20 km depth. Beneath this region, a singular conductive zone in the mid-crust is shown, inferred to represent a complex of relatively mafic magma, consistent with petrologic and other geophysical data. Discrete low-resistivity (conductive) plumes rise from the margins of this magmatic complex and connect 1-to-1 with the locations of high-temperature geothermal fields at the surface. These plumes are inferred to represent pathways of heat transfer through the crust via magmatic intrusion and exsolved saline fluids, and inform connections with other geothermal areas located nearby to the major fields. An additional conductive plume rises beneath the area where the most voluminous recent volcanism has occurred in the Okataina Volcanic Centre; the summit region (e.g. Makatiti dome) in the Haroharo Volcanic Complex that has erupted 39 km3 (magma equivalent) in the past 9 ka. While hot springs locate around this conductive plume at the base of the rhyolite domes, the plume itself is arrested at 1–2 km depth beneath the resistive (inferred impermeable) young (5 ka old) rhyolite lava domes, and is considered to represent a deep ‘blind’ geothermal system in the TVZ.

Magnetotelluric Investigations of the Kīlauea Volcano, Hawaii

AbstractIn 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kīlauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kīlauea and the East and Southwest Rift zones, which has implications for understanding Kīlauea's plumbing system. An improved understanding of the rift zones has implications in understanding large‐scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi‐station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three‐dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kīlauea. Melt fractions between 0.096 and 0.117 are calculated for the Kīlauea and Puʻuʻōʻō low resistivity anomalies, respectively.