Fact-checked by Grok 2 weeks ago

Magnetotellurics

Magnetotellurics (MT) is a passive electromagnetic geophysical method that utilizes naturally occurring variations in the Earth's electric and magnetic fields to map the subsurface electrical resistivity structure, providing insights into geological formations from the shallow crust to the upper mantle.^[1] The technique measures orthogonal components of the electric field (E_x and E_y) and magnetic field (H_x, H_y, and sometimes H_z) at the surface, deriving the impedance tensor that relates these fields to subsurface conductivity variations across a broad frequency range from approximately 10^{-4} to 10^4 Hz.^[2] The natural electromagnetic signals employed in MT originate from two primary sources: global lightning discharges, which generate higher-frequency fields (>1 Hz) that propagate as plane waves in the Earth-ionosphere waveguide, and low-frequency variations (<1 Hz) induced by magnetospheric currents driven by solar wind interactions with Earth's magnetic field.^[3] These fields penetrate the Earth, inducing secondary currents whose amplitudes and phases are influenced by subsurface resistivity contrasts, allowing for the computation of apparent resistivity and phase angles that inform 1D, 2D, or 3D inversion models of the subsurface.^[1] Developed independently in the early 1950s by Louis Cagniard in France, and through theoretical work by Andrey Tikhonov and Tsuneji Rikitake, MT has evolved from analog recordings to modern digital systems, with significant advancements including the introduction of geomagnetic depth sounding (GDS) in 1959 and biennial international EM induction workshops since 1972.^[4] Key technological progress in the 1990s included digital data storage and 2D inversion software, followed by 3D inversion capabilities in the 2010s, enhancing resolution for complex structures.^[4] MT is widely applied in resource exploration, including mineral deposits, hydrocarbons, and geothermal reservoirs, as well as in academic studies of lithospheric architecture, tectonic processes, and hazards like geomagnetically induced currents.^[1] Notable large-scale implementations include Australia's AusLAMP program, launched in 2013, which has deployed over 1,500 long-period MT sites as of 2023 to create continent-scale 3D resistivity models revealing mantle metasomatism and crustal anisotropy.^[5]^[4] Despite its advantages in deep penetration without artificial sources, MT faces challenges from cultural noise interference and source variability, often requiring remote reference techniques for data quality improvement.^[3]

Overview

Definition and Principles

Magnetotellurics (MT) is a passive geophysical exploration technique that measures the natural variations in Earth's electromagnetic fields to infer the subsurface electrical resistivity structure. It relies on time-varying magnetic fields from external sources, such as global thunderstorms and solar-induced ionospheric currents, which induce secondary electric fields and currents within the conductive Earth. These measurements allow imaging of resistivity contrasts associated with geological features like faults, fluids, or mineral deposits, providing insights into subsurface properties without the need for artificial energy sources.^[6]^[7] The core physical principle of MT is electromagnetic induction governed by Faraday's law, which states that a time-varying magnetic field induces an electric field, leading to diffusive propagation of electromagnetic energy into the Earth. Under the plane wave assumption, the incident fields are treated as vertically propagating plane waves with horizontal polarization, where the primary magnetic field penetrates the surface and induces horizontal electric fields that interact with subsurface conductivity variations. This diffusion process attenuates the fields with depth, enabling the mapping of resistivity as the primary parameter, in contrast to seismic methods that primarily image acoustic velocity variations.^[6]^[8]^[7] A key concept in MT is the skin depth, which quantifies the penetration distance of electromagnetic fields into the Earth and depends on the frequency of the variations and the material's conductivity. The skin depth \delta is given by

\delta = \sqrt{\frac{2}{\omega \mu \sigma}},

where \omega is the angular frequency, \mu is the magnetic permeability, and \sigma is the electrical conductivity (or inversely, resistivity \rho = 1/\sigma). Lower frequencies (longer periods) yield greater skin depths, allowing deeper exploration, typically from tens of meters to hundreds of kilometers. This frequency-dependent penetration is fundamental to resolving vertical resistivity structure.^[8]^[7]^[6]

Historical Development

The theoretical foundations of magnetotellurics trace back to the late 19th century, when British physicist Arthur Schuster recognized that geomagnetic variations could induce electric currents in the Earth's crust, linking diurnal magnetic fluctuations to telluric potentials observed at the surface.^[9] This insight laid the groundwork for understanding electromagnetic induction in the subsurface, though practical applications remained undeveloped for decades. The method's practical development occurred in the mid-20th century, with Soviet mathematician Andrey Tikhonov and Japanese seismologist Tsuneji Rikitake proposing in 1950 techniques to infer deep crustal electrical properties from natural electromagnetic fields, independently paralleled by French geophysicist Louis Cagniard's 1953 formulation of the magneto-telluric approach, which enabled the first soundings by measuring orthogonal electric and magnetic field components.^[10]^[11]^[12] These contributions shifted focus from theoretical speculation to geophysical prospecting, with early field tests in the Soviet Union and France demonstrating the potential for subsurface resistivity mapping without artificial sources.^[10] By the 1960s, magnetotellurics expanded internationally, particularly in the United States through adoption by the U.S. Geological Survey (USGS) for geothermal and mineral exploration, alongside the establishment of processing standards and the acquisition of initial continental-scale profiles across North America and Europe.^[13] This era marked the transition to routine surveys, fostering global collaboration and refining data interpretation for crustal studies.^[10] Commercialization accelerated in the 1970s and 1980s, driven by companies like Phoenix Geophysics, founded in 1975, which pioneered robust instrumentation for resource exploration, coinciding with the shift from analog recording devices to digital systems that improved signal fidelity and reduced noise in field operations.^[14] From the 1990s onward, magnetotellurics integrated with seismic and gravity methods for multidimensional imaging, bolstered by open-source software such as ModEM in the 2010s for three-dimensional inversions, and large-scale datasets from initiatives like EarthScope's USArray (2006–2018), which deployed nearly 1,800 stations across the contiguous United States.^[15]

Theoretical Foundations

Electromagnetic Induction Basics

Electromagnetic induction forms the foundation of magnetotellurics (MT), where time-varying natural magnetic fields penetrate the Earth and induce secondary electric currents that reveal subsurface conductivity structures. In the MT frequency range (typically 10^{-4} to 10^3 Hz), the interaction is governed by low-frequency electromagnetic phenomena, allowing the use of the quasi-static approximation, which neglects displacement currents relative to conduction currents because \sigma \gg \omega \epsilon, where \sigma is conductivity, \omega is angular frequency, and \epsilon is permittivity.^[16] This approximation simplifies Maxwell's equations to Faraday's law, \nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}, and Ampère's law, \nabla \times \mathbf{H} = \mathbf{J}, with the constitutive relations \mathbf{J} = \sigma \mathbf{E} and \mathbf{B} = \mu \mathbf{H}, assuming constant magnetic permeability \mu = \mu_0.^[17] Taking the curl of Faraday's law and substituting Ampère's law yields the diffusion equation for the magnetic field, \frac{\partial \mathbf{B}}{\partial t} = \frac{1}{\sigma \mu} \nabla^2 \mathbf{B}, which describes the diffusive propagation of electromagnetic fields in conducting media rather than wave-like behavior.^[17] A similar equation holds for the electric field. This diffusive nature implies that fields attenuate and phase-shift with depth, contrasting with free-space propagation. For a plane-wave solution in a uniform half-space, assume a horizontal magnetic field H_y e^{i\omega t} incident normally from above; it induces a horizontal electric field E_x = Z H_y, where Z is the complex impedance relating the fields at the surface.^[16] The fields penetrate to a characteristic skin depth \delta = \sqrt{\frac{2}{\omega \mu \sigma}}, beyond which they decay exponentially as e^{-z/\delta} while phase-shifting by z/\delta radians.^[16] Lower frequencies correspond to larger \delta, enabling deeper probing (e.g., periods of seconds probe kilometers, while higher frequencies resolve shallow structures). For example, in a medium with \sigma = 0.01 S/m and \mu = 4\pi \times 10^{-7} H/m, \delta \approx 5000 m at 1 Hz. This frequency-dependent attenuation establishes the scale for MT sounding.^[17] In real Earth settings, electrical anisotropy—where conductivity varies with direction due to aligned minerals or structures—alters field polarization and can produce phase differences exceeding 90° between orthogonal components, complicating 1D interpretations.^[18] Similarly, lateral heterogeneity introduces 3D effects, such as current channeling or galvanic distortion, which scatter fields and deviate from plane-wave assumptions, requiring advanced modeling to account for non-uniform conductivity distributions.^[19]

MT Response Functions and Modeling

