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.     

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.     

基于GPS的中、高轨道航天器定轨研究

利用GPS确定空间航天器轨道，尤其是GPS在中、高空间航天器定轨中的应用，针对中、高轨道中可见星少及GPS信号微弱等问题，提出了有效的微弱信号检测算法，提高了GPS接收机的接收能力，同时提出了利用伪加速度测量和轨道约束进行定轨解算的新方法，使得在HEO(高偏心轨道)中利用两颗可见星可以解算定位解，对GEO(静止地球轨道)只用一颗星就可求解，仿真结果表明该方法具有很好的效果。     

粗差拟准检定法在星载GPS低轨卫星定轨中的应用

利用自适应卡尔曼滤波进行星载GPS低轨卫星定轨时,必须解决量测方程中经常存在的粗差问题．在分析以往方法的优缺点后,用拟准检定法来探测和修正量测方程中存在的粗差．该法的优点是辨识粗差准确率高,能同时定位多个粗差．另外,为了克服星载GPS低轨卫星定轨的滤波器可能出现的数值不稳定性及发散现象,还采用了UD分解算法及Sage自适应滤波器．最后用一个CHAMP卫星的模拟算例验证本方法的可行性和有效性．     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

基于GPS/北斗组合系统的高轨卫星定位技术研究

针对目前高轨GPS信号可用性差及定位精度低的特点, 对GPS/北斗组合系统的 高轨卫星定位技术进行研究, 对比分析了单GPS系统与GPS/北斗组合系统的卫 星可见性和几何精度因子. 结果表明, GPS/北斗组合系统比单GPS系统的卫星可 见性好, 且定位精度高. 同时通过提出在星载接收机上采用高精度原子钟, 可实现三星定位, 降低对接收机的技术要求.      

Near real-time BDS GEO satellite orbit determination and maneuver analysis with reversed point positioning

The Geostationary Earth Orbit (GEO) satellite is a crucial part of the BeiDou Navigation Satellite System (BDS) constellation. However, due to various perturbation forces acting on the GEO satellite, it drifts gradually over time. Thus, frequent orbit maneuvers are required to maintain the satellite at its designed position. During the orbit maneuver and recovery periods, the orbit quality of the maneuvered satellite computed with broadcast navigation ephemeris will be significantly degraded. Furthermore, the conventional dynamic Precise Orbit Determination (POD) approach may not work well, because of a lack of publicly available satellite information for modeling the thrust forces. In this paper, a near real-time approach free of thrust forces modeling is proposed for BDS GEO satellite orbit determination and maneuver analysis based on the Reversed Point Positioning (RPP). First, the station coordinates and receiver clock offsets are estimated by GPS/BDS combined Single Point Positioning (SPP) with single-frequency phase-smoothed pseudorange observations. Then, with the fixed station coordinates and receiver clock offsets, the RPP method can be conducted to determine the GEO satellite orbits. When no orbit maneuvers occur, the proposed method can obtain orbit accuracies of 0.92, 2.74, and 8.30?m in the radial, along-track, and cross-track directions, respectively. The average orbit-only Signal-In-Space Range Error (SISRE) is 1.23?m, which is slightly poorer than that of the broadcast navigation ephemeris. Using four days of GEO maneuvered datasets, it is further demonstrated that the derived orbits can be employed to characterize the behaviors of GEO satellite maneuvers, such as the time span of the maneuver as well as the satellite thrusting accelerations. These results prove the efficiency of the proposed method for near real-time GEO satellite orbit determination during maneuvers.     

Estimation of differential code biases with Jason-2/3 onboard GPS observations

