Spatial structural characteristics of the Deda ancient landslide in the eastern Tibetan Plateau:Insights from Audio-frequency Magnetotellurics and the Microtremor Survey Method

It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the landslides and indicate the potential for future reactivation.This study examines the Deda ancient landslide,situated in the Chalong-ranbu fault zone,where creep deformation suggests a complex underground structure.By integrating remote sensing,field surveys,Audio-frequency Magnetotellurics(AMT),and Microtremor Survey Method(MSM)techniques,along with engineering geological drilling for validation,to uncover the landslide’s spatial feature s.The research indicates that a fault is developed in the upper part of the Deda ancient landslide,and the gully divides it into Deda landslide accumulation zoneⅠand Deda landslide accumulation zoneⅡin space.The distinctive geological characteristics detectable by MSM in the shallow subsurface and by AMT in deeper layers.The findings include the identification of two sliding zones in the Deda I landslide,the shallow sliding zone(DD-I-S1)depth is approximately 20 m,and the deep sliding zone(DD-I-S2)depth is 36.2-49.9 m.The sliding zone(DD-Ⅱ-S1)depth of the DedaⅡlandslide is 37.6-43.1 m.A novel MSM-based method for sliding zone identification is proposed,achieving less than 5%discrepancy in depth determination when compared with drilling data.These results provide a valuable reference for the spatial structural analysis of large-deepseated landslides in geologically complex regions like the eastern Tibetan Plateau.

3D Step-by-step inversion strategy for audio magnetotellurics data based on unstructured mesh

A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the inversion region is discretized with unstructured tetrahedral elements.The inversion proceeds step-by-step from a coarse mesh to a fine mesh.In the inversion iteration process,a mesh is adaptively optimized according to the spatial gradient information about the model resistivity to fine delineate the boundaries of abnormal bodies.In the early stage of inversion execution,a coarse mesh is used for inversion,and the inversion stability is improved by reducing the number of inversion elements.In addition,mesh refinement is performed in the iterative inversion process.The inversion results obtained from the previous mesh are used as the reference and initial models for the next mesh iterative inversion.The step-by-step inversion strategy can ensure that the inversion is performed in the correct direction,improving the inversion stability and results gradually.Synthetic results show that the step-by-step inversion strategy with a Gauss-Newton method for 3D AMT inversion is stable and reliable,which lays a foundation for further practical 3D AMT data inversion.

Proterozoic tectonothermal processes imaged with magnetotellurics and seismic reflection in southern Australia

Over the last two decades,co-located seismic and magnetotelluric(MT) profiles provided fundamental geophysical data sets to image the Australian crust.Despite their complimentary nature,the data are processed and often interpreted separately without common processes in mind.We here qualitatively compare 2 D resistivity inversion models derived from MT and seismic reflection profiles across a region of Archean-Proterozoic Australia to address the causes of variations in seismic response and anomalous conductivity in the crust.We find that there exists a spatial association between regions of low reflectivity in seismic sections and low resistivity in co-located2 D MT modelled sections.These relationships elucidate possible signatures of past magmatic and fluid-related events.Depending on their diffuse or discrete character,we hypothesize these signatures signify fossil melting of the crust due to mafic underplating,magma movement or hydrothermal fluid flow through the crust.The approach discussed herein is a process-oriented approach to interpretation of geophysical images and a significant extension to traditional geophysical methods which are primarily sensitive to a singular bulk rock property or state.

3D model of Svecofennian Accretionary Orogen and Karelia Craton based on geology,reflection seismics,magnetotellurics and density modelling:Geodynamic speculations

