地球及其形状非对称性的重力学研究

依据卫星重力、大地水准面观测结果,对表示地球形状的大地水准面异常进行了不同阶次的球谐函数计算和分析,发现2~6阶大地水准面异常表现出了地球双重非对称的基本形状,而高于6阶的大地水准面异常只表示了地球形态的局部变化特征。另外,采用阻尼最小二乘方法,利用全球地震层析成像资料和大地水准面异常资料联合反演了三维全球地幔密度异常。结果显示了地幔密度异常不仅在横向上,而且在纵向上也存在着明显的不均匀性。对比分析低阶大地水准面异常和地幔密度异常结果表明,地球形状非对称性主要是由下地幔的物质密度不均匀引起的。     

全球地幔密度异常及其构造意义

利用扣除地形、莫霍面和核幔边界起伏影响的中长波大地水准面异常和全球地震层析成像资料,采用阻尼最小二乘方法反演计算了全球地幔6个不同层面上的密度异常分布。分析了全球密度异常与板块构造的关系,探讨了全球密度异常分布对板块运动的作用。全球地幔密度异常结果表明存在两个主要的密度异常中心：一个位于东经80°,北纬0°;另一个位于东经240°,北纬10°附近。     

利用谱分析法估计GRACE Follow-On地球重力场的空间分辨率

利用轨道扰动引力谱和大地水准面累计误差谱分析的方法估计未来GRACE(gravity recovery and climateexperiment)Follow-On卫星反演地球重力场的空间分辨率。基于GRACE Follow-On卫星的轨道特性,计算其在高空所受到的径向扰动引力,并根据谱特性及星载加速度计的测量噪声水平分析该卫星能反演重力场的阶数。利用EGM96重力场模型分别计算200 km和250 km轨道高度处的扰动引力谱。分析其特性表明:在两个轨道高度处分别能反演281阶和242阶的地球重力场模型。给出大地水准面累计误差谱模型,并计算200 km和250 km轨道高度处大地水准面累计误差谱。分析其谱特性表明:在两个轨道高度处分别能反演至286阶和228阶的地球重力场模型。     

扰动位的综合确定

利用地球重力场任意一种有关信息都可以描述地球重力场一定的情况。根据卫星轨道摄动观测求定的引力位球谐系数只能表示地球重力场的长波分量。大地水准面起伏是地球扰动引力场越来越丰富的有用信息，但目前用其计算的引力位系数也只是在低阶较准确。重力异常、垂线偏差、单层密度和纯重力异常都利于求定高阶位系数，其中与大地水准面起伏有关的量，如纯重力异常和单层密度，用它们计算位系数等于联合应用大地水准面和重力异常，故用其计算的位系数在低阶次精度也较好。重力异常垂直梯度是描述扰动引力场细部最有利的信息。本文给出利用各种类型观测资料计算位系数精度估计式，提出综合利用各种资料求定位系数依资料类型的谱特性赋权的方法。     

Estimation of the Resolution of Earth＇s Gravity Field for GRACE Follow-On Using the Spectrum Method

利用轨道扰动引力谱和大地水准面累计误差谱分析的方法估计未来GRACE（gravity recovery and climate experimenl）Follow—On卫星反演地球重力场的空间分辨率。基于GRACEFollow—On卫星的轨道特性,计算其在高空所受到的径向扰动引力,并根据谱特性及星载加速度计的测量噪声水平分析该卫星能反演重力场的阶数。利用EGM96重力场模型分别计算200km和250km轨道高度处的扰动引力谱。分析其特性表明：在两个轨道高度处分别能反演281阶和242阶的地球重力场模型。给出大地水准面累计误差谱模型,并计算200km和250km轨道高度处大地水准面累计误差谱。分析其谱特性表明：在两个轨道高度处分别能反演至286阶和228阶的地球重力场模型。     

下地幔侧向不均匀及绝对板块运动对大地水准面的影响

本文提出并研究了板块上大地水准面受下地幔侧向不均匀性及板块绝对运动的影响而随时间变化的问题。利用已有的下地幔侧向不均匀性模型、板块绝对运动速度模型进行了计算,并分核幔边界不存在起伏和存在起伏两种情况进行了讨论。结果表明,在相当长时间内,其垂直变化绝对幅度可达 200米。确定古大地水准面应考虑这一变化。     

基于遗传Elman神经网络进行矿区GPS高程拟合

目前,城市、平原地区的似大地水准面建立精度已经达到厘米级,但在矿区进行高程拟合时,由于地面高低起伏没有规则,其似大地水准面的拟合精度并不理想。针对此问题,本文提出利用遗传算法优化Elman神经网络的方法精化似大地水准面,采用移去-恢复法对残差进行建模,使用EGM 2008地球重力场模型和地形起伏信息来精化求解似大地水准面和参考椭球面之间的高程异常,同时着重分析了地球重力场模型以及地形变化信息对高程异常求解的重要性,并使用某矿区实测数据(GPS、水准)对所提方法进行验证,实验结果表明:文中所提方法的精度要优于二次曲面拟合模型和单一Elman模型,其外符合精度达到了1.14 cm,可以代替四等水准测量。     

广东地区地球重力场模型的比较与选择

为了优化广东地区的似大地水准面,需要选择最适合该地区的地球重力场模型。本文对该地区6个具有代表性的地球重力场模型及简单组合模型,从高程异常、重力异常和垂线偏差3个方面的精度进行比较分析。结果表明,各模型所代表的似大地水准面与该地区的大地水准面之间均存在着一定的系统偏差,而EGM2008/EIGEN-6C简单组合模型精度则优于其他模型,此模型可作为该地区似大地水准面优化的参考重力场首选模型。     

基于地球物理信息构建全球重力异常模型方法的探讨

