Magnetotelluric Data Research Articles

Geophysical transects in the Abitibi greenstone belt of Canada from the mineral-exploration-oriented Metal Earth project

The Metal Earth project integrates geophysics, geology, geochemistry, and geochronology to improve the understanding of metal endowment in Precambrian terranes. Magnetics (airborne), gravity, magnetotellurics, and reflection seismic methods are the primary geophysical tools employed. Data were collected along 13 transects in the initial phase of the project. All geophysical tools are crucial for understanding the structure of the shallow, middle, and deeper crust and identifying pathways along which the constituents of critical minerals might have migrated from a source to a deposit. The magnetic data are used predominantly to help map the geology away from the transects, and the gravity data are useful for extending the near-surface geology to depths up to 8 km. The magnetotelluric data show the upper Archean crust to about 10 km as highly resistive, except for some conductive subvertical zones that correspond to major deformation zones, many of which are known to be metalliferous. This suggests that these conductive zones could have been hydrothermal fluid pathways feeding the mineral deposits. These zones can be traced to larger horizontal conductive zones in the midcrust. The seismic reflection data are consistent with and complement this: the upper crust is primarily nonreflective; however, the midcrust shows many horizontal reflectors, usually with a consistent dip to the north. Processing crooked-line seismic data is problematic, and techniques have been developed to improve the imaging, including multifocusing, 3D processing, full-waveform inversion, and cross-dip moveout methods. Passive seismic data have also been collected. Ambient-noise surface-wave tomography can be used to infer broad zones of similar seismic velocity between major reflectors, while receiver function analysis has been used to identify deeper structures such as horizontal features at or below the Moho and a dipping structure evident to about 70 km depth.

High-Speed Railways Interference Signal Characteristics and Multiple Remote References Denoising of Magnetotelluric Data in Jizhong Depression, China

In the economically developed Beijing–Tianjin–Hebei region, magnetotelluric data are susceptible to contamination from cultural noise, which can be caused, for example, by urban stray currents, high-speed railways, or high-voltage lines. The multiple remote references method is an effective tool that can be used to suppress interference and improve signal-to-noise ratios. Therefore, this paper first introduces the basic principles of multiple remote references and then takes high-speed railway noise as an example. The characteristics of the time domain and frequency domain of the high-speed railway noise signals are analyzed. Then, we use two remote reference stations (with a single remote reference and multiple remote references) to process the data interfering with the high-speed railway and compare the results. Finally, the multiple remote references method is used to process the data for the entire section. Coupled with the known geological and seismic data, the inversion results well-reflect the deep underground geological structure.

Retrieval of Subsurface Resistivity from Magnetotelluric Data Using a Deep-Learning-Based Inversion Technique

Inversion is a fundamental step in magnetotelluric (MT) data routine analysis to retrieve a subsurface geoelectrical model that can be used to inform geological interpretations. To reduce the effect of non-uniqueness and local minimum trapping problems and improve calculation speeds, a data-driven mathematical method with a deep neural network was developed to estimate the subsurface resistivity. In this study, a deep learning (DL) inversion technique using a revised multi-head convolutional neural network (CNN) architecture was investigated for MT data analysis. We created synthetic datasets consisting of 100,000 random samples of resistivity layers to train the network’s parameters. The trained model was validated with independent noised datasets, and the predicted results displayed reasonable accuracy and reliability, which demonstrates the potential application of DL inversion for real-world MT data. The trained model was used to analyze MT data collected in the southwestern Athabasca Basin, Canada. The calculated results from the DL method displayed a detailed subsurface resistivity distribution compared to traditional iterative inversion. Since this approach can predict a resistivity model without multiple forward modeling operations after the CNN model is created, this framework is suitable to speed up the computation of multidimensional MT inversion for subsurface resistivity.

Three-dimensional electrical structure and seismogenic environment of the crust–mantle in the Lushan earthquake region, China