With the continuous deployment of Low Earth Orbit (LEO) satellites, the estimation of differential code biases (DCBs) based on GNSS observations from LEO has gained increasing attention. Previous studies on LEO-based DCB estimation are usually using the spherical symmetry ionosphere assumption (SSIA), in which a uniform electron density is assumed in a thick shell. In this study, we propose an approach (named the SHLEO method) to simultaneously estimate the satellite and LEO onboard receiver DCBs by modeling the distribution of the global plasmaspheric total electron content (PTEC) above the satellite orbit with a spherical harmonic (SH) function. Compared to the commonly used SSIA method, the SHLEO model improves the GPS satellite DCB estimation accuracy by 13.46% and the stability by 22.34%, respectively. Compared to the GPS satellite DCBs estimated based on the Jason-3-only observations, the accuracy and monthly stability of the satellite DCBs can be improved by 14.42% and 26.8% when both Jason-2 and Jason-3 onboard observations are jointly processed. Compared with the Jason-2 solutions, the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations have an improved consistency of better than 18.26% and 9.71% with the products provided by the Center for Orbit Determination in Europe (CODE) and Chinese Academy of Sciences (CAS). Taking the DCB products provided by the German Aerospace Center (DLR) as references, there is no improvement in accuracy of the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations than the Jason-2 solutions alone. A periodic variation is found in the time series of both the Jason-3 and Jason-2 onboard receiver DCB estimates. Preliminary analysis of the PTEC distribution based on the estimated SH coefficients are also presented.     

卫星与测站钟差支持条件下的GEO卫星精密定轨

在基于伪距的GEO卫星精密定轨中, GEO卫星的静地特性导致定轨解算无法对星地组合钟差进行有效估计, 需要独立的时间同步支持. 本文讨论了卫星和测站钟差支持条件下的GEO卫星定轨原理, 利用仿真数据系统地分析了中国区域网跟踪条件下GEO卫星的定轨精度, 从定性和定量角度分析了钟差二次项、星地时间同步精度、站间时间同步精度及系统差等因素对定轨精度的影响.      

Particle filter-based real-time phase line bias estimation for GNSS-based attitude determination with common-clock receivers

Global navigation satellite system (GNSS)-based attitude determination has been widely adopted in a wide variety of terrestrial, sea, air, and space applications. Recently, the emergence of commercial multi-GNSS common-clock receivers has brought new opportunities for high-precision GNSS-based attitude determination with single-differenced (SD) model. However, the key requirement of using this approach is the accurate estimation of the troublesome line bias (LB) in real-time. In this contribution, we propose a particle filter-based real-time phase LB estimation approach that apply to SD model with single-system single-frequency observations from common-clock receiver. We first analyzed the relationship between the integer ambiguity ratio value and the phase LB. It is proved that the accuracy of a given phase LB value can be qualified by the related ambiguity resolution ratio value, and the normalized ratio value can therefore be used to represent the likelihood function of observations. Then, we presented the particle filter-based real-time phase LB estimation procedure, and assessed its performance using GPS L1/BDS B1I observations from two datasets collected with different types of common-clock receivers in terms of the accuracy and convergence time of phase LB estimation, the computation load, and the positioning and attitude determination accuracy with respect to the double-differenced (DD) model. Experimental results demonstrated that the phase LB could be accurately estimated with short convergence time (generally within 15 epochs). Moreover, compared with the classical DD approach, the particle filter-based SD approach delivers comparable positioning root-mean-square (RMS) errors in the North and East components but significantly smaller RMS errors in the Up component. Accordingly, the achievable yaw accuracy is comparable whereas the pitch accuracy is remarkably improved. The improvements of positioning accuracy in the Up component and pitch accuracy are approximately 35.7 % to 63.7 %, and 33.3 % to 63.1 %, respectively. Additionally, the single-epoch computation time with our particle filter-based SD approach is generally 0.08 s, which is obviously larger than the DD approach but could still meet the requirements of real-time applications below 10 Hz sampling.     

Performance improvement of real-time PPP ambiguity resolution using a regional integer clock

