AIUB-CHAMP02S: The influence of GNSS model changes on gravity field recovery using spaceborne GPS

The gravity field model AIUB-CHAMP02S, which is based on six years of CHAMP GPS data, is presented here. The gravity field parameters were derived using a two step procedure: In a first step a kinematic trajectory of a low Earth orbiting (LEO) satellite is computed using the GPS data from the on-board receiver. In this step the orbits and clock corrections of the GPS satellites as well as the Earth rotation parameters (ERPs) are introduced as known. In the second step this kinematic orbit is represented by a gravitational force model and orbit parameters.     

GPS-only gravity field recovery with GOCE,CHAMP, and GRACE

Gravity missions such as the Gravity field and steady-state Ocean Circulation Explorer (GOCE) are equipped with onboard Global Positioning System (GPS) receivers for precise orbit determination (POD), instrument time-tagging, and the extraction of the long wavelength part of the Earth’s gravity field. The very low orbital altitude of the GOCE satellite and the availability of dense 1 s GPS tracking data are ideal characteristics to exploit the contribution of GPS high-low Satellite-to-Satellite Tracking (hl-SST) to gravity field determination. We present gravity field solutions based on about 8 months of GOCE GPS hl-SST data from 2009 and compare the results with those obtained from the CHAllenging Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) missions. The very low orbital altitude of GOCE significantly improves gravity field recovery from GPS hl-SST data above degree 20, but not for the degrees below 20, where the quality of the spherical harmonic coefficients remains essentially unchanged. Despite the limited time span of GOCE data used, the gravity field of the Earth can be resolved up to about degree 115 using GPS data only. Empirically determined phase center variations (PCVs) of the GOCE onboard GPS helix antenna are, however, mandatory to achieve this performance.     

Onboard orbit determination using GPS observations based on the unscented Kalman filter

Spaceborne GPS receivers are used for real-time navigation by most low Earth orbit (LEO) satellites. In general, the position and velocity accuracy of GPS navigation solutions without a dynamic filter are 25 m (1σ) and 0.5 m/s (1σ), respectively. However, GPS navigation solutions, which consist of position, velocity, and GPS receiver clock bias, have many abnormal excursions from the normal error range for space operation. These excursions lessen the accuracy of attitude control and onboard time synchronization. In this research, a new onboard orbit determination algorithm designed with the unscented Kalman filter (UKF) was developed to improve the performance. Because the UKF is able to obtain the posterior mean and covariance accurately by using the second-order Taylor series expansion through the sampled sigma points that are propagated by using the true nonlinear system, its performance can be better than that of the extended Kalman filter (EKF), which uses the linearized state transition matrix to predict the covariance. The dynamic models for orbit propagation applied perturbations due to the 40 × 40 geo-potential, the gravity of the Sun and Moon, solar radiation pressure, and atmospheric drag. The 7(8)th-order Runge–Kutta numerical integration was applied for orbit propagation. Two types of observations, navigation solutions and C/A code pseudorange, can be used at the user’s discretion. The performances of the onboard orbit determination were verified using real GPS data of the CHAMP and KOMPSAT-2 satellites. The results of the orbit determination were compared with the precision orbit ephemeris (POE) of the CHAMP and KOMPSAT-2 satellites.     

Gravity field models from kinematic orbits of CHAMP,GRACE and GOCE satellites

The aim of our work is to generate Earth’s gravity field models from GPS positions of low Earth orbiters. Our inversion method is based on Newton’s second law, which relates the observed acceleration of the satellite with forces acting on it. The observed acceleration is obtained as numerical second derivative of kinematic positions. Observation equations are formulated using the gradient of the spherical harmonic expansion of the geopotential. Other forces are either modelled (lunisolar perturbations, tides) or provided by onboard measurements (nongravitational perturbations). From this linear regression model the geopotential harmonic coefficients are obtained.     

