A VGGNet-based correction for satellite altimetry-derived gravity anomalies to improve the accuracy of bathymetry to depths of 6500 m

Understanding the topographic patterns of the seafloor is a very important part of understanding our planet.Although the science involved in bathymetric surveying has advanced much over the decades,less than 20%of the seafloor has been precisely modeled to date,and there is an urgent need to improve the accuracy and reduce the uncertainty of underwater survey data.In this study,we introduce a pretrained visual geometry group network(VGGNet)method based on deep learning.To apply this method,we input gravity anomaly data derived from ship measurements and satellite altimetry into the model and correct the latter,which has a larger spatial coverage,based on the former,which is considered the true value and is more accurate.After obtaining the corrected high-precision gravity model,it is inverted to the corresponding bathymetric model by applying the gravity-depth correlation.We choose four data pairs collected from different environments,i.e.,the Southern Ocean,Pacific Ocean,Atlantic Ocean and Caribbean Sea,to evaluate the topographic correction results of the model.The experiments show that the coefficient of determination(R~2)reaches 0.834 among the results of the four experimental groups,signifying a high correlation.The standard deviation and normalized root mean square error are also evaluated,and the accuracy of their performance improved by up to 24.2%compared with similar research done in recent years.The evaluation of the R^(2) values at different water depths shows that our model can achieve performance results above 0.90 at certain water depths and can also significantly improve results from mid-water depths when compared to previous research.Finally,the bathymetry corrected by our model is able to show an accuracy improvement level of more than 21%within 1%of the total water depths,which is sufficient to prove that the VGGNet-based method has the ability to perform a gravity-bathymetry correction and achieve outstanding results.

Characteristics of gravity anomalies and tectonic analysis of Enderby Land in East Antarctica and its adjacent areas

Enderby Land in East Antarctica and its adjacent areas,which are closely related to the Indian Plate in their geological evolution,have become one of the key zones for studies on how the Antarctic continent evolves.Based on the isostasy and flexure theories of the lithosphere and using the CRUST1.0 model as the depth constraint,this paper uses the gravity field model EIGEN-6C4 and topographic data to calculate the isostatic gravity anomalies of Enderby Land and its adjacent areas.Then,the crustal thickness of the study area is calculated,and three comprehensive geophysical interpretation profiles that vertically span the study area are plotted.The results show that the flexural isostatic gravity anomalies in Enderby Land and its adjacent areas are closely related to the regional tectonic setting,and the anomalies in different regions differ substantially,ranging from−50×10^(−5)m/s^(2)to 85×10^(−5)m/s^(2).A zone of high isostatic gravity anomalies(30×10^(−5)−80×10^(−5)m/s^(2))is distributed outside the Cooperation Sea and Queen Maud Land,which may be plate remnants generated by early rifting.Except for the Kerguelen Plateau,which was formed by a hotspot and has a crustal thickness of 15 km,the thickness of the oceanic crust in other parts of the study area changes slightly by approximately 4–9 km,with the thinnest part being in Enderby Basin.The thickness of the inland crust along the coastline increases with the elevation,with the maximum thickness reaching 34 km.The isostatic gravity anomalies corresponding to the zone of high magnetic anomalies along the continental margin of Queen Maud Land are negative and small,with an isostatic adjustment trend indicating Moho surface uplift,and those on the edge of central Enderby Land are near zero,approaching the isostatic state,which may be caused by the magmatism at the early stage of rifting.The continental-oceanic boundary should be close to the contour line of the crustal thickness 10–12 km on the outer edge of the coastline.

Constraining the crustal structure under the central and western Tian Shan based on teleseismic receiver functions and gravity anomalies

