Crustal structure of the Qiangtang and Songpan-Ganzi terranes(eastern Tibet) from the 2-D normalized full gradient of gravity anomaly

Numerous geophysical studies have revealed the lithospheric structure of the Qiangtang and the Songpan-Ganzi terranes in the eastern Tibetan Plateau.However,crust-mantle evolution and crustal response to the Indian lithospheric subduction are still controversial.Answering these questions requires additional information regarding crustal structure.In this study,the 2-D normalized full gradient(NFG)of the Bouguer gravity anomaly was used to investigate anomalous sources and interpret the crustal structure underneath the Qiangtang and Songpan-Ganzi terranes.The NFG-derived structures with loworder harmonic numbers(N=33 and N=43)showed that an anomalous source beneath the southern Qiangtang terrane had a characteristic northeastward-dipping shape,suggesting the northeastward motion of the crustal material induced by underthrusting Indian lithospheric mantle.The NFG images with harmonic number N=53 showed a large-scale anomalous source in the lower crust of the transformational zone from the Qiangtang terrane to the Songpan-Ganzi terrane,consistent with thickening crust and resistance of lower crustal flow.The anomalous source demonstrated by the NFG results with harmonic number N=71,located in the upper crust underneath the Ganzi-Yushu fault,suggested a seismogenic body of the 2010 M6.9 Yushu event.

Estimation of Depth to Salt Domes from Normalized Full Gradient of Gravity Anomaly and Examples from the USA and Denmark

We present an estimation of depth of anomalous bodies using normalized full gradient （NFG） of gravity anomaly. Maxima in the NFG map can locate the bodies and indicate their depth. Model calculation using a sphere and application of the NFG method to gravity anomalies over salt domes in the USA and Denmark shows effectiveness of the method. However, the accuracy of depth estimation strongly depends on the number of term N in the Fourier series used to calculate the NFG. An optimum N for the calculation can be given from a test.

Distribution law and susceptibility of geohazards across a gradient belt of the Western Sichuan Plateau

Across a gradient belt of the Western Sichuan Plateau,geohazards have seriously limited economic and social development.According to incomplete statistics,15,673 geohazards have been recorded in the study area.In order to mitigate the threat of geohazards to human engineering activities in the region,an overall understanding of the distribution pattern of geohazards and susceptibility assessment are necessary.In this paper,a gradient belt of the Western Sichuan Plateau and its zoning criteria were defined.Subsequently,on the basis of relief amplitude,distance to faults,rainfall,and human activities,three indicators of endogenic process were introduced:Bouguer gravity anomaly gradient,vertical deformation gradient,and horizontal deformation gradient.Thereafter,the distribution patterns of geohazards were investigated through mathematical statistics and ArcGIS software.By randomly selecting 10,449 hazards,a geohazard susceptibility map was generated using the Information Value(IV)model.Finally,the IV model was validated against 5224 hazards using the Area Under Curve(AUC)method.The results show that 47.6%of the geohazards were distributed in the zone of steep slope.Geohazards showed strong responses to distance to faults,human activities,and annual rainfall.The distribution of geohazards in the gradient belt of the Western Sichuan Plateau is more sensitive to vertical internal dynamics factors(such as vertical deformation gradient and Bouguer gravity anomaly gradient)without any apparent sensitivity to horizontal internal dynamics factors.The areas of high and very-high risk account for up to 32.22%,mainly distributed in the Longmenshan and Anning River faults.According to the AUC plot,the success rate of the IV model for generating the susceptibility map is 76%.This susceptibility map and geohazard distribution pattern can provide a reference for geological disaster monitoring,preparation of post-disaster emergency measures,and town planning.

Refining geoid and vertical gradient of gravity anomaly

We have derived and tested several relations between geoid （N） and quasi-geoid （~） with model validation. The elevation correction consists of the first-term （Bouguer anomaly） and second-term （vertical gradient of gravity anomaly）. The vertical gradient was obtained from direct measurement and terrain calcula- tion. The test results demonstrated that the precision of geoid can reach centimeter-level in mountains less than 5000 meters high.

