月球Airy均衡状态与月壳厚度估计

月球水准面异常和表面地形变化是其内部密度不均匀和各个界面的起伏变化的体现,因此利用水准面和地形之比(geoid to topography ratio,GTR)可估计月球均衡和月壳厚度。本文基于月球重力场模型SGM100h和地形模型STM359_grid-02,经过去除表面玄武岩填充和深层异常质量影响,并结合理论Airy均衡模型中GTR与参考月壳厚度的关系,计算得到了新的月壳厚度模型。该模型的月壳平均厚度为36.9 km,背面比正面平均厚约13.5 km,Apollo12/14登陆点的月壳厚度分别是28.3 km和29.1 km。在各月海盆地存在着中央较薄、四周逐渐增厚的趋势。     

利用低阶大地水准面异常反演大尺度核幔边界起伏

核幔边界（core-mantle boundary,CMB）是地球内部最重要的物理化学界面之一,地核和地幔通过核幔边界发生多种相互作用,这对地球重力场、地球自转及地磁场等都能产生重要影响。大地水准面异常是地球重力场的重要观测量,反映了地球内部的物质密度异常及界面变化等重要信息。推导了通过大地水准面异常反演核幔边界起伏的公式,利用2~4阶大地水准面异常反演了大尺度核幔边界起伏形态。结果显示,核幔边界起伏的径向幅度达±5 km、与Morelli的地震层析成像结果的幅度接近,但在形态上略有差异。以高为5 km、底边长为1 000 km的棱柱体模型模拟计算了核幔边界密度异常引起的大地水准面异常响应,结果与观测大地水准面异常比较接近。     

Minimization and estimation of geoid undulation errors

The objective of this paper is to minimize the geoid undulation errors by focusing on the contribution of the global geopotential model and regional gravity anomalies, and to estimate the accuracy of the predicted gravimetric geoid.The geopotential model's contribution is improved by (a) tailoring it using the regional gravity anomalies and (b) introducing a weighting function to the geopotential coefficients. The tailoring and the weighting function reduced the difference (1) between the geopotential model and the GPS/levelling-derived geoid undulations in British Columbia by about 55% and more than 10%, respectively.Geoid undulations computed in an area of 40° by 120° by Stokes' integral with different kernel functions are analyzed. The use of the approximated kernels results in about 25 cm () and 190 cm (maximum) geoid errors. As compared with the geoid derived by GPS/levelling, the gravimetric geoid gives relative differences of about 0.3 to 1.4 ppm in flat areas, and 1 to 2.5 ppm in mountainous areas for distances of 30 to 200 km, while the absolute difference (1) is about 5 cm and 20 cm, respectively.A optimal Wiener filter is introduced for filtering of the gravity anomaly noise, and the performance is investigated by numerical examples. The internal accuracy of the gravimetric geoid is studied by propagating the errors of the gravity anomalies and the geopotential coefficients into the geoid undulations. Numerical computations indicate that the propagated geoid errors can reasonably reflect the differences between the gravimetric and GPS/levelling-derived geoid undulations in flat areas, such as Alberta, and is over optimistic in the Rocky Mountains of British Columbia.Paper presented at the IAG General Meeting, Beijing, China, August 8–13, 1993.     

An isostatically compensated gravity model using spherical layer distributions

A new isostatic model for the Earths gravity field is presented based on a simple hypothesis of layers approximating constant density contrasts. The spherical layer distribution used to describe the hydrostatic equilibrium of the Earths masses leads to a new set of spherical harmonic coefficients for the gravitational potential. First attempts to quantify the information content of these coefficients led to the outcome that they seem to explain the observed gravity field for a certain wavelength band, while they are insufficient for short and very long wavelengths. A synthesis of the derived coefficients over specific degree ranges provided a computation of band-limited geoid undulations on a global scale. The association of these potential quantities with known tectonic structures, such as the topography of the core–mantle boundary, strengthens the belief that the interpretation of Earth gravity models, especially those arising from global digital elevation models, should be considered in close relation with deep-Earth structure.     

基于实测重力数据的秦岭地区均衡异常与垂向构造应力分析