The Tian Shan is a vast range that spans several countries in Asia.Understanding its evolutionary history may provide valuable insights into intracontinental orogenic dynamics.In this study,we explored the crustal characteristics of the Tian Shan and their relationships to the tectonic evolution of the region.A new H-stacking method that combines the P receiver function and gravity anomalies was used to estimate the thickness and ratio of P-to S-wave velocities(Vp/Vs)for 91 broadband seismic stations in the central and western Tian Shan.Our results revealed significant lateral variations in crustal thickness and Vp/Vs.A—45-km-thick crust and an intermediate-high Vp/Vs(-1.74-1.84)were found in the Kazakh Shield and Tarim Basin,which we interpreted to indicate a mafic crystalline basement and lower crust.The central Tian Shan varied greatly in crustal thickness(40-64 km)and Vp/Vs ratio(1.65-2.00).which may be due to crustal shortening,mafic underplating,and crustal melting.In contrast,we observed a relatively thin crust(42-50 km)with an intermediate Vp/Vs ratio(-1.78)in the western Tian Shan.The differences in the crustal structures between the western and central Tian Shan imply that the Talas-Fergana Fault may be trans-lithospheric.

Gravity anomalies determined from mean sea surface model data over the Gulf of Mexico

With the improvements in the density and quality of satellite altimetry data,a high-precision and high-resolution mean sea surface model containing abundant information regarding a marine gravity field can be calculated from long-time series multi-satellite altimeter data.Therefore,in this study,a method was proposed for determining marine gravity anomalies from a mean sea surface model.Taking the Gulf of Mexico(15°–32°N,80°–100°W)as the study area and using a removal-recovery method,the residual gridded deflections of the vertical(DOVs)are calculated by combining the mean sea surface,mean dynamic topography,and XGM2019e_2159 geoid,and then using the inverse Vening-Meinesz method to determine the residual marine gravity anomalies from the residual gridded DOVs.Finally,residual gravity anomalies are added to the XGM2019e_2159 gravity anomalies to derive marine gravity anomaly models.In this study,the marine gravity anomalies were estimated with mean sea surface models CNES_CLS15MSS,DTU21MSS,and SDUST2020MSS and the mean dynamic topography models CNES_CLS18MDT and DTU22MDT.The accuracy of the marine gravity anomalies derived by the mean sea surface model was assessed based on ship-borne gravity data.The results show that the difference between the gravity anomalies derived by DTU21MSS and CNES_CLS18MDT and those of the ship-borne gravity data is optimal.With an increase in the distance from the coast,the difference between the gravity anomalies derived by mean sea surface models and ship-borne gravity data gradually decreases.The accuracy of the difference between the gravity anomalies derived by mean sea surface models and those from ship-borne gravity data are optimal at a depth of 3–4 km.The accuracy of the gravity anomalies derived by the mean sea surface model is high.

Crustal Uplift in the Longmen Shan Mountains Revealed by Isostatic Gravity Anomalies along the Eastern Margin of the Tibetan Plateau

This study examines the relationship between high positive isostatic gravity anomalies （IGA）, steep topography and lower crustal extrusion at the eastern margin of the Tibetan Plateau. IGA data has revealed uplift and extrusion of lower crustal flow in the Longmen Shan Mountains （the LMS）. Firstly, The high positive IGA zone corresponds to the LMS orogenic belt. It is shown that abrupt changes in IGA correspond to zones of abrupt change of topography, crustal thickness and rock density along the LMS. Secondly, on the basis of the Airy isostasy theory, simulations and inversions of the positive IGA were conducted using three-dimensional bodies. The results indicated that the LMS lacks a mountain root, and that the top surface of the lower crust has been elevated by 11 km, leading to positive IGA, tectonic load and density load. Thirdly, according to Watts＇s flexural isostasy model, elastic deflection occurs, suggesting that the limited （i.e. narrow） tectonic and density load driven by lower crustal flow in the LMS have led to asymmetric flexural subsidence in the foreland basin and lifting of the forebulge. Finally, based on the correspondence between zones of extremely high positive IGA and the presence of the Precambrian Pengguan-Baoxing complexes in the LMS, the first appearance of erosion gravels from the complexes in the Dayi Conglomerate layer of the Chengdu Basin suggest that positive IGA and lower crustal flow in the LMS took place at 3.6 Ma or slightly earlier.

Calculation of geoid undulations and gravity anomaliesin the Northwest Pacific by using the Topex/Poseidonand Geosat altimeter data