Full gravity gradient tensors from vertical gravity by cosine transform

We present a method to calculate the full gravity gradient tensors from pre-existing vertical gravity data using the cosine transform technique and discuss the calculated tensor accuracy when the gravity anomalies are contaminated by noise. Gravity gradient tensors computation on 2D infinite horizontal cylinder and 3D ＂Y＂ type dyke models show that the results computed with the DCT technique are more accurate than the FFT technique regardless if the gravity anomalies are contaminated by noise or not. The DCT precision has increased 2 to 3 times from the standard deviation. In application, the gravity gradient tensors of the Hulin basin calculated by DCT and FFT show that the two results are consistent with each other. However, the DCT results are smoother than results computed with FFT. This shows that the proposed method is less affected by noise and can better reflect the fault distribution.

Edge enhancement of gravity anomalies and gravity gradient tensors using an improved small sub-domain filtering method

In order to enhance geological body boundary visual effects in images and improve interpretation accuracy using gravity and magnetic field data, we propose an improved small sub-domain filtering method to enhance gravity anomalies and gravity gradient tensors. We discuss the effect of Gaussian white noise on the improved small sub-domain filtering method, as well as analyze the effect of window size on geological body edge recognition at different extension directions. Model experiments show that the improved small sub-domain filtering method is less affected by noise, filter window size, and geological body edge direction so it can more accurately depict geological body edges than the conventional small sub-domain filtering method. It also shows that deeply buried body edges can be well delineated through increasing the filter window size. In application, the enhanced gravity anomalies and calculated gravity gradient tensors of the Hulin basin show that the improved small sub-domain filtering can recognize more horizontal fault locations than the conventional method.

Theoretical simulation of the effect of deformation on local gravity in a density gradient zone

We modeled the effect of the deformation of a Density Gradient Zone （DGZ） on a local gravity field using a cubical model and introduced a new method to simulate a complex DGZ （CDGZ）. Then, we analyzed the features of the model for the influence of the deformation of the DGZ on the local gravity field. We concluded that land-based gravity is not sensitive to the thickness of the DGZ and that the magnitude of the contribution of the DGZ is one order less than that of the volume strain with the same displacement.

Interpretation of Gravity Anomalies by Multi-Scale Evaluation of Maxima of Gradients and 3D Modelling in Bipindi Region (South-West Cameroon)

The gravity maps display, in Bipindi zone, local oval culminations of low anomalies indicative of a presence of intrusive light body in a subsurface but the nature, the form and the position of this body are still unknown. The analyses of established gravimetric anomaly maps, the multi-scale evaluation of maxima of gradients and the quantitative interpretation of residual anomalies by 3D modelling permit characterizing the intrusive light body situated at Bipindi. The multi-scale evaluation of maxima of gradients shows that the modelling of the intrusive light body of Bipindi can be done without the problem of interference of anomalies from different sources. The 3D model of Bipindi zone shows two dissymmetrical blocks of the same type of rock with a density contrast of -0.095 g·cm-3 in comparison with the density of the surrounding metamorphic rocks. The two blocks are at a distance about 3 km from one to another. The upper surfaces of these blocks lie at a depth between 1 and 2 km. Their lower surfaces have two landings;one lies at a depth of about 8 km and another at a depth about 14 km. A consideration of the density of the modelled body, of the ranges of densities of specific rocks present in the general region indicates that the body may be composed of nepheline syenites. The intrusive body of Bipindi is situated in a senestral shear zone. The area situated between the two blocks of this intrusive body may be indicated for a detail study in the domain of mineral research.

Combined analysis of gravity and magnetic anomalies using normalized source strength