确定秦岭地区准确的均衡异常与垂向构造应力对认识该地区深部构造特征、动力学机制等具有重要意义。目前,高精度重力数据是获取均衡异常与垂向构造应力的重要手段之一。采用高精度重力/GPS联测数据,得到秦岭地区的自由空气重力异常与布格重力异常,反演均衡异常与垂向构造应力。研究结果表明,均衡面深度在40~49 km之间,莫霍面深度范围为39~48 km,垂向构造应力大小在-28~24 MPa之间。秦岭北侧的渭河盆地处于不均衡状态;四川盆地北部的均衡异常与垂向构造应力几乎为零,地壳处于均衡状态;在秦岭出现局部的负均衡异常,表明存在一定的地壳运动。     

Spherical harmonic analysis of the CRUST 2.0 global crustal model

The spherical harmonic analysis of the 2×2 density and stratification information contained in the global crustal model CRUST 2.0 is presented from the viewpoint of gravity field recovery and interpretation. Using a standard Airy/Heiskanen (A/H) isostatic hypothesis and a radially distributed compensation mechanism, two models of topographic/isostatic (t/i) potential harmonic coefficients are obtained up to degree and order 90. The CRUST-derived coefficients are compared with the spectrum of uncompensated topography, with the EGM96 observed gravity field, and with the t/i spectrum based on an A/H hypothesis with a constant compensation depth of 30 km. The signal degree variances of both CRUST models decrease quite smoothly towards degree 90, while the seven-layer model approaches the EGM96 spectrum for degrees 80–90. The significant deviation of the CRUST spectra from the A/H combined spectrum may prove of relevance to local and regional applications investigating the validity of current isostatic hypotheses.Acknowledgments. Sincere thanks go to Nikolaos Pavlis and three unknown referees for their thoughtful comments. Figure 1 has been produced using the mapping package m_map by R. Pawlowicz, which is a MATLAB toolbox that can be freely downloaded from http://www2.ocgy.ubc.ca/~rich/map.html     

The exterior Airy/Heiskanen topographic–isostatic gravity potential, anomaly and the effect of analytical continuation in Stokes' formula

This study deals with the external type of topographic–isostatic potential and gravity anomaly and its vertical derivatives, derived from the Airy/Heiskanen model for isostatic compensation. From the first and the second radial derivatives of the gravity anomaly the effect on the geoid is estimated for the downward continuation of gravity to sea level in the application of Stokes' formula. The major and regional effect is shown to be of order H ³ of the topography, and it is estimated to be negligible at sea level and modest for most mountains, but of the order of several metres for the highest and most extended mountain belts. Another, global, effect is of order H but much less significant Received: 3 October 1997 / Accepted: 30 June 1998     

Geoid determination with density hypotheses from isostatic models and geological information

 Geoid determination by Stokes's formula requires a complete knowledge of the topographical mass density distribution in order to perform gravity reductions to the geoid boundary. However, deeper masses are also of interest, in order to produce a smooth field of gravity anomalies which will improve results from interpolation procedures. Until now, in most cases a constant mass density has been considered, which is a very rough approximation of reality. The influence on the geoid height coming from different mass density hypotheses given by the isostatic models of Pratt/Hayford, Airy/Heiskanen and Vening Meinesz is studied. Apart from a constant mass density value, additional density information deduced from geological maps and thick sedimentary layers is considered. An overview of how mass density distributions act within Stokes's theory is given. The isostatic models are considered in spherical and planar approximation, as well as with constant and lateral variable mass density of the topographical and deeper masses. Numerical results in a test area in south-west Germany show that the differences in the geoid height due to different density hypotheses can reach a magnitude of more than 1 decimetre, which is not negligible in a precise geoid determination with centimetre accuracy. Received: 7 January 2002 / Accepted: 20 September 2002 M. Kuhn now at: Western Australian Centre for Geodesy, Curtin University of Technology, GPO Box U1987, Perth, WA 6845, Australia Acknowledgements. The author would gratefully thank Prof. Dr.-Ing. B. Heck, who was the supervisor of my PhD thesis, and the second examiner Prof. Dr.-Ing. K.H. Ilk, as well as all other colleagues for their support of this work. Particular thanks go to the Landesvermessungsamt Baden–Württemberg (Survey Department of Baden–Württemberg), Bureau Gravimetrique International (BGI, France) for providing the gravity data and the Geologisches Landesamt Baden–Württemberg (Geological Department of Baden–Württemberg) for providing data and maps of the sediment layers within the Rhine Valley. Grateful thanks goes to Prof. W.E. Featherstone and the reviewers Prof. S.D. Pagiatakis, Dr. U. Marti as well as an unknown reviewer for their helpful comments on this paper.     

