Search

The controlled source audio-frequency magnetotelluric method (CSAMT) stands out for its economic efficiency and widespread application in geophysical monitoring. However, the separate inversion of time-lapse monitoring data encounters challenges in comparing and identifying abnormal changes due to variations in data fitting. Furthermore, the utilization of a method akin to Cagniard apparent resistivity for inversion necessitates the simultaneous observation of at least two components of the electromagnetic field, making it unsuitable for extensive three-dimensional observations. This paper proposes a 3D time-lapse electric field inversion algorithm for CSAMT, addressing the complexities in geophysical monitoring. The algorithm introduces two regularization factors and defines an objective function with both temporal and spatial constraints. Synthetic testing reveals the stability of the 3D time-lapse electric field inversion algorithm, demonstrating its effectiveness in delineating underground variations. This solution resolves the challenges posed by the independent inversion of time-lapse monitoring data. [ABSTRACT FROM AUTHOR]

To understand the complex crustal structure beneath the central Tarim Basin, we perform joint inversions of P-wave receiver functions, ambient seismic noise Rayleigh wave dispersion, and magnetotelluric data along a nearly N-S oriented linear array. In this contribution, we present a multistep flowchart to enhance the efficiency and reliability of the joint inversion. Synthetic examples indicate that this method is effective and appropriate for realistic conditions. By using these three complementary datasets, we can construct accurate 1-D models, particularly in the presence of thick overlying low-velocity layers. We apply this approach to data recorded at 45 stations within the Tarim Basin and generate quasi-2-D images that reveal distinct lateral variations along the N-S profile. In the near-surface region, we observe that the conductive and low-velocity sedimentary layer is thinner beneath uplift units and thicker beneath structural depressions, which correlates well with surface geology and previous studies. Additionally, we interpret several major faults, concealed crustal detachment layers, and discuss the geological implications. A conspicuous feature of the integrated images is the identification of a localized and deep-rooted doming structure beneath the Bachu uplift, characterized by high-velocity and high-resistivity anomalies. This structure also exhibits elevated density, intense magnetism and the thinnest overlying sedimentary layer. Consequently, we infer that these anomalies are probably linked to a hybrid basement consisting of plume-related intrusions and mafic dyke swarms. • We present a multistep flowchart to jointly invert receiver functions, ambient seismic noise, and magnetotelluric data. • Examination of the applicability on our joint inversion method for imaging the crust of Tarim basin. • Lateral heterogeneities in the near-surface structure and a localized doming structure beneath the Bachu Uplift. [ABSTRACT FROM AUTHOR]

Electrical anisotropy is increasingly recognized as an important aspect of the resistivity models required to explain magnetotelluric (MT) observations. However, a limited number of practical MT inversion algorithms that can consider anisotropy have been published to date. To address this problem, we have developed a three‐dimensional (3‐D) MT inversion algorithm that recovers a 3‐D resistivity model that considers arbitrary electrical anisotropy. The inversion uses the same inversion algorithm as the widely used ModEM inversion algorithm, and a novel forward modeling algorithm to consider the anisotropic Earth. The algorithm was tested on both synthetic and field MT data. Inversions considered both a completely general anisotropy tensor with six components and approximations with less parameters. Synthetic inversions show that the two horizontal components of resistivity and the anisotropy strike can be well recovered, while the vertical component of resistivity is poorly resolved, primarily because current flow in MTs is dominantly horizontal. The synthetic examples confirm the limitation of the axial anisotropic inversion technique when applied to MT data produced by a resistivity model with arbitrary anisotropy. The synthetic inversions also showed that inversion of data from an isotropic model will not result in an artificially anisotropic model. Compared to the isotropic inversion model of the real MT data, the anisotropic model clearly shows some features that are consistent with the mapped geology. As expected, the results showed that a given data set can be fit by a range of models, with an inherent trade‐off from 3‐D heterogeneity to 3‐D anisotropy. This uncertainty can be reduced with the use of prior information in the inversion. Plain Language Summary: The resistivity model of real Earth generally contains anisotropic features. Field magnetotelluric (MT) data may be significantly influenced by these features. To interpret such data, inversion algorithms which can consider electrical anisotropy are highly demanded. However, a limited number of practical MT inversion algorithms that can consider anisotropy have been published to date. To address this problem, we have developed a three‐dimensional (3‐D) MT inversion algorithm that recovers a 3‐D resistivity model that considers arbitrary electrical anisotropy. The algorithm was tested to be effective on both synthetic and field MT data. Compared to the isotropic inversion models of the tested MT data, the resolved anisotropic models are more reliable. Key Points: A three‐dimensional non‐linear conjugate gradient inversion algorithm is implemented for the inversion of magnetotelluric (MT) data to produce an electrical resistivity model with arbitrary anisotropyInversion examples of synthetic and real MT data validate the usefulness of the algorithmReal MT data inversion obtains an anisotropic model which is more consistent with the mapped geology than the corresponding isotropic model [ABSTRACT FROM AUTHOR]