Obtaining reliable GNSS uncalibrated phase delay (UPD) or integer clock products is the key to achieving ambiguity-fixed solutions for real-time (RT) precise point positioning (PPP) users. However, due to the influence of RT orbit errors, the quality of UPD/integer clock products estimated with a globally distributed GNSS network is difficult to ensure, thereby affecting the ambiguity resolution (AR) performance of RT-PPP. In this contribution, by fully utilising the consistency of orbital errors in regional GNSS network coverage areas, a method is proposed for estimating regional integer clock products to compensate for RT orbit errors. Based on Centre National d’Études Spatiales (CNES) RT precise products, the regional GPS/BDS integer clock was estimated with a CORS network in the west of China. Results showed that the difference between the estimated regional and CNES global integer clock/bias products was generally less than 5 cm for GPS, whereas clock differences of greater than 10 cm were observed for BDS due to the large RT orbit error. Compared with PPP using global integer clock/bias products, the AR performance of PPP using the regional integer clock was obviously improved for four rover stations. For single GPS, the horizontal and vertical accuracies of ambiguity-fixed PPP solutions were improved by 56.2% and 45.3% on average, respectively, whereas improvements of 67.5% and 50.5% in the horizontal and vertical directions, respectively, were observed for the combined GPS/BDS situation. Based on a regional integer clock, the RMS error of a kinematic ambiguity-fixed PPP solution in the horizontal direction could reach 0.5 cm. In terms of initialisation time, the average time to first fix (TTFF) in combined GPS/BDS PPP was shortened from 18.2 min to 12.7 min. In view of the high AR performance realised with the use of regional integer clocks, this method can be applied to scenarios that require high positioning accuracy, such as deformation monitoring.     

基于实时SSR数据的LEO卫星精密单点定位研究

对目前低轨卫星实时定位的方法进行了研究,现在通常采用GPS定位,使用广播星历和普通晶振,实时定位精度一般在10m以内,不能满足高精度实时定位的需求。IGS组织在全球范围内对GPS跟踪分析,生成精密星历、精密钟差产品、按SSR格式的广播星历和钟差修正产品并在网上发布。对这些IGS产品进行了调查,提出在现有测控支持情况下,可以通过高密度上注SSR信息流实现在轨高精度定位。以某型号低轨微小卫星在轨导航增强载荷为应用背景,用IGS03产品中的1057和1058数据对双频GPS接收机的星历和钟差进行修正,采用递推最小二乘估计和LAMDA模糊度固定过对载波相位和伪距信息进行处理,在短时间内获得亚米级定位结果。     

VLBI observations to the APOD satellite

The APOD (Atmospheric density detection and Precise Orbit Determination) is the first LEO (Low Earth Orbit) satellite in orbit co-located with a dual-frequency GNSS (GPS/BD) receiver, an SLR reflector, and a VLBI X/S dual band beacon. From the overlap statistics between consecutive solution arcs and the independent validation by SLR measurements, the orbit position deviation was below 10?cm before the on-board GNSS receiver got partially operational. In this paper, the focus is on the VLBI observations to the LEO satellite from multiple geodetic VLBI radio telescopes, since this is the first implementation of a dedicated VLBI transmitter in low Earth orbit. The practical problems of tracking a fast moving spacecraft with current VLBI ground infrastructure were solved and strong interferometric fringes were obtained by cross-correlation of APOD carrier and DOR (Differential One-way Ranging) signals. The precision in X-band time delay derived from 0.1?s integration time of the correlator output is on the level of 0.1?ns. The APOD observations demonstrate encouraging prospects of co-location of multiple space geodetic techniques in space, as a first prototype.     

Beidou wide-area augmentation system clock error correction and performance verification

On December 27, 2018, the Beidou-3 System (BDS-3) has completed the deployment of 18 Medium Earth Orbit (MEO) satellites combined as a space constellation. In addition to the augmentation information for the new system signals B1C and B2a, the BDS-3 is compatible with the three augmentation information broadcast by the BDS-2 system for B1I and B3I signals: equivalent clock error correction, User Differential Ranging Error (UDRE) and Grid Ionosphere Vertical Error (GIVE). In this paper, the observation data of Beidou monitoring network are used to analyze the pseudo-range observation quality of the smooth transition signals B1I and B3I of BDS-2 and BDS-3. At the same time, the relationship between the equivalent clock error correction and the prediction error of the satellite clock is analyzed by using the Two-Way Satellite Time and Frequency Transfer (TWSTFT) data. The results show that the correlation between the equivalent clock error correction and the monitored clock error by using the TWSTFT data is greater than 60%. We calculate the UDRE by using the equivalent clock error correction. The results show that the satellite equivalent clock error correction can improve the accuracy of User Equivalent Range Error (UERE) by about 50%. This paper also compares the positioning accuracy of the BDS-2 satellites with the BDS-2 satellites combined BDS-3 satellites. The results show that the three-dimensional positioning accuracy is improved by about 30% after the BDS-3 satellites are added.     