Characteristics of GOCE orbits based on Satellite Laser Ranging

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first European Space Agency’s (ESA) Earth Explorer core mission. Through its extremely low, about 260?km above the Earth, circular, sun-synchronous orbit, the satellite gained high spatial resolution and accuracy gravity gradient, and ocean circulation data. Global Positioning System (GPS) receivers, mounted on the spacecraft, allowed the determination of reduced-dynamic and kinematic GOCE orbits, whereas Laser Retroreflector Array (LRA) dedicated to Satellite Laser Ranging (SLR) allowed an independent validation of GPS-derived orbits. In this paper, residuals between different GPS-based orbit types and SLR observations are used to investigate the sensitivity and the influence of solar, geomagnetic, and ionospheric activities on the quality of kinematic and reduced-dynamic GOCE orbits. We also analyze the quality of data provided by individual SLR sites, by detecting time biases using ascending and descending sun-synchronous GOCE orbit passes, and the residual analysis of the measurement characteristics, i.e., the dependency of SLR residuals as a function of nadir and horizontal angles. Results show a substantial vulnerability of kinematic orbit solutions to the solar F10.7 index and the ionospheric activity measured by the variations of the Total Electron Content (TEC) values. The sensitivity of kinematic orbits to the three-hour-range KP index is rather minor. The reduced-dynamic orbits are almost insensitive to indices describing ionospheric, solar, and geomagnetic activities. The investigation of individual SLR sites shows that some of them are affected by time bias errors, whereas other demonstrate systematics, such as a dependency between observation residuals and the satellite nadir angle or the horizontal azimuth angle from the SLR station to the direction of the satellite.     

GOCE orbit analysis: Long-wavelength gravity field determination using the acceleration approach

The restricted sensitivity of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gradiometer instrument requires satellite gravity gradiometry to be supplemented by orbit analysis in order to resolve long-wavelength features of the geopotential. For the hitherto published releases of the GOCE time-wise (TIM) and GOCE space-wise gravity field series—two of the official ESA products—the energy conservation method has been adopted to exploit GPS-based satellite-to-satellite tracking information. On the other hand, gravity field recovery from data collected by the CHAllenging Mini-satellite Payload (CHAMP) satellite showed the energy conservation principle to be a sub-optimal choice. For this reason, we propose to estimate the low-frequency part of the gravity field by the point-wise solution of Newton’s equation of motion, also known as the acceleration approach. This approach balances the gravitational vector with satellite accelerations, and hence is characterized by (second-order) numerical differentiation of the kinematic orbit. In order to apply the method to GOCE, we present tailored processing strategies with regard to low-pass filtering, variance–covariance information handling, and robust parameter estimation. By comparison of our GIWF solutions (initials GI for “Geodätisches Institut” and IWF for “Institut für WeltraumForschung”) and the GOCE-TIM estimates with a state-of-the-art gravity field solution derived from GRACE (Gravity Recovery And Climate Experiment), we conclude that the acceleration approach is better suited for GOCE-only gravity field determination as opposed to the energy conservation method.     

A simulation of Martian gravity field recovery by using a near equatorial orbiter

An improvement to the Martian gravity field may be achieved by means of future orbiting spacecraft with small eccentricity and low altitude exemplified through a newly proposed mission design that may be tested in upcoming reconnaissance of Mars. Here, the near equatorial orbital character (with an inclination approximating 10°, eccentricity as 0.01 and semi-major axis as 4000 km) is considered, and its contribution to Martian gravity field solution is analyzed by comparing it with a hypothetical polar circular orbiter. The solution models are evaluated in terms of the following viewpoints: power spectra of gravity field coefficients, correlations of low degree zonal coefficients, precise orbit determination, and error distribution of both Mars free air gravity anomaly and areoid. At the same time, the contributions of the near equatorial orbiters in low degree zonal coefficients time variations are also considered. The present results show that the near equatorial orbiter allows us to improve the accuracy of the Martian gravity field solution, decrease correlation of low degree zonal coefficients, retrieve much better time variable information of low degree zonal coefficients, improve precise orbit determination, and provide more accurate Mars free air gravity anomaly and areoid around the equatorial region.     

Proposal of application of same beam VLBI measurements in precision orbit determination of lunar orbiter and return capsule and lunar gravity field simulation in Chang’E-3 mission

