GLONASS precise orbit determination with identification of malfunctioning spacecraft

Due to the continued development of the GLONASS satellites, precise orbit determination (POD) still poses a series of challenges. This study examines the impact of introducing the analytical tube-wing model for GLONASS-M and the box-wing model for GLONASS-K in a series of hybrid POD strategies that consider both the analytical model and a series of empirical parameters. We assess the perturbing accelerations acting on GLONASS spacecraft based on the analytical model. All GLONASS satellites are equipped with laser retroreflectors for satellite laser ranging (SLR). We apply the SLR observations for the GLONASS POD in a series of GNSS + SLR combined solutions. The application of the box-wing model significantly improves GLONASS orbits, especially for GLONASS-K, reducing the STD of SLR residuals from 92.6 to 27.6 mm. Although the metadata for all GLONASS-M satellites reveal similar construction characteristics, we found differences in empirical accelerations and SLR offsets not only between GLONASS-M and GLONASS-M+ but also within the GLONASS-M+ series. Moreover, we identify satellites with inferior orbit solutions and check if we can improve them using the analytical model and SLR observations. For GLONASS-M SVN730, the STD of the SLR residuals for orbits determined using the empirical solution is 48.7 mm. The STD diminishes to 41.2 and 37.8 mm when introducing the tube-wing model and SLR observations, respectively. As a result, both the application of the SLR observations and the analytical model significantly improve the orbit solution as well as reduce systematic errors affecting orbits of GLONASS satellites.

Stability of CubeSat Clocks and Their Impacts on GNSS Radio Occultation

Global Navigation Satellite Systems’ radio occultation (GNSS-RO) provides the upper troposphere-lower stratosphere (UTLS) vertical atmospheric profiles that are complementing radiosonde and reanalysis data. Such data are employed in the numerical weather prediction (NWP) models used to forecast global weather as well as in climate change studies. Typically, GNSS-RO operates by remotely sensing the bending angles of an occulting GNSS signal measured by larger low Earth orbit (LEO) satellites. However, these satellites are faced with complexities in their design and costs. CubeSats, on the other hand, are emerging small and cheap satellites; the low prices of building them and the advancements in their components make them favorable for the GNSS-RO. In order to be compatible with GNSS-RO requirements, the clocks of the onboard receivers that are estimated through the precise orbit determination (POD) should have short-term stabilities. This is essential to correctly time tag the excess phase observations used in the derivation of the GNSS-RO UTLS atmospheric profiles. In this study, the stabilities of estimated clocks of a set of CubeSats launched for GNSS-RO in the Spire Global constellation are rigorously analysed and evaluated in comparison to the ultra-stable oscillators (USOs) onboard the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC-2) satellites. Methods for improving their clock stabilities are proposed and tested. The results (i) show improvement of the estimated clocks at the level of several microseconds, which increases their short-term stabilities, (ii) indicate that the quality of the frequency oscillator plays a dominant role in CubeSats’ clock instabilities, and (iii) show that CubeSats’ derived UTLS (i.e., tropopause) atmospheric profiles are comparable to those of COSMIC-2 products and in situ radiosonde observations, which provided external validation products. Different comparisons confirm that CubeSats, even those with unstable onboard clocks, provide high-quality RO profiles, comparable to those of COSMIC-2. The proposed remedies in POD and the advancements of the COTS components, such as chip-scale atomic clocks and better onboard processing units, also present a brighter future for real-time applications that require precise orbits and stable clocks.

Multi-GNSS Combined Orbit and Clock Solutions at iGMAS

Global navigation services from the quad-constellation of GPS, GLONASS, BDS, and Galileo are now available. The international GNSS monitoring and assessment system (iGMAS) aims to evaluate the navigation performance of the current quad systems under a unified framework. In order to assess impact of orbit and clock errors on the positioning accuracy, the user range error (URE) is always taken as a metric by comparison with the precise products. Compared with the solutions from a single analysis center, the combined solutions derived from multiple analysis centers are characterized with robustness and reliability and preferred to be used as references to assess the performance of broadcast ephemerides. In this paper, the combination method of iGMAS orbit and clock products is described, and the performance of the combined solutions is evaluated by various means. There are different internal precisions of the combined orbit and clock for different constellations, which indicates that consistent weights should be assigned for individual constellations and analysis centers included in the combination. For BDS-3, Galileo, and GLONASS combined orbits of iGMAS, the root-mean-square error (RMSE) of 5 cm is achieved by satellite laser ranging (SLR) observations. Meanwhile, the SLR residuals are characterized with a linear pattern with respect to the position of the sun, which indicates that the solar radiation pressure (SRP) model adopted in precise orbit determination needs further improvement. The consistency between combined orbit and clock of quad-constellation is validated by precise point positioning (PPP), and the accuracies of simulated kinematic tests are 1.4, 1.2, and 2.9 cm for east, north, and up components, respectively.