In magnetotellurics, the primary observable is the impedance tensor \mathbf{Z}, a 2×2 complex matrix that linearly relates the horizontal electric field \mathbf{E} to the horizontal magnetic field \mathbf{H} at the Earth's surface according to \mathbf{E} = \mathbf{Z} \mathbf{H}, where the components are E_x = Z_{xx} H_x + Z_{xy} H_y and E_y = Z_{yx} H_x + Z_{yy} H_y.^[20] This tensor encapsulates the subsurface electrical conductivity structure and is estimated from time-series measurements of \mathbf{E} and \mathbf{H} via Fourier transforms at discrete frequencies \omega.^[7] In a one-dimensional (1D) Earth, \mathbf{Z} simplifies such that off-diagonal elements dominate, with Z_{xy} = -Z_{yx} and diagonal elements near zero, reflecting plane-wave propagation.^[21] From the impedance tensor, key derived quantities include the apparent resistivity \rho_a(\omega) and phase \phi(\omega), typically computed from the off-diagonal element Z_{xy} for transverse electric (TE) mode in 1D or 2D settings. The apparent resistivity is given by

\rho_a(\omega) = \frac{1}{\mu_0 \omega} |Z_{xy}|^2,

where \mu_0 is the magnetic permeability of free space, providing a frequency-dependent measure of effective subsurface resistivity that increases with depth as \sqrt{\rho / \omega} due to electromagnetic diffusion.^[6] The phase is \phi = \arg(Z_{xy}), which for a homogeneous conductor is 45° and shifts to indicate transitions between resistive and conductive layers, with values between 0° and 90° signifying increasing conductivity with depth.^[7] These parameters form the basis for sounding curves used in interpretation.^[21] Another important response function is the tipper vector \mathbf{T}, also known as the induction vector, which relates the vertical magnetic field component H_z to the horizontal fields via H_z = T_x H_x + T_y H_y, where \mathbf{T} = (T_x, T_y) is complex-valued.^[6] In terms of the impedance tensor, \mathbf{T} can be expressed as T_x + i T_y = \frac{Z_{xy} - Z_{xx}}{Z_{yy} - Z_{yx}} for certain approximations, though it is directly computed from field data.^[22] The real part of \mathbf{T}, often plotted as induction arrows, points toward regions of higher conductivity, aiding in mapping lateral heterogeneities, while its magnitude diminishes in 1D structures where H_z = 0. Dimensionality of the subsurface is assessed using indicators derived from \mathbf{Z}, such as the skew angle \beta = \frac{|Z_{xy} - Z_{yx}|}{|Z_{xy}| + |Z_{yx}|}, which quantifies deviation from 2D symmetry.^[23] Values of \beta < 0.1 suggest 1D or 2D conditions suitable for simplified modeling, while higher \beta indicates 3D complexity, influencing the choice of inversion approach.^[22] The phase tensor, formed from the imaginary parts of \mathbf{Z} elements, provides a distortion-free alternative for dimensionality analysis, with its skew related to \beta.^[24] Forward modeling in magnetotellurics computes \mathbf{Z} and \mathbf{T} for a given conductivity model \sigma(\mathbf{r}) by solving Maxwell's equations under plane-wave excitation. For 1D layered Earth models, analytical solutions employ Hankel transforms to integrate radial wave contributions, yielding explicit expressions for \rho_a and \phi as functions of layer resistivities and thicknesses.^[25] These are computationally efficient and form the basis for initial interpretations, with fast Hankel transform algorithms accelerating evaluation for arbitrary layer sequences.^[25] In two-dimensional (2D) cases, assuming strike-parallel invariance, numerical methods like finite differences discretize the domain on a staggered grid to solve for TE and transverse magnetic (TM) modes separately, capturing galvanic and inductive effects.^[26] For three-dimensional (3D) structures, finite-difference or finite-element methods on unstructured grids are required, often incorporating absorbing boundaries to simulate infinite domains, though they demand significant computational resources.^[27] Sensitivity kernels describe how perturbations in subsurface resistivity \delta \rho affect the observed \mathbf{Z}, formalized as Fréchet derivatives \frac{\partial Z_{ij}}{\partial \rho(\mathbf{r})}, which peak at depths scaling with the skin depth \delta = \sqrt{2 / (\mu_0 \omega \sigma)} and broaden with frequency.^[28] These kernels reveal trade-offs in resolution, with horizontal smearing due to diffusive propagation, emphasizing the need for multi-frequency data to constrain structure. Inversion of MT data to recover \sigma(\mathbf{r}) from \mathbf{Z} suffers from non-uniqueness, as multiple conductivity models can produce fitting responses within data errors, exacerbated by the ill-posed nature of the inverse problem and limited depth sensitivity at high frequencies.^[29] Regularization, such as minimum-structure norms, mitigates this by favoring smooth models consistent with prior geologic information.^[30]

Field Methods

Survey Design and Site Selection

Site selection for magnetotelluric (MT) surveys prioritizes locations that minimize cultural electromagnetic noise, such as those distant from power lines, urban areas, railroads, pipelines, and other anthropogenic sources that can contaminate natural signals.^[31] Ideal sites are chosen for their geological relevance to the target structures, ensuring the array covers key features like faults or basins, while also considering accessibility for field crews and equipment transport in remote or rugged terrains.^[31]^[32] Flat or gently sloping terrain facilitates electrode and sensor deployment, whereas steep or vegetated areas may require additional logistical planning to avoid unstable ground or excessive vegetation interference.^[33] Survey geometries are designed to match the expected dimensionality of the subsurface targets, with linear profile lines typically used for two-dimensional (2D) imaging along strike directions and grid arrays for three-dimensional (3D) structures in complex geology.^[31] Station spacing is determined by the desired resolution and target depth, where broader arrays spanning 5-50 km provide adequate coverage for crustal features at 1-10 km depths, while tighter spacings of 100 m or less delineate shallow targets.^[31]^[32] For regional surveys, 20-100 stations are common to balance spatial resolution with fieldwork efficiency.^[32] Frequency range considerations guide survey planning, with broadband MT setups targeting 10^{-4} to 10^4 Hz to image from shallow sediments to deep crustal layers, though lower frequencies below 0.001 Hz often necessitate extended overnight recordings.^[32]^[31] Trade-offs involve coverage versus cost, as wider arrays for low-frequency penetration increase logistical demands, while audio-magnetotelluric (AMT) subsets above 1 Hz suit shallower investigations but face a "deadband" gap around 1–5 kHz where signal strength diminishes.^[31]^[34] Remote reference strategies enhance data quality by deploying additional stations at distant sites, typically 20-1000 km away depending on noise coherence, to average out local biases and improve signal-to-noise ratios through cross-correlation of magnetic fields.^[31]^[32] At minimum, one remote site is used, but multiple references are preferred for robust noise rejection in culturally contaminated regions.^[31] Environmental permits and safety protocols are essential for MT fieldwork in remote areas, requiring approvals for land access, vegetation disturbance, and cultural heritage protection, often coordinated with local authorities or indigenous groups.^[35] Safety measures include monitoring weather to halt operations during thunderstorms, securing cables against wildlife, and implementing health plans for isolation risks like medical evacuations in rugged terrains.^[31]^[35]

Instrumentation and Data Acquisition