The Longmenshan (LMS) southern section is the front of the eastwards material migration from the Qinghai–Tibet Plateau. To understand whether the 2022 Ms6.1 Lushan earthquake was a strong aftershock of the 2013 Ms7.0 Lushan earthquake and why the aftershocks of the two earthquakes were concentrated relative to the 2008 Wenchuan earthquake, a three-dimensional electrical structural model of the Lushan earthquake region is established based on magnetotelluric (MT) array data. The model shows high-conductor layers (HCL) in the Songpan–Ganzi terrane (SGT) middle crust, Sichuan Basin (SCB) shallow surface, western Sichuan foreland depression (WSCFD) and LMS Moho. The analysis suggests that the HCL in the study area at a depth of 10 km is caused by saline fluids in the fracture zone or porous sedimentary rock voids, while the high-conductor body in the middle and lower crust and upper mantle is caused by partial melting of rocks. The Lushan earthquake area is subject to compressional deformation horizontally as the eastwards migration of Qinghai–Tibet Plateau material is blocked by the SCB and vertically by the upwards extension of high-conductor plastic material at depth. The aftershocks of the Lushan earthquake area are mainly confined within the Baoxing–Zhonglin–Daxing–Lushan high-conductor region, which has noticeable electrical structural differences compared to the overall high-resistivity structure of the LMS southern section; these differences explain why the aftershocks of the two earthquakes did not expand in the northeast and southwest directions of LMS. Because the 2022 and 2013 Lushan earthquakes both occurred in the same high-conductor region in the deep part of the Dachuan–ShuangShi fault to the Danyi blind fault, the 2022 Lushan earthquake was a strong aftershock of the 2013 Lushan earthquake in terms of electrical structural characteristics.

Comparison study of different forward modeling approaches by using the finite-element method for 3D magnetotelluric data in anisotropic media

We have examined three approaches for analyzing 3D magnetotelluric (MT) data with general anisotropy. These included the direct electromagnetic (EM) approach, which solved the electric field Helmholtz equation directly, and the Coulomb-gauged and Ungauged approaches. The three approaches are discretized using finite-element (FE) methods. First, the correctness of the three approaches is validated using a typical isotropic 3D model and a 2D anisotropic model. The FE solutions of each approach obtained by the direct solver MUMPS are compared with the analytical solutions to demonstrate that all approaches generated accurate results. Subsequently, an anisotropic 3D model is designed to investigate the forward characteristics of the different approaches. The FE solutions are obtained using the quasi-minimum residual (QMR) solver with a symmetric successive overrelaxation (SSOR) preconditioner in this model. The forward modeling results of each approach were compared in terms of stiffness matrices, the convergence of the iterative solver, and apparent resistivity and phase. We find that the FE solution of the Coulomb-gauged and Ungauged approaches produces larger stiffness matrices, but the proportion of non-zero elements in the Coulomb-gauged approach is the least. The Ungauged approach can yield faster solutions over a wide range of frequencies. Finally, a complex numerical example is provided to demonstrate the good convergence behavior of the Ungauged approach.

A novel method for processing noisy magnetotelluric data based on independence of signal sources and continuity of response functions

Various schemes for processing magnetotelluric (MT) data have been reported aiming at suppressing the strong effect of artificial electromagnetic noises, especially coherent noise that is correlated between electric and magnetic time series. Many of the recent denoising schemes are based on decomposing MT data into the responses of the natural signals and noises. However, a steady differentiation of the natural signals from noises independent of any empirical choice of parameter setting is critically important. In addition, improper subtraction of values from the separated signals can lead to losing useful values of the natural signals or missing noise-affected values, which may result in failure in deriving the true MT responses. We propose a novel data-processing method that applies frequency-domain independent component analysis (FDICA) to both the local MT data and the reference magnetic data. Among the separated signals, the proposed method can quantitatively distinguish the natural signals from the noise-affected components by calculating the ratio of the cross-power spectrum with the reference data to the auto-power spectrum for each component. When determining which values to subtract from the separated signals, we introduce an evaluation index with respect to two characteristics of the MT response function: stationarity in the time domain and smoothness in the frequency domain. We conducted two types of experiments: with MT time series severely contaminated by synthetic coherent noises and with MT field data interfered with direct-current (DC) railways. Consequently, we confirmed the noise-suppression performance superiority of the proposed method over the conventional methods of MT data processing.

Three-dimensional anisotropic inversion and electrostratigraphic imaging of marine magnetotelluric data to understand the control of crustal deformations by pre-existing lithospheric structures in the Mexican Ridges Fold belt, Southwestern Gulf of Mexico