The geoid undulations in the Northwest Pacific were calculated by usingtheTopex/Poseidon and Geosataltimeter data. Firstly the bias between two types of the altimeter data was removed and the geoid undulations in theNorthwest Pacific were acquired by a long wave bias diminishing model with a resolution of 30 km and precision of 14cm. Then an algorithm of inversion of gravity anomalies was derived , and the gravity anomalies in the East China Seawere calculated by using the algorithm and the geoid undulations. The rms of difference between the in situ measure-ments the gravity anomalies from altimeter data was 3 .8× 10-5 m/s2. A method to colculate the gravity anomaliesin a larger area was developed which combined gravity anomalies in four subregions overlapping each other into one data set in a larger region. The error analysys shaws that the model and result of the inversion of gravity anomalies were reliable.

Bathymetry predicting using the altimetry gravity anomalies in South China Sea

In South China Sea(112°E-119°E, 12°N-20°N), 81159 ship soundings published by NGDC(National Geophysics Data Center) and the altimetry gravity anomalies published by SIO(Scripps Institute of Oceanography) were used to predict bathymetry by GGM(gravity-geologic method) and SAS(Smith and Sandwell) method respectively. The residual 40576 ship soundings were used to estimate precisions of the predicted bathymetry models. Results showed that: the standard deviation of difference between the GGM model and ship soundings was 59.75 m and the relative accuracy was 1.86%. The SAS model is60.07 m and 1.87%. The power spectral densities of the ETOPO1, SIO, GGM and SAS models were also compared and analyzed. At last, we presented an integrated bathymetry model by weighted averaging method, the weighted factors were determined by precisions of the ETOPO1, SIO, GGM, and SAS model respectively.

Active Depths of Main Faults in the Ying-Qiong Basin Investigated by Multi-Scale Wavelet Decomposition of Bouguer Gravity Anomalies and Power Spectral Methods

The Ying-Qiong Basin is located on the northwestern margin of the South China Sea and at the junction of the South China Block and the Indochina Block.It is characterized by complex geological structures.The existing seismic data in the study area is sparse due to the lack of earthquake activities.Because of the limited source energy and poor coverage of seismic data,the knowledge of deep structures in the area,including the spatial distribution of deep faults,is incomplete.Contrarily,satellite gravity data cover the entire study area and can reveal the spatial distribution of faults.Based on the wavelet multi-scale decomposition method,the Bouguer gravity field in the Ying-Qiong Basin was decomposed and reconstructed to obtain the detailed images of the first-to sixth-order gravitational fields.By incorporating the known geological features,the gravitational field responses of the main faults in the Ying-Qiong Basin were identified in the detailed fields,and the power spectrum analysis yielded the depths of 1.4,8,15,26.5,and 39 km for the average burial depths of the bottom surfaces from the first-to fifth-order detailed fields,respectively.The four main faults in the Yinggehai Basin all have a large active depth range:fault A(No.1)is between 5 and 39 km,fault B is between 26.5 and 39 km,and faults C and D are between 15 and 39 km.However,the depth of active faults in the Qiongdongnan Basin is relatively shallow,mainly between 8 and 26.5 km.

Pre-seismic gravity anomalies before Linkou Ms6.4 earthquake by continuous gravity observation of Crustal Movement Observation Network of China

A Ms6.4 earthquake occurred at Linkou country, Heilongjiang Province （44.8°N, 129.9°E） on January 2, 2016 at a depth of 580 km. Pre-seismic graviW anomalies obtained at a 1 Hz sampling rate from Crustal Movement Observation Network of China （CMONOC） are analyzed after the earthquake. The results show that： （1） different from previous studies, both pre-seismic amplitude perturbation and co-seismic amplitude perturbation are not critical inversely proportional to epicentral distance; （2） unlike shallow earthquake, the pre-seismic and co-seismic amplitude perturbation of gravity illustrate syn- chronous spatial variation characters with decrease of epicentral distance for Linkou earthquake. This may because Linkou earthquake is a deep earthquake and occurred in Pacific Plate subduction zone; （3） compared to basement and semi-basement, cave can provide a better observation environment for gPhone gravimeter to detect pre-seismic gravity anomalies.

Spatial Analysis of Gravity Data in the Basement of the Yaoundé-Yoko Area from the Global Gravity Model: Implication on the Sanaga Fault (South-Cameroon)

