多类地球重力场模型的高程异常精度比较

分析了最新的多类地球重力场模型解算GNSS点高程异常的精度,利用中国某测区内的GNSS/水准数据与地球重力场模型解算值相比较,选择出较优精度的重力场模型,在测区范围内解算精度达到21.09 cm,消除系统偏差后的精度最高为3.94 cm.结果表明在本区域地球重力场模型确定的重力似大地水准面与GNSS/水准的似大地水准面存在系统偏差,消除系统偏差后的结果精度有所提升,可在铁路工程测量中发挥一定的作用.     

利用LP多普勒数据解算月球重力场模型的分析

利用GSFC/NASA提供的GEODYNII/SOLVE软件处理月球探测器LP(Lunar Prospector)扩展任务阶段最后3个月的测速和测距跟踪数据,解算月球重力场模型,比较分析全部位系数解算和部分中高阶次位系数解算2种方法,并对2种方法解算得到的模型进行阶方差分析.同时通过轨道残差、重复轨道误差以及模型大地水准面等几个方面对基于部分中高阶次住系数解算的有效性和得到的模型进行精度评估.结果表明基于部分中高阶次位系数解算的方法进行恢复月球重力场是可行的.该方法将用于"嫦娥工程"中月球重力场的解析.     

西藏墨脱地区工程精化似大地水准面的研究

采用地球重力场模型 ,考虑区域重力异常和地形改正 ,建立了西藏墨脱地区 1′× 1′似大地水准面模型。实际检测结果表明 ,该似大地水准面模型相对精度达到了± 0 .0 2 5m ,可代替几何水准测量 ,满足该地区工程建设的要求。     

利用GOCE模拟观测反演重力场的Torus法

在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。     

联合GRACE/GOCE重力场模型和GPS/水准数据确定我国85高程基准重力位

GRACE、GOCE卫星重力计划的实施,对确定高精度重力场模型具有重要贡献。联合GRACE、GOCE卫星数据建立的重力场模型和我国均匀分布的649个GPS/水准数据可以确定我国高程基准重力位,但我国高程基准对应的参考面为似大地水准面,是非等位面,将似大地水准面转化为大地水准面后确定的大地水准面重力位为62 636 854.395 3m~2s~(-2),为提高高阶项对确定大地水准面的贡献,利用高分辨率重力场模型EGM2008扩展GRACE/GOCE模型至2190阶,同时将重力场模型和GPS/水准数据统一到同一参考框架和潮汐系统,最后利用扩展后的模型确定的我国大地水准面重力位为62 636 852.751 8m~2s~(-2)。其中组合模型TIM_R4+EGM2008确定的我国85高程基准重力位值62 636 852.704 5m~2s~(-2)精度最高。重力场模型截断误差对确定我国大地水准面的影响约16cm,潮汐系统影响约4~6cm。     

由CHAMP星载GPS相位双差数据解算地球引力场模型

利用7d的CHAMP星载GPS相位观测数据和48个IGS跟踪站的观测数据，构造星地双差相位观测量，进行GPS数据预处理；利用Cowell Ⅱ数值法进行轨道积分和分块Bayes最小二乘参数估计，解算了地球引力场位系数。该模型与EGM96相比（70阶次），大地水准面起伏差异最大为2．872m，差弄精度为0．522m，平均差异为-0．003m，这说明本文解算的地球重力场模型与EGM96没有系统性差异。     

测高-重力边值问题的局部剪裁解

给出了基于最小方差原理的测高-重力超定边值问题解的局部建模剪裁解法,并采用由卫星测高资料计算的6′×6′海洋重力异常和30′×30′大地水准面数据,计算了展开到1 800阶次的中国南海及邻近海域的测高-重力试验模型.与其他模型比较表明,这种由超定剪裁法解得到的局部模型能够较好地改善重力场模型中波长部分的精度.     

积分法恢复地球重力场模拟分析