The traditional combined gravity and magnetic analysis uses the linear regression of the first order vertical derivative of the gravity anomaly and the reduction to the pole(RTP) magnetic anomaly,and provides the quantitative or semi-quantitative interpretation by calculating the correlation coefficient,slope,and intercept.In the calculation process,due to the remanent magnetization,the RTP anomaly still contains the effect of oblique magnetization,as a result,the homologous gravity and magnetic anomalies may display irrelevant results in the linear regression calculation.To solve this problem,we present a new combined analysis using normalized source strength(NSS).Based on the Poisson's relation,the gravity field can be transformed into the pseudomagnetic field of the direction of geomagnetic field magnetization under the homologous condition.The NSS of the pseudomagnetic field and that of the original magnetic field are calculated,which are insensitive to the remanence,and then the linear regression analysis is carried out.The approach is tested using synthetic model under complex magnetization,the results show that it can still identify the gravity and magnetic anomalies from same source under strong remanence,and can establish the Poisson's ratio.Finally,this approach is applied in Wudalianchi in China.The results demonstrated that this approach is feasible and can provide the reference for further data processing and interpretation.

Polar 3D Transformation of the Full Gradient of Attractive Potential

The method of 3D polar transformation of full gravity potential gradient vectors is based on the geometric properties of the crossing points of complete gradient of the potential to localize the source region that causes the observed anomaly. The cross-points—poles—are defined for rectangular polygons of different sizes where the full gradient vector is defined at every vertex. The polygon size range could be specified. The set of poles, positive and negative, is then represented on the 3D chart in the form of clusters of dots or cubes and can be considered as a model image of the sources, intended for visual analysis and further interpretation.

Gravity pattern in southeast margin of Tibetan Plateau and its implications to tectonics and large earthquakes

There are many active faults in the southeast margin of Tibetan Plateau,where three large active faults zones,the Longmenshan,Xianshuihe and Anninghe,merge to form a"Y"shape.Strong crustal deformation and a complicated fault distribution accompany strong earthquake activity in this zone.In this paper,we investigate a multi-scale gravity anomaly in the southeastern margin of the Tibetan Plateau using the wavelet transform;we find that the pattern of the gravity field is closely related to the fault system in the study area.Analyzing the characteristics of this Bouguer gravity anomaly at different orders indicates that the eastern Himalayan syntaxis has produced a strong eastward push during its northward movement,resulting in a shortening of the crust from west to east and a rapid uplift of the Tibetan Plateau.The Songpan–Garzêand Sichuan–Yunnan blocks have been forced to slip and extrude southward and eastward laterally.The distributions of seven large earthquakes from 1970 to 2018 reflects the relationship between large earthquakes and characteristics of the gravity anomaly.Comparing the tectonic backgrounds of several earthquakes reveals that the large earthquakes occur usually in the high gravity anomaly gradient zone,which corresponds in general to the boundary zones of the blocks.We infer that large earthquakes occur primarily in high Bouguer gravity anomaly zones in the upper crust,while low Bouguer gravity anomalies encompass the lower crust and the uppermost mantle.

Simultaneous Modelling of Gravity and Magnetic Data in a Measured Heat Flux Area to Characterize Geothermal Heat Sources: A Case for Eburru Geothermal Complex, Kenya

Forward modelling of gravity and magnetic data was done simultaneously to show the correlation between gravity and magnetic anomalies on a measured heat flux region. The results were used to characterize the heat source structures in Eburru area. Modelling was done using Oasis montaj geosoft software which is an iteration process where the gravity and magnetic anomalies were calculated and compared to the observed residual anomaly until there was a fit. The start model was constructed based on depths from Euler deconvolution and models constrained using stratigraphy data from the existing wells in the study area. Forward modelling of gravity and magnetic data revealed intrusions within the Earth’s subsurface with depth to the top of the sources ranging from 739 m to 5811 m. The density of the sources ranges between 3.0 g/cm3 and 3.2 g/cm3 while their magnetic susceptibility was zero. This implies that intrusions from the mantle with a magnetic susceptibility of zero have temperatures exceeding the curie temperature of rocks. The density of the intrusions modelled was higher than 2.67 g/cm3, the average crustal density, hence it explains the observed positive gravity anomaly. The results also revealed that areas with high heat flux have shallow heat sources and if the heat sources are deep, then there must be a good heat transfer mechanism to the surface.