In this work, gravity anomalies from the XGM2016 global gravity model are used to study the basement of the Yaounde, Yoko area. The aim is to locate the characteristic tectonic faults and to characterize the geometry of the basement of these localities in order to improve the knowledge of the structural and tectonic basement of the study area. Numerical filters (vertical gradient, horizontal gradient, upward continuation) and Euler deconvolution were applied to the gravity anomalies respectively for qualitative and quantitative analysis. The results of the qualitative analysis allowed us to establish the lineament map of the study area;ranging from 0 to 35 km depth. For the quantitative analysis, the work is done in two parts: 1) highlighting the distribution of depths of geological structures in the basement of the study area;2) 2D1/2 modeling of geological structures to highlight the geometry of the basement of Yaounde, Yoko area. Thus, from five suitably selected profiles, the established models reveal the presence of eight blocks of geological structures of different densities and analyze their implications on the Sanaga Fault. Moreover, the models show that the positive anomalies characteristics for the Sanaga Fault reflect the anomalous character due to the strong dominance of the shale intrusion in the basement.

Implications on Gravity Anomaly Measurements Associated with Different Lithologies in Turkana South Subcounty

The use of gravity data has demonstrated capability for monitoring lithological changes on a large scale as a consequence of differentiating basement and sedimentary of buried valleys. Gravity anomalies are associated with lateral contrasts in density and therefore deformation by faulting or folding will be manifested if accompanied by lateral density changes, otherwise, the vice versa is true. The study’s objective is to evaluate the effectiveness of gravity method in establishing different lithologies in an area. The study has revealed that regional anomaly gravity map presents high anomalies in the Northern region in the NW-SE trend and low anomalies in the southern trend in NW-SE, while the residual anomaly gravity map shows different trends for the low and high gravity anomalies. The gravity anomalies are well interpreted in line with the lithologies of the study area rather than the deformation of the same lithologies. There are observed high values of gravity anomaly values (ranging from -880.2 to -501.2 g.u.) where there are eolian unconsolidated rocks overlying the basement compared to low gravity anomaly values (ranging from -1338.9 to -1088.7 g.u.) where the andesites, trachytes and phonolites overly the basement. The different regional gravity anomalies relate well with different rock densities in the study area along the line profile for radially averaged power spectrum. The gravity highs are noted in the eastern point and are associated with andesites, trachytes, basalts and igneous rocks, while the gravity lows are associated with sandstone, greywacke, arkose, and eolian unconsolidated rock. The utilization of the information from the Power spectrum analysis demonstrates that the depth to the deepest basement rock is 12.8 km which is in the eastern flank, while the shallowest to the basement of 1.1 km to the western flank.

Characteristic of Gravity and Magnetic Anomalies in the Daba Shan and the Sichuan basin, China: Implication for Architecture of the Daba Shan

The Dabashan nappe structural belt links the Hannan block to the west with the Huangling block to the east between Yangxian and Xiangfan. The Dabashan arc-shaped fold belt formed during late Jurassic and was superposed on earlier Triassic folds. To achieve an improved understanding of the deep tectonics of the Dabashan nappe structural belt, we processed and interpreted the gravity and magnetic data for this area using new deep reflection seismic and other geophysical data as constraints. The results show that the Sichuan basin and Daba Mountains lie between the Longmenshan and Wulingshan gravity gradient belts. The positive magnetic anomalies around Nanchong-Tongjiang-Wanyuan-Langao and around Shizhu result from the crystalline basement. Modeling of the gravity and magnetic anomalies in the Daba Mountains and the Sichuan basin shows that the crystalline basement around Nanchong-Tongjiang-Wanyuan-Langao extends to the northeast underneath the Wafangdian fault near Ziyang. The magnetic field boundary in the Zhenba-Wanyuan-Chengkou-Zhenping area is the major boundary of the Dabashan nappe thrusting above the Sichuan Basin. This boundary might be the demarcation between the south Dabashan and the north Dabashan structural elements. The low gravity anomaly between Tongjiang and Chengkou might be partly caused by thickened lower crust. The local low gravity anomaly to the south of Chengkou-Wanyuan might result from Mesozoic strata of low density in the Dabashan foreland depression area.

Application of Singular Value Decomposition(SVD)to the Extraction of Gravity Anomalies Associated with Ag-Pb-Zn-W Polymetallic Mineralization in the Bozhushan Ore Field,Southwestern China

