重力卫星五年运行对求定地球重力场模型的进展和展望

对近年升空的CHAMP和GRACE和将于2007年升空的GOCE卫星在测定地球重力场方面的技术特点和初步成果进行了回顾、比较和评估。并对它们今后在静态和动态的地球重力场构模方面可能的进展作一展望。现在只用一颗重力卫星的轨道摄动数据,就可以以前所未有的可靠性和精确性来求定地球重力场的长波和部分中波。如CHAMP重力卫星的33个月数据所求定的地球重力场模型,相对于曾利用多颗卫星资料所推算的GRlM5-S1重力场模型,在长波方面的精度和可靠性都有很大改善。而GRACE重力卫星的110天数据所导出成果的空间分辨率,又优于CHAMP的33个月的数据成果。GRACE卫星还有一个重要任务,就是测定重力场非潮汐的短期性或准实时的变化。还介绍了新发表的一个联合地球重力场模型EIGEN-CG03C,360完全阶次,分辨率约30′。CG03C同CHAMP／GRACE以前的重力场模型比较,在400km波长的精度方面改善了一个量级,大地水准面的精度改善了3cm,重力异常的精度改善了0．4mgal。     

基于卫-卫跟踪观测技术利用能量守恒法恢复地球重力场的数值模拟研究

根据卫－卫跟踪观测技术的测量原理,基于能量守恒法建立了一种新的双星相互跟踪和三星相互跟踪的卫星观测方程. 通过数值模拟,采用预处理共轭梯度法恢复120阶地球重力场. 模拟结果表明：第一,双星相互跟踪恢复地球重力场的精度和美国喷气动力实验室公布的EIGEN GRACE02S的结果相符合;第二,三星相互跟踪恢复地球重力场的精度较双星提高约2倍.     

全球水储量变化的GRACE卫星检测

利用GRACE月尺度变化的地球重力场反演了全球水储量变化,并与陆地水文资料、卫星测高资料及海洋模式得到的结果进行了比对.通过对SOURE台站重力变化的陆地水储量变化计算结果和GRACE重力场系数截断为15阶得到的结果比较,发现两者比较接近,且年周期变化特征明显.对于亚马逊流域,当重力场系数截断为15阶且平滑半径使用106 m时,GRACE反演的区域平均水储量厚度的周年变化振幅为15.6×10-2m,小于使用平滑半径为4×105m的23.7×10-2m.在研究长江流域时,本文对水文资料做球谐系数展开,并与GRACE数据做同样的截断和平滑处理,结果发现GRACE反演的水厚度变化与水文资料结果基本上符合.对于纬度±66°之间的海洋区域,GRACE反演的海水质量变化接近于结合卫星测高和海洋模式得到的结果,但对于2°×2°网格,则在一些区域差异明显,最大超过了0.2 m,中误差为3.8×10-2m.可见,当前GRACE卫星时变重力场只能确定出上千公里及以上尺度区域的水储量变化.     

基于多源数据分析维多利亚湖流域水储量变化

本文利用CSR发布的GRACE RL06时变重力场模型,结合两种水文模式、卫星测高、降雨和蒸散等多源数据,从多个角度综合系统地分析维多利亚湖流域2003-01-2017-06的陆地水储量变化.比较了正向建模方法和单一尺度因子对泄漏误差的改正效果,经对比采用正向建模方法在此流域效果更好.基于多源数据得出以下三点与此前研究...     

GRACE时变信号恢复方法对比分析

针对利用GRACE重力数据反演地表质量变化过程中产生的时变信号削弱和泄露现象,分别采用尺度因子法和迭代恢复法对信号进行恢复,并选用CSR Mascon产品从时间序列变化和空间分布两个角度对两种恢复方法的效果进行对比分析。结果表明:1)尺度因子法虽实现上较为简单,但其计算结果完全依赖于所选的先验模型,当所选模型的可靠性较差时,会影响到GRACE反演结果的可靠性。2)迭代恢复法的恢复过程仅受GRACE原始观测值的约束,能较好地恢复研究区域的长期趋势信号和周期信号,但在部分区域可能存在信号过量恢复的现象。建议采用迭代恢复法进行GRACE时变信号恢复。     