GPS授时校频方法研究与试验结果

为了解决多目标综合测量系统各测站之间时间同步和频率校准问题,提出了利用GPS(Global Positioning System)单星或多星共视方法进行站间时间同步与校频,给出了这两种方法的计算公式,分析了星历误差、星钟误差、电离层折射误差、对流层折射误差、多径效应和接收机硬件延迟对时间同步精度的影响.为了验证GPS授时校频精度,进行了相关试验.通过与铯原子钟比对,表明利用GPS可实现纳秒级时间同步,校频精度也优于5.0×10-11,多星共视具有更高的同步校频精度.      

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.     

Autonomous broadcast ephemeris improvement for GNSS using inter-satellite ranging measurements

This paper discusses the concept of using inter-satellite ranging (ISR) measurements of the satellites of a Global Navigation Satellite System (GNSS) for autonomous broadcast ephemeris improvement. Firstly the inter-satellite ranging is modeled to obtain the clock and orbit error observables. The orbit error observable is analyzed and its observation equation is provided. Both least-squares estimation and Kalman Filter approach are proposed to estimate satellite clock errors, while solely the Kalman Filter is used to estimate the orbit errors. All these algorithms are validated using true broadcast ephemeris and precise ephemeris of GPS along with the simulated ranging noise. Based on the settings adopted in the test, the result shows that the orbit accuracy of the precise ephemeris using the proposed method is around 20–50 cm and the accuracy of the satellite clock can reach 20 cm while ranging noise is assumed to be 0.45 m (1σ). The User Ranging Error (URE) is improved from 1.05 m to 0.34 m, which is comparable to other sources of precise ephemeris or even better, while the proposed approach has many advantages such as compatibility and accessibility. It is also noted that the proposed method may provide useful functions in determining inter-constellation coordinate and time differences autonomously for better interoperability and interchangeability in the multi GNSS operation era.     

星载GPS几何法实时定轨有关问题的研究

首先讨论了星载GPS几何法实时定轨的绝对定位方法和各种差分技术。由于伪距差技术能克服GPS卫星的星历误差、卫星钟误差，特别是SA误差的影响，而且实现难度不大，所以应用它来实时定轨。实测数据的处理表明，它能明显提高定轨的精度。然后分析了星载GPS所受扰动影响的情况，对应用抗差估计削弱GPS卫星信号扰动的影响进行了试验，试验的结果说明抗差估计能进一步提高星载GPS几何法定轨的精度。     

Impact of GPS antenna phase center variations on precise orbits of the GOCE satellite

The first European Space Agency Earth explorer core mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) has been launched on March 17, 2009. The 12-channel dual-frequency Global Positioning System receiver delivers 1 Hz data and provides the basis for precise orbit determination (POD) on the few cm-level for such a very low orbiting satellite (254.9 km). As a member of the European GOCE Gravity Consortium, which is responsible for the GOCE High-level Processing Facility (HPF), the Astronomical Institute of the University of Bern (AIUB) provides the Precise Science Orbit (PSO) product for the GOCE satellite. The mission requirement for 1-dimensional POD accuracy is 2 cm. The use of in-flight determined antenna phase center variations (PCVs) is necessary to meet this requirement. The PCVs are determined from 154 days of data and the magnitude is up to 3-4 cm. The impact of the PCVs on the orbit determination is significant. The cross-track direction benefits most of the PCVs. The improvement is clearly seen in the orbit overlap analysis and in the validation with independent Satellite Laser Ranging (SLR) measurements. It is the first time that SLR could validate the cross-track component of a LEO orbit.