Gravity anomaly and crustal density structure in Jilantai rift zone and its adjacent region

This paper deals with the interpretation of Bouguer gravity anomalies measured along a 250 km long Suhaitu-Etuokeqi gravity profile located at the transitional zone of the Alxa and Ordos blocks where geophysical characteristics are very complex. The analysis is carried out in terms of the ratio of elevation and Bouguer gravity anomaly, the normalized full gradient of a section of the Bouguer gravity anomaly （Gh） and the crustal density structure reveal that （1） the ratio of highs and lows of elevation and Bouguer gravity anomaly is large between Zhengyiguan fault （F4） and Helandonglu fault （F6）, which can be explained due to crustal inhomogeneities related to the uplift of the Qinghai-Tibet block in the northeast; （2） the main active faults correspond to the Gh contour strip or cut the local region, and generally show strong deformation characteristics, for example the Bayanwulashan mountain front fault （F1） or the southeast boundary of Alxa block is in accord with the western change belt of Gh, a belt about 10 km wide that extends to about 30 km; （3） Yinchuan- Pingluo fault （FS） is the seismogenic structure of the Pin- gluo M earthquake, and its focal depth is about 15 km; （4） the Moho depth trend and Bouguer gravity anomaly vari- ation indicates that the regional gravity field is strongly correlated with the Moho discontinuity.

融合多源海洋大地测量数据的南海海底地形多层感知机反演

本文融合SIO(Scripps Institution of Oceanography)发布的垂线偏差、重力异常和垂直重力梯度数据及NCEI(National Centers for Environmental Information)发布的船载测深数据,利用多层感知机神经网络(Multi-Layer Perceptron,MLP)建立南海海域(108°E—121°E,6°N—23°N)分辨率为1'×1'的海底地形模型(MLP_Depth).首先,将642716个船载测深控制点的位置信息与周围4'×4'格网点处的地球重力信息(垂线偏差、重力异常、垂直重力梯度)作为输入数据,将船载测深控制点处实测水深值作为输出数据,训练MLP神经网络模型,训练结束时决定系数R2为99%,平均绝对误差MAE为39.33 m.然后,将研究区域内1'×1'格网正中心点处的输入数据输入于MLP模型中,可得格网正中心点处的预测海深值.最后,根据预测海深值建立研究区域范围内分辨率为1'×1'的MLP_Depth模型.将MLP_Depth模型预测水深与160679个检核点处实测水深对比,其差值的标准差STD(75.38 m)、平均绝对百分比误差MAPE(5.89%)与平均绝对误差MAE(42.91 m)皆优于GEBCO_2021模型、topo_23.1模型、ETOPO1模型与检核点实测水深差值的STD(108.88 m、113.41 m、229.67 m)、MAPE(6.11%、6.94%、18.37%)与MAE(47.33 m、52.24 m、130.08 m).同时,为了研究不同区域内利用该方法建立的海底地形模型的精度,本文在研究区域内分别建立了A、B区域的海底地形模型(MLP_Depth_A、MLP_Depth_B).经过验证得:MLP_Depth_A、MLP_Depth_B相比于MLP_Depth模型具有更高的精度,更能反应海底地形的变化趋势.

重力数据探测地下目标的贝叶斯方法研究

针对地下目标空间位置和形状均未知的情况,基于重力数据,分别以区域测量和线测量2种场景进行建模分析,研究贝叶斯推断用于地下目标探测的效果。结果表明,将贝叶斯推断用于重力梯度探测,比微重力探测更具优势,能有效反演出目标的未知参数。

基于CEP的重力自适应并行EKF匹配算法

