The geoid-to-quasigeoid difference using an arbitrary gravity reduction model

Recently it was proved that the classical formula for computing the geoid to quasigeoid separation (GQS) by the Bouguer gravity anomaly needs a topographic correction. Here we generalize the modelling of the GQS not only to Bouguer types of anomalies, but also to arbitrary reductions of topographic gravity. Of particular interest for practical applications should be isostatic and Helmert types of reductions, which provide smaller and smoother components, more suitable for interpolation and calculation, than the Bouguer reduction.     

中国海及邻域重力场特征及其构造解释

利用中国海及邻域 2’× 2’由Seaset,Geosat,ERS 1及TOPEX/POSEIDON等卫星测高资料解算的自由空气重力异常 ,结合ETOPO5全球海底地形资料计算了中国海及邻域布格重力异常 ,并反演了得到了中国海及邻域地壳厚度 .经对重力异常特征进行分析 ,讨论了重力异常与大地构造及其活动的相关性 .     

TOPOGRAPHIC AND ISOSTATIC CORRECTIONS TO GRAVITY SURVEYS IN MOUNTAINOUS AREAS*

In mountainous areas it is an inadequate procedure to reduce gravity observations by merely subtracting the effect of an infinite flat slab of material between the station and sea-level, and adding a terrain correction. A programme is described which directly computes the effect of masses above sea-level, and mass-deficiencies below it; as well as the effect of compensating masses under the Pratt and Airy isostatic systems. As an example, the method has been applied to a regional gravity survey of Papua and New Guinea where it is seen to remove the usually high correlation of the Bouguer anomaly with local topography.     

基于EGM2008重力场模型的GPS高程转换方法及精度分析

研究利用重力场模型直接进行GPS高程转换的方法,即模型高程异常法和模型重力位法.根据误差传播定律,导出了由重力场模型位系数方差计算模型高程异常以及模型正常高的误差公式,并对影响两种GPS高程转换方法精度的因素进行分析.利用具有代表性的EGM96、EIGEN-51C和EGM2008重力场模型及某线路GPS/水准数据对两种方法进行实验研究与精度分析,结果表明:理论公式及实验结果都显示两种高程转换方法是等价的,转换精度相当;两种方法使用过程中,均可能需要系统偏差校正;EGM2008重力场模型具有非常好的分辨率及精度,在局部区域,基于EGM2008模型的GPS高程转换精度有望达到cm级.     

An estimation of the height system bias parameter <Emphasis Type="Italic">N</Emphasis><Subscript>0</Subscript> using least squares collocation from observed gravity and GPS-levelling data

This paper deals with the analysis of gravity anomaly and precise levelling in conjunction with GPS-Levelling data for the computation of a gravimetric geoid and an estimate of the height system bias parameter N_o for the vertical datum in Pakistan by means of least squares collocation technique. The long term objective is to obtain a regional geoid (or quasi-geoid) modeling using a combination of local data with a high degree and order Earth gravity model (EGM) and to determine a bias (if there is one) with respect to a global mean sea surface. An application of collocation with the optimal covariance parameters has facilitated to achieve gravimetric height anomalies in a global geocentric datum. Residual terrain modeling (RTM) technique has been used in combination with the EGM96 for the reduction and smoothing of the gravity data. A value for the bias parameter N_o has been estimated with reference to the local GPS-Levelling datum that appears to be 0.705 m with 0.07 m mean square error. The gravimetric height anomalies were compared with height anomalies obtained from GPS-Levelling stations using least square collocation with and without bias adjustment. The bias adjustment minimizes the difference between the gravimetric height anomalies with respect to residual GPS-Levelling data and the standard deviation of the differences drops from 35 cm to 2.6 cm. The results of this study suggest that N_o adjustment may be a good alternative for the fitting of the final gravimetric geoid as is generally done when using FFT methods.     

A regional strategy for geothermal exploration with emphasis on gravity and magnetotellurics