本文用积分法模拟恢复了地球重力场，同真实重力场进行比较，分析了结果。结果表明，本文使用模型及程序算法正确，可以用于恢复重力场，精度在10^-12的量级上，对应的大地水准面阶误差小于0．2mm，大地水准面累积误差30阶处小于0．5mm。     

利用谱分析法估计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阶的地球重力场模型。     

全国及部分省市地区高精度高分辨率似大地水准面的研究和实施

结合我国重力和地形资料及国内外较优的重力场模型,研制适合我国重力场特征的360阶重力场模型WDM94;建立中国新一代分米级似大地水准面CQG2000,包括建立新的以GPS/水准为基础的高程异常控制网、利用海洋卫星测高数据计算海洋大地水准面、陆地重力似大地水准面的研制及陆海似大地水准面的拼接等;研制江苏省、海南省、深圳市、大连市、南京市及"南水北调"西线工程具有厘米级精度的局域似大地水准面模型;结合GPS技术和高精度似大地水准面模型,研制GPS测图软硬件一体化系统.本研究项目获得2004年度国家科技进步二等奖.     

利用卫星测高数据反演海洋重力异常研究

全面研究了利用卫得测高数据反演海洋重力异常3种主要方法（即Stokes数据解析反解以及逆Vening-Meinesz公式）的技术特点，建立了3种算法的数学模型及其谱计算式，在以1440阶次位模型定义的标准场中完成了3种算法的数值比较和内部检核，通过仿真试验实现了3种算法的可靠性和稳定性检验，最后，本文利用卫得测高实测对南中国海地区的海洋重力异常进行了实际反演，并将反演结果同船测数据进行了比较。     

联合使用位模型和地形信息的陆区航空重力向下延拓方法

为了规避传统逆Poisson积分向下延拓解算过程的不适定性问题,借鉴导航定位中的"差分"概念,利用超高阶位模型直接计算海域航空重力测量向下延拓改正数的方法。本文在此基础上提出联合使用重力位模型和地形高数据,计算陆部航空重力向下延拓总改正数的改进方案,以飞行高度面与地面对应点的位模型差分信息表征总改正数的中长波分量,以相对应的局部地形改正差分修正量表征总改正数的中高频成分,从而实现航空重力数据点对点向地面的全频段延拓。在地形变化不同区域,联合使用EGM2008位模型、地面实测重力和高分辨率高程数据进行了实际数值计算和精度评估,验证了该方法的有效性。     

Improving gravity anomalies over China marginal sea from retracked Geosat and ERS-1 data

The quality of altimeter data and ocean tide model is critical to the recovery of coastal gravity anomalies. In this contribution, three retracking methods (threshold, improved threshold and Beta-5) are investigated with the aim of improving the altimeter data over a shallow water area. Comparison indicates that the improved threshold is the best retracking method over China Sea. Two ocean tide models, NAO99b and CSR4.0, are analyzed. Results show that different tide models used in the processing of altimeter data may result in differences more than 10 mGal in recovered coastal gravity anomalies. Also, NAO99b is more suitable than CSR4.0 over the shallow water area of China Sea. Finally, gravity anomalies over China Sea are calculated from retracked Geosat/GM and ERS-1/GM data by least squares collocation. Comparison with shipborne gravimetry data demonstrates that gravity anomalies from retracked data are significantly superior to those from non-retracked data. Our results have the same order as the other two altimeter-derived gravity models: Sandwell&Smith(V16) and DNSC08.     

Long-wavelength global gravity field models: GRIM4-S4, GRIM4-C4