SUMMARY3-D imaging of the lithosphere in the Mexican Ridges fold belt is important for understanding how the crustal deformations in this basin relate to deep tectonic processes and structures inherited from extinct Jurassic seafloor spreading. Here, we use broad-band (0.0001–0.4642 Hz) marine magnetotelluric data from the basin to reconstruct the 3-D anisotropic resistivities of the lithosphere and their spatial gradients. The resistivity gradients maxima enabled independent definition of important geological boundaries (seen on collocated seismic reflection data) and estimation of crustal thickness. We found anomalous layered zones of low resistivity and high electrical anisotropy at 5–8 km depth (coinciding with the regional detachment zone in Eocene shales in 3-D seismic data) and in the upper mantle which we interpret as indicating intense deformation and/or recent magmatic influence. We also found a banded crystalline basement structure across the fossil spreading centre comprising WSW–ENE trending, 6–10 km wide, electrically resistive subvertical sheets with conductive and anisotropic borders, which merge into a basal resistive stock-like body at 15–20 km depth. These are cut or bounded by later NNW trending major faults. These WSW and NNW structural trends correlate with the previously interpreted transform and normal faults that formed during the Late Jurassic opening of the Gulf of Mexico only if rotated clockwise by 25–30°. Surprisingly, the rugged thrust-related seabed is offset at the projected positions of the steep resistive-conductive basement sheets (which also have spatially coincident high magnetic intensity and seismicity) enabling us to infer they represent magmatic intrusions facilitated by pre-existing faults. Their conductive borders spatially coincide with possibly fluid-filled vertical fracture-sets in the overlying sediments seen in seismic data which we interpret as hydrothermal fluid pathways. We infer that a magmatic body recently intruded our study area, its ascent controlled by pre-existing basement structures, and influenced the deformation of the Neogene sequences and the seafloor topography.

Geoelectric Model of the Central Part of the North Caucasus: Three-Dimensional Inversion

Abstract ––Interpretation of magnetotelluric data within the North Caucasus marginal massif, the uplift of the Main Range, and the Stavropol Arch, performed on the basis of one- and two-dimensional inversions of the impedance tensor components and three-dimensional mathematical modeling, makes it possible to compile the starting and test models that are necessary for a three-dimensional inversion. It is shown by the three-dimensional inversion of all the impedance tensor components, carried out on a test three-dimensional model, that it is possible to estimate the parameters of conducting model blocks. The three-dimensional inversion of all the experimental impedance matrix components, performed with account for the results obtained on the test three-dimensional model, significantly corrects the parameters of the conducting blocks identified at the stages of one- and two-dimensional inversions, as well as the three-dimensional mathematical modeling of magnetotelluric fields in the central part of the North Caucasus. In the resulting three-dimensional geoelectrical model of the region, the position of low-resistance blocks correlates with the location of suture zones, deep faults, volcanic chambers, and domains with converted earthquake waves, whose velocities are reduced and whose absorption is increased. This can be explained by the dependence of the conductivity of crustal blocks on water saturation. Earthquake hypocenters are grouped near low-resistance anomalies.

Geoelectric Model of the Northwestern Caucasus: Three-Dimensional Inversion

Interpretation of magnetotelluric observations made within the Northwestern Caucasus, by means of one-dimensional, two-dimensional inversion and three-dimensional mathematical modeling programs, has helped build testing and starting geoelectric models that are necessary to test and adapt the program of three-dimensional inversion of the impedance tensor components. The performed three-dimensional inversion of the experimental magnetotelluric data has significantly changed the parameters of conducting blocks within the folded structures of the region, identified at the previous stages of MT data interpretation. In the resulting three-dimensional geoelectric model, the position of low-resistance blocks correlates with: suture zones, deep faults, mud volcanoes, and domains characterized by an increased absorption of transverse and longitudinal seismic waves. The electrical resistivity of the most low-resistance anomalies is explained by the degree of their saturation with the water fraction of the fluid.

Undulating electrical lithosphere–asthenosphere boundary beneath Northeast China; as revealed by long–period magnetotelluric data