As part of the resource evaluation and exploration program conducted by Los Alamos Scientific Laboratory for the national Hot Dry Rock (HDR) Geothermal Program, a regional magnetotelluric (MT) survey of New Mexico and Arizona is being performed. The MT lines are being located in areas where the results of analysis of residual gravity anomaly maps of Arizona and New Mexico, integrated with other geologic and geophysical studies indicate the greatest potential for HDR resources.The residual gravity anomalies are derived by applying the concept of predicting gravity anomalies from topography. This can be accomplished by employing reductions similar to those used in some isostatic investigations, in which a regional topographic surface is used as the Bouguer reduction datum. The datum is derived by comparison of various harmonics of Bouguer anomalies and elevations of stations. Topography can be used to predict Bouguer anomalies because of isostatic compensation; the resultant anomalies can be considered high frequency residual anomalies or isostatic anomalies corrected for regional compensation. Such maps have been produced for Arizona, New Mexico, west Texas, and Chihuahua, Mexico.The main objective of the MT project is to produce a regional geoelectric contour map of the pervasive deep electrical conductor within the crust and/or upper mantle beneath the Colorado Plateau and the adjacent Basin and Range Province and Rio Grande Rift. The MT survey consists of 200 sites along several long profiles with site spacing of 15–20 km. Pre-existing available MT data are being integrated with the new data. After the data are processed, a one-dimensional inversion is applied to the sounding curve and used as a starting point for 2-D modeling. Such a project and ultimate map will be of major value in studying the regional geophysics and tectonics of the southwest United States as they now apply to HDR resources in particular and geothermal resources in general.Electrical conductivity anomalies of large areal extent are of particular interest in geothermal exploration. Correlation analysis of large conductive anomalies with other geophysical, geological, and geotectonic data is being performed. Preliminary analysis of the data has suggested several major regions of anomalously shallow high electrical conductivity. Among these is the Aquarius area of northwest Arizona which is the site of a longwavelength residual anomaly low, which when modeled and correlated with other geophysical data can be shown to be possibly related to low density and high temperature in the crust at depths of 20 km or less. Preliminary analysis of MT data indicates the possible existence of a mid-crustal high electrical conductivity anomaly in this same region.     

ESTIMATION OF THE DEGREE OF POLYNOMIAL FITTED TO GRAVITY ANOMALIES AND ITS APPLICATION1

A method for estimating the degree of polynomial fitted to gravity anomalies to evaluate the regional anomaly is presented. The anomaly can be fitted by polynomials of different degrees with the least-squares method and the optimum degree of the polynomial evaluating this anomaly can be estimated from the point of discontinuity of the gradient on a graph of variance against the polynomial degree. The Bouguer gravity is initially separated by upward continuation to a proper height and then the degree of regional polynomial to fit the Bouguer anomaly can be estimated. Theoretical and field examples show the effectiveness of the method.     

Crustal structure and tectonic setting of the south central Andes from gravimetric analysis

Based on terrestrial gravity data, in this paper we prepared a map of Bouguer anomalies, which was filtered to separate shallow and deep gravity sources. Based on a density model and gravimetric inversion techniques, the discontinuous crust-mantle boundary and the top of crystalline basement were modeled. Subsequently, the equivalent elastic thickness (Te) was evaluated, considering information from the crust-mantle discontinuity and topographic load, finding high Te values in the eastern Andean foothills and west of the Velasco range. These results are consistent with the positive isostatic and residual Bouguer anomaly values, which suggest the presence of high-density rocks in the mid-to upper crust. In addition, petrographic and geochemical analysis conducted in surface outcrops suggest a mantle origin.     

渭河盆地及邻区重力异常小波多尺度分解与解释

基于EGM2008重力场模型计算获得了渭河盆地及邻区布格重力异常。采用小波多尺度分解方法对布格重力异常进行了4阶小波逼近和小波细节分解,同时基于平均径向对数功率谱方法定量化地计算出1~4阶小波细节和小波逼近所对应的场源平均埋深。结合区域地质和地震资料,对获得的重力场结果进行分析,得到如下结论:①鄂尔多斯地块、渭河盆地、秦岭造山带3个一级构造单元的布格重力异常之间存在明显差异;构造区内部重力异常也存在横向的显著差异。布格重力异常的走向、规模、分布特征与二级构造区及主要的断裂具有一定的对应关系。②渭河盆地及邻区布格重力异常1~4阶细节对应4~23 km不同深度的场源信息,鄂尔多斯地块南缘东、西部的地壳结构存在明显的差异;渭河盆地凹陷、凸起构造区边界清晰,断裂边界与重力异常边界具有较好的一致性;秦岭造山带重力异常连贯性不好,东、西部重力异常变化特征表现出明显的差异。③渭河盆地及邻区布格重力异常分布与莫霍面埋深具有非常明显的镜像关系。渭河盆地及邻区地震主要分布在六盘山—陇县—宝鸡断裂带、渭河断裂与渭南塬前断裂交汇处、韩城断裂与双泉—临猗断裂交汇处。渭河盆地及邻区重力异常主要由中上地壳剩余密度体所影响,这可能是该区地震以浅源地震为主的主要原因。     