针对现有惯性/重力/重力梯度组合导航EKF匹配算法滤波状态方程存在模型误差以及惯性导航系统定位存在累积误差而造成的滤波失准乃至发散问题,提出一种基于CEP(Circular Error Probable,圆概率误差)的自适应并行EKF匹配算法.该算法首先通过自适应因子调节状态预测信息的权重,削弱预设动力学模型不准确产生的误差;同时利用惯性导航系统圆概率误差半径构建基于CEP的移动窗口分层模型,然后根据滤波量测值与窗口坡度等信息,对移动窗口范围进行约束;最后组建分层窗口并行滤波器,得到最优匹配结果.南海海域实验结果表明,基于CEP的自适应并行EKF匹配算法相较于传统EKF算法和自适应EKF算法的水下重力匹配导航定位精度分别提升了74.0%和49.8%.该算法能够在一定程度上克服惯性导航系统由于时间推移误差积累的缺陷,提高系统导航定位精度,增加匹配算法的鲁棒性.

基于地形单元分区的BP神经网络反演海底地形

海底地形和海洋重力场在一定波段频率域上存在相关性,传统海底地形反演方法主要考虑海底地形与重力数据之间的线性关系,而忽略非线性项作用对反演结果的影响。基于此,提出了基于海底地形单元分区的BP神经网络反演方法,改善海底地形反演精度。选取东印度洋北部东经90°海岭附近海域(83°E~92°E,10°S~10°N)为实验区,利用船测水深、卫星测高重力等数据,根据地形特征将实验区分为海盆区和海岭区,使用BP神经网络方法分别进行海底地形反演,构建了实验区1′×1′局部海底地形模型。结果表明,基于地形单元分区的反演结果模型平均相对误差精度可达1.45%,比未分区的反演结果均方根误差降低了42 m,验证了本文方法的有效性和可行性。

似大地水准面与大地水准面转换分析

传统的似大地水准面与大地水准面近似转换在地形复杂的山区难以满足精度需求。本文基于似大地水准面与大地水准面严密转换方法,计算珠峰地区转换改正数值大小,分析似大地水准面与大地水准面转换的空间变化。结果表明,严密方法计算的转换值变化范围为-3.270~-0.119 m,其中地形位差改正极值达到1.272 m,重力梯度改正极值达到-0.138 m,这两项改正影响幅度较大,在山区需顾及;严密方法计算和近似计算结果都与地形高度具有一定的相关性,但前者结果相较于后者更为平滑,在布格重力异常与地形变化的共同影响下,严密方法局部特征与地形高度负相关。

基于重力水平总梯度倾斜角的加强处理方法在断裂识别中的应用

断裂两侧因密度差异而产生重力异常带,对水平总梯度方法处理的重力异常,通过极值点的追踪,进行断裂识别。对于埋深较大的深大断裂和中等埋深的次一级断裂,要在重力水平总梯度基础上,利用倾斜角的概念,对重力水平总梯度异常进行“倾斜角”加强处理,实现对弱异常的增强提取,进而识别出更多断裂信息。利用重力水平总梯度倾斜角的加强处理方法,在博格达山前复杂构造带米泉地区预测了区域性深大断裂及次一级断裂的发育特征,综合电法预测成果,研究区发育逆冲或逆掩断层,断裂识别结果对断裂发育模式的认识具有重要的辅助作用。

Determination of vertical gradient of gravity anomaly with topographic data

The role of vertical gradient of gravity anomaly in the solution and interpretation of geological and geophysical problems on superficial layer is increasingly evident. However, it is difficult to directly observe the vertical gradient of gravity anomaly in mountainous area. Moreover,aerial and satellite gravity gradient measurments have not been put into practice up to date yet.The vertical gradient of gravity anomaly is calculated by use of topographic data and several problems concerned are investigated and discussed. It turns out in theory and practice that the relative error of the vertical gradient determined with the method is only 5%, and when the vertical gradient of gravity anomaly reaches the maximum 1 000 ns-2 the error is only 50 ns-2. In addition, how to divide ring zones in calculation and the errors caused by the inaccuracy of landform elevation and rock density has been quantitatively analysed.