A 3D model of deep crustal structure of the Archaean Karelia Craton and late Palaeoproterozoic Svecofennian Accretionary Orogen including the boundary zone is presented.The model is based on the combination of data from geological mapping and reflection seismic studies,along profiles 1-EU,4 B,FIRE-1-2 a-2 and FIRE-3-3 a,and uses results of magnetotelluric soundings in southern Finland and northern Karelia.A seismogeological model of the crust and crust-mantle boundary is compared with a model of subhorizontal velocity-density layering of the crust.The TTG-type crust of the Palaeoarchaean and Mesoarchaean microcontinents within the Karelia Craton and the Belomorian Province are separated by gently dipping greenstone belts,at least some of which are palaeosutures.The structure of the crust was determined mainly by Palaeoproterozoic tectonism in the intracontinental settings modified by a strong collisional compression at the end of the Palaeoproterozoic.New insights into structure,origin and evolution of the Svecofennian Orogen are provided.The accretionary complex is characterized by inclined tectonic layering:the tectonic sheets,~15 km thick,are composed of volcanic-sedimentary rocks,including electro-conductive graphite-bearing sedimentary rocks,and electro-resistive granitoids,which plunge monotonously and consecutively eastward.Upon reaching the level of the lower crust,the tectonic sheets of the accretionary complex lose their distinct outlines.In the seismic reflection pattern they are replaced by a uniform acoustically translucent medium,where separate sheets can only be traced fragmentarily.The crust-mantle boundary bears a diffuse character:the transition from crust to mantle is recorded by the disappearance of the vaguely drawn boundaries of the tectonic sheets and in the gradual transition of acoustically homogeneous and translucent lower crust into transparent mantle.Under the effect of endogenic heat flow,the accretionary complex underwent high-temperature metamorphism and partial melting.Blurring of the rock contacts,which in the initial state created contrasts of acoustic impedance,was caused by partial melting and mixing of melts.The 3 D model is used as a starting point for the evolutionary model of the Svecofennian Accretionary Orogen and for determination of its place in the history of the Palaeoproterozoic Lauro-Russian intracontinental orogeny,which encompassed a predominant part of the territory of Lauroscandia,a palaeocontinent combining North American and East European cratons.The model includes three stages in the evolution of the Lauro-Russian Orogen(-2.5,2.2-2.1 and 1.95-1.87 Ga).The main feature of the Palaeoproterozoic evolution of the accretionary Svecofennian Orogen and Lauroscandia as a whole lay in the causal link with evolution of a superplume,which initiated plate-tectonic events.The Svecofennian-Pre-Labradorian palaeo-ocean originated in the superplume axial zone;the accretionary orogens were formed along both continental margins due to closure of the palaeo-ocean.

Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes

Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.

Imaging the Architecture of Mineral Systems and the Pathways of Ore-forming Fluids across Mongolia with Magnetotellurics

In the framework of a mineral system approach,a combination of components is required to develop a mineral system.This includes the whole-lithosphere architecture,which controls the transport of ore-forming fluids,and favorable tectonic and geodynamic processes,occurring at various spatial and temporal scales,that influence the genesis and evolution of ore-forming fluids(Huston et al.,2016;Groves et al.,2018;Davies et al.,2020).Knowledge of the deep structural framework can advance the understanding of the development of a mineral system and the emplacement of mineral deposits.Deep geophysical exploration carried out with this aim is increasingly important for targeting new ore deposits in unexplored and underexplored regions(Dentith et al.,2018;Dentith,2019).

Application of integrated geophysical techniques in geothermal exploration in Binhai County, Jiangsu Province

Integrated geophysical technology is a necessary and effective means for geothermal exploration.However,integration of geophysical technology for large‐scale surveys with those for geothermal reservoir localization is still in development.This study used the controlled source audio‐frequency magnetotelluric method technology for large‐scale exploration to obtain underground electrical structure information and micromotion detection technology to obtain underground wave velocity structure information.The combination of two detection technologies was used for local identification of geothermal reservoirs.Further,auxiliary correction and inversion constraint were implemented through the audio magnetotelluric sounding technology for maximum authenticity restoration of the near‐and transition‐field data.Through these technology improvements,a geothermal geological model was established for the Binhai County of Jiangsu Province in China and potential geothermal well locations were identified.On this basis,a geothermal well was drilled nearly 3000m deep,with a daily water volume of over 2000m3/day and a geothermal water temperature of 51°C at the well head.It is found that predictions using the above integrated geophysical exploration technology are in good agreement with the well geological formation data.This integrated geophysical technology can be effectively applied for geothermal exploration with high precision and reliability.

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