High accuracy differenced phase delay can be obtained by observing multiple point frequencies of two spacecraft using the same beam Very Long Baseline Interferometry (VLBI) technology. Its contribution in lunar spacecraft precision orbit determination has been performed during the Japanese lunar exploration mission SELENE. In consideration that there will be an orbiter and a return capsule flying around the moon during the Chinese lunar exploration future mission Chang’E-3, the contributions of the same beam VLBI in spacecraft precision orbit determination and lunar gravity field solution have been investigated. Our results show that the accuracy of precision orbit determination can be improved more than one order of magnitude after including the same beam VLBI measurements. There are significant improvements in accuracy of low and medium degree coefficients of lunar gravity field model obtained from combination of two way range and Doppler and the same beam VLBI measurements than the one that only uses two way range and Doppler data, and the accuracy of precision orbit determination can reach meter level.     

Modeling and analysis of Earth’s gravity field measurement performance by inner-formation flying system

The Earth’s gravity field can be measured with high precision by constructing the purely gravitational orbit of the inner-satellite in Inner-formation Flying System (IFS), which is independently proposed by Chinese scholars and offers a new way to carry out gravity field measurement by satellite without accelerometers. In IFS, for the purpose of quickly evaluating the highest degree of recovered gravity field model and geoid error as well as analyzing the influence of system parameters on gravity field measurement, an analytical formula was established by spectral analysis method. The formula can reflect the analytical relationship between gravity field measurement performance and system parameters such as orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and total measurement time. This analytical formula was then corrected by four factors introduced from numerical simulation of IFS gravity field measurement. By comparing computation results from corrected analytical formula and the actual gravity field measurement performance by CHAMP, the correctness and rationality of this analytical formula were verified. Based on this analytical formula, the influences of system parameters on IFS gravity field measurement were analyzed. It is known that gravity field measurement performance is a monotone decreasing function of orbit altitude, the inner-satellite orbit determination error, the inner-satellite residual disturbances, data sampling interval and the reciprocal of total measurement time. There is a match relationship between the inner-satellite orbit determination error and residual disturbances, in other words, the change rate of gravity field measurement performance with one of them is seriously restricted by their relative size. The analytical formula can be used to quantitatively evaluate gravity field measurement performance fast and design IFS parameters optimally. It is noted that the analytical formula and corresponding conclusions are applied to any gravity satellite which measures gravity field by satellite perturbation orbit.     

近地卫星的GPS自主定轨算法研究

为了用GPS实现近地卫星的高精度自主定轨,提出了采用序贯UD协方差分解滤波算法作为定轨算法.完整地给出了序贯UD分解算法的计算公式,详细地推导了系统的状态方程和观测方程,建立了由9维状态向量和8维观测向量组成的导航系统模型.最后作了数字仿真试验,试验结果表明,文中的定轨算法既有较高的精度又便于星上计算机的实现.      

Enhancing the kinematic precise orbit determination of low earth orbiters using GPS receiver clock modelling

Clock error estimation has been the focus of a great deal of research because of the extensive usage of clocks in GPS positioning applications. The receiver clock error in the spacecraft orbit determination is commonly estimated on an epoch-by-epoch basis, along with the spacecraft’s position. However, due to the high correlation between the spacecraft orbit altitude and the receiver clock parameters, estimates of the radial component are degraded in the kinematic approach. Using clocks with high stability, the predictable behaviour of the receiver oscillator can be exploited to improve the positioning accuracy, especially for the radial component. This paper introduces two GPS receiver clock models to describe the deterministic and stochastic property of the receiver clock, both of which can improve the accuracy of kinematic orbit determination for spacecraft in low earth orbit. In particular, the clock parameters are estimated as time offset and frequency offset in the two-state model. The frequency drift is also estimated as an unknown parameter in the three-state model. Additionally, residual non-deterministic random errors such as frequency white noise, frequency random walk noise and frequency random run noise are modelled. Test results indicate that the positioning accuracy could be improved significantly using one day of GRACE flight data. In particular, the error of the radial component was reduced by over 40.0% in the real-time scenario.     