Magnetotelluric instrumentation relies on sensitive sensors to measure the natural variations in Earth's electric and magnetic fields. For the magnetic field components (Hx, Hy, and optionally Hz), fluxgate magnetometers or induction coil magnetometers are commonly used. Fluxgate magnetometers provide stable measurements at low frequencies with a noise level of approximately 4 × 10^{-12} T/√Hz at 1 Hz, making them suitable for long-period (>100 s) observations.^[36] Induction coil magnetometers, preferred for broadband applications, achieve noise levels below 10^{-8} nT²/Hz at 1 s periods, equivalent to sensitivities around 10^{-12} T/√Hz in the audio-frequency range.^[37] These sensors are typically oriented north-south and east-west using GPS for precise alignment.^[37] Electric field components (Ex, Ey) are recorded using non-polarizing electrodes, such as Ag/AgCl or Pb/PbCl types, which offer low noise performance critical for detecting weak telluric currents. Ag/AgCl electrodes exhibit noise levels below 1 μV/√Hz, with power spectral densities as low as 10^{-14} V²/Hz at 1 s periods, ensuring reliable measurements over electrode separations of 50–100 m.^[37]^[38] Electrodes are buried in moist soil or backfilled with bentonite to minimize contact resistance and cultural noise interference.^[37] Recording systems consist of multi-channel digitizers with 24- to 32-bit analog-to-digital converters (ADCs) and dynamic ranges exceeding 130 dB, enabling capture of signals from DC to several kHz.^[37] Sampling rates typically reach 100 Hz or higher for broadband surveys, with GPS timing modules ensuring synchronization across sites to within milliseconds. Autonomous recorders, powered by lithium batteries (often 120 Ah or more), support deployments of 4–6 weeks for long-period data, incorporating solid-state storage for terabytes of time-series data.^[37] These systems facilitate simultaneous recording of Ex, Ey, Hx, and Hy (with Hz optional for 3D studies), producing continuous time-series essential for frequency-domain analysis.^[39] Acquisition protocols emphasize continuous, simultaneous logging to preserve phase relationships between fields. For broadband magnetotellurics, time-series are segmented into overlapping windows, typically 15–30 minutes per frequency decade, to provide sufficient signal cycles for robust transfer function estimation across 10^{-4} to 10^3 Hz.^[40] Remote reference stations are often co-recorded to enhance signal-to-noise ratio by mitigating local noise.^[39] Calibration procedures include in-situ tests to verify sensor orientation, gain, and response. Magnetic sensors undergo frequency-dependent calibration using secondary coils to generate known fields, while electrodes are checked for DC offset and impedance stability. Special attention is given to handling dead bands in audio frequencies (e.g., 50–60 Hz cultural noise), where alternative processing or site relocation may be required.^[37] Since the 2010s, modern trends have shifted toward miniaturized, low-power nodal systems for deploying dense arrays (hundreds of stations over tens of km²). These devices integrate GPS synchronization, solid-state memory, and ultra-low-noise amplifiers in compact, battery-efficient packages, enabling high-resolution 3D imaging with reduced logistical demands.^[41]

Data Processing and Interpretation

Preprocessing and Noise Handling

Preprocessing of magnetotelluric (MT) data begins with time-series editing to ensure the integrity of the recorded electric and magnetic field measurements. This step involves identifying and removing anomalous spikes caused by lightning strikes or instrument malfunctions, which can distort the natural electromagnetic signals. Common techniques include threshold-based detection, where data points exceeding several standard deviations from the local mean are flagged and interpolated or excised, often using a running median filter to preserve signal continuity. For uneven sampling due to gaps in recording, spectral methods such as decimation with anti-aliasing filters are applied to produce evenly spaced segments suitable for Fourier transform analysis.^[42]^[43] MT time series are susceptible to various noise sources that must be characterized and mitigated during preprocessing. Local cultural noise, such as harmonics from power lines at 50/60 Hz and their multiples, introduces coherent interference that contaminates higher-frequency bands. In contrast, regional noise from ionospheric and magnetospheric variations affects broader scales but is generally more stationary. A notable challenge is the "dead band" around 1-10 Hz, where natural source energy is minimal, leading to poor signal-to-noise ratios and requiring specialized enhancement techniques.^[43]^[44] To achieve robust estimates of MT response functions, preprocessing employs advanced statistical methods that account for noise contamination. Remote reference processing, which uses simultaneous measurements from a distant "quiet" site, reduces local noise bias by computing cross-power spectra between local electric fields and remote magnetic fields, weighted by coherence to emphasize reliable segments. This approach, pioneered in the late 1970s, significantly improves impedance tensor accuracy in noisy environments. Outlier rejection is further enhanced through inter-station ratios, where discrepancies in phase or amplitude between nearby sites flag contaminated data for exclusion, ensuring consistency across the survey array.^[45]^[46] Quality control in MT preprocessing relies on quantitative metrics to validate data reliability before inversion. Coherence between orthogonal electric (E) and magnetic (H) field components serves as a primary indicator, with thresholds typically exceeding 0.8 considered indicative of high-quality data suitable for interpretation. Error bars on response functions are estimated using jackknife resampling, which systematically omits subsets of the time series to assess variability, providing a measure of statistical robustness without assuming Gaussian noise. These metrics guide the selection of usable frequency bands and inform subsequent analysis steps.^[47]^[48] Several software tools facilitate these preprocessing tasks, offering standardized workflows for filtering and editing. The BIRRP (Bounded Influence Remote Reference Processing) package implements robust remote reference techniques with built-in spike detection and coherence-based weighting. Similarly, implementations of the EMTF (Electromagnetic Transfer Functions) framework, such as the open-source Aurora tool, support initial time-series cleaning and spectral analysis for broad compatibility with MT data formats. These tools emphasize non-endorsed, community-driven solutions for reproducible noise handling.^[49]^[50]

Dimensionality Assessment and Inversion

Dimensionality assessment in magnetotellurics involves analyzing the impedance tensor derived from preprocessed data to determine whether the subsurface conductivity structure can be approximated as one-dimensional (1D), two-dimensional (2D), or requires a three-dimensional (3D) model.^[22] The strike angle, which indicates the principal direction of geoelectric strike, is estimated by rotating the impedance tensor \mathbf{Z} to minimize the difference between its off-diagonal elements, often using the angle \alpha that diagonalizes the tensor such that |Z_{xy}| \approx |Z_{yx}| and |Z_{xx}|, |Z_{yy}| are minimized.^[51] This rotation aligns the measurement coordinates with the regional geologic strike, facilitating dimensionality classification.^[52] The Swift skew, defined as \text{skew} = \frac{|\mathbf{Z}_{xx} + \mathbf{Z}_{yy}|}{|\mathbf{Z}_{xy} - \mathbf{Z}_{yx}|}, quantifies deviations from 1D or 2D symmetry; values below 0.1–0.2 typically indicate 1D/2D structures, while higher values suggest 3D effects due to galvanic distortions or complex geometry.^[53] Introduced in early theoretical work, this invariant parameter helps identify sites suitable for simpler inversions.^[54] For rotation-invariant analysis immune to galvanic distortions, the phase tensor \boldsymbol{\Phi} = [\mathrm{Re}(\mathbf{Z})]^{-1} \mathrm{Im}(\mathbf{Z}) (where \mathrm{Re}(\mathbf{Z}) and \mathrm{Im}(\mathbf{Z}) are the real and imaginary parts of \mathbf{Z}) is employed, with its ellipticity \beta = \frac{|\phi_2 - \phi_1|}{\phi_2 + \phi_1} (phases \phi_1, \phi_2) indicating 3D complexity when exceeding 0.2; in 1D cases, \boldsymbol{\Phi} is circular with equal phases.^[22]^[55] Inversion frameworks reconstruct resistivity models from these assessed data. For 1D layered models, Occam's inversion minimizes model roughness subject to data fit, producing smooth conductivity-depth profiles via a linear least-squares solution with a smoothness constraint \int (\frac{d \log \rho}{dz})^2 dz.^[56] Alternatively, the Marquardt-Levenberg algorithm iteratively updates a layered starting model using damped non-linear least-squares to fit apparent resistivity and phase curves.^[57] For 2D and 3D cases, non-linear least-squares methods like WSINV3D solve the forward problem on finite-difference grids, minimizing the data misfit \phi_d = \sum | \frac{d_{obs} - d_{pred}}{\sigma_d} |^2 through conjugate-gradient optimization. Regularization stabilizes these ill-posed inversions by incorporating prior information. Smoothness constraints, such as minimum structure norms, penalize rapid lateral or vertical resistivity changes to favor geologically plausible models, often using \phi_m = \int |\nabla \log \rho|^2 dV.^[58] Minimum gradient support regularization promotes compact anomalies by minimizing \int |\nabla \log \rho| dV, enhancing resolution of sharp features like faults.^[59] Joint inversions with gravity or seismic data incorporate priors by coupling model parameters, reducing non-uniqueness through shared structures.^[60] Uncertainty quantification relies on the posterior covariance matrix approximated from the Hessian \mathbf{H} = \mathbf{J}^T \mathbf{W}_d \mathbf{J} + \lambda \mathbf{W}_m, where \mathbf{J} is the Jacobian and \mathbf{W}_d, \mathbf{W}_m are data and model weighting matrices; standard deviations scale with \sigma_m \approx \sqrt{\text{diag}(\mathbf{H}^{-1})}.^[61] Resolution matrices, derived as \mathbf{R} = (\mathbf{J}^T \mathbf{W}_d \mathbf{J} + \lambda \mathbf{W}_m)^{-1} \mathbf{J}^T \mathbf{W}_d \mathbf{J}, reveal parameter trade-offs, with diagonal elements near 1 indicating well-resolved depths.^[62] In practice, inversions iteratively refine models until the normalized data misfit falls below an error floor, such as 5% of observed values, ensuring geophysical consistency without overfitting noise.

Variants

Natural Source Variants