Despite 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.

Three-dimensional magnetotelluric modeling in the spherical and Cartesian coordinate systems:A comparative study

With the increase in the coverage area of magnetotelluric data,three-dimensional magnetotelluric modeling in spherical coordinates and its differences with respect to traditional Cartesian modeling have gradually attracted attention.To fully understand the influence of the Earth’s curvature and map projection deformations on Cartesian modeling,qualitative and quantitative analyses based on realistic three-dimensional models need to be examined.Combined with five representative map projections,a type of model conversion method that transforms the original spherical electrical conductivity model to Cartesian coordinates is described in this study.The apparent resistivity differences between the spherical western United States electrical conductivity model and the corresponding five Cartesian models are then compared.The results show that the cylindrical equal distance map projection has the smallest error.A meridian convergence correction resulting from the deformation of the map projection is introduced to rotate the Cartesian impedance tensor from grid north to geographic north,which reduces differences from the spherical results.On the basis of the magnetotelluric field data,the applicability of the Cartesian coordinate system to western and contiguous United States models is quantitatively evaluated.Precise interpretations of the contiguous United States model were found to require spherical coordinates.

MT探测居里面深度的可行性研究

居里面是地球内部铁磁性物质向顺磁性物质转换的界面,在温度略低于居里点时物质磁化率会快速升高,这被称为Hopkinsin峰.对于地球内部而言Hopkinsin峰是只有几百米至几公里厚的薄层,由于其与居里温度的关系,因此其底界面的深度可以作为居里面深度的估计.传统的居里面深度探测方法包括谱分析方法、等层模型方法和温度-深度剖面法.这些方法是人们研究地球内部热结构和居里面深度的重要手段,但是谱分析和等层模型的结果均有一些固有的缺点,如横向分辨率太低等;而地热方程的结果则受地表因素影响十分严重,并且地球内部的热源分布也不是十分清楚,这导致了其结果是不可靠的.本文提出用MT方法探测居里面深度,通过对几种简单一维模型进行的正反演数值试验,论证了该方法的可行性.结果表明,用MT方法研究居里面性质,不但可以得到居里面深度,还可以得到居里面顶部Hopkinsin峰所对应介质的电学和磁学性质,但必须同时反演岩层的电阻率和磁导率,才能获得较为可靠的居里面深度估计.

Parallel rapid relaxation inversion of 3D magnetotelluric data

We implement a parallel algorithm with the advantage of MPI （Message Passing Interface） to speed up the rapid relaxation inversion for 3D magnetotelluric data. We test the parallel rapid relaxation algorithm with synthetic and real data. The execution efficiency of the algorithm for several different situations is also compared. The results indicate that the parallel rapid relaxation algorithm for 3D magnetotelluric inversion is effective. This parallel algorithm implemented on a common PC promotes the practical application of 3D magnetotelluric inversion and can be suitable for the other geophysical 3D modeling and inversion.

Three-dimensional conjugate gradient inversion of magnetotelluric sounding data

Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.

Nonlinear inversion for magnetotelluric sounding based on deep belief network

To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network inputs are the apparent resistivities of known models,and the outputs are the model parameters.The optimal network structure is achieved by determining the numbers of hidden layers and network nodes.Secondly,the learning process of the DBN is implemented to obtain the optimal solution of network connection weights for known geoelectric models.Finally,the trained DBN is verified through inversion tests,in which the network inputs are the apparent resistivities of unknown models,and the outputs are the corresponding model parameters.The experiment results show that the DBN can make full use of the global searching capability of the restricted Boltzmann machine(RBM)unsupervised learning and the local optimization of the back propagation(BP)neural network supervised learning.Comparing to the traditional neural network inversion,the calculation accuracy and stability of the DBN for MT data inversion are improved significantly.And the tests on synthetic data reveal that this method can be applied to MT data inversion and achieve good results compared with the least-square regularization inversion.