The Bozhushan Ore Field,located at the western margin of the South China Block,is an important area for Ag-Pb-Zn-W polymetallic mineralization which may be associated with the Late Cretaceous granitic magmaism.In this paper,the singular value decomposition(SVD)was effectively applied to decompose gravity data at scale of 1:50000 within the Bozhushan Ore Field to extract deep ore-finding information.Two gravity anomaly images displaying different scales of the ore-controlling factors were obtained.(1)The low-pass filtered image may reflect the deeply buried geological structures,hidden intrusions and concealed ore bodies.The negative gravity anomaly may reflect the overall distribution of granite bodies in the Bozhushan Ore Field.One negative gravity anomaly area may correspond to the exposed part of the Baozhushan granitic intrusion and the other corresponds to the concealed part of the granitic intrusion.The granitic intrusions are the main ore-controlling factors in this ore district.(2)The band-pass filtered image depicts the shallow concealed geological structures and geological bodies within this study area.There are two obvious negative gravity anomalies,which may be created by the hidden granites at different depths at both northwestern and southeastern sides of the exposed granitic intrusion.Thus the two negative gravity anomalies are favorable prospecting areas for various type of polymetallic ore deposits at depth.The gravity anomalies extracted by using the SVD exactly reflect the distribution of the ore deposits,structures and intrusions,which will give new insights for further mineral exploration in the study area.

Robust inversion analysis of local gravity anomalies caused by geological dislocation model of faults

Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.

Using WGM2012 to Compute Gravity Anomaly Corrections of Leveling Observations in China

Gravity Anomaly Correction(GAC)is a very important term in leveling data processing.In most cases,it is troublesome for field surveyors to measure gravity when leveling.In this paper,based on the complete Bouguer Gravity Anomaly(BGA)map of WGM2012,the feasibility of replacing in-situ gravity surveying in China is investigated.For leveling application,that is to evaluate the accuracy of WGM2012 in China.Because WGM2012 is organized with a standard rectangle grid,two interpolation methods,bilinear interpolating and Inverse Distance Weighted(IDW)interpolating,are proposed.Four sample areas in China,i.e.,Hanzhong,Chengdu,Linzhi and Shantou,are selected to evaluate the systems bias and precision of WGM2012.Numerical results show the average system bias of WGM2012 BGA in west China is about-100.1 mGal(1 mGal=10^(-5) m/s^(2))and the standard deviation is about 30.7 mGal.Tests in Shantou indicate the system bias in plain areas is about-130.4 mGal and standard deviation is about 6.8 mGal.All these experiments means the accuracy of WGM2012 is limited in high mountain areas of western China,but in plain areas,such as Shantou,WGM2012 BGA map is quite good for most leveling applications after calibrating the system bias.

Gravity and Magnetic Evidence for the Geological Setting of Major Mineral Systems of the Main Metallogenic Belts in South China: A Qualitative Analysis

South China is characterized by large-area multistage magmatism.It boasts a huge number of polymetallic deposits such as W-Sn,Cu-Au,rare earth deposits,thus serving as a"giant granary"of metal mineral resources in China(Lüet al.,2021).

Methodology for local correction of the heights of global geoid models to improve the accuracy of GNSS leveling

At present,one of the methods used to determine the height of points on the Earth’s surface is Global Navigation Satellite System(GNSS)leveling.It is possible to determine the orthometric or normal height by this method only if there is a geoid or quasi-geoid height model available.This paper proposes the methodology for local correction of the heights of high-order global geoid models such as EGM08,EIGEN-6C4,GECO,and XGM2019e_2159.This methodology was tested in different areas of the research field,covering various relief forms.The dependence of the change in corrected height accuracy on the input data was analyzed,and the correction was also conducted for model heights in three tidal systems:"tide free","mean tide",and"zero tide".The results show that the heights of EIGEN-6C4 model can be corrected with an accuracy of up to 1 cm for flat and foothill terrains with the dimensionality of 1°×1°,2°×2°,and 3°×3°.The EGM08 model presents an almost identical result.The EIGEN-6C4 model is best suited for mountainous relief and provides an accuracy of 1.5 cm on the 1°×1°area.The height correction accuracy of GECO and XGM2019e_2159 models is slightly poor,which has fuzziness in terms of numerical fluctuation.