Broadband magnetotellurics (MT) employs natural electromagnetic fields spanning frequencies from approximately $10^{-3} to $10^{3} Hz to image subsurface electrical resistivity structures at depths ranging from 1 to 100 km. These signals primarily originate from global lightning discharges, which dominate the higher-frequency band (roughly 1–1000 Hz), and from magnetospheric currents induced by solar wind interactions, which provide the lower-frequency components (below 1 Hz). This broad spectral coverage allows for comprehensive crustal imaging without artificial signal generation, distinguishing it from controlled-source techniques.^[63]^[1]^[34] Long-period MT extends this approach to ultra-low frequencies of $10^{-4}–0.1 Hz, targeting deeper mantle structures at depths of 100–1000 km, where electrical conductivity anomalies reveal insights into asthenospheric properties and lithospheric boundaries. The signals at these periods are inherently weak due to distant ionospheric and magnetospheric sources, necessitating extended recording durations—often spanning months—to accumulate sufficient data for reliable impedance estimates and to overcome low signal-to-noise ratios. Such prolonged observations are essential for resolving subtle conductivity variations in the upper mantle, though they demand robust, low-noise instrumentation.^[64]^[65]^[66] Geomagnetic depth sounding (GDS) represents a specialized subset of natural-source MT that analyzes only the horizontal magnetic field (H-field) components to infer large-scale, global electrical conductivity distributions, often focusing on regional to continental anomalies. Developed historically in auroral zones, where geomagnetic variations from substorm activity and polar currents are amplified, GDS originated from early 20th-century observatory measurements that exploited these enhanced signals for initial conductivity profiling. This approach simplifies data requirements compared to full MT but limits resolution to broader structures, complementing traditional MT for synoptic studies.^[67]^[68]^[69] These natural-source variants provide key advantages, including non-invasive deployment across expansive areas for large-scale subsurface imaging, leveraging ubiquitous global electromagnetic fluctuations without logistical burdens of energy sources. However, limitations arise in shallow-depth resolution (less than 1 km) due to gaps in natural signal availability, such as diurnal or seasonal variations in lightning activity and cultural noise interference in populated regions, which can obscure high-frequency data essential for near-surface details.^[2]^[6]^[70] Recent adaptations incorporate space weather observations from post-2015 satellite missions, such as ESA's Swarm constellation, to refine source field modeling and mitigate external field distortions in MT data processing. These datasets enable more accurate separation of local induction responses from global geomagnetic perturbations, enhancing inversion reliability for deep structures amid varying solar activity.^[71]^[72]

Controlled Source Methods

Controlled source methods in magnetotellurics employ artificial electromagnetic transmitters to generate controlled signals, supplementing or replacing natural sources for enhanced data quality in the audio-frequency range. These techniques, inspired by traditional magnetotelluric principles, allow for targeted investigations by producing predictable fields that penetrate the subsurface, particularly useful where natural signals are weak or contaminated.^[73] One prominent variant is Controlled Source Audiomagnetotellurics (CSAMT), which utilizes low-frequency transmitters operating from 1 to 10 kHz to image depths of 0.1 to 2 km. The primary source in CSAMT is a grounded horizontal electric dipole, typically 1 to 2 km in length, positioned several skin depths away from the receiver array to minimize near-field distortions. This setup generates a quasi-static electromagnetic field that induces telluric currents, enabling resistivity mapping of shallow geological structures such as faults and bedrock interfaces.^[73]^[74] Semi-airborne magnetotellurics represents another controlled source approach, featuring ground-based high-power transmitters paired with airborne receivers towed by aircraft or unmanned aerial systems (UAS). In this configuration, transmitters—often horizontal electric dipoles 800 m to 1.5 km long—emit signals that excite subsurface currents, while receivers mounted on multicopters or helicopters at low altitudes (e.g., 40 m) measure the induced magnetic fields. This method enhances resolution in rugged or inaccessible terrains by combining the deep penetration of ground sources with the broad coverage and efficiency of aerial surveys.^[75] Source effects in these methods arise from the finite geometry of the transmitter, leading to distinctions between near-field, transition, and far-field zones. In the near-field (within approximately one skin depth), magnetic fields decay as 1/r² and electric fields as 1/r³, resulting in apparent resistivities proportional to r² and near-zero phases, which deviate significantly from plane-wave assumptions. The far-field, beyond about three skin depths, approximates plane-wave behavior suitable for standard magnetotelluric inversion. Corrections for these effects, including 3D source geometry and topographic influences, involve forward modeling with finite element methods to transform data into equivalent plane-wave responses, though challenges persist due to unknown subsurface conductivity.^[76] Compared to natural source magnetotellurics, controlled source methods offer higher signal-to-noise ratios in noisy environments, as the known artificial signal can be isolated from cultural or natural interference, enabling reliable data acquisition even at low frequencies. However, they suffer from limited range—constrained by transmitter power and distance (typically up to several kilometers)—and higher operational costs due to the need for deploying and powering sources, along with added logistical complexities.^[73]^[77]

Applications

Commercial Exploration

Magnetotellurics (MT) plays a significant role in hydrocarbon exploration by imaging subsurface structures such as salt domes and underlying reservoirs, particularly in complex geological settings like the Gulf of Mexico. In marine MT surveys at the Gemini prospect, sea-floor data revealed a resistive salt body with surrounding sediments exhibiting resistivities around 1-100 ohm-m, enabling accurate mapping of the salt base to within 5-10% of its burial depth and aiding identification of potential hydrocarbon traps beneath the salt.^[78] These surveys are often integrated with seismic methods to enhance resolution, as joint inversions improve delineation of subsalt reservoirs where seismic imaging alone is limited by salt-induced distortions.^[79] In mineral exploration, MT has been applied to detect conductive targets associated with ore deposits, such as porphyry copper systems and kimberlites, particularly in Australia during the 2000s. For instance, MT imaging in western Victoria identified lower crustal conductors at 20-30 km depths linked to orogenic gold mineralization, guiding targeted drilling in gold mine districts.^[80] Similarly, broadband MT surveys have delineated conductive kimberlite pipes, which exhibit low-resistivity signatures due to serpentinization, supporting exploration in prospective terranes. For geothermal resource mapping, MT is used to characterize reservoirs and overlying clay caps, as demonstrated in Iceland's high-enthalpy fields. In the Krafla geothermal area, MT inversions revealed low-resistivity clay caps (typically 10-50 ohm-m) at depths corresponding to temperatures of 100-300°C, overlying hotter reservoirs up to 500°C, which informed the Iceland Deep Drilling Project (IDDP) well targeting supercritical conditions.^[81] These models helped identify drilling opportunities by distinguishing permeable reservoirs from impermeable caps, contributing to successful resource delineation. Recent applications include a 2021-2023 MT survey in Watson Lake, Canada, to assess geothermal potential in the region.^[82] The cost-effectiveness of MT surveys provides economic benefits by reducing the risk of dry wells in exploration campaigns. By improving subsurface imaging, MT lowers overall exploration costs compared to drilling without prior geophysical constraint. In environmental applications, MT supports groundwater assessment and CO2 storage site evaluation. MT soundings in sedimentary basins have mapped aquifer extents and salinities, aiding sustainable water resource management in arid regions.^[83] For CO2 sequestration, post-2015 monitoring at saline aquifer sites like Aquistore, Canada, utilized MT to establish baseline resistivity models and detect plume migration through changes in electrical conductivity.^[84] These efforts ensure site integrity by verifying containment in deep formations.^[85]

Scientific Research