GPS天线相位中心校正对GRACE卫星精密定轨的影响

使用GRACE卫星星载GPS观测数据,研究地面获取的先验PCV模型和利用残差法估计的在轨PCV模型对低轨卫星精密定轨的影响,并采用GFZ精密轨道对比和SLR检核手段对其进行评估。结果表明,使用先验PCV模型会降低GRACE卫星定轨精度,相反利用在轨PCV模型可以提高定轨精度,提升数量级可达mm级。     

GRACE与SLR J_2项时间序列的小波分析

利用小波功率谱分析和离散小波变换,对CNES/GRGS、CSR、GFZ等研究机构发布的GRACE时变重力位模型的J2项时间序列进行了周期分析和特征提取,并与利用SLR跟踪数据获得的J2项时间序列的分析结果进行了比较,结果显示:SLR跟踪数据获得的J2项相较GRACEJ2项具有更明显的季节性变化特征.GRACE 时变重力位模型的J2项中混淆进了161天的周期信号,是由海洋潮汐模型中存在的太阳半日分潮S2项的误差引起的.     

On the use of iterative re-weighting least-squares and outlier detection for empirically modelling rates of vertical displacement

The proper identification and removal of outliers in the combination of rates of vertical displacements derived from GPS, tide gauges/satellite altimetry, and GRACE observations is presented. Outlier detection is a necessary pre-screening procedure in order to ensure reliable estimates of stochastic properties of the observations in the combined least-squares adjustment (via rescaling of covariance matrices) and to ensure that the final vertical motion model is not corrupted and/or distorted by erroneous data. Results from this study indicate that typical data snooping methods are inadequate in dealing with these heterogeneous data sets and their stochastic properties. Using simulated vertical displacement rates, it is demonstrated that a large variety of outliers (random scattered and adjacent, as well as jointly influential) can be dealt with if an iterative re-weighting least-squares adjustment is combined with a robust median estimator. Moreover, robust estimators are efficient in areas weakly constrained by the data, where even high quality observations may appear to be erroneous if their estimates are largely influenced by outliers. Four combined models for the vertical motion in the region of the Great Lakes are presented. The computed vertical displacements vary between  − 2 mm/year (subsidence) along the southern shores and 3 mm/year (uplift) along the northern shores. The derived models provide reliable empirical constraints and error bounds for postglacial rebound models in the region.     

A study on the reference frame consistency in recent Earth gravitational models

All gravity field functionals obtained from an Earth gravitational model (EGM) depend on the underlying terrestrial reference frame (TRF), with respect to which the EGM’s spherical harmonic coefficients refer to. In order to maintain a coherent framework for the comparison of current and future EGMs, it is thus important to investigate the consistency of their inherent TRFs, especially when their use is intended for high precision studies. Following the methodology described in an earlier paper by Kleusberg (1980), the similarity transformation parameters between the associated reference frames for several EGMs (including the most recent CHAMP/GRACE models at the time of writing this paper) are estimated in the present study. Specifically, the differences between the spherical harmonic coefficients for various pairs of EGMs are parameterized through a 3D-similarity spatial transformation model that relates their underlying TRFs. From the least-squares adjustment of such a parametric model, the origin, orientation and scale stability between the EGMs’ reference frames can be identified by estimating their corresponding translation, rotation and scale factor parameters. Various aspects of the estimation procedure and its results are highlighted in the paper, including data weighting schemes, the sensitivity of the results with respect to the selected harmonic spectral band, the correlation structure and precision level of the estimated transformation parameters, the effect of the estimated differences of the EGMs’ reference frames on their height anomaly signal, and the overall feasibility of Kleusberg’s formulae for the assessment of TRF inconsistencies among global geopotential models.     

基于能量守恒的星间距离变率与地球重力场频谱关系的建立与分析(英文)

The spectral relationship between range-rate measurements and the gravity potential for low-low satellite-to-satellite tracking mission was established based on the energy conservation theory. Then the performances of satellite separation, the orbital altitude, and the accuracy of range-rate measurements in recovering the earth’s gravity field were simulated and analyzed by this method. Finally, the cumulative geoid errors of the reference mode were obtained by using the configuration parameters of the GRAC...