Three‐dimensional inversion of magnetotelluric data for a resistivity model with arbitrary anisotropy

Author(s):  
Wenxin Kong ◽  
Handong Tan ◽  
Changhong Lin ◽  
Martyn Unsworth ◽  
Benjamin Lee ◽  
...  
2009 ◽  
Vol 46 (2) ◽  
pp. 139-154 ◽  
Author(s):  
Erşan Türkoğlu ◽  
Martyn Unsworth ◽  
Dinu Pana

Geophysical studies of upper mantle structure can provide constraints on diamond formation. Teleseismic and magnetotelluric data can be used in diamond exploration by mapping the depth of the lithosphere–asthenosphere boundary. Studies in the central Slave Craton and at Fort-à-la-Corne have detected conductors in the lithospheric mantle close to, or beneath, diamondiferous kimberlites. Graphite can potentially explain the enhanced conductivity and may imply the presence of diamonds at greater depth. Petrologic arguments suggest that the shallow lithospheric mantle may be too oxidized to contain graphite. Other diamond-bearing regions show no upper mantle conductor suggesting that the correlation with diamondiferous kimberlites is not universal. The Buffalo Head Hills in Alberta host diamondiferous kimberlites in a Proterozoic terrane and may have formed in a subduction zone setting. Long period magnetotelluric data were used to investigate the upper mantle resistivity structure of this region. Magnetotelluric (MT) data were recorded at 23 locations on a north–south profile extending from Fort Vermilion to Utikuma Lake and an east–west profile at 57.2°N. The data were combined with Lithoprobe MT data and inverted to produce a three-dimensional (3-D) resistivity model with the asthenosphere at 180–220 km depth. This model did not contain an upper mantle conductor beneath the Buffalo Head Hills kimberlites. The 3-D inversion exhibited an eastward dipping conductor in the crust beneath the Kiskatinaw terrane that could represent the fossil subduction zone that supplied the carbon for diamond formation. The low resistivity at crustal depths in this structure is likely due to graphite derived from subducted organic material.


2020 ◽  
Author(s):  
Qinghua Huang ◽  
Tao Ye ◽  
Xiaobin Chen

<p>Influenced by the extrusion of Tibetan blocks and Indo-Burmese collision, the region in western Yunnan is associated with active seismicity and Quaternary volcanoes. Based on broadband magnetotelluric data collected in western Yunnan, we obtain a three-dimensional crustal electrical resistivity model after various data processing and three-dimensional inversion test. The above resistivity model reveals the seismogenic structures of the moderate and strong earthquakes in this tectonic region. We investigate the possible relationship between the seismicity and the electrical structure in western Yunnan region. The results indicate that earthquakes in this region tend to occur in the transition zone between the resistive and conductive structures. Our results also show that one resistive body imaged at the mid-lower crust may have blocked the previously proposed crustal channel flow along this intra-continental block boundary to the east of Tibetan Plateau. Our resistivity model suggests a bifurcation of the crustal flow in western Yunnan. This bifurcated crustal flow structure may play an important dynamical role in the seismogenesis of the earthquakes in western Yunnan.</p>


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Hiroshi Ichihara ◽  
Toru Mogi ◽  
Hideyuki Satoh ◽  
Yusuke Yamaya

Abstract The Hidaka collision zone, the collision boundary between the NE Japan and Kurile arcs, is known to be an ideal region to study the evolution of island arcs. The hypocenter of the 2018 Hokkaido Eastern Iburi earthquake (Mw 6.6) in the western part of the Hidaka collision zone was unusually deep for an inland earthquake, and the reverse fault that caused the earthquake has an uncharacteristically steep dip. In this study, we used three-dimensional inversion to reanalyze broadband magnetotelluric data acquired in the collision zone. The inverted resistivity model showed a significant area of high resistivity around the center of the collision boundary. We also identified a conductive zone beneath an area of serpentinite mélange in a zone of high P–T metamorphic rocks west of the high-resistivity zone. The conductive zone possibly reflects areas rich in pore fluids related to the formation and elevation of the serpentinites. Sensitivity tests indicated the need for additional magnetotelluric survey data to delineate the resistivity distribution around the epicentral area of the 2018 earthquake although the resistivity model showed a conductive zone in this area.


2020 ◽  
Vol 17 (3) ◽  
pp. 349-360
Author(s):  
Xiang-Zhong Chen ◽  
Yun-He Liu ◽  
Chang-Chun Yin ◽  
Chang-Kai Qiu ◽  
Jie Zhang ◽  
...  

2021 ◽  
Vol 73 (1) ◽  
Author(s):  
Matthew J. Comeau ◽  
Michael Becken ◽  
Alexey V. Kuvshinov ◽  
Sodnomsambuu Demberel