Experimental Verification and Research for the Distortion in the Integrated Frequency Responses of the High-Pressure Sealed Cabin and Magnetic Field Sensor

Although magnetotelluric sounding method applied to the land is advanced, there are many difficulties when it is applied to marine environment, one of which is how to lay magnetic field sensors down to the seafloor to complete measurements. To protect the magnetic field sensors from intense erosion and high pressure, suitable high-pressure sealed cabins must be designed to load them. For the consideration of magnetic measurement and marine operation, the sealed pressure cabin should be nonmagnetic and transportable. Among all optional materials, LC4 super.hard aluminum alloy has the highest performance of price/quality ratio to make the sealed pressure cabin. However, it does not mean that the high-pressure sealed cabin made using LC4 will be perfect in performance. In fact, because of its weak magnetism, the pressure cabin made using LC4 has distorting effect on frequency responses of the magnetic field sensors sealed in it. This distorting effect does not affect the use of the magnetic field sensor, but if we want to eliminate its effect, we should study it by experimental measurements. In our experiment tests, frequency sweep magnetic field as excitation signal was used, and then responses of the magnetic field sensor before and after being loaded into the high-pressure sealed cabin were measured. Finally, normalized abnormal curves for the frequency responses were obtained, through which we could show how the high-pressure sealed cabin produces effects on the responses of the magnetic field sensor. Experimental results suggest that the response distortion induced by the sealed pressure cabin appears on mid- and high-frequency areas. Using experimental results as standardization data, the frequency responses collected from seafloor magnetotelluric measurements can be corrected to restore real information about the seafloor field source.

Geophysical signatures of fluids in a reactivated Precambrian collisional suture in central India

The Central India Tectonic Zone （CITZ） marks the trace of a major suture zone along which the south Indian and the north Indian continental blocks were assembled through subduction-accretion- collision tectonics in the Mesoproterozoic. The CITZ also witnessed the major, plume-related, late Cretaceous Deccan volcanic activity, covering substantial parts of the region with continental flood basalts and associated magmatic provinces. A number of major fault zones dissect the region, some of which are seismically active. Here we present results from gravity modeling along five regional profiles in the CITZ, and combine these results with magnetotelluric （MT） modeling results to explain the crustal architecture. The models show a resistive （more than 2000 Ω. m） and a normal density （2.70 g/cm3） upper crust suggesting dominant tonalite-trondhjemite-granodiorite （TTG） composition. There is a marked correlation between both high-density （2.95 g/cm3） and low-density （2.65 g/cm3） regions with high conductive zones （〈80 Ω. m） in the deep crust. We infer the presence of an interconnected grain boundary network of fluids or fluid-hosted structures, where the conductors are associated with gravity lows. Based on the conductive nature, we propose that the lower crustal rocks are fluid reservoirs, where the fluids occur as trapped phase within minerals, fluid-filled porosity, or as fluid-rich structural conduits. We envisage that substantial volume of fluids were transferred from mantle into the lower crust through the younger plume-related Deccan volcanism, as well as the reactivation, fracturing and expulsion of fluids transported to depth during the Mesoproterozoic subduction tectonics. Migration of the fluids into brittle fault zones such as the Narmada North Fault and the Narmada South Fault resulted in generating high pore pressures and weakening of the faults, as reflected in the seismicity. This inference is also supported by the presence of broad gravity lows near these faults, as well as the low velocity in the lower crust beneath regions of recent major earthquakes within the CITZ.

High resolution tomography of the Tanlu fault zone near Hefei with passive seismic and magnetotelluric linear array data