Magnetotellurics has significantly advanced the imaging of crustal and lithospheric structures, particularly in revealing conductive layers associated with the electrical asthenosphere beneath continental regions. Studies have identified conductors with resistivities around 100 ohm-m at depths of 100-200 km, often interpreted as zones of partial melt or hydrated minerals that facilitate electrical conduction. For instance, in the western United States, magnetotelluric data from the EarthScope program delineated a low-resistivity layer at approximately 100-150 km depth, attributed to interconnected fluids or melt in the asthenosphere, providing insights into lithospheric deformation and mantle upwelling.^[86] These findings highlight MT's role in mapping the lithosphere-asthenosphere boundary (LAB), where resistivity contrasts mark transitions from rigid, resistive lithosphere (>1000 ohm-m) to more ductile, conductive asthenosphere, influencing plate tectonics and continental stability.^[87] In mantle studies, MT has elucidated differences in conductivity profiles between oceanic and continental settings, with the 1990s Mantle Electromagnetic and Tomography (MELT) experiment serving as a seminal effort along the East Pacific Rise. The experiment revealed anisotropic melt zones in the oceanic upper mantle, with conductivities elevated to ~0.1-1 S/m (1-10 ohm-m) due to aligned partial melts beneath the ridge axis, contrasting with the more isotropic, higher-resistivity continental mantle.^[88] These oceanic profiles demonstrate focused melt accumulation at depths of 50-100 km, driven by mantle upwelling, while continental counterparts often show broader, fluid-influenced conductors, underscoring MT's utility in comparing mantle dynamics across tectonic regimes.^[89] MT investigations of tectonics and subduction zones have linked arc conductors to fluid migration from dehydrating slabs. In Japan, extensive MT arrays across the northeastern subduction zone imaged conductive features at 10-50 km depths, corresponding to dehydration zones where slab-released fluids enhance crustal conductivity to ~10-50 ohm-m.^[90] These conductors, often aligned parallel to the arc, trace fluid pathways that weaken the forearc and influence seismicity, as seen in three-dimensional models revealing upward-migrating fluids from the Pacific slab at 80-100 km depth.^[91] Such imaging supports models of subduction-driven volatile cycling, where MT-detected anomalies correlate with seismic reflectors and geochemical signatures of fluid alteration. For volcano and earthquake monitoring, MT has detected precursor resistivity changes linked to stress-induced fluid mobilization. Around the 2011 Tohoku earthquake, post-event MT surveys revealed resistivity decreases in the crust near the rupture zone, from >1000 ohm-m to ~100 ohm-m, interpreted as fluid influx along faults following the megathrust event.^[92] Temporal analyses of crustal structures before and after the quake indicated subtle pre-event anomalies in apparent resistivity, potentially signaling pore pressure buildup in conductive layers at 10-20 km depth.^[93] These observations extend to volcanic settings, where MT monitors shallow conductors (<10 km) for unrest indicators, such as resistivity drops preceding eruptions due to hydrothermal activity. Global initiatives have leveraged large-scale MT deployments to probe plate boundaries and lithospheric evolution. The USArray Transportable Array, operational since the 2000s, provided continental-scale MT data across the United States, imaging lithospheric conductors associated with ancient sutures and modern tectonics, such as a ~200 km thick resistive lithosphere thinning beneath the western cordillera.^[94] Similarly, the SinoProbe project in the 2010s and 2020s, focusing on South China, deployed MT arrays to delineate plate boundaries between the Yangtze and Cathaysia blocks, revealing NE-SW oriented conductive interfaces at 50-100 km depth indicative of multi-microcontinental assembly and subduction remnants.^[95] These collaborative efforts integrate MT with seismic data to model deep Earth dynamics, fostering interdisciplinary understanding of continental growth and intraplate deformation.^[96]

Challenges and Future Directions

Method Limitations

Magnetotellurics (MT) suffers from inherent resolution trade-offs, particularly in lateral and vertical directions, due to the diffusive nature of electromagnetic wave propagation in conductive media. Lateral resolution decreases with depth and is generally on the order of the target depth, making it challenging to delineate fine-scale horizontal structures at greater depths. In one-dimensional (1D) interpretations, equivalence problems arise where thin, highly conductive layers can produce similar apparent resistivity responses to thicker, less conductive layers with equivalent conductance, leading to ambiguity in layer thickness and resistivity estimates.^[97]^[98] Noise vulnerabilities further constrain MT applicability, especially in anthropogenically influenced environments. Cultural interference from power lines, pipelines, and other man-made sources generates strong, coherent electromagnetic signals that overwhelm natural MT fields, particularly at frequencies above 10 Hz, thereby limiting reliable data acquisition in urban or industrialized areas. Additionally, static shift distortions, caused by near-surface heterogeneities acting as galvanic capacitors, introduce frequency-independent multiplicative offsets to apparent resistivity curves, which can be partially mitigated through telluric field corrections but often require complementary data for full decoupling.^[99]^[100] In three-dimensional (3D) settings, MT faces increased complexities from galvanic distortions induced by topographic variations and near-surface conductivity contrasts, which perturb the electric field and complicate impedance tensor interpretations. Full 3D inversions, necessary for resolving such effects, incur high computational costs, often requiring hours to days of processing on multi-processor clusters for large datasets, due to the need for iterative forward modeling on discretized grids.^[101]^[102] Frequency coverage in MT surveys is hampered by inherent gaps in the natural source spectrum, notably the "dead band" between 0.1 and 10 Hz where signal strength is minimal, reducing continuity in resistivity-depth profiles and increasing uncertainty in mid-crustal imaging. During solar minima, the predictability and intensity of low-frequency sources diminish due to reduced geomagnetic activity, further degrading data quality for deep-sounding applications.^[44]^[103] Compared to active-source electromagnetic methods, MT is less effective for very shallow targets (<100 m) where higher-frequency signals are prone to cultural noise and near-surface attenuation, and for high-resistivity structures, as the method's sensitivity favors conductive anomalies over resistive ones.^[104]^[105]

Recent Advances

Recent advances in magnetotellurics (MT) have leveraged machine learning to enhance data processing and inversion efficiency. Neural network-based inversions, such as those using Bayesian Neural Networks (BNNs), have enabled the recovery of subsurface resistivity profiles with integrated uncertainty quantification, addressing the ill-posed nature of traditional methods.^[106] For instance, physics-informed deep learning models incorporating MT wave propagation laws have achieved high accuracy in 3D inversions while reducing computational demands; one study using a Mix Transformer network (MT-MitNet) for forward modeling reported a reduction in computation time from over 2800 seconds to under 5 seconds per model, facilitating faster iterative inversions.^[107] Additionally, machine learning techniques like convolutional neural networks (CNNs) and long short-term memory (LSTM) networks have been applied for automatic noise classification and denoising, improving signal-to-noise ratios in contaminated datasets by learning features from time-series data without manual intervention.^[108] Dense array technologies have expanded MT surveys into challenging environments, particularly urban areas prone to cultural noise. Recent deployments of broadband MT instruments in dense configurations, with station spacings as fine as 100-500 meters, have enabled high-resolution imaging of shallow structures in seismically active urban zones, such as the Huaibei plain in China, where 3D MT arrays revealed deep seismogenic environments amid infrastructure interference.^[109] These advancements have improved resolution of lateral heterogeneities, with examples demonstrating resistivity contrasts at depths up to 5 km in city-scale profiles. Joint inversion techniques integrating MT with seismic and gravity data have gained traction for robust subsurface imaging, particularly in complex geological settings like desert basins. Such methods enhance model reliability by coupling complementary sensitivities.^[110]^[111] Space-based enhancements from satellite missions like ESA's Swarm have improved the characterization of long-period MT source fields since 2015, aiding global-scale conductivity modeling. By providing high-precision geomagnetic data, Swarm observations help separate ionospheric and magnetospheric signals from local MT responses, enhancing accuracy in low-frequency (periods >100 s) inversions for mantle studies.^[112] Recent analyses (2020-2025) have utilized Swarm magnetic field measurements to refine source field predictions, reducing bias in deep Earth resistivity estimates by incorporating tidal and core-generated variations.^[113] Open-source software developments have democratized advanced MT analysis. Updates to the ModEMpy package in 2021 introduced Bayesian frameworks for 3D inversions, enabling probabilistic uncertainty estimation through Markov chain Monte Carlo sampling integrated with deterministic solvers.^[114] The IRIS Consortium's MT data hub, launched in the early 2020s, provides community-accessible datasets from global arrays, supporting reproducible research and machine learning training on diverse real-world time-series.^[115] Emerging fiber-optic MT sensing prototypes promise distributed measurements for real-time monitoring. Distributed electric field sensing via fiber optics in borehole environments converts strain from electromagnetic fields into measurable signals, offering continuous profiling along cable lengths up to kilometers with resolutions better than 1 mV/km. As of 2025, initial prototypes combining magnetostrictive coatings on fibers for magnetic sensing have demonstrated feasibility for hybrid MT-seismic arrays, potentially revolutionizing permanent installations in geothermal and volcanic sites by enabling dense, low-maintenance networks.^[116]

References