The geometry of continental lithosphere–asthenosphere boundary (LAB) significantly affects the movements of tectonic plates and flow of the upper mantle. However, the properties and geometry of this boundary beneath tectonically active regions in Phanerozoic continents are mysterious and elusive. Ninety–three long–period magnetotelluric (LMT) observations across Northeast China (NEC) provide a unique opportunity to image the shape of the LAB. Probabilistic resistivity–depth model inversion was employed to obtain the one–dimensional resistivity structure beneath each LMT site. The depth of the electrical LAB is constrained by the most probable interface at which the resistivity changes rapidly from high to low. The inversion results show that the electrical LAB is undulating beneath NEC. The North–South Gravity Lineament marks the transition from deep to shallow LAB. The thickest lithosphere occurs beneath the Great Xing'an Range (GXAR), where the LAB is ∼200 km deep. The electrical lithosphere is relatively thin in the east, with the thinnest lithosphere (∼50 km) being found beneath the Changbai Mountain volcano. We speculate that melt depletion, together with possible mantle downwelling are the likely causes for the formation of the thick GXAR lithosphere. The unresolved lithosphere beneath the northern Songliao Basin could be ascribed to the strong attenuation of electromagnetic waves in highly conductive sedimentary layers or the weak difference in the resistivity of the lithosphere and asthenosphere. The lithosphere steps that can form the edge-driven convection in the asthenosphere may play key roles in forming the widely-spread intraplate magmatism in NEC.

New constraints on the structure and composition of the lithospheric mantle beneath the Slave craton, NW Canada from 3-D magnetotelluric data – Origin of the Central Slave Mantle Conductor and possible evidence for lithospheric scale fluid flow

Previously collected magnetotelluric data were used to produce an improved 3-D resistivity model of the lithosphere beneath the Slave craton. A region of low resistivity (<10 Ωm) at ∼100 km depth was imaged beneath the central Slave craton, and identified as the Central Slave Mantle Conductor (CSMC) reported in previous studies. Earlier studies interpreted the CSMC as graphite films but new analysis from laboratory experiments indicates that it is unlikely that graphite films are stable/continuous at this depth. Using recent laboratory experiments we review alternative conduction mechanisms for the CSMC including (a) metasomatic phlogopite, (b) grain boundary sulphides and (c) high density fluids including hydrous carbonatite melts or brines. None of these hypotheses are without flaws requiring either unreasonably high-volume fractions or have petrological challenges. However, brines are the preferred explanation as they (1) can explain geophysical observations, (2) are observed in fibrous diamonds coincident with the CSMC and (3) require small volume fractions (<1%). We hypothesize their emplacement was related to a Mesozoic subduction event. This implies that upper mantle conductors should not be used as an indicator of carbon concentration in the mantle to prospect for diamondiferous regions. Our resistivity model indicates a lithospheric thickness of 210 ± 10 km beneath the Slave craton, consistent with seismic and xenolith studies. A water content of 10–150 ppm at 100–170 km depth in the lithospheric mantle is inferred from resistivity, broadly agreeing with estimates from mantle xenoliths from the central Slave craton lithosphere. Below 170 km the bulk lithosphere is "dry" (<10 ppm), which may explain why the cratonic root has been resistant to erosion by the underlying asthenosphere. Low resistivity fingers extend to surface in multiple locations that may represent paleo-fluid flow pathways through the lithosphere, broadly resembling those observed beneath Olympic Dam, Australia.

Three-dimensional audio magnetotelluric imaging of the Yangyi geothermal field in Tibet, China

The Yangyi geothermal field in southern Tibet has enormous potential for geothermal energy development. It is the high-temperature geothermal field at the highest elevation in the world, and the second largest high-temperature geothermal field in China. The total planned installed capacity of the geothermal power station is 32 MW. At present, no detailed three-dimensional (3D) geophysical survey has been carried out in the Yangyi geothermal field, and in-depth research is urgently needed. The main goal of our geophysical exploration is to study the geothermal structure of this region. In this study, the audio magnetotelluric (AMT) data of 85 stations were collected in the Yangyi geothermal field. The average distance between the stations was about 500 m, and the dimensions of the exploration area were 6 × 7 km. First, we analyzed the magnetotelluric data using the phase tensor method to qualitatively estimate the dimensionality of the subsurface resistivity structure. Then, by performing 3D inversion of 85 AMT data with relatively high quality, the 3D resistivity model from the surface to a depth of 2 km was determined, the high-temperature hydrothermal system under the Yangyi geothermal field was displayed for the first time. In addition, the main characteristics of the geothermal system were identified, and the genetic mechanism of the geothermal system was discussed based on the geological information and borehole data.

Integrated subsurface investigation for magmatic sulfide mineralization in Betul Fold Belt, central India