As the largest fault trending NNE-SSW to NE-SW in the eastern Eurasia Continent,the Tanlu fault zone(TLFZ)extends over 2,400 km within China,roughly from Wuxue,Hubei Province,to Russia.Since the Quaternary period,the TLFZ has been an earthquake-prone area in eastern China where several major earthquakes resulted by tectonic compression occurred,causing tremendous casualties and significant economic losses.Many studies on different segments of the TLFZ have been carried out in the past few decades.However,numerous key questions regarding the fault zone remain unanswered due to a lack of clear subsurface characterization and fault delineation.In this study,we present high-resolution tomographic results across the TLFZ to the east of Hefei,where one 22-km-long passive seismic array with densely spaced short-period nodes,and a 24-km-long magnetotelluric array were deployed adjacent to each other.We find the velocity and resistivity variations are highly consistent with known surface geology.Sharp property contrasts in both the seismic shear wave velocity and electrical resistivity profiles clearly delineate the Tanlu F1 fault(TLF-1)near Hefei.More interestingly,an upwelling with distinct high velocity is imaged within the Hefei Basin to the west of the TLF-1,whereas a slanted block with lowvelocity and low-resistivity seems to cut into or thrust upon the high-grade to low-grade middle-pressure rocks in the Zhangbaling uplift right below the Tanlu F2 ductile shear fault(TLF-2).The presented results show a new approach to characterize deep subsurface structure of the TLFZ beyond 2-km depths using passive data,which it is often difficult for active seismic surveys with refracted and reflected waves to image.

Analysis of prospecting polymetallic metallogenic belts by comprehensive geophysical method

This paper is based on the analysis and research on the silver-lead-zinc polymetallic ore in New Ballyhoo Banner in southern Manzhouli of Inner Mongolia.Because metal mineralization brings rock formations,the geophysical features such as low resistivity,high polarization rate and uneven distribution of magnetization,the comprehensive geophysical methods are adopted including high-precision magnetic measurement,high-power induced polarization,IP field middle gradient and controlled source audio-frequency magnetotellurics.In the survey work of multi-metal ore deposits,from surface sweeping to single point measurement,and from single point to section going deeper layer by layer,the resolution of measurement is continuously improved,and various geophysical methods support and complement each other,so explorers can successfully predict the direction,scale and volume of the metallogenic belts in conjunction with geochemical exploration,geological survey and drilling.It has provided a strong basis for completing the exploration task of predicting the reserve volume of ore bodies.The research conclusions of this exploration case have thus a high reference value in the same type of exploration work.

An Asthenospheric Upwelling Beneath Central Mongolia——Implications for Intraplate Surface Uplift and Volcanism

Intraplate processes,such as continental surface uplift and intraplate volcanism,are enigmatic and the underlying mechanisms responsible are not fully understood.Central Mongolia is an ideal natural laboratory for studying such processes because of its location in the continental interior far from tectonic plate boundaries.

Crustal Electrical Structure and Deep Metallogenic Potential in Northern Wuyi Area(South China),based on Magnetotelluric Data

The northern Wuyi area,which is located in the northern Wuyi metallogenic belt,has superior mineralization conditions.The Pingxiang-Guangfeng-Jiangshan-Shaoxing fault(PSF)extends across the whole region regardless of whether or how the PSF relates to the near-surface mineralization.We carried out an MT survey in the region and obtained a reliable 2D model of the crustal electrical structure to a depth of 30 km.In the resistivity model,we inferred that a continuous high conductivity belt that ranges from the shallow to deep crust is a part of the PSF.Then,we estimated the fluid content and pressure gradient to identify the deep sources of fluid as well as its pattern of motion pattern.Finally,we proposed a model for the deep metallogenic migration processes that combines geological data,fluid content data,pressure gradient data,and the subsurface resistivity model.The model analysis showed that the Jiangnan orogenic belt and the Cathaysia block formed the PSF during the process of com.The deep fluid migrated upward through the PSF to the shallow crust.Therefore,we believe that the PSF is an ore-forming fluid migration channel and that it laid the material basis for large-scale mineralization in the shallow crust.

Three-dimensional magnetotelluric regularized inversion based on smoothness-constrained model

How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.