The challenges in long-term altimetry calibration for addressing the problem of global sea level change

Long-term change of the global sea level resulting from climate change has become an issue of great societal interest. The advent of the technology of satellite altimetry has modernized the study of sea level on both global and regional scales. In combination with in situ observations of the ocean density and space observations of Earth’s gravity variations, satellite altimetry has become an essential component of a global observing system for monitoring and understanding sea level change. The challenge of making sea level measurements with sufficient accuracy to discern long-term trends and allow the patterns of natural variability to be distinguished from those linked to anthropogenic forcing rests largely on the long-term efforts of altimeter calibration and validation. The issues of long-term calibration for the various components of the altimeter measurement system are reviewed in the paper. The topics include radar altimetry, the effects of tropospheric water vapor, orbit determination, gravity field, tide gauges, and the terrestrial reference frame. The necessity for maintaining a complete calibration effort and the challenges of sustaining it into the future are discussed.     

基于GPS的地球资源卫星精密星历外推新方法

随着技术的发展,通过星载GPS接收机直接确定卫星星历成为卫星定位的一个重要手段.GPS接收机获取的卫星星历数据是某一时刻的瞬时状态,要获取连续的卫星星历数据还需要进一步处理.常用的处理方法有几何法与动力学法.在GPS接收机给定瞬时星历频率较低的情况下,几何法的计算误差比较大,特别是只有一组瞬时星历时,无法用几何法进行轨道的外推.在分析地球资源卫星轨道特点的基础上,提出一种新的轨道缩减动力模型,该模型将卫星运动在直角坐标系中分解为简谐运动,利用模型实现了轨道外推的算法.通过试验验证,该算法可以达到较高的精度.      

基于自适应滤波的编队卫星实时相对定轨

提出基于自适应滤波的编队卫星实时相对定轨算法,利用2005-12-09—10两颗GRACE（Gravity Recovery and Climate Experiment）卫星的GPS（Global Positioning System）实测数据进行实时相对定轨试验计算,采用JPL（Jet Propulsion Laboratory）轨道对试验结果外部检核,结果表明:①自适应滤波相对定轨通过自适应因子,可以较好地平衡编队卫星的观测信息和相对动力学信息,其相对定轨结果精度优于Kalman滤波相对定轨结果;②自适应滤波相对定轨结果随着星间基线缩短而精度提高;③两颗GRACE卫星采用单频伪距和广播星历进行自适应滤波相对定轨,可以得到精度优于6cm的星间基线。     

GPS在卫星精密测轨中的应用

自从美国全球定位系统(GPS)试验系统建立以来,发展了各种各样的方法,利用GPS 测量来精密确定卫星轨道。在1984年,LANDSAT-5上装载一个GPSPAC 的仪器进行了飞行,以验证用GPS 定轨的精确度。1982年在美国喷气推进试验室((?)PL),为TOPEX 卫星研究了一种GPS 距离和距离变化率的双差分法,以满足其厘米级的精度要求。除GPS 对低轨地球卫星定轨应用外,有人也提出了用GPS 测量对高轨地球卫星甚至同步卫星的定轨方法及精度分析。文章对上述各种方法给以简要的介绍和评述。     

Yarkovsky-Schach effect on space debris motion

The Yarkovsky-Schach effect is a small perturbation affecting Earth satellites and space debris illuminated by the Sun. It was first applied to the orbit of LAGEOS satellites as an explanation of the residuals in orbital elements. In this work, we carry out several numerical integration tests taking into consideration various orbit and rotation parameters, in order to analyse this effect in a broader context. The semi-major axis variations remain small and depend on the spin axis attitude with respect to the Sun. We show that the force amplitude is maximised for orbits inclined with $i\approx$?20–30°. We also observe the influence on other orbital elements, notably on the orbit inclination. However, these effects are clearly observed only on long timescales; in our simulations, we propagated the orbits for 200?y. The Yarkovsky-Schach effect is thus confirmed to have a minuscule magnitude. It should be taken into account in studies requiring high-precision orbit determination, or on expanded timescales.     

基于经验加速度的低轨卫星轨道预报新方法