Summary.  GFZ Potsdam and GRGS Toulouse/Grasse jointly developed a new pair of global models of the Earth's gravity field to satisfy the requirements of the recent and future geodetic and altimeter satellite missions. A precise gravity model is a prerequisite for precise satellite orbit restitution, tracking station positioning and altimeter data reduction. According to different applications envisaged, the new model exists in two parallel versions: the first one being derived exclusively from satellite tracking data acquired on 34 satellites, the second one further incorporating satellite altimeter data over the oceans and terrestrial gravity data. The most recent “satellite-only” gravity model is labelled GRIM4-S4 and the “combined” gravity model GRIM4-C4. The models are solutions in spherical harmonics and have a resolution up to degree and order 60 plus a few resonance terms in the case of GRIM4-S4, and up to degree/order 72 in the case of GRIM4-C4, corresponding to a spatial resolution of 555 km at the Earth's surface. The gravitational coefficients were estimated in a rigorous least squares adjustment simultaneously with ocean tidal terms and tracking station position parameters, so that each gravity model is associated with a consistent ocean tide model and a terrestrial reference frame built up by over 300 optical, laser and Doppler tracking stations. Comprehensive quality tests with external data and models, and test arc computations over a wide range of satellites have demonstrated the state-of-the-art capabilities of both solutions in long-wavelength geoid representation and in precise orbit computation. Received 1 February 1996; Accepted 17 July 1996     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Regional geopotential model improvement for the Iranian geoid determination

Spherical harmonic expansions of the geopotential are frequently used for modelling the earth’s gravity field. Degree and order of recently available models go up to 360, corresponding to a resolution of about50 km. Thus, the high degree potential coefficients can be verified nowadays even by locally distributed sets of terrestrial gravity anomalies. These verifications are important when combining the short wavelength model impact, e.g. for regional geoid determinations by means of collocation solutions. A method based on integral formulae is presented, enabling the improvement of geopotential models with respect to non-global distributed gravity anomalies. To illustrate the foregoing, geoid computations are carried out for the area of Iran, introducing theGPM2 geopotential model in combination with available regional gravity data. The accuracy of the geoid determination is estimated from a comparison with Doppler and levelling data to ±1.4m.     

Gravity-field improvement in the Mediterranean Sea by estimating the bottom topography using collocation

The contribution of bathymetry to the prediction of quantities related to the gravity field (e.g., gravity anomalies, geoid heights) is discussed in an extended test area of the central Mediterranean Sea. Sea gravity anomalies and a priori statistical characteristics of depths are used in a least-squares collocation procedure in order to produce new depths, giving a better smoothing of the gravity field when using a remove-restore procedure. The effect of the bottom topography on gravity-field modeling is studied using both the original and the new depths through a residual terrain modeling reduction. The numerical tests show a considerable smoothing of the sea gravity anomalies and the available altimeter heights when the new depth information is taken into account according to the covariance analysis performed. Moreover, geoid heights are computed by combining the sea gravity anomalies either with the original depths or with the new ones, using as a reference surface the OSU91A geopotential model. Comparing the computed geoid heights with adjusted altimeter sea-surface heights (SSHs), better results are obtained when subtracting the attraction of the new depth information. Similar results are obtained when predicting gravity anomalies from altimeter SSHs where the terrain effect on altimetry is based on the new bottom topography. Received: 10 September 1996 / Accepted: 4 August 1997     

用T/P测高数据反演海域重力异常

以反解 Stokes公式为数学模型 ,应用由 T/ P测高数据计算的大地水准面高反演了海域平均重力异常 ,并与船测平均重力异常和 OSU91A位模型计算的平均重力异常进行了比对分析 ,得出了一些有益的结论。     

重力三层点质量模型的构造与分析

点质量模型理论是研究区域重力场的一个非常重要的方法,本文简要介绍了点质量模型逼近区域重力场的原理,计算分析了构造点质量模型过程中系数矩阵元素的特性。以32~34N和103~105E为计算中心区域,利用EGM2008的720阶次的位系数计算出的重力异常作为观测数据,在36阶次位系数模型的基础上,构造了四层点质量组分频段从低到高来逼近该区域重力场。数值试验的结果表明：三层点质量模型效果较好,点质量模型计算的扰动重力在径向上的截断误差优于2 .