CRUST 2.0模型在均衡位模型构建中的应用研究

基于均衡理论构制地球重力场模型中一个关键问题是扰动质量的确定,不同的地壳密度分布将得到不同的扰动质量。本文立足于面凝聚模型和Airy模型两种均衡补偿机制,研究了CRUST 2.0全球地壳模型在构制高分辨率地球重力场模型中的应用,推导了顾及地球物理信息的两种均衡重力场模型构建公式,分别讨论了CRUST 2.0模型和补偿深度在两种补偿机制中构建重力场模型的贡献。数值分析表明,地壳模型中Moho面深度的算术平均值22.97 km不是最佳补偿深度,而40 km相对最优,补偿深度的对模型超高阶部分影响较小;CRUST 2.0模型能够在361-1080频段内较好地改善原模型;相同补偿深度的面凝聚模型和Airy模型在不同频段的优劣性不一致。     

模拟大地水准面应用于GPS水准内插的研究

采用移去-恢复技术和多面函数法对已知GPS水准数据进行了内插处理。试验结果表明．利用实际地球物理资料得到的模拟大地水准面精度最高。     

Fourier geoid determination with irregular data

The fast Fourier transform (FFT) and, recently, the fast Hartley transform (FHT) have been extensively used by geodesists for efficient geoid determination. For this kind of efficiency, data must be given on a regular grid and, consequently, a pre-processing step of interpolation is required when only point measurements are available. This paper presents a way of computing a grid of geoid undulations N without explicitly gridding the data. The method is applicable to all FFT or FHT techniques of geoid or terrain effects determination, and it works with planar as well as spherical formulas. This method can be used not only for, e.g., computing a grid of undulations from irregular gravity anomalies g but it also lends itself to other applications, such as the gridding of gravity anomalies and, since the contribution of each data point is computed individually, the update of N- or g-grids as soon as new point measurements become available. In the case that there are grid cells which contain no measurements, the results of gravity interpolation or geoid estimation can be drastically improved by incorporating into the procedure a frequency-domain interpolating function. In addition to numerical results obtained using a few simple interpolating functions, the paper presents briefly the mathematical formulas for recovering missing grid values and for transforming values from one grid to another which might be rotated and/or scaled with respect to the first one. The geodetic problems where these techniques may find applications are pointed out throughout the paper.Presented at theIAG General Meeting, Beijing, P.R. China, Aug. 6–13, 1993     

Degree variances of the Earth's potential,topography and its isostatic compensation

A spherical harmonic expansion of the earth's gravitational potential and equivalent rock topography to degree and order 180 is described. The potential implied by the topography considered as uncompensated and with isostatic compensation has been computed. Good agreement with the observed potential field is found when the depth of compensation in the Airy theory is assumed to be 50 km. At the higher degrees the correlation coefficient between the potential expansion and the equivalent rock topography is about 0.5. The Lachapelle equations for the topographic isostatic potential were tested using 1^ox1^o equivalent rock topography. The degree variances agree at the lower degrees but at degree 36 the Lachapelle results using 5^o data underestimate the potential degree variances by about one-third.     