物理大地测量的主要任务是研究地球形状及其外部重力场。而无论是通过求解司托克斯反问题来确定大地水准面及其重力位,还是通过莫洛琴斯基问题求解似大地水准面和实际地球重力位的过程中都需要用到一个重要的量：空间重力异常。在实际计算中往往需要格网化的重力异常,因此需要建立重力异常模型。本文主要在地球物理信息的基础上,重点论述了两种推求重力异常的方法,并探讨了利用全球重力场信息源诊断与融合技术推估重力异常模型的方法。     

月球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。在各月海盆地存在着中央较薄、四周逐渐增厚的趋势。     

A new isostatic model of the lithosphere and gravity field

Based on the analysis of various factors controlling isostatic gravity anomalies and geoid undulations, it is concluded that it is essential to model the lithospheric density structure as accurately as possible. Otherwise, if computed in the classical way (i.e. based on the surface topography and the simple Airy compensation scheme), isostatic anomalies mostly reflect differences of the real lithosphere structure from the simplified compensation model, and not necessarily the deviations from isostatic equilibrium. Starting with global gravity, topography and crustal density models, isostatic gravity anomalies and geoid undulations have been determined. The initial crust and upper-mantle density structure has been corrected in a least squares adjustment using gravity. To model the long-wavelength (>2000 km) features in the gravity field, the isostatic condition (i.e. equal mass for all columns above the compensation level) is applied in the adjustment to uncover the signals from the deep-Earth interior, including dynamic deformations of the Earths surface. The isostatic gravity anomalies and geoid undulations, rather than the observed fields, then represent the signals from mantle convection and deep density inhomogeneities including remnants of subducted slabs. The long-wavelength non-isostatic (i.e. the dynamic) topography was estimated to range from –0.4 to 0.5 km. For shorter wavelengths (<2000 km), the isostatic condition is not applied in the adjustment in order to obtain the non-isostatic topography due to regional deviations from classical Airy isostasy. The maximum deviations from Airy isostasy (–1.5 to 1 km) occur at currently active plate boundaries. As another result, a new global model of the lithosphere density distribution is generated. The most pronounced negative density anomalies in the upper mantle are found near large plume provinces, such as Iceland and East Africa, and in the vicinity of the mid-ocean ridge axes. Positive density anomalies in the upper mantle under the continents are not correlated with the cold and thick lithosphere of cratons, indicating a compensation mechanism due to thermal and compositional density.     

Inverse Vening Meinesz formula and deflection-geoid formula: applications to the predictions of gravity and geoid over the South China Sea

Using the spherical harmonic representations of the earth's disturbing potential and its functionals, we derive the inverse Vening Meinesz formula, which converts deflection of the vertical to gravity anomaly using the gradient of the H function. The deflection-geoid formula is also derived that converts deflection to geoidal undulation using the gradient of the C function. The two formulae are implemented by the 1D FFT and the 2D FFT methods. The innermost zone effect is derived. The inverse Vening Meinesz formula is employed to compute gravity anomalies and geoidal undulations over the South China Sea using deflections from Seasat, Geosat, ERS-1 and TOPEX//POSEIDON satellite altimetry. The 1D FFT yields the best result of 9.9-mgal rms difference with the shipborne gravity anomalies. Using the simulated deflections from EGM96, the deflection-geoid formula yields a 4-cm rms difference with the EGM96-generated geoid. The predicted gravity anomalies and geoidal undulations can be used to study the tectonic structure and the ocean circulations of the South China Sea. Received: 7 April 1997 / Accepted: 7 January 1998     

利用CHAMP科学轨道数据和星载加速度计数据反演地球重力场

基于卫星动力学理论,采用德国地球科学中心GFZ提供的CHAMP精密轨道数据和星载加速度计数据,反演了36阶地球重力场模型CDS01S。用不同模型之间的位系数差比较模型CDS01S、EIGEN3P、EIGEN1S及EGM96,表明CDS01S模型的位系数最接近于EIGEN3P;比较上述几种模型的位系数精度,表明CDS01S模型的位系数精度高于EGM96;用CDS01S和GGM01C的前30阶位系数分别计算全球2°×2°网格的大地水准面起伏,两者之间的标准偏差为4.7 cm。     

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.     

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

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

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.     

局部重力异常向上延拓的实用算法

针对局部重力异常向上延拓计算复杂、耗时长的问题,该文基于泊松积分离散化的基本原理,提出一种快速的局部格网重力异常向上延拓的实用算法;并结合中国东北和青藏高原地区大地水准面的重力异常格网数据,采用该延拓方法分别计算了空中10、50、100km处的重力异常,将其与等高度的EIGEN-6C4模型结果对比分析。实验结果表明:在顾及边界效应影响的情况下,相对于EIGEN-6C4模型,中国东北和青藏高原地区重力异常向上延拓的最大均方根误差分别优于1.5和3.5mGal;在保证精度可用的前提下,计算效率可以有大幅度提高,证明了该方法解算局部重力异常向上延拓的适用性。     

利用重力异常数据构建重力梯度场的方法比较

为实现大范围、高精度基准重力梯度数据库的构建,考虑到重力梯度场对地形质量的敏感效应,一般利用恒密度数字高程模型来求取重力梯度值,从而忽略了地形密度变化以及水准面以下密度异常对重力梯度的影响。根据重力位理论中求解边值问题的数值应用方法,直接利用重力异常数据求取重力梯度场,弥补了密度变化和密度异常在重力梯度上的反映。根据模型算例和实测重力异常数据求取了剖面重力梯度值,结果表明,限于重力数据空间分辨率的影响,利用重力异常数据可恢复中长波段重力梯度场。该方法与地形数据求取重力梯度和卫星重力梯度测量等方法技术相结合,对重力梯度数据库的建设具有实际应用价值。