Magnetotelluric (MT), Electrical Resistivity Tomography (ERT), Time-Domain Induced Polarization (TDIP), Geochemical and Geological studies followed by drilling and down-hole logging were undertaken with in the Betul Fold Belt (BFB) in Central India, to demarcate zones of magmatic Ni-Cu-PGE sulfide mineralization. The BFB is predominantly composed of circular to elongate gabbro bodies of the Padhar Mafic-Ultramafic Complex, intruded into a sequence of bimodal volcanic rocks and quartzites. Near-surface samples of ultramafic rocks were subjected to precise geochemical analysis and scanned by an electron microscopy with an energy dispersive spectrometer (SEM-EDS). This work indicated the presence of pyrite, pyrrhotite, chalcopyrite, pentlandite, and minor amounts of W-Cd bearing boweiite and palladenite assemblage. These minerals are regarded as favorable to the occurrence of Ni-Cu-PGE sulfide mineralization. MT data derived from two profiles were analyzed and modeled using 2D and 3D inversion algorithms. The robust conductivity anomalies obtained from the MT model have been interpreted coupled with electrical tomography, geology, and geochemistry data. The near-surface shallow depth conductors observed in the ERT sections are interpreted as the sulfide mineralized zones. They corroborate the MT results. These conductive zones reflect the occurrence of the magmatic Ni-Cu-PGE bearing sulfide mineralization associated with rocks of the mantle-derived Padhar Mafic Ultramafic Complex. This geophysical data, in conjunction with petrological and geochemical analysis of drill core samples have allowed the identification of the origin and paragenesis of sulphide mineralization in the study area. Geochemical studies suggest that the parental magma was generated from a subduction modified, metasomatized and an enriched mantle source that was subsequently emplaced in a magmatic continental arc setting. The interpreted conductors, observed at shallow depths (∼200-300 m), have been generated by secondary hydrothermal fluid circulation leading to vein formation in the host Padhar Mafic-Ultramafic Complex. The MT and electrical tomography models delineate the geological boundaries of the sulfide-bearing mineralized deposits in the BFB.

Nature-inspired stochastic hybrid technique for joint and individual inversion of DC and MT data

Here, a new naturally-inspired stochastic nonlinear joint and individual inversion technique for integrating direct current (DC) and magnetotelluric (MT) data interpretation-based simulation of a swarm intelligence combo with specific capabilities for exploitation of the variable weight particle swarm optimizer (vPSO) and exploration of the grey wolf optimizer (GWO), vPSOGWO, is used. They are particularly notable for their capacity for information exchange while hunting for food. Through synthetic MT and DC data contaminated with various sets of random noise, the applicability of the anticipated vPSOGWO algorithm based joint and individual inversion algorithm was assessed. The field examples, collected from diversified different geological terrains, including Digha (West Bengal), India; Sundar Pahari (Jharkhand), India; Puga Valley (Ladakh), India; New Brunswick, Canada; and South Central Australia, have shown the practical application of the proposed algorithm. Further, a Bayesian probability density function (bpdf) for estimating a mean global model and uncertainty assessment in posterior; and a histogram for model resolution assessment have also been created using 1000 inverted models. We examined the inverted outcomes and compared them with results from other cutting-edge methodologies, including the GWO, PSO, genetic algorithm (GA), Levenberg–Marquardt (LM), and ridge-regression (RR). Our findings showed that the current methodology is more effective than the GWO, PSO, GA, LM, and RR techniques at consistently improving the convergence of the global minimum. In contrast to earlier approaches, the current cutting-edge strategy vPSOGWO offers an improved resolution of an additional significant crustal thickness of about 65.68 ± 1.96 km over the Puga Valley, in which the inverted crustal thickness determined by vPSOGWO agrees well with the published crustal thickness over the Puga Valley. The new technology brings simulations closer to genuine models by significantly reducing uncertainty and enhancing model resolution.

Geophysical Characterisation and Oil–Gas Resource Analysis of the Southern Huaying Mountain Fault Zone, Sichuan Basin, China