Geoid determination using adapted reference field, seismic Moho depths and variable density contrast

 The traditional remove-restore technique for geoid computation suffers from two main drawbacks. The first is the assumption of an isostatic hypothesis to compute the compensation masses. The second is the double consideration of the effect of the topographic–isostatic masses within the data window through removing the reference field and the terrain reduction process. To overcome the first disadvantage, the seismic Moho depths, representing, more or less, the actual compensating masses, have been used with variable density anomalies computed by employing the topographic–isostatic mass balance principle. In order to avoid the double consideration of the effect of the topographic–isostatic masses within the data window, the effect of these masses for the used fixed data window, in terms of potential coefficients, has been subtracted from the reference field, yielding an adapted reference field. This adapted reference field has been used for the remove–restore technique. The necessary harmonic analysis of the topographic–isostatic potential using seismic Moho depths with variable density anomalies is given. A wide comparison among geoids computed by the adapted reference field with both the Airy–Heiskanen isostatic model and seismic Moho depths with variable density anomaly and a geoid computed by the traditional remove–restore technique is made. The results show that using seismic Moho depths with variable density anomaly along with the adapted reference field gives the best relative geoid accuracy compared to the GPS/levelling geoid. Received: 3 October 2001 / Accepted: 20 September 2002 Correspondence to: H.A. Abd-Elmotaal     

A numerical investigation on height anomaly prediction in mountainous areas

This paper provides numerical examples for the prediction of height anomalies by the solution of Molodensky's boundary value problem. Computations are done within two areas in the Canadian Rockies. The data used are on a grid with various grid spacings from 100 m to 5 arc-minutes. Numerical results indicate that the Bouguer or the topographicisostatic gravity anomalies should be used in gravity interpolation. It is feasible to predict height anomalies in mountainous areas with an accuracy of 10 cm (1) if sufficiently dense data grids are used. After removing the systematic bias, the differences between the geoid undulations converted from height anomalies and those derived from GPS/levelling on 50 benchmarks is 12 cm (1) when the grid spacing is 1km, and 50 cm (1) when the grid spacing is 5. It is not necessary, in most cases, to require a grid spacing finer than 1 km, because the height anomaly changes only by 3 cm (1) when the grid spacing is increased from 100 m to 1000 m. Numerical results also indicate that, only the first two terms of the Molodensky series have to be evaluated in all but the extreme cases, since the contributions of the higher order terms are negligible compared to the objective accuracy.     

全球1°×1°海洋岩石圈有效弹性厚度模型

利用岩石圈挠曲均衡原理,联合海底地形模型和测高重力异常数据,采用滑动窗口导纳技术（moving window admittance technique,MWAT）计算了全球1°×1°海洋岩石圈有效弹性厚度模型。结果表明,海洋岩石圈有效弹性厚度总体较小,10 km以下区域约占70%,15 km以下约占85.4%,均值和标准差分别约为10 km和6.7 km。岩石圈有效弹性厚度20 km以上的区域主要位于海沟外隆地区,洋中脊岩石圈有效弹性厚度一般小于5 km;被动大陆边缘,如澳洲大陆南缘,岩石圈有效弹性厚度一般也较小;太平洋的海山密集分布地区,岩石圈有效弹性厚度一般为10~20 km。     

The effect on the geoid of lateral topographic density variations

Gravimetric geoid determination by Stokes formula requires that the effects of topographic masses be removed prior to Stokes integration. This step includes the direct topographic and the downward continuation (DWC) effects on gravity anomaly, and the computations yield the co-geoid height. By adding the effect of restoration of the topography, the indirect effect on the geoid, the geoid height is obtained. Unfortunately, the computations of all these topographic effects are hampered by the uncertainty of the density distribution of the topography. Usually the computations are limited to a constant topographic density, but recently the effects of lateral density variations have been studied for their direct and indirect effects on the geoid. It is emphasised that the DWC effect might also be significantly affected by a lateral density variation. However, instead of computing separate effects of lateral density variation for direct, DWC and indirect effects, it is shown in two independent ways that the total geoid effect due to the lateral density anomaly can be represented as a simple correction proportional to the lateral density anomaly and the elevation squared of the computation point. This simple formula stems from the fact that the significant long-wavelength contributions to the various topographic effects cancel in their sum. Assuming that the lateral density anomaly is within 20% of the standard topographic density, the derived formula implies that the total effect on the geoid is significant at the centimetre level for topographic elevations above 0.66 km. For elevations of 1000, 2000 and 5000 m the effect is within ± 2.2, ± 8.8 and ± 56.8 cm, respectively. For the elevation of Mt. Everest the effect is within ± 1.78 m.     

