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 lan...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.展开更多
Different geophysical exploration methods have significant differences in terms of exploration depth,especially in frequency domain electromagnetic(EM)exploration.According to the definition of skin depth,this differe...Different geophysical exploration methods have significant differences in terms of exploration depth,especially in frequency domain electromagnetic(EM)exploration.According to the definition of skin depth,this difference will increase with the effective detection frequency of the method.As a result,when performing three-dimensional inversion on single type of EM data,it is not possible to effectively distinguish the subsurface geoelectric structure at the full scale.Therefore,it is necessary to perform joint inversion on different type of EM data.In this paper we combine the magnetotelluric method(MT)with the controlled-source audio-magnetotelluric method(CSAMT)to study the frequency-domain three-dimensional(3D)joint inversions,and we use the unstructured finite-element method to do the forward modeling for them,so that the numerical simulation accuracies of different electromagnetic methods can be satisfied.By combining the two sets of data,we can obtain the sensitivity of the electrical structure at different depths,and depict the full-scale subsurface geoelectric structures.In actual mineral exploration,the 3D joint inversion is more useful for identifying subsurface veins in the shallow part and blind mines in the deep part.It can delineate the morphological distribution of ore bodies more completely and provide reliable EM interpretations to guide the mining of minerals.展开更多
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 ofte...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.展开更多
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 i...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.展开更多
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...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.展开更多
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 fr...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.展开更多
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 attr...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.展开更多
基金supported by the National Natural Science Foundation of China(42372339)the China Geological Survey Project(DD20221816,DD20190319)。
文摘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.
基金Supported by the National Natural Science Foundation of China(No.42074120).
文摘Different geophysical exploration methods have significant differences in terms of exploration depth,especially in frequency domain electromagnetic(EM)exploration.According to the definition of skin depth,this difference will increase with the effective detection frequency of the method.As a result,when performing three-dimensional inversion on single type of EM data,it is not possible to effectively distinguish the subsurface geoelectric structure at the full scale.Therefore,it is necessary to perform joint inversion on different type of EM data.In this paper we combine the magnetotelluric method(MT)with the controlled-source audio-magnetotelluric method(CSAMT)to study the frequency-domain three-dimensional(3D)joint inversions,and we use the unstructured finite-element method to do the forward modeling for them,so that the numerical simulation accuracies of different electromagnetic methods can be satisfied.By combining the two sets of data,we can obtain the sensitivity of the electrical structure at different depths,and depict the full-scale subsurface geoelectric structures.In actual mineral exploration,the 3D joint inversion is more useful for identifying subsurface veins in the shallow part and blind mines in the deep part.It can delineate the morphological distribution of ore bodies more completely and provide reliable EM interpretations to guide the mining of minerals.
文摘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.
文摘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.
基金funded by the Anhui Province Science Program(2018-g-1-4)the National Science Program(Nos.41630320,42174087)the National Key R&D Program of China(No.2016YFC0600201)。
文摘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.
基金the RF state project(No.0135-2019-0036)partly supported by the project(No.0135-20180039)of the Program of the Presidium of the RAS“Fundamental problems of geological and geophysical study of lithospheric processes”the Russian Foundation for Basic Research for supporting project(No.19-05-00336)。
文摘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.
基金the National Natural Science Foundation of China(Nos.42220104002,42104073,and 41630317).
文摘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.
基金supported by the Construction S&T Project of Department of Transportation of Sichuan Province(Grant No.2023A02,No.2024A04,No.2020A01)the Sichuan Science and Technology Program(Grant No.2022YFG0141)+3 种基金the Research Project of Sichuan Highway Planning,Survey,Design,and Research Institute Ltd.(Grant No.KYXM2021000049,No.KYXM2022000038,No.KYXM202300056)the National Natural Science Foundation of China(41630640)the National Science Foundation of Innovation Research Group(41521002)the National Natural Science Foundation of China(41790445).
基金supported by the Construction S&T Project of Department of Transportation of Sichuan Province (Grant No.2023A02,No.2024A04,No.2020A01)the Sichuan Science and Technology Program (Grant No.2022YFG0141)+3 种基金the Research Project of Sichuan Highway Planning,Survey,Design,and Research Institute Ltd. (Grant No.KYXM2021000049,KYXM2022000038,No.KYXM2023000056)the National Natural Science Foundation of China (41630640)the National Natural Science Foundation of China (41790445)the National Science Foundation of Innovation Research Group (41521002)。