[1]
Magnetotelluric Method - an overview | ScienceDirect Topics
The MT method utilizes naturally occurring, broadband electromagnetic waves over the Earth's surface to image subsurface resistivity structure.
[2]
8. The Magnetotelluric Method | Electromagnetic Methods in Applied ...
In the magnetotelluric (MT) method, natural electromagnetic fields are used to investigate the electrical conductivity structure of the earth.<|control11|><|separator|>
[3]
Magnetotellurics — Electromagnetic Geophysics - EM GeoSci
The magnetotelluric (MT) method is a passive electromagnetic (EM) exploration method that measures orthogonal components of the electric and magnetic fields ...
[4]
Introduction to the magnetotellurics special issue of exploration ...
Sep 18, 2024 · The geophysical technique of magnetotellurics (MT) was first described in 1953 and the associated technique of geomagnetic depth sounding (GDS) ...
[5]
THE MAGNETOTELLURIC METHOD | GeoScienceWorld Books
SUMMARY. In the magnetotelluric (MT) method, natural electromagnetic fields are used to investigate the electrical conductivity structure of the earth.
[6]
[PDF] Basics of the magnetotelluric method
This relation is crucial as Equation (1.15) states that the skin depth is only due to the bulk conductivity of the overlying material and the period range used.Missing: physical | Show results with:physical
[7]
[PDF] Electromagnetic Methods. - Marine EM Laboratory
In uniform conductors, the rate at which the primary field is absorbed is exponential, and the characteristic length scale is given by the skin depth, which we ...
[8]
XV. The diurnal variation of terrestrial magnetism - Journals
Schuster Arthur. 1889XV. The diurnal variation of terrestrial magnetismPhilosophical Transactions of the Royal Society of London. (A.)180467–518http://doi ...
[9]
Magnetotellurics: Principles and practice | Journal of Earth System ...
Tikhonov A N 1950 On determining electrical characteristics of the deep ... Cite this article. Vozoff, K. Magnetotellurics: Principles and practice.
[10]
[PDF] BASIC THEORY OF THE MAGNETO-TELLURIC METHOD ... - MTNet
This “principle of similitude” greatly simplifies the representation of a master set of curves, such as is given for use in geologic interpretation. In addition ...
[11]
A practical, low-noise coil system for magnetotellurics
Abstract. Magnetotellurics is a geophysical technique which was developed by Cagnaird (1953) and Tikhonov (1950) and later refined by other scientists worldwide ...
[12]
Phoenix Geophysics Ltd. - MSTA CANADA
Since 1975, Phoenix Geophysics has provided instrumentation and services to clients in more than 75 countries and has carried out surveys for many of the world ...
[13]
Magnetotelluric Array - USArray
From 2006-2018, the EarthScope USArray magnetotelluric (MT) program collected observations of the Earth's naturally occurring electric and magnetic fields.
[14]
BASIC THEORY OF THE MAGNETO‐TELLURIC METHOD OF ...
The magneto-telluric method uses the ratio of electric and magnetic fields to determine subsurface structure and resistivities, and is ideal for large ...
[15]
Basic theoretical concepts (Chapter 2) - Practical Magnetotellurics
Practical Magnetotellurics - February 2005. ... Starting from Maxwell's equations, we derive the diffusion equation for the electric field.
[16]
Anomalous phases exceeding 90° in magnetotellurics: anisotropic ...
... doi.org/10.1046/j.1365-246X.2003.02050.x. Published: 01 October 2003 ... anisotropy, electrical conductivity, magnetotellurics, phases. Issue Section ...
[17]
Anisotropy and phase splits in magnetotellurics - ScienceDirect.com
Magnetotellurics in south western Sweden: evidence for electrical anisotropy in the lower crust? J. Geophys. Res., 93 (1988), pp. 7897-7907. View in Scopus ...
[18]
The Magnetotelluric Impedance Tensor and its Properties
Aug 6, 2025 · A concept of the frequency-normalized impedance (FNI) function is introduced for the visualization of all components of the magnetotelluric ...
[19]
Dictionary:Cagniard apparent resistivity - SEG Wiki
Oct 14, 2024 · In the magnetotelluric method, an apparent resistivity in ohm-meters at period T obtained from 0.2T |Z|2, where Z is the Cagniard impedance ( ...Missing: formula | Show results with:formula
[20]
magnetotelluric phase tensor | Geophysical Journal International
4. These arrows are graphical representations of the real part of the vertical magnetic field transfer function or tipper vector (K) defined by the equation. ...
[21]
[PDF] Dimensionality and directionality analysis of magnetotelluric data by ...
The skew angle b measures the deviation from two-dimensionality and it represents 'quasi-2D' (if —b— \ 38). In case of 3D, b is non-zero and the direction ...
[22]
Magnetotellurics: the CBB or phase tensor and Bahr's 1988 analysis
May 13, 2020 · The two major axes fail orthogonality by 1 7 . 0 ∘ , which CBB refers to as the “skew” angle. (d) ...
[23]
A hybrid fast Hankel transform algorithm for electromagnetic modeling
A hybrid fast Hankel transform algorithm has been developed that uses several complementary features of two existing algorithms: Anderson's digital ...
[24]
[PDF] ANALYSIS OF MAGNETOTELLURIC DATA USING THE REILLY ...
Thus the measured properties (contained in the impedance tensor) are presented as an apparent conductivity and an apparent permittivity tensor. This means of ...
[25]
Interpretation of 1D Vector Controlled-Source Audio-Magnetotelluric ...
Aug 8, 2025 · This paper presents 1D vector CSAMT interpretation. A full solution 1D CSAMT forward modeling has been developed and used to interpret both ...
[26]
[PDF] Joint Time-Lapse Acquisition and Inversion of Passive Seismic and ...
Because the wavelength of the P-wave is longer than that of the S-wave, their sensitivity kernels overlap even more. If in addition, only those P-phases and ...
[27]
Fast calculation of spatial sensitivity kernels for scattered waves in ...
Feb 24, 2022 · Here we present a method for rapidly calculating the P- and S-wave receiver function sensitivity kernels, based upon the shortest path method and Dijkstra's ...
[28]
https://www.energy.ca.gov/sites/default/files/2024-07/CEC-500-2024-075_0.pdf
[29]
[PDF] MT Acquisition Best Practices: Survey Design and Field Procedures
• Cow fences – analogue worse than digital. • Radar. Page 44. Noise Solutions - Cultural. • Avoid or minimize. • Filters remove much of 60Hz/harmonics. • Record ...
[30]
[PDF] Magnetotellurics an Introduction for Geologists
Magnetotellurics (MT for short) is an electromagnetic geophysical technique that uses natural electric fields from lightning sources, solar.
[31]
[PDF] MAGNETOTELLURIC SURVEY ON THE MT. SPURR ...
Dec 27, 2011 · Vegetation is estimated to be the single biggest factor in determining the location of the MT stations. Very steep topography, unstable ...
[32]
MT Geophysical Survey Methods - Zonge International
The MT geophysical survey method combines measurements of the earth's electric field and magnetic field over a wide band of frequencies.
[33]
Developing a Health and Safety Plan for Hazardous Field Work in ...
Aug 6, 2025 · This article details development of a HASP to cover field researchers and ship personnel conducting scientific research in a remote area of the ...
[34]
[PDF] Instrumentation for marine magnetotelluric and controlled source ...
The first sea-floor measurements using the MT method were carried out in the 1960s by Jean Filloux and colleagues. ... standards of amplifiers and so the ...
[35]
[PDF] “Nuts and bolts” of magnetotelluric instrumentation and field-work
Magnetic induction coil sensors are usually calibrated using a secondary coil that creates an accurately known magnetic field. Calibration of coils for ...Missing: digitizers | Show results with:digitizers
[36]
[PDF] Low-Noise-Ag_AgCl-Electric-Field-Sensor-System-for ... - EMGS -
We have demonstrated Ag/AgCl E field sensors with noise level less than 1 Hz nV down to 10mHz in the laboratory, and down to ~0.1Hz on the seabed. We have also ...
[37]
[PDF] Magnetotellurics: Principles and practice
The magnetotelluric method was proposed by Tikhonov (1950), by Kato and Kikuchi. (1950), by Rikitake (1950) and by Cagniard (1953), although the principle was.
[38]
[PDF] New possibilities of audiomagnetotellurics in mineral exploration
Another advantage is high productivity (15-30 minutes of recording per site) and as a result relatively low cost. That's why AMT is widely used for mineral ...
[39]
Ultra-low power high precision magnetotelluric receiver array based ...
Dense 3D magnetotelluric (MT) data acquisition owns the benefit of suppressing the static shift and topography effect, can achieve high precision and high ...Missing: nodal devices
[40]
Robust processing of magnetotelluric data in the AMT dead band ...
The energy sources for magnetotellurics (MT) at frequencies above 8 Hz are electromagnetic waves generated by distant lightning storms propagating globally ...
[41]
[PDF] Processing of noisy magnetotelluric data using digital filters and ...
Various noise-reduction techniques based on digital filters are discussed with special reference to persistent noise signals, e.g. from power lines, DC-operated ...
[42]
Survey — Electromagnetic Geophysics - EM GeoSci
Little energy is emitted between 0.1 and 10 Hz, creating a dead-band in the frequency spectrum for MT methods. Receivers: A typical MT station measures the ...
[43]
https://scispace.com/pdf/processing-of-noisy-magnetotelluric-data-using-digital-heou7dcafn.pdf
[44]
[PDF] BIRRP: Bounded Influence, Remote Reference Processing
Oct 1, 2001 · The BIRRP program computes magnetotelluric and geomagnetic depth sounding response. Page 2. 2 functions using a bounded influence, remote ...
[45]
[PDF] Remote reference magnetotelluric processing algorithm based on ...
Jul 20, 2017 · The magnetic coherence detects the data segment (gray dots in Fig. 10c) where the polarization direction of the electric and magnetic fields ...
[46]
On the stability of magnetotelluric transfer function estimates and the ...
Fixed-weight jackknife error bars were larger than those computed in the conventional fashion, but still significantly underestimated the true variability in ...
[47]
BIRRP -Bounded Influence Remote Reference Processing
The BIRRP program computes magnetotelluric and geomagnetic depth sounding response functions using a bounded influence, remote reference method.Missing: EMTF | Show results with:EMTF
[48]
Aurora: An open-source Python implementation of the EMTF ...
The Aurora software package robustly estimates single station and remote reference electromagnetic transfer functions (TFs) from magnetotelluric (MT) time ...Missing: Birrp | Show results with:Birrp
[49]
[PDF] Magnetotelluric strike rules - MTNet
Next we consider the form of the total impedance tensor when measured in the direction of the regional strike. For. 0 = 8,. Z(%) = L. -P2 sin 2AOZj,. ( ...
[50]
Stable 2D magnetotelluric strikes and impedances via the phase ...
The strike direction is provided by the phase tensor, which is known to provide unstable estimates for noisy data. We obtain stable directions in three steps.
[51]
Magnetotelluric measurements across the Sorgenfrei-Tornquist ...
The simple Swift skew (Swift 1967) generally lies between 0 and 1 for real data, indicating deviations from 1-D/2-D conductivity distribution. However, it is ...
[52]
Magnetotelluric Imaging of the Lithosphere Across the Variscan ...
Nov 27, 2017 · The Swift's skew is defined as the ration of the diagonal components to nondiagonal components of the impedance tensor (Swift, 1967, 1986).
[53]
The magnetotelluric phase tensor - Caldwell - 2004
Jul 15, 2004 · This tensor expresses how the phase relationships change with polarization in the general case where the conductivity structure is 3-D.
[54]
[PDF] Occam's inversion to generate smooth, two-dimensional models ...
Magnetotelluric (MT) data are inverted for smooth. 2-D models using an extension of the existing I-D algorithm, Occam's inversion. Since an MT data set.
[55]
Comparison of 1D magnetotelluric inversion using Levenberg ...
Apr 16, 2015 · One-dimensional (1D) inversion problem for magnetotelluric (MT) data are addressed using two different schemes: The Levenberg-Marquardt (LM) inversion and ...
[56]
Regularized magnetotelluric inversion based on a minimum support ...
Nov 21, 2017 · Among a number of possible stabilizing functionals, smoothing constraints are most commonly used, which produce spatially smooth inversion ...Missing: smoothness | Show results with:smoothness
[57]
Two‐dimensional magnetotelluric inversion using reflection seismic ...
Mar 29, 2017 · We present a novel 2-D magnetotelluric (MT) inversion scheme, in which the local weights of the regularizing smoothness constraints are ...
[58]
General formalism for the efficient calculation of the Hessian matrix ...
It is known that one can gain valuable information about model uncertainties from the analysis of Hessian matrix. ... magnetotellurics and controlled-source ...
[59]
Uncertainty and Resolution Analysis of 2D and 3D Inversion Models ...
Sep 24, 2019 · The model resolution matrix measures to what extent the unknown true solution maps into the preferred solution. However, most current ...Missing: Hessian | Show results with:Hessian
[60]
Magnetotellurics Point to Serpentinization as a Potential Source of ...
Sep 26, 2025 · This study helps clarify the conditions required for natural hydrogen production and highlights the Bulqizë area as significant for future ...Missing: seminal | Show results with:seminal
[61]
Magnetotellurics - Nord-West Ltd.
The first one (“AMT dead band”) is from 1 to 5 kHz, and the second one (“MT dead band”) lies in the gap between the two main types of MT sources (1 – 10 s) ...
[62]
Effect of metasomatism on the electrical resistivity of the lithospheric ...
Long period MT data are widely applied for LAB depth determination because of the contrast between the more resistive lithospheric mantle (generally few tens to ...
[63]
Influence of geomagnetic storms on the quality of magnetotelluric ...
Jul 20, 2022 · We combined the coherence between the electric and magnetic fields and the result of the MT sounding curve to evaluate the MT impedance quality ...Missing: EH | Show results with:EH
[64]
[PDF] 5.07 Geomagnetism - Marine EM Laboratory
Arthur. Schuster (1889) observed the relationship between the diurnal magnetic variation and Earth potentials and Sydney Chapman (1919) inferred that Earth's.
[65]
[PDF] Geomagnetism applications - USGS Publications Warehouse
Magnetotelluric Sounding of the Earth's Crust ... observatories, called INTERMAGNET (international real-time magnetic observatory network).
[66]
[PDF] Geomagnetic Induction Studies in Scandinavia - MTNet
Other papers will be concerned with the more traditional geomagnetic depth sounding (GDS) and magneto-telluric (MT) techniques, modelling studies, and a.
[67]
Field Testing of a Controlled-Source Wide Frequency Range ... - MDPI
Magnetotellurics suffers from limitations such as unstable natural sources and low signal-to-noise ratios. Controlled-Source magnetotellurics (CSMT) is a ...
[68]
Developing a New Ground Electric Field Model for Geomagnetically ...
Aug 9, 2025 · Nowcasting and validating Earth's electric field response to extreme space weather events using magnetotelluric data: Application to the ...
[69]
Modelling Internal and External Geomagnetic Fields Using Satellite ...
The aim of this article is to provide a brief review of current geomagnetic field modelling techniques, which are based on measurements of Earth's magnetic ...Missing: magnetotellurics | Show results with:magnetotellurics
[70]
CONTROLLED SOURCE AUDIO-FREQUENCY ... - GeoScienceWorld
Controlled source audio-frequency magnetotellurics (CSAMT) is a frequency-domain electromagnetic sounding technique which uses a fixed grounded dipole or ...Missing: commercialization | Show results with:commercialization
[71]
Controlled-Source Audio-Magnetotellurics (CSAMT)
Jan 24, 2025 · Transmitter frequencies traditionally range from 5 Hz to 10 kHz. The CSAMT survey system utilizes the versatile DIAS32-Gen2 receiver and ...Missing: explanation depths sources
[72]
Evaluation of a Semi-Airborne Electromagnetic Survey Based on a ...
Jan 7, 2022 · The concept of the method is to deploy high-power transmitters on the ground, which excite subsurface currents and induce strong magnetic fields ...
[73]
Rapid relaxation inversion of CSAMT data - Oxford Academic
In the far field, CSAMT data are indistinguishable from plane wave data and an obvious strategy is to use magnetotelluric (MT) inversion to interpret them. This ...
[74]
[PDF] Processing of Noisy Controlled Source Audio Magnetotelluric ...
To remove outliers in the form of transient noise spikes in the time series, a method based on comparing long-term and short-term averages (Bastani, 2001) has ...
[75]
[PDF] An overview of the ZTEM and AirMt airborne electromagnetic systems
Jun 2, 2012 · For both. ZTEM and AirMt, data are typically measured from 30 Hz to. 720 Hz, giving detection depths to 1 km or more, depending on the terrain ...
[76]
Marine magnetotellurics for base‐of‐salt mapping: Gulf of Mexico ...
A sea‐floor magnetotelluric (MT) survey was conducted over the Gemini subsalt petroleum prospect in the Gulf of Mexico (GOM) to demonstrate that the base of ...
[77]
Two-dimensional determinant inversion of marine magnetotelluric ...
The marine magnetotelluric (MT) method has been used to investigate seafloor geology for hydrocarbon exploration (Constable et al., 1998; Hoversten et al., 2000 ...
[78]
Lower crustal resistivity signature of an orogenic gold system
Aug 4, 2021 · In this paper we compiled thirty years of long-period magnetotelluric (MT) and geomagnetic depth sounding (GDS) data across western Victoria and ...
[79]
None
Nothing is retrieved...<|control11|><|separator|>
[80]
‪Knútur Árnason‬ - ‪Google Scholar‬
Imaging the magmatic system beneath the Krafla geothermal field, Iceland: A new 3-D electrical resistivity model from inversion of magnetotelluric data. B Lee, ...
[81]
[PDF] annual report 2016 - SEG.org
Jun 7, 2023 · The high cost of deepwater exploration motivated the development of commercial marine magnetotelluric. (MT) exploration in 1995, but it wasn ...
[82]
Controlled source electromagnetic technology and hydrocarbon ...
Aug 6, 2025 · (2010) but view the data in terms of commercial success rates. ... EM technology offers several benefits such as it avoids drilling dry wells ...
[83]
Groundwater exploration by magnetotelluric method within the ...
Feb 7, 2022 · The Magnetotelluric (MT) geophysical method has successfully been applied in groundwater exploration within the Birimian rocks of Mankessim.
[84]
Pre-injection magnetotelluric surveys at the Aquistore CO2 ...
Magnetotelluric (MT) soundings were conducted in a 4 × 4 km area at the Aquistore CO 2 sequestration site at Estevan, Saskatchewan, Canada in 2013, 2014 and ...Missing: post- | Show results with:post-
[85]
[PDF] The feasibility of MT tipper data to monitor CO2 storage sites
The EM method we focus on in this paper is magnetotellurics. (MT), also called natural source electromagnetics. MT uses natural fields as a source and records ...
[86]
A Lithosphere‐Asthenosphere Boundary and Partial Melt Estimated ...
Aug 20, 2020 · At a broad scale, we observe a resistive lid (>100 Ohm m) that ranges in thickness between 20 to 80 km at shallow depths. At greater depths ...
[87]
Constraints from magnetotelluric measurements on magmatic ...
Jun 1, 2013 · We next examine sensitivity in order to provide bounds on the resistivities of features in the mantle (Fig. 8). While the MT method is sensitive ...Missing: challenges | Show results with:challenges
[88]
Mantle dynamics beneath the East Pacific Rise at 17°S: Insights ...
Feb 17, 2006 · [2] The Mantle Electromagnetic and Tomography (MELT) experiment was undertaken to map the distribution of melt from its sources deep within the ...
[89]
[PDF] 98-031 Oceanus F/W 97 final
The goal of the MELT Experiment was to map basaltic melt, from its source within the mantle to the base of the oceanic crust at the mid-ocean ridge crest ...
[90]
Electrical image of subduction zone beneath northeastern Japan
Oct 8, 2015 · This indicates that fluids released from the subducting slab at depths of 80–100 km move upward toward the back-arc direction and that ...Abstract · Introduction · Experiment · Discussion
[91]
Three-dimensional magnetotelluric imaging of crustal fluids and ...
Dec 4, 2014 · Northeastern (NE) Japan can be classified as a typical subduction zone and is therefore a suitable test site for the study of fluid ...
[92]
2D resistivity model around the rupture area of the 2011 Tohoku-oki ...
May 15, 2023 · In this study, we modeled the two-dimensional resistivity distribution across the slip area based on ocean-bottom electromagnetic measurements.
[93]
[PDF] Temporal Resistivity Change of Crustal Resistivity Structure Before ...
Temporal Resistivity Change of Crustal Resistivity Structure Before and After the 2011. Tohoku Earthquake ... Keywords: resistivity, temporal variation, ...
[94]
Structure and tectonics of the northwestern United States from ...
Sep 15, 2014 · The magnetotelluric component of the EarthScope USArray program has covered over 35% of the continental United States.
[95]
Magnetotelluric evidence for the multi-microcontinental composition ...
Aug 4, 2020 · We present three profiles that reveal the geoelectrical structure of eastern South China. Distinct conductive interfaces oriented NE–SW are identified.
[96]
Progress of Deep Geological Survey Project under the China ...
Magnetotelluric array of South China. 5.2. Determination of the boundary of Jiangnan Orogen. The gravity-based multiscale edge detection is mainly used to ...Missing: SINOPADE | Show results with:SINOPADE
[97]
[PDF] Simple Equivalence Analysis for Magnetotelluric (MT) Depth ...
We studied the resolving capability of MT by using simple synthetic 1D models where a thin conductive layer is embedded in a homogeneous more resistive medium.Missing: lateral | Show results with:lateral
[98]
Exploring nonlinear inversions: A 1D magnetotelluric example
A classic example that demonstrates the nonuniqueness of MT data is the equivalence of the conductivity-thickness product (conductance) of a thin layer. If we ...Missing: sheet | Show results with:sheet
[99]
Magnetotelluric Noise Attenuation Using a Deep Residual ... - MDPI
The MT signal is non-stationary and broadband. It can be easily affected by cultural noises, from power lines, railway lines, road vibration, and vehicles [6].
[100]
Correcting for static shift of magnetotelluric data with airborne ... - SE
Galvanic distortions of magnetotelluric (MT) data, such as the static-shift effect, are a known problem that can lead to incorrect estimation of ...Missing: noise vulnerabilities cultural interference
[101]
[PDF] Measurements and Modeling of Topography and Galvanic Distortion
Here, we present a joint inversion for 3D resistivity and site-specific distortion parameters, from MT impedances and vertical magnetic transfer function (Tz) ...Missing: complexities | Show results with:complexities
[102]
[PDF] 3D Inversion of Magnetotelluric (MT) Resistivity Data from Krýsuvík ...
Feb 11, 2013 · The inversion is a very heavy computational task and each iteration took about 5 hours and the total computing time was about 150 hours.
[103]
[PDF] Geomagnetic Activity and its Influence on Magnetotelluric Transfer ...
Quiet geomagnetic periods are characterized by a low signal level which can lead to the degradation of response quality due to the increased uncertainty on ...<|control11|><|separator|>
[104]
Two-Dimensional Magnetotelluric Modelling of Ore Deposits
Dec 15, 2017 · However, for 2D and 3D models, there is limited sensitivity to the normal (layered background) resistivity distribution, mostly in the form of ...
[105]
The advantages of complementing MT profiles in 3-D environments ...
Sep 23, 2016 · However, issues in recovering high electrical resistivity anomalies are somewhat of a limitation for the MT method in certain resource ...
[106]
Inversion of Magnetotelluric Data using Bayesian Neural Networks
Aug 5, 2025 · This paper explores using Bayesian Neural Networks (BNNs) for Magnetotelluric (MT) inversion, a geophysics technique, to recover resistivity ...Missing: 2021 | Show results with:2021
[107]
Magnetotelluric Deep Learning Forward Modeling and Its ... - MDPI
Jul 23, 2023 · Using the trained MT-MitNet network can effectively reduce the forward computation time. However, the Occam inversion program runs in the ...
[108]
Magnetotelluric data denoising method combining two deep ...
Jan 4, 2023 · In this paper, two deep learning methods, CNN and LSTM, are combined to denoise MT data. On the one hand, the use of CNN for signal-noise ...
[109]
Utilizing dense magnetotelluric array to analyze deep seismogenic ...
The study area is primarily located within urban zones, characterized by dense transportation infrastructure, including roads and railways. Residential ...
[110]
Drone-Enabled 3D Magnetometric Resistivity Imaging for Geological ...
This study proposes a novel drone-based semi-airborne total-field magnetometric resistivity (SA-TFMMR) system for high-resolution detection of conductive ...
[111]
Enhancing electrical structure in magnetotelluric inversion by the ...
Dec 2, 2024 · A framework for 3-D joint inversion of MT, gravity and seismic refraction data. ... Received: 05 October 2024; Accepted: 11 November 2024;
[112]
Seismically constrained gravity inversions reveal magmatic and ...
Sep 20, 2023 · Seismically constrained gravity inversions reveal magmatic and metamorphic processes at a major lithospheric boundary in northwestern Australia.
[113]
Core field changes from eleven years of Swarm satellite observations
The Swarm satellite trio has provided global vector magnetic field measurements, with high precision and absolute accuracy, for the past eleven years.
[114]
Satellite monitoring of long period ocean-induced magnetic field ...
Dec 2, 2024 · We present the expected spatio-temporal characteristics of the ocean-induced magnetic field (OIMF) at satellite altitude on periods of months to decades.
[115]
A Reduced Order Approach for Probabilistic Inversions of 3D ...
Nov 13, 2021 · We present a novel strategy to invert 3D magnetotelluric (MT) data together with other data sets in a fully probabilistic manner Perform the ...
[116]
Magnetotelluric Arrays | SAGE - IRIS
The 2003-2018 EarthScope program MT observations were used to study different geologic features and processes across North America but are critical to ...
[117]
https://ui.adsabs.harvard.edu/abs/2023AGUFM.S41C0318Y/abstract