高原山区局部大地水准面精化方法研究

局部大地水准面精化的实质是精确计算出大地水准面的起伏变化情况。一般情况下,需要密度足够的重力数据,依重力异常密集计算大地水准面差距或高程异常。但是在大陆西部高原山区重力点密度是不够的,无法达到大地水准面精化的目的。本文从理论上证实了用地形和岩石密度数据进行局部大地水准面精化的可行性。     

Applications of downward-continuation in gravimetric geoid modeling: case studies in Western Canada

The objective of this study is to evaluate two approaches, which use different representations of the Earth’s gravity field for downward continuation (DC), for determining Helmert gravity anomalies on the geoid. The accuracy of these anomalies is validated by 1) analyzing conformity of the two approaches; and 2) converting them to geoid heights and comparing the resulting values to GPS-leveling data. The first approach (A) consists of evaluating Helmert anomalies at the topography and downward-continuing them to the geoid. The second approach (B) downward-continues refined Bouguer anomalies to the geoid and transforms them to Helmert anomalies by adding the condensed topographical effect. Approach A is sensitive to the DC because of the roughness of the Helmert gravity field. The DC effect on the geoid can reach up to 2 m in Western Canada when the Stokes kernel is used to convert gravity anomalies to geoid heights. Furthermore, Poisson’s equation for DC provides better numerical results than Moritz’s equation when the resulting geoid models are validated against the GPS-leveling. On the contrary, approach B is significantly less sensitive to the DC because of the smoothness of the refined Bouguer gravity field. In this case, the DC (Poisson’s and Moritz’s) contributes only at the decimeter level to the geoid model in Western Canada. The maximum difference between the geoid models from approaches A and B is about 5 cm in the region of interest. The differences may result from errors in the DC such as numerical instability. The standard deviations of the h−H−N for both approaches are about 8 cm at the 664 GPS-leveling validation stations in Western Canada.     

Gravity field convolutions without windowing and edge effects

A new set of formulas has been developed for the computation of geoid undulations and terrain corrections by FFT when the input gravity anomalies and heights are mean gridded values. The effects of the analytical and the discrete spectra of kernel functions and that of zero-padding on the computation of geoid undulations and terrain corrections are studied in detail.Numerical examples show that the discrete spectrum is superior to the analytically-defined one. By using the discrete spectrum and 100% zero-padding, the RMS differences are 0.000 m for the FFT geoid undulations and 0.200 to 0.000 mGal for the FFT terrain corrections compared with results obtained by numerical integration.     

Density interface topography recovered by inversion of satellite gravity gradiometry observations

A radial integration of spherical mass elements (i.e. tesseroids) is presented for evaluating the six components of the second-order gravity gradient (i.e. second derivatives of the Newtonian mass integral for the gravitational potential) created by an uneven spherical topography consisting of juxtaposed vertical prisms. The method uses Legendre polynomial series and takes elastic compensation of the topography by the Earth’s surface into account. The speed of computation of the polynomial series increases logically with the observing altitude from the source of anomaly. Such a forward modelling can be easily applied for reduction of observed gravity gradient anomalies by the effects of any spherical interface of density. An iterative least-squares inversion of measured gravity gradient coefficients is also proposed to estimate a regional set of juxtaposed topographic heights. Several tests of recovery have been made by considering simulated gradients created by idealistic conical and irregular Great Meteor seamount topographies, and for varying satellite altitudes and testing different levels of uncertainty. In the case of gravity gradients measured at a GOCE-type altitude of \(\sim \)300 km, the search converges down to a stable but smooth topography after 10–15 iterations, while the final root-mean-square error is \(\sim \)100 m that represents only 2 % of the seamount amplitude. This recovery error decreases with the altitude of the gravity gradient observations by revealing more topographic details in the region of survey.