Magnetotelluric (MT) and wide‐angle seismic reflection/refraction surveys play a fundamental role in understanding the crustal rheology and lithospheric structure of the Earth. In recent years, the integration of the two methods in order to improve the robustness of the inversion has started to gain attention. We present a new approach for joint 3‐D inversion of MT and wide‐angle seismic reflection/refraction data to accurately determine crustal structures and Moho depth. Based on H‐κ stacking of teleseismic receiver functions, we estimate an initial reference Moho. This is used as input for the subsequent MT/seismic joint inversion, where the Moho interface is updated and crustal structures are added to the model. During the joint inversion process, structural similarity is facilitated through the cross‐gradient constraint. Synthetic model tests show an improvement of the inversion results over separate inversions. In particular, the tests based on two geologically realistic models demonstrate that the crustal structure and even the trade‐off between velocity and Moho interface can be sufficiently resolved by combined MT and seismic data sets when using the estimates from analysis of receiver functions. These results show that the new method can provide useful constraints on crustal structures including their geophysical properties and discontinuities. Key Points: A new coupled algorithm is developed for 3‐D imaging of magnetotelluric and wide‐angle seismic reflection/refraction dataTwo strategies are presented to mitigate the problem of accurately determining crustal structures and Moho interface simultaneouslyThe inclusion of teleseismic receiver functions is proved to be effective in the new algorithm [ABSTRACT FROM AUTHOR]

Controlled-source audio-frequency magnetotellurics (CSAMT) has developed rapidly in recent years and are widely used in the area of mineral and oil resource exploration as well as other fields. The current theory, numerical simulation, and inversion research are based on the assumption that the underground media have resistivity isotropy. However a large number of rock and mineral physical property tests show the resistivity of underground media is generally anisotropic. With the increasing application of CSAMT, the demand for probe accuracy of practical exploration to complex targets continues to increase. The question of how to evaluate the influence of anisotropic resistivity to CSAMT response is becoming important. To meet the demand for CSAMT response research of resistivity anisotropic media, this paper examines the CSAMT electric equations, derives and realizes a three-dimensional (3D) staggered-grid finite difference numerical simulation method of CSAMT resistivity axial anisotropy. Through building a two-dimensional (2D) resistivity anisotropy geoelectric model, we validate the 3D computation result by comparing it to the result of controlled-source electromagnetic method (CSEM) resistivity anisotropy 2D finite element program. Through simulating a 3D resistivity axial anisotropy geoelectric model, we compare and analyze the responses of equatorial configuration, axial configuration, two oblique sources and tensor source. The research shows that the tensor source is suitable for CSAMT to recognize the anisotropic effect of underground structure. [ABSTRACT FROM AUTHOR]

The eastern terminations of the Himalayan orogeny, named Namche Barwa, are considered a vital natural laboratory in the Tibetan plateau for geodynamics due to its distinctive geological and geomorphological characteristics. Magnetotelluric (MT) data measured at 83 sites around the Namche Barwa are imaged by three-dimensional (3-D) inversion to better reveal the crustal structure of the eastern Himalaya. The results show a complex and heterogeneous electrical structure beneath the Namche Barwa. The electrical conductors distributed in the middle and lower crust around the Namche Barwa provide additional evidence for the 'crustal flow' model if they are considered as some parts of the flow in a relatively large-scale region. The near-surface resistivity model beneath the inner part of Namche Barwa conforms with the locations of hot spring and fluid inclusions, the brittle-ductile transition, and the 300°C-400°C isotherm from previous hydrothermal studies. Relatively resistive upper crust (>800 Ωm) is underlain by a more conductive middle to lower crust (<80 Ωm). The electrical characteristics of the thermal structure at shallow depth indicate an accumulation of hydrous melting, a localized conductive steep dipping zone for decompression melting consistent with the 'tectonic aneurysm' model for explaining the exhumation mechanism of metamorphic rocks at Namche Barwa. The results also imply that both surface processes and local tectonic responses play a vital role in the evolution of Namche Barwa. An alternative hypothesis that the primary sustained heat source accounts for the local thermal-rheological structure beneath Namche Barwa is also discussed. [ABSTRACT FROM AUTHOR]

The Luobusa chromite mine in Tibet, the largest chromite deposit discovered in China to date, is hosted in the Luobusa ophiolite in the eastern section of the Yarlung Zangbo suture zone. A combined seismic reflection and magnetotelluric survey was carried out on the mantle peridotites in which the chromitite deposits occur in order to investigate their origin and to prospect for other deposits in the area. The survey indicates that the Luobusa chromitite-bearing ultramafic block is a tectonic slice thrust onto the southern margin of the Gangdese batholith by multiple south-dipping faults. The maximum thickness of the ultramafic rock mass is 4 km, and it is underlain by the upper crust of the Indian Plate. Low-resistivity bodies in the upper crust of the Luobusa area lie at a depth of up to 5 km, and strata at a depth of 7 km cross the Yarlung Zangbo suture zone at a relatively small angle. Between depths of 20–50 km, there are a number of oval low-resistivity bodies, the lower parts of which might have been close to the Moho at a depth of approximately 65 km near the Yarlung Zangbo River. It is inferred that such low-resistivity bodies arise from partially molten bodies in the middle and lower crust and are related to deep tectonic activity and mantle material upwelling in the Yarlung Zangbo suture zone. [ABSTRACT FROM AUTHOR]