Precise GRACE baseline determination using GPS

Precision relative navigation is an essential aspect of spacecraft formation flying missions, both from an operational and a scientific point of view. When using GPS as a relative distance sensor, dual-frequency receivers are required for high accuracy at large inter-satellite separations. This allows for a correction of the relative ionospheric path delay and enables double difference integer ambiguity resolution. Although kinematic relative positioning techniques demonstrate promising results for hardware-in-the-loop simulations, they were found to lack an adequate robustness in real-world applications. To overcome this limitation, an extended Kalman Filter modeling the relative spacecraft dynamics has been developed. The filter processes single difference GPS pseudorange and carrier phase observations to estimate the relative position and velocity along with empirical accelerations and carrier phase ambiguities. In parallel, double difference carrier phase ambiguities are resolved on both frequencies using the least square ambiguity decorrelation adjustment (LAMBDA) method in order to fully exploit the inherent measurement accuracy. The combination of reduced dynamic filtering with the LAMBDA method results in smooth relative position estimates as well as fast and reliable ambiguity resolution. The proposed method has been validated with data from the gravity recovery and climate experiment (GRACE) mission. For an 11-day data arc, the resulting solution matches the GRACE K-Band Ranging System measurements with an accuracy of 1 mm, whereby 83% of the double difference ambiguities are resolved.     

卫星重力测量数据处理软件系统的设计

介绍了自主开发的卫星重力测量数据处理软件GRASTAR,给出了该软件的整体设计框架和功能。该软件主要采用动力学法实现,应用CHAMP卫星和GRACE卫星的观测数据反演地球重力场模型。利用模拟方法验证了该软件的正确性,并利用GRASTAR处理了126 d的CHAMP卫星数据,解算出直到40阶次重力场模型的初步结果。     

Short-arc analysis of intersatellite tracking data in a gravity mapping mission

 A technique for the analysis of low–low intersatellite range-rate data in a gravity mapping mission is explored. The technique is based on standard tracking data analysis for orbit determination but uses a spherical coordinate representation of the 12 epoch state parameters describing the baseline between the two satellites. This representation of the state parameters is exploited to allow the intersatellite range-rate analysis to benefit from information provided by other tracking data types without large simultaneous multiple-data-type solutions. The technique appears especially valuable for estimating gravity from short arcs (e.g. less than 15 minutes) of data. Gravity recovery simulations which use short arcs are compared with those using arcs a day in length. For a high-inclination orbit, the short-arc analysis recovers low-order gravity coefficients remarkably well, although higher-order terms, especially sectorial terms, are less accurate. Simulations suggest that either long or short arcs of the Gravity Recovery and Climate Experiment (GRACE) data are likely to improve parts of the geopotential spectrum by orders of magnitude. Received: 26 June 2001 / Accepted: 21 January 2002     

Comparison of present SST gravity field models

IntroductionSince 2000 , with the launch of CHAMP satellite,several series of high-accuracy and high-resolutionstatic Earth’s gravityfield models have been createdbased on aboundent SST data. With these models ,the research in solid geophysics ,oceangraphy,andgeodesy can be promoted greatly[1].In this paper ,the SSTgravity models’accura-cyin various frequently domainis studied.First-ly,the difference among these models is compu-ted and compared,and then their accuracyis an-alyzed. Fina…     

利用GRACE重力和重力梯度变化估计2012年苏门答腊地震断层参数

2012年4月11日苏门答腊北部附近海域发生M_W=8.6地震,国际上采用不同数据获得了该地震的断层滑动模型,但断层走滑性质存在较大差异.卫星重力GRACE（Gravity Recovery and Climate Experiment）观测覆盖震中区域,可以提供很好的断层参数估计约束.本文采用GRACE月重力场模型数据提取了此次地震同震重力和重力梯度变化,表明了北向分量的优越性.并与断层模型理论计算结果进行了比较,分析了T_xx（北-北向重力梯度）分量对断层参数的敏感性,以及基于粒子群算法以及Okubo位错模型反演了该地震断层参数.结果表明GRACE观测到的同震重力梯度变化在空间形态分布上与断层模型模拟结果符合较好,但振幅差异较大.重力和重力梯度的北向分量可以很好地压制条带误差,其中g_N（北向重力变化）和T_xx的误差最小,其次是T_xz（北-上向重力梯度）,误差最大和对条带误差抑制效果最差的是T_xy（北-东向重力梯度）.T_xx对断层的深度不敏感,对其余断层参数较敏感.基于GRACE反演得到的断层的走向角为113.63°,倾角为89.99°,滑移角为175.26°,平均滑移量为28.18 m,相应的矩震为8.71级,且此次地震的走滑性质为右旋走滑.     