As the west rim of an ejective fold zone, the Huaying Mountain fault zone (HMFZ) in the eastern Sichuan Basin (SB) plays an important role in the tectonic evolution of the SB. The distribution and characteristics of HMFZ are strongly associated with tectonic activities and have greatly impacted the distribution of oil and gas reservoirs. However, its distribution and characteristics have remained poorly understood due to a lack of geophysical data, especially areal gravity survey and magnetotellurics (MT) survey, which are admittedly advantageous for detecting the edges of geological structures. Therefore, we carried out the ground geophysical surveys that areal gravity survey and MT survey, and acquired 1:250,000-scale real gravity data and MT data for the first time in this area. Optimized edge-detection methods were adopted to process the areal gravity data, allowing us to characterize the planar distribution of faults more reliably and convincingly. We found that the southern HMFZ is well developed and primarily trends in NNE and NE, whereas the subordinate faults trend in N-S and W-E. Vertical information for the faults extracted using the improved depth from the extreme points method revealed that the fault dominantly dipped to the SE, which was consistent with the results of MT inversion. Based on the spatial distribution of the faults, we further discussed the gravity anomaly, fault distribution, Luzhou palaeo uplift, and the distribution and characteristics of oil-gas resources, and found the convincing evidence to analysis the distribution of oil and gas resources in this region.

The Western Segment of the Precambrian Suture Between the Yangtze and Cathaysia Blocks: Constraints From Magnetotelluric Data in Southwest China

AbstractThe western segment of the suture zone between the Yangtze and Cathaysia blocks, which is the most important tectonic boundary related to the formation and rifting of south China, is enigmatic and not fully understood due to the sporadic exposure of Precambrian strata and ophiolites. Three‐dimensional electrical resistivity models derived from inversion of magnetotelluric data identified a lithospheric‐scale conductive zone extending northeastwards beneath the Youjiang basin, which was interpreted as the western segment of the suture zone. The high conductivity and coincident high magnetic anomalies closely match the location of Carlin‐type gold deposits, which can be explained by fluids and gold‐bearing sulfide minerals in a fossil suture zone. Inconsistent with the southeast‐dip resolved at the eastern segment of the suture zone (the Jiangshan‐Shaoxing fault) in the earlier study, the slightly north‐dipping geometry at the western suture zone implies the reactivation by northward subduction and closure of the Paleo‐Tethys Ocean.

Salair–Gornaya Shoria Junction (Northwestern Central Asian Orogenic Belt): Deep Structure and Tectonics from Magnetotelluric Data

Abstract —The Salair fold-thrust orogenic belt (Salair orogen, Salair) is located in the northwestern Altai–Sayan fold area within the Central Asian Orogenic Belt. The Salair orogen is an allochthon overriding the Kuznetsk Basin on a system of imbricate thrusts. The southern flank of the Salair thrust system is tectonically juxtaposed against the Gornaya Shoria terrane which differs markedly from Salair in its geological setting. The Salair and Gornaya Shoria terranes are separated by the Nenya-Chumysh Basin, a deep Mesozoic trough. The Salair orogen is composed of Cambrian–Early Ordovician island arc volcanic and sedimentary rocks, widespread garnet amphibolites and gneisses of the Angurep complex in its southern flank, and the Shalap subduction-related melange in the Alambai ophiolite suture. The southern Salair orogen and its junction with Gornaya Shoria have been imaged down to the lower crust by magnetotelluric (MT) soundings, which is an efficient tool for investigating the deep structure and tectonic history of orogenic areas. The MT surveys were performed at 25 stations on a 120 km long profile. MT data revealed an up to 70 km wide low-resistivity zone (a conductor) traceable till a depth of 20 km between the Salair and Gornaya Shoria terranes. The low-resistivity zone has a complex structure with its outer and interior boundaries dipping almost vertically. The conductor lies under several major geological structures: the Shalap melange, the Nenya-Chumysh Basin, and the NE trending Altai–Salair right-lateral strike-slip fault. The Altai–Salair fault, along which the Salair allochthon was displaced relative to Gorny Altai and Gornaya Shoria, joins the Salair system of imbricate thrusts. The Nenya-Chumysh Basin at the Salair–Gornaya Shoria junction is a deep trough having an asymmetric transversal profile with a steep western side and a shallower-dipping stepped eastern side. The southeastern flank of the basin is a wide area of thin sediments over the Paleozoic basement dipping gently in the northwestern direction. The revealed deep structure of the Nenya-Chumysh trough is consistent with its tectonic model implying an Early Cretaceous basin superposed over an early Jurassic pull-apart basin. Early Mesozoic motions on major faults is a regional-scale phenomenon known from many areas of southern West Siberia.