Real-Time Precise Orbit Determination for LEO between Kinematic and Reduced-Dynamic with Ambiguity Resolution

The real-time integer-ambiguity resolution of the carrier-phase observation is one of the most effective approaches to enhance the accuracy of real-time precise point positioning (PPP), kinematic precise orbit determination (KPOD), and reduced-dynamic precise orbit determination (RPOD) for low earth orbit (LEO) satellites. In this study, the integer phase clock (IPC) and wide-lane satellite bias (WSB) products from CNES (Centre National d’Etudes Spatiales) are used to fix ambiguity in real time. Meanwhile, the three models of real-time PPP, KPOD, and RPOD are applied to validate the contribution of ambiguity resolution. Experimental results show that (1) the average positioning accuracy of IGS stations for ambiguity-fixed solutions is improved from about 7.14 to 5.91 cm, with an improvement of around 17% compared to the real-time float PPP solutions, with enhancement in the east-west direction particularly significant, with an improvement of about 29%; (2) the average accuracy of the estimated LEO orbit with ambiguity-fixed solutions in the real-time KPOD and RPOD mode is improved by about 16% and 10%, respectively, with respect to the corresponding mode with the ambiguity-float solutions; (3) the performance of real-time LEO RPOD is better than that of the corresponding KPOD, regardless of fixed- or float-ambiguity solutions. Moreover, the average ambiguity-fixed ratio can reach more than 90% in real-time PPP, KPOD, and RPOD.

Integrity monitoring for precise orbit determination of LEO satellites

Due to an increasing requirement for high accuracy orbital information for low Earth orbit (LEO) satellites, precise orbit determination (POD) of LEO satellites is a topic of growing interest. To assure the safety and reliability of the applications requiring high accuracy LEO orbits in near-real-time, integrity monitoring (IM) is an essential operation of the POD process. In this contribution, the IM strategy for LEO POD in both the kinematic and reduced-dynamic modes is investigated. The overbounding parameters of the signal-in-space range error are investigated for the GPS products provided by the International GNSS Service’s Real-Time Service and the Multi-GNSS Advanced Demonstration of Orbit and Clock Analysis service. Benefiting from the dynamic models used and the improved model strength, the test results based on the data of the LEO satellite GRACE FO-1 show that the average-case mean protection levels (PLs) can be reduced from about 3–4 m in the kinematic mode to about 1 m in the reduced-dynamic mode in the radial, along-track and cross-track directions. The overbounding mean values of the SISRE play the dominant role in the final PLs. In the reduced-dynamic mode and average-case projection, the IM availabilities reach above 99% in the radial, along-track and cross-track directions with the alert limit (AL) set to 2 m. The values are still above 98% with the AL set to 4 m, when the duty cycle of tracking is reduced to 40%, e.g., in the case of power limits for miniature satellites such as CubeSats.

Precise Orbit Determination and Maneuver Assessment for TH-2 Satellites Using Spaceborne GPS and BDS2 Observations

The TH-2 satellite system, including the TH-2A and TH-2B, is the first distributed interferometric synthetic aperture radar (InSAR) satellite system in China. During the in-orbit operation, the TH-2A satellite should perform three maneuvers per day to keep the formation flying geometry. We estimate those maneuvers in the precise orbit determination (POD) by the GPS and BDS2 measurements on board, respectively. The residuals of the POD show that the effects caused by orbital maneuvers can be well eliminated for both the GPS and BDS2 data. The precision of the BDS2-based POD is better than 8.0 cm in the three-dimensional direction (3D) compared with the orbit derived from the GPS observations. Such a precision level of the satellite orbit satisfies the InSAR mission requirement of the TH-2. In addition, the relative error of velocity changes is employed to evaluate the maneuver estimations by the POD using the regional navigation system of BDS2. The results show that the relative error of velocity changes between the GPS- and BDS2-based POD is less than 7.0%, which indicates that the maneuver performance extracted from the regional BDS2 data is as good as that extracted from the global GPS data. In the GNSS fused processing, we found that the independent receiver clock offsets should be taken into account, since the time tag corrections for the GPS and BDS2 observations collected on the TH-2 spaceborne receivers were different. The precision of the GPS and BDS2 (GC) combined single point positioning (SPP) can be improved by 12–14% compared with the GPS-only solution when the position dilution of precision (PDOP) of GPS exceeds three. The overlap comparisons of the GC combined orbits show that the internal orbit precision of the TH-2 satellites is better than 0.7 cm. However, the improvement of the GC combined POD result is only 3–4% with respect to the GPS-only solution, which is limited to the precision of the precise orbit and clock products of BDS2 at the present stage.