利用独立成分分析检测2004年和2012年印度洋地震的重力变化

2004年12月苏门答腊发生M_W9.3地震,造成巨大的质量重新分布.利用GRACE卫星月重力场数据计算了研究区域地面1°×1°网格点上的重力变化时间序列,采用主成分分析和独立成分分析两种方法,提取了重力变化的空间与时间特征,结果显示震中两侧区域的重力变化呈两极分布,其中东侧重力下降,西侧重力增加.相较于传统的主成分分析方法,独立成分分析能更好地从原始信号中提取地震的信息,能分解出具有显著阶跃变化的独立成分.除了2004年的苏门答腊大地震外,独立成分分析还分解得到了2012年Andaman地震的特征,与该次地震的空间特征与时间序列基本一致.相较于适合定量分析的多项式拟合方法,独立成分分析更适合大范围区域的定性分析,建议将两种方法相结合,取长补短,从而为GRACE地震监测提供一种更为客观、有效的方法.     

联合星载GPS和KBRR星间测速数据反演GRACE时变重力场模型

高精度GRACE卫星时变重力场反演一直是卫星重力测量中的难题.为了恢复高精度的时变地球重力场模型,本文联合GRACE卫星的星载GPS和KBR星间测速观测数据,在对GRACE卫星进行精密定轨的同时,解算出60阶月平均地球重力场模型.通过对GRACE卫星的定轨精度、星载GPS相位和KBR星间测速数据的拟合残差以及时变地球重力场模型解算精度等分析,表明:(1)与美国宇航局喷气推进实验室(JPL)发布的约化动力学精密轨道相比,本文确定GRACE卫星轨道三维位置误差小于5 cm.(2)星载GPS相位数据拟合残差为5~8 mm,KBR星间测速数据拟合残差为0.18~0.30μm·s~(-1).(3)解算的月平均重力场模型与美国德克萨斯大学空间研究中心(CSR)、德国地学研究中心(GFZ)和JPL发布的RL05模型精度接近,时变信号在全球范围内具有很好的空间分布一致性.通过计算亚马逊流域和长江流域的水储量变化,本文与上述三个机构的计算结果无明显差异,且相关系数均达0.9以上.可见,本文建立的卫星轨道与重力场同解算法具有反演高精度GRACE时变重力场能力,为我国卫星重力场反演提供了重要的技术支持.     

不同Mascon模型解比较分析

简要介绍不同机构发布的Mascon模型解数据,并在全球和区域尺度上对Mascon模型解、球谐系数法以及GLDAS水文模型陆地水储量变化计算结果进行对比分析。结果表明,不同Mascon模型解具有较好的一致性;与球谐系数法反演结果相比,Mascon模型解能提高反演结果信噪比值,空间分辨率更高,保留更多有效信号,与GLDAS水文模型结果相关性更好;不同Mascon模型解中,JPL Mascon模型解信噪比值最大,与GLDAS水文模型时间序列相关性最好。     

基于经验正交函数分解的GRACE时变地球重力场滤波处理

针对传统GRACE去相关滤波方法很难同时完美实现真实信号的保留与条带误差的削弱等问题,提出将经验正交函数滤波与去相关滤波相结合的改进策略。分别对真实GRACE数据和模拟数据进行滤波处理,结果表明,改进的滤波方法继承了经验正交函数滤波方法与去相关滤波方法的优点,能够在尽可能地削弱南北条带误差影响的同时,最大程度地保留真实地球物理信号。     

中国西部区域GRACE卫星重力变化

针对研究中国西部长期重力变化的问题,该文利用GRACE月重力场模型,基于Slepian变换构建中国西部区域局部重力场,并采用GIA和GLDAS等模型扣除相应误差,通过时间序列分析方法得到该区域2003—2013年卫星重力变化。结果表明,东天山呈现负重力变化,可能由区域冰川消融引起的;西天山的正重力变化,应与冰川小幅增加和构造运动引起的壳幔物质积累有关;青藏高原内陆正重力变化可能由区域地壳抬升、粘弹性地壳的构造应变以及壳幔物质的质量迁移与积累引起;青藏高原边界区域的负重力变化,应与冰川加速消融以及地下水抽取有关。Slepian方法较好地克服了滤波平滑处理带来的重力变化信号压制及细节平滑,可为中国西部的地壳运动、地震活动和气候变化等研究提供数据支撑。