Crustal Structures of the Qimantagh Metallogenic Belt in the Northern Tibetan Plateau from Magnetotelluric Data and Their Correlation to the Distribution of Mineral Deposits

Crustal structure and fluid or melt originating in the deep crust and mantle are critical in regional magmatic mineral systems. However, the crustal structure and the processes that entrain and focus fluids from a deep-source region to a metallogenic belt remain relatively undisclosed. We present a magnetotelluric (MT) study of the eastern Qimantagh Metallogenic Belt (QMB) in the northern Tibetan Plateau. Data from 33 MT stations in two sections and 7 dispersed stations are acquired using a surface electromagnetic prospecting (SEP) system in frequency band ranges from 320 Hz to 0.00034 Hz. Data are converted by Bostick conversion and two-dimensional (2D) nonlinear conjugate gradient inversion. Our MT results reveal the geoelectrical crustal structure of the QMB, which consists of a southern low-resistivity domain that reflects the Kumukuri rift, a high-resistivity middle domain that represents the southern QMB in the central Kunlun belt, and a northern low-resistivity domain that covers the northern QMB and southwestern Qaidam block. We present a comprehensive tectonic and geophysical model of QMB based on the MT interpretation and geological analysis. We infer the high-resistivity domain as a reflection of a rigid crust and detached lithospheric mantle, this belt separate the QMB into northern and southern QMB. Most of the mineral deposits are found in the northern low-resistivity domain of QMB. Our study and findings provide an understanding of the tectonic evolution of the northern Tibetan Plateau, the crustal structure that controls the temporal and spatial distribution of magmatic rocks, and the geological signature associated with mineral deposits.

Smooth magnetotelluric impedance estimation by optimization

In urban and industrialised areas, magnetotelluric (MT) recordings often contain a significant amount of high-power near-field noise which, constituting non-Gaussian noise for most of the recording length, may lead to the failure of the common transfer functions (TF) estimation methods. The most evident consequence of this is the appearance of sharp points in the TF curves versus frequency, deteriorating their inherent smoothness. The TF's smoothness is in fact one of the main criteria adopted by MT practitioners to evaluate the data quality and effectiveness of the estimation procedure; when not reached by the estimation, it may be a-posteriori introduced by splines or filtering procedures, but these approaches lack physical soundness. We propose a new heuristic algorithm aimed at obtaining the maximum TF versus frequency smoothness, through an optimization scheme rejecting EM events in frequency-domain. The algorithm searches for frequency-dependent power thresholds to be applied to the events, rejecting those making the TF sharp. The TF's smoothness is the objective function to be minimized, and the model space is constituted by the infinite set of threshold vectors for each analysed frequency. After the process, the Ordinaly Least Squares (OLS) residual distribution is much closer to a Rayleigh one, allowing for an unbiased estimation. We applied the procedure both to synthetic and real MT data. Measured data has been collected for a geothermal exploration project in the East Africa Rift, where strong coherent noise affected the recordings. The combination of the technique with the remote-reference technique effectively reduced the bias due to the artificial signal. The physical consistency of the resulting TF has been tested by dispersion relationships.

Intra‐continental Orogeny: Insights from Magnetotelluric Data into the Mesozoic Uplift History of the Eastern Jiangnan Orogen in South China

AbstractDespite extensive efforts to understand the tectonic evolution of the Jiangnan Orogen in South China, the orogenic process and its mechanism remain a matter of dispute. Previous geodynamic studies have mostly focused on collisional orogeny, which is commonly invoked to explain the Jiangnan Orogen. However, it is difficult for such hypotheses to reconcile all the geological and geophysical data, especially the absence of ultrahigh‐pressure metamorphic rocks. Based on the magnetotelluric data, we present a group of resistivity models produced through the combination of two ‐dimensional and three‐dimensional inversions, revealing the geo‐electrical structures of Jiangnan and a typical collisional orogen. In our models, the resistive crust is separated into three parts by a prominent conductive layer with opposite dipping directions on both sides. A special thrust‐nappe system, which is different from that developed in a typical collisional process, is revealed in the Jiangnan Orogen. This structure suggests a process different from the simple collisional orogeny. To interpret our observations, an ‘intra‐continental orogeny’ is proposed to address the development of the Jiangnan Orogen in the Mesozoic. Furthermore, this ‘reworked’ process may contain at least two stages caused by the decoupling of the lithosphere, which is revealed by an extra conductive layer beneath Jiangnan.