青藏高原东南缘多尺度重力场及构造含义

本文基于EGM2008重力场模型研究了青藏高原东南缘均衡重力异常和多尺度的布格重力异常特征,以鲁甸和景谷地震为例,认识其深部构造环境和动力学过程,为该区域的构造运动和地震孕育环境研究提供依据.结果表明,研究区布格重力异常和均衡重力异常与地质构造格局相关性较好,川滇地块剧烈的区域布格重力异常和非均衡状态与其强烈的地壳变形...     

应用多种来源重力异常编制中国海陆及邻区空间重力异常图及重力场解读

本文通过分析陆地实测空间重力异常数据、海洋船载测量空间重力异常数据、卫星测高重力异常,布格重力异常数据、EGM2008地球重力模型数据等多种来源数据的性质和精度,并对相关数据进行对比,研究了编制1:500万中国海陆空间重力异常图的数据使用方案和技术方法.在地形较为平坦、实测数据分布均匀的陆区,使用实测数据,在地形复杂,实测数据稀少以及没有实测数据的陆区或岛屿,利用布格重力异常反推空间异常的方法合成平均空间重力数据,西藏地区的数据对比实验证明合成平均空间重力异常数据是一种有效的数据补充.利用三观测列方差分解法在南海地区对船载测量空间重力数据和美国SS系列及丹麦DNSC08GRA卫星重力数据进行了方差分解计算,结果表明不同来源的卫星测高重力数据具有很大的一致性,数据精度较以往有了很大的提高.海区空间重力数据使用原则是在船载重力测量数据校准下,全面使用卫星测高重力数据进行编图.海陆过渡区的异常处理应以EGM2008地球重力模型重力场为基准参考场,实现海陆异常平缓过渡,无缝连接.对中国海陆空间重力异常场进行了小波变换处理,对空间重力异常场进行了解读,勾画出三横四竖的一级重力梯级带及其所围限的8个一级重力异常区,并划分了二级重力异常区和梯级带,为块体构造学体系中大地构造格架的建立提供了地球物理证据.     

A REGIONAL GRAVITY STUDY OF THE NORTHERN NORTH SEA (56–62° N)*

A Bouguer gravity anomaly map is presented of the North Sea and adjacent land areas in Norway and Denmark, covering an area situated between 56° and 62°N, 1°W and 10°E. The gravity data from the UK sector of the North Sea, the land and offshore areas of Denmark, and the land areas of Norway have been published before. However, the gravity data from the Norwegian sector of the North Sea are new. A large number (about 60) of individual gravity features can be defined in the mapped area. Most of those situated in the UK sector of the North Sea and on land in Norway have been discussed earlier; however, most of the anomalies found elsewhere which are qualitatively interpreted here have not been discussed before. An interpreted Bouguer anomaly map is presented which identifies all these features. The majority of the gravity anomalies encountered in the mapped area can be shown to be associated with one of the following geological features: (i) basement highs, (ii) large bodies of heavy basic or ultrabasic rock in the crystalline basement, (iii) large igneous intrusions within the sedimentary column and thick accumulations of volcanic rocks or their associated eruption centers, (iv) major basement faults. Large-scale geological structures such as the Central, Viking and Sogn Grabens and the East Shetland, Stord, Forth Approaches and Norwegian-Danish Basins are essentially in isostatic equilibrium and are only locally marked by relatively weak gravity minima. A residual gravity anomaly map has been produced by subtracting from the observed Bouguer anomalies the estimated gravity effect of an assumed thinned crust. This residual gravity anomaly map shows a number of features of the Bouguer anomaly field with greater clarity.     

A geoid solution for airborne gravity data

Airborne gravity data is usually attached with satellite positioning of data points, which allow for the direct determination of the gravity disturbance at flight level. Assuming a suitable gridding of such data, Hotine’s modified integral formula can be combined with an Earth Gravity Model for the computation of the disturbing potential (T) at flight level. Based on T and the gravity disturbance data, we directly downward continue T to the geoid, and we present the final solution for the geoid height, including topographic corrections. It can be proved that the Taylor expansion of T converges if the flight level is at least twice the height of the topography, and the terrain potential will not contribute to the topographic correction. Hence, the simple topographic bias of the Bouguer shell yields the only topographic correction. Some numerical results demonstrate the technique used for downward continuation and topographic correction.     