研究将定轨过程中的经验加速度应用于地球低轨卫星轨道预报的新方法. 利用GPS伪距观测数据和简化动力学最小二乘批处理方法对地球低轨卫星定 轨, 其中卫星位置、速度及大气阻力系数和辐射光压系数可以直接用于轨道预报. 作为简化动力学最重要特征的经验加速度呈现准周期、余弦曲线特点, 可通过 傅里叶级数拟合建模. 确定性动力学模型与补偿大气阻力模型误差的切向经验 加速度级数拟合模型组成增强型动力学模型用于提高轨道预报精度. 应用 GRACE-A星载GPS伪距观测数据和IGS超快星历定轨并进行轨道预报, 结果表明 轨道预报初值位置精度达到0.2m, 速度精度达到1×10-4m·s-1, 预报3天位置精度优于60m, 比只利用确定性动力学模型进行预报精度平 均提高2.3倍. 先定轨后预报的模式可用在星上自主精确导航系统中.      

GPS定位数据在海洋一号卫星中的应用与比较

简要介绍了海洋一号卫星星载GPS接收机的定位原理、流程和应用，探讨了一种在轨定位结果的确认和互验方法．借助于卫星工具包软件(STK)，利用NASA网站公布的HY-1卫星两行根数(TLE)进行卫星轨道推算，生成星下点位置，并与相应时刻星载GPS接收机实测数据得到的星下点位置进行比较，由此得到两种方法定位结果之间的偏差，用实际在轨数据验证了两者的位置符合程度．     

基于双频GPS观测的简化动力学最小二乘批处理精密定轨

星载双频GPS载波相位和伪距观测量已成为低轨卫星获取精确三维位置和速度信息的主要方式. 本文以非差消电离载波相位和伪距组合作为观测量,应用简化动力学最小二乘批处理方法进行地球低轨卫星的精密定轨,并给出完整定轨流程. 采用逐段常量的经验加速度对动力学模型误差进行补偿,描述了经验加速度敏感矩阵及稀疏带状矩阵求逆的有效计算方法. 利用GRACE-A卫星GPS观测数据对定轨位置精度进行分析,结果显示,三维位置定轨精度优于5cm,经验加速度在径向、切向和法向上的补偿水平不超过40nm·s-2,大气阻力系数和辐射光压系数的估计值符合物理实际,星载接收机钟差大致呈线性并具有短周期小波动.      

The GLONASS-M satellite yaw-attitude model

The proper modeling of the satellites’ yaw-attitude is a prerequisite for high-precision Global Navigation Satellite System (GNSS) positioning and poses a particular challenge during periods when the satellite orbital planes are partially eclipsed. Whereas a lot of effort has been put in to examine the yaw-attitude control of GPS satellites that are in eclipsing orbits, hardly anything is known about the yaw-attitude behavior of eclipsing GLONASS-M satellites. However, systematic variations of the carrier phase observation residuals in the vicinity of the orbit’s noon and midnight points of up to ±27 cm indicate significant attitude-related modeling issues. In order to explore the GLONASS-M attitude laws during eclipse seasons, we have studied the evolution of the horizontal satellite antenna offset estimates during orbit noon and orbit midnight using a technique that we refer to as “reverse kinematic precise point positioning”. In this approach, we keep all relevant global geodetic parameters fixed and estimate the satellite clock and antenna phase center positions epoch-by-epoch using 30-second observation and clock data from a global multi-GNSS ground station network. The estimated horizontal antenna phase center offsets implicitly provide the spacecraft’s yaw-attitude. The insights gained from studying the yaw angle behavior have led to the development of the very first yaw-attitude model for eclipsing GLONASS-M satellites. The derived yaw-attitude model proves to be much better than the nominal yaw-attitude model commonly being used by today’s GLONASS-capable GNSS software packages as it reduces the observation residuals of eclipsing satellites down to the normal level of non-eclipsing satellites and thereby prevents a multitude of measurements from being incorrectly identified as outliers. It facilitates continuous satellite clock estimation during eclipse and improves in particular the results of kinematic precise point positioning of ground-based receivers.