AbstractCrustal architecture strongly influences the development and emplacement of mineral zones. In this study, we image the crustal structure beneath a metallogenic belt and its surroundings in the Bayankhongor area of central Mongolia. In this region, an ophiolite belt marks the location of an ancient suture zone, which is presently associated with a reactivated fault system. Nearby, metamorphic and volcanic belts host important mineralization zones and constitute a significant metallogenic belt that includes sources of copper and gold. However, the crustal structure of these features, and their relationships, are poorly studied. We analyze magnetotelluric data acquired across this region and generate three-dimensional electrical resistivity models of the crustal structure, which is found to be locally highly heterogeneous. Because the upper crust (< 25 km) is found to be generally highly resistive (> 1000 Ωm), low-resistivity (< 50 Ωm) features are conspicuous. Anomalous low-resistivity zones are congruent with the suture zone, and ophiolite belt, which is revealed to be a major crustal-scale feature. Furthermore, broadening low-resistivity zones located down-dip from the suture zone suggest that the narrow deformation zone observed at the surface transforms to a wide area in the deeper crust. Other low-resistivity anomalies are spatially associated with the surface expressions of known mineralization zones; thus, their links to deeper crustal structures are imaged. Considering the available evidence, we determine that, in both cases, the low resistivity can be explained by hydrothermal alteration along fossil fluid pathways. This illustrates the pivotal role that crustal fluids play in diverse geological processes, and highlights their inherent link in a unified system, which has implications for models of mineral genesis and emplacement. The results demonstrate that the crustal architecture—including the major crustal boundary—acts as a first‐order control on the location of the metallogenic belt.


2019 ◽  
Vol 219 (1) ◽  
pp. 129-147 ◽  
Author(s):  
M Lajaunie ◽  
J Gance ◽  
P Nevers ◽  
J-P Malet ◽  
C Bertrand ◽  
...  

SUMMARY This work presents a 3-D resistivity model of the Séchilienne unstable slope acquired with a network of portable resistivimeters in summer 2017. The instrumentation consisted in distributed measuring systems (IRIS Instruments FullWaver) to measure the spatial variations of electrical potential. 23 V-FullWaver receivers with two 50 m dipoles have been deployed over an area of circa 2 km2; the current was injected between a fixed remote electrode and a mobile electrode grounded successively at 30 locations. The data uncertainty has been evaluated in relation to the accuracy of electrodes positioning. The software package BERT (Boundless Electrical Resistivity Tomography) is used to invert the apparent resistivity and model the complex data set providing the first 3-D resistivity model of the slope. Stability tests and synthetic tests are realized to assess the interpretability of the inverted models. The 3-D resistivity model is interpreted up to a depth of 500 m; it allows identifying resistive and conductive anomalies related to the main geological and hydrogeological structures shaping the slope. The high fracturation of the rock in the most active zone of the landslide appears as a resistive anomaly where the highest resistivity values are located close to the faults. A major drain formed by a fault in the unaltered micaschist is identified through the discharge of a perched aquifer along the conductive zone producing an important conductive anomaly contrasting with the unaltered micaschist.


Geophysics ◽  
1985 ◽  
Vol 50 (10) ◽  
pp. 1610-1617 ◽  
Author(s):  
Simon Spitz

A serious limitation to conventional data analysis is that the data refer mainly to elongated bodies. When three‐dimensional distortions are present, quantitative interpretation based only on the off‐diagonal elements of the conventionally rotated impedance tensor is inadequate, because these off‐diagonal elements are insensitive to the tensor trace. The impedance tensor eigenstate formulation proposed in the literature defines a complete set of parameters suitable for recognition of three‐dimensionality. Generally, though, the eigenvalues do not stand for the off‐diagonal elements of an impedance tensor measured in a physical coordinate system. It is shown how the eigenvalues are modified when the relationship between coordinate system rotations and the eigenstate formulation is clarified. A generalization of the conventional analysis results, but the rotation angle obtained is neither unique nor complete To improve the situation, two new analytical rotation angles are proposed. These angles define two complete intrinsic coordinate systems suitable for magnetotelluric data analysis when a general three‐dimensional structure is involved.


Geophysics ◽  
2012 ◽  
Vol 77 (1) ◽  
pp. E33-E42 ◽  
Author(s):  
Arild Buland ◽  
Odd Kolbjørnsen

We have developed a Bayesian methodology for inversion of controlled source electromagnetic (CSEM) data and magnetotelluric (MT) data. The inversion method provided optimal solutions and also the associated uncertainty for any sets of electric and magnetic components and frequencies from CSEM and MT data. The method is based on a 1D forward modeling method for the electromagnetic (EM) response for a plane-layered anisotropic earth model. The inversion method was also designed to invert common midpoint (CMP)-sorted data along a 2D earth profile assuming locally horizontal models in each CMP position. The inversion procedure simulates from the posterior distribution using a Markov chain Monte Carlo (McMC) approach based on the Metropolis-Hastings algorithm. The method that we use integrates available geologic prior knowledge with the information in the electromagnetic data such that the prior model stabilizes and constrains the inversion according to the described knowledge. The synthetic examples demonstrated that inclusion of more data generally improves the inversion results. Compared to inversion of the inline electric component only, inclusion of broadside and magnetic components and an extended set of frequency components moderately decreased the uncertainty of the inversion. The results were strongly dependent on the prior knowledge imposed by the prior distribution. The prior knowledge about the background resistivity model surrounding the target was highly important for a successful and reliable inversion result.


Sign in / Sign up

Export Citation Format

Share Document