Allowance for the vertical gradient in an airborne gravity survey

When processing the results of an airborne gravity survey flown in the Arctic at flight heights of 2000 to 3000 m, the authors encountered conditions, when the values of the free-air anomaly of gravity exceeded the values taken from the anomaly map, with a trend towards “heavier” values. The formula for calculation of the vertical gradient as a function of the latitude of the locality has been refined. The tentative calculations reported in the paper also indicate that a free-air gravity anomaly above the mountainous areas will contain additional positive systematic components, since only the effect of the vertical gradient of the normal gravity field is taken into account here, and the flight’s height during an airborne gravity survey, above a mountainous area, should be held constant, otherwise it would be necessary to take into account additional systematic biases, obtained in different flights, which are caused by the different values of the vertical gradient at different heights.     

华北北部均衡重力异常的初步研究

本文利用地面实测的重力资料和地形高程资料,采用爱里-海斯于宁(Airy-Heiskenen)均衡理论模式,构制了华北北部N37°—41°,E115°—120°范围内的均衡重力异常图和1°×1°均衡重力异常图。并对该图进行了分析、解释,同时结合本区地质构造、现代构造运动以及地震活动性特点分别进行了讨论     

应用卫星重力信息对横断山系地区布格重力异常特异分布的纠正

在云南省西部,跨越中、缅两国交界的横断山系地区（97°E~102°E,24°N~30°N）有近一半的面积尚没有重力测点、即重力数据空白区和重力测点稀少的普查级测区.以前的有关文献、图集中所给出对此地区的重力场都是十分模糊的结果与图件.因此应用这些资料无法详细地研究该地区重力场特征与深部地壳结构（构造）.本文应用卫星重力异常资料作为“近似空间重力异常”,经中间层改正后给出的“计算布格重力异常”,其分布特征与该地区的地形高程呈很好的镜像相关.对相应山脉、河谷以及断裂构造都有所反映.特别是在横断山系地区该布格重力异常呈现为近南北的走向.为此,据该“计算布格重力异常”,并选定对该区有代表性的一条重力异常剖面作正反演计算,以得到其地壳深部结构剖面.结果表明,在横断山脉地区的地壳厚度在51~56 km间起伏变化;滇西北云岭山系以及玉龙山区的地壳厚度约在60 km以上. 最后,对所得结果与图件进行了讨论,并提出了几点认识和纠正的建议.     

Detailed gravity study of the offshore structure of Piton de la Fournaise volcano,Reunion Island

We present an interpretation of gravity data acquired in 1984 by the French R/V Jean Charcot on the submarine part of the eastern flank of Piton de la Fournaise volcano. We comment on the Bouguer anomaly map and give a quantitative interpretation of three gravity profiles. The main results are that a gravity high over Grand Brûlé, the lower subaerial part of the eastern flank, does not extend far offshore and that an anomalous topographic feature, discovered in 1982 on the submarine eastern flank, is characterized by a large negative anomaly. We propose three hypotheses to explain the origin of this anomaly, i. e., it marks the site of a new volcano, or it is a consequence of lateral volcanism from a volcano older than Piton de la Fournaise, or more probably, it represents a great landslide deposit.     

我国重力场新的统计特征参数的计算分析

基于地形信息,构造了完全布格异常和完全空间异常的协方差和代表误差模型参数等新统计量,利用某试验区密集的重力和地形数据进行了统计分析,并将计算方案推广到了全国大范围、多区域内的统计计算,给出了6种不同地形类别区域（平原、丘陵、小山区、中山区、大山区、特大山区）的完全空间异常和完全布格异常的方差、协方差以及代表误差模型参数.试验结果表明：依据本文的统计模型、算法与思路,在实际测量数据的支撑下,可以给出全国范围的统计参数的网格值、等值线图等,提供后续重力场相关研究工作使用.     

广义线性高程拟合法在物探中的应用

地球物理勘探工作需要三维坐标为其定位,尤其是重力勘探,对高程有较高的要求。作为一种先进的现代化测量技术,GPS能实时提供高精度的平面坐标,但因其使用W G S-84坐标系统,高程为大地高,没有实际的地球物理意义。地球物理勘探用的是正常高,故GPS测量的大地高无法直接使用。若联测高程,不仅增加工作难度,时间成本和经济成本也将大大增加。本文以长江中下游某些地区地球物理矿产勘查为例,介绍利用广义线性拟合法来获取高程。该模型不仅能对拟合结果进行评价,亦能评价模型参数的优劣。