On the Accuracy Assessment of the Latest Releases of GOCE Satellite-Based Geopotential Models with EGM2008 and Terrestrial GPS/Levelling and Gravity Data over Egypt

The Global Geopotential Models (GGMs) of GOCE (Gravity Recovery and steady- state Ocean Circulation Explorer) differ globally as well as regionally in their accuracy and resolution based on the maximum degree and order (d/o) of the fully normalized spherical harmonic (SH) coefficients, which express each GGM. The main idea of this study is to compare the free-air gravity anomalies and quasi geoid heights determined from several recent GOCE-based GGMs with the corresponding ones from the Earth Gravitational Model 2008 (EGM2008) over Egypt on the one hand and with ground-based measurements on the other hand. The results regarding to the comparison of GOCE-based GGMs with terrestrial gravity and GPS/levelling data provide better improvement with respect to EGM2008. The 4th release GOCE-based GGM developed with the use of space-wise solution strategy (SPW_R4) approximates the gravity field well over the Egyptian region. The SPW_R4 model is accordingly suggested as a reference model for recovering the long wavelength (up to SH d/o 200) components of quasi geoid heights when modelling the gravimetric quasi-geoid over the Egypt. Finally, three types of transformation models: Four-, Five- and Seven-parameter transformations have been applied to reduce the data biases and to provide a better fitting of quasi geoid heights obtained from the studied GOCE-based GGMs to those from GPS/levelling data. These models reveal that the standard deviation of vertical datum over Egypt is at the level of about 32 cm.

Consideration of the Daily Variation of Gravity on the Manifestation of Gravitational Shielding

The result of mathematical and physical analysis of the daily change in gravity is presented. The subject of consideration was the manifestation of semi-daily factors in diurnal variations of gravity. The assumption is investigated, according to which the cause of the half-day factors is the gravitational shielding of the planet Earth. Gravitational shielding is considered as a function of the size and thickness of celestial bodies and growing with distance from their poles. It is concluded that the planet Earth has the property of partial gravitational shielding, and the Moon does not have enough thickness to exhibit a tangible gravitational shielding. The obtained mathematical results correspond to the existing experimental data. It is suggested that gravitational shielding is the cause of the precession of the perihelion of Mercury and the peculiarities of the manifestation of tidal processes. It is assumed that gravitational shielding is one of the main reasons for the presence of hidden substances in the Universe. It is concluded that the physical picture with mutual shielding of interaction elements corresponds to the classical ideas of Fatio and Lesage. This approach is proposed as an alternative point of view to the existing theory on the description of tidal processes. It is shown that the existing basic approach to the description of tidal forces is unsatisfactory: the factors underlying the existing approaches have values 10 times less than those observed and cannot be considered as the reason for the manifestation of half-day manifestations in the daily change in gravity. The work is a continuation of the implementation by the author of P. Dirac’s ideas about accounting for the size of microparticles in physical theory.

Frozen subduction in the Yangtze block:insights from the deep seismic profiling and gravity anomaly in east Sichuan fold belt

The Sichuan basin is the main part of the middle-upper Yangtze block, which has been experienced a long-term tectonic evolution since Archean. The Yangtze block was regarded as a stable block until the collision with the Cathaysia block in late Neoproterozoic. A new deep seismic reflection profile conducted in the eastern Sichuan fold belt（ESFB） discovered a serials of south-dipping reflectors shown from lower crust to the mantle imply a frozen subduction zone within the Yangtze block. In order to prove the speculation, we also obtain the middle-lower crustal gravity anomalies by removing the gravity anomalies induced by the sedimentary rocks and the mantle beneath the Moho, which shows the mid-lower crustal structure of the Sichuan basin can be divided into eastern and western parts. Combined with the geochronology and Aeromagnetic anomalies, we speculated the Yangtze block was amalgamated by the West Sichuan and East Sichuan blocks separated by the Huayin-Chongqing line. The frozen subduction zone subsequently shifted to a shear zone accommodated the lower crustal shortening when the decollement at the base of the Nanhua system functioned in the upper plate.