Satellite Laser Ranging Residuals Research Articles

SLR Validation and Evaluation of BDS-3 MEO Satellite Precise Orbits

Starting from February 2023, the International Laser Ranging Service (ILRS) began releasing satellite laser ranging (SLR) data for all BeiDou global navigation satellite system (BDS-3) medium earth orbit (MEO) satellites. SLR data serve as the best external reference for validating satellite orbits, providing a basis for comprehensive evaluation of the BDS-3 satellite orbit. We utilized the SLR data from February to May 2023 to comprehensively evaluate the orbits of BDS-3 MEO satellites from different analysis centers (ACs). The results show that, whether during the eclipse season or the yaw maneuver season, the accuracy was not significantly decreased in the BDS-3 MEO orbit products released from the Center for Orbit Determination in Europe (CODE), Wuhan University (WHU), and the Deutsches GeoForschungsZentrum (GFZ) ACs, and the STD (Standard Deviation) of SLR residuals of those three ACs are all less than 5 cm. Among these, CODE had the smallest SLR residuals, with 9% and 12% improvement over WHU and GFZ, respectively. Moreover, the WHU precise orbits exhibit the smallest systematic biases, whether during non-eclipse seasons, eclipse seasons, or satellite yaw maneuver seasons. Additionally, we found some BDS-3 satellites (C32, C33, C34, C35, C45, and C46) exhibit orbit errors related to the Sun elongation angle, which indicates that continued effort for the refinement of the non-conservative force model further to improve the orbit accuracy of BDS-3 MEO satellites are in need.

Impacts of inter-satellite links on the ECOM model performance for BDS-3 MEO satellites

Inter-satellite link (ISL) plays an essential role in current and future Global Navigation Satellite System (GNSS). In this study, we investigate the impact of ISL observations on precise orbit determination for BeiDou-3 Navigation Satellite System (BDS-3) Medium Earth Orbit (MEO) satellites based on different Extended CODE Orbit Models (ECOM). Thanks to the better observation geometry of the Ka-band ISL data compared to the L-band data for BDS-3 MEO satellites, the ISL solution substantially reduces Orbit Boundary Discontinuity (OBD) errors, except for C30, which suffers from unstable Ka-band hardware delay. From the external quality analysis, ISL significantly enhances the reliability of the orbit of MEO satellites manufactured by the China Academy of Space Technology (CAST). The standard deviation (STD) of the satellite laser ranging (SLR) residuals is approximately 2.5 cm, and the root mean square (RMS) is reduced by 10–23% compared to L-band solutions. Besides, the Sun-elongation angle dependent systematic error in SLR residuals nearly vanishes based on the reduced 5-parameter ECOM (ECOM1) or extended 7-parameter ECOM (ECOM2) with ISL data. This is because the ISL reduces the correlation between state parameters and solar radiation pressure (SRP) parameters as well as those among SRP parameters, leading to a more accurate estimation of both orbit and SRP perturbations, particularly those along B direction. This confirms that the deficiency of the SRP models for BDS-3 CAST satellites can be compensated by using better observation geometry from ISL data. On the other hand, for the satellite manufactured by Shanghai Engineering Center for Microsatellites (SECM), the ISL allows for a more accurate estimation of the Bc1 parameter in the ECOM1 model. This only reduces linear systematic error, possibly because the impact generated by the satellite bus cannot be entirely absorbed by the B-direction parameters.

Analysis and diagnosis of abnormal SLR validation results for BeiDou-3 SECM-B MEO C225 and C226 satellite orbits

Satellite laser ranging (SLR) is the only stand-alone and external tool for verifying the accuracy of microwave-based GNSS orbits. On February 2, 2023, the International Laser Ranging Service initiated the SLR global tracking campaign for the complete BeiDou-3 medium Earth orbit (MEO) constellation. This initiative provides an opportunity to unveil the potential issue for BeiDou-3 MEO orbit modeling. An examination of the SLR validation results of precise BeiDou-3 MEO orbits across different IGS/Multi-GNSS Experiment (MGEX) analysis centres (ACs) reveals that the SLR residual standard deviation for BeiDou-3 SECM-B MEOs C225 and C226 orbits reaches 18 cm. This figure is approximately four to nine times those for other BeiDou-3 MEOs. This contradicts the fact that the orbit quality of C225 and C226 is comparable to other satellites. By analyzing the correlation between SLR residuals and de-trend clock residuals, detector-specific bias performance, and the distribution characteristics of SLR residuals relative to satellite azimuth and nadir angles, we diagnose that abnormal SLR validation results for C225 and C226 originate from a large deviation in the official horizontal offsets for the laser retroreflector arrays (LRAs) in the satellite reference frame. LRA horizontal offsets for C225 and C226 are estimated based on SLR residuals. With our adjusted LRA offset estimates, the SLR residual metrics for C225 and C226 orbits (across various IGS/MGEX ACs) are comparable to the results obtained for other BeiDou-3 MEOs. Additionally, our adjustments enable the identification of their orbit modeling errors across different IGS/MGEX ACs according to the pattern of SLR residual variation. From the study, spacecraft manufacturers should be encouraged to conduct a comprehensive review of ground calibrations for BeiDou LRA offsets. This is crucial for optimizing BeiDou orbit model and enhancing its applications in high-precision geosciences.

Earth Rotation Parameters Derived from BDS-3 New Signals B1C/B2a Dual-Frequency Combination Observations

The Earth rotation parameters (ERP) play a crucial role in defining the global reference frame and the Global Navigation Satellite System (GNSS) is one of the important tools used to obtain ERP, including polar motion (PM), its rates, and length of day (LOD). The latest IGS Repro3 ERP products, which provided the IGS contribution to the latest ITRF2020, were generated without consideration of the Beidou Navigation Satellite System (BDS) observations. The global BDS, namely the BDS-3 constellation, has been completely constructed from July 2020 and the observing stations are evenly distributed globally now. Two couple dual-frequency combinations, including the B1I/B3I and B1C/B2a combinations, are commonly used for BDS-3 ionosphere-free combination usage. With the goal of identifying the optimal dual-frequency combination for BDS-3 ERP estimates for the future ITRF definition with a consideration of BDS-3, this research evaluated the performance of ERP estimation using B1I/B3I and B1C/B2a combinations. Firstly, we conducted a comparison of the ambiguity resolutions. The mean percentage of successfully resolved ambiguities for the BDS-3 B1C/B2a combination is 86.5%, surpassing that of B1I/B3I. The GNSS satellite orbits and ERP were estimated simultaneously, thus the accuracy of orbits could also reflect the performance of the ERP estimates. Subsequently, we validated the orbits of 22 BDS-3 Medium Earth Orbit (MEO) satellites using Satellite Laser Ranging (SLR), and the root mean square error (RMS) of the SLR residuals for the 3-day arc orbit with B1C/B2a signals was 5.72 cm, indicating superior accuracy compared with the B1I/B3I combination. Finally, we compared the performance of ERP estimation, considering both internal and external accuracy. For the internal accuracy, B1C/B2a-based solutions demonstrated a reduction in mean formal errors of approximately 17% for PM, 22% for LOD, and 21% for PM rates compared with B1I/B3I-based solutions. In terms of external accuracy, we compared BDS-3-derived ERP estimates with the IERS 20C04 products. The B1C/B2a combination exhibited a slightly better standard deviation performance and a significant reduction in mean bias by 56%, 54%, 39%, 64%, and 23% for X, Y polar motion, X, Y polar motion rates, and LOD, respectively, compared with B1I/B3I solutions. In conclusion, the results highlight the excellent signal quality for BDS-3 B1C/B2a and its superiority in ERP estimation when compared with the B1I/B3I combination.

Satellite laser ranging to BeiDou-3 satellites: initial performance and contribution to orbit model improvement

In January 2023, the International Laser Ranging Service (ILRS) approved the tracking of 20 additional BeiDou-3 Medium Earth Orbit (BDS-3 MEO) satellites, integrating them into the ILRS tracking network. Before that, only 4 BDS-3 MEO satellites had been tracked. BDS satellites employ highly advanced GNSS components and technological solutions; however, microwave-based orbits still contain systematic errors. Satellite Laser Ranging (SLR) tracking is thus crucial for better identification and understanding of orbit modeling issues. Orbit improvements are necessary to consider BDS in future realizations of terrestrial reference frames, supporting the determination of global geodetic parameters and utilizing them for the co-location of GNSS and SLR in space. In this study, we summarize the first 6 months of SLR tracking 24 BDS-3 MEO satellites. The study indicates that the ILRS network effectively executed the request to track the entire BDS-3 MEO constellation. The number of observations is approximately 1300 and 450 for high- and low-priority BDS-3 satellites, respectively, over the 6 months. More than half of the SLR observations to BDS-3 MEO satellites were provided by 5 out of the 24 laser stations, which actively measured GNSS targets. For 14 out of 24 BDS-3 MEO satellites, the standard deviation of SLR residuals is at the level of 19–20 mm, which is comparable with the quality of the state-of-the-art Galileo orbit solutions. However, the SLR validation of the individual satellites revealed that the BDS-3 MEO constellation consists of more ambiguous groups of satellites than originally reported in the official metadata files distributed by the BDS operators. For 8 BDS-3 satellites, the quality of the orbits is noticeably inferior with a standard deviation of SLR residuals above 100 mm. Therefore, improving orbit modeling for BDS-3 MEO satellites remains an urgent challenge for the GNSS community.

Enhanced precise orbit determination for GPS and BDS-2/3 with real LEO onboard and ground observations

Due to the current geographical situation in China that it is difficult to achieve globally well distributed ground monitoring stations for Chinese BeiDou Navigation System (BDS), the BDS Medium Earth Orbit (MEO) satellite orbital accuracy can be seriously affected by the insufficient length of continuous observation arc based on regional ground network. What is more, the highly static characteristics of BDS Geosynchronous Orbiting (GEO) satellites resulting in meter-level orbit determination accuracy is far from satisfactory for precise positioning services. In order to overcome these problems, the low Earth orbit (LEO) satellites serving as fast-moving space monitoring stations can be used as a good supplement for the current BDS ground monitoring stations layout. In this study, real measured spaceborne data from two LEO satellites, namely Fengyun-3D (FY-3D) and Tianjin University No. 1 (TJU-01), are collected and studied for the augmentation of BDS precise orbit determination (POD). The multipath errors of spaceborne code observations are comprehensively analyzed. It is found that not only BDS-2 relevant satellite-induced code bias but also near-field multipath errors, especially distinct in GPS, exist in a spaceborne environment. The piece-wise linear model is therefore established and applied to correct the code observations to reduce the multipath errors in the integrated POD processing. The GPS and BDS-2/3 POD solutions are compared with orbital overlapping differences and Satellite Laser Ranging (SLR) residuals validation. It is demonstrated that including the two LEO satellites can improve the orbit determination accuracy of GPS and BDS-2/3 significantly under a regional ground network, especially in the along component. The average overlapping orbital Root Mean Square (RMS) is decreased by 18.4 %, 34.3 %, 24.9 %, 21.4 %, and 7.8 % for the GPS, GEO, Inclined Geosynchronous Orbit (IGSO), BDS-2 MEO, and BDS-3 MEO satellites, respectively, in the along component. The RMS of SLR residuals is decreased from 37.1 cm to 17.8 cm for the GEO satellites, and the overall improvement is about 8 %-12 % for the BDS MEO satellites. There has also been a slight improvement under a global network, with an average improvement of 2–3 cm in the SLR residuals RMS. However, significant systematic errors are also observed in the SLR residuals in BDS GEO, IGSO, and MEO satellites, and are expected to be improved with refined observation models based on more LEO satellites, i.e., the phase center variation (PCV).

Combination and SLR validation of IGS Repro3 orbits for ITRF2020

In preparation for the International Terrestrial Reference Frame 2020, the International GNSS Service analysis centers released the results of the third reprocessing campaign (IGS Repro 3) of all the GNSS network solutions backwards starting from 1994. For the first time, the IGS reprocessing products included not just GPS and GLONASS, but also the Galileo constellation. In this study, we describe the methodology and results of the orbit combination provided by the IGS Analysis Center Coordinator (IGS ACC) at Geoscience Australia. The quality of the combined orbit products was cross-checked with the individual IGS Repro3 Analysis Center (AC) contributions. The internal consistency of the individual Analysis Center (AC) solutions with the combined orbits was assessed based on the root mean square of the 3D orbit differences. In 2020, the mean consistency of the combination is at the level of 9, 23, and 15 mm for GPS, GLONASS, and Galileo, respectively. The external validation of the orbits was performed using Satellite Laser Ranging (SLR). We proposed a novel approach to handling detector-specific biases in the results of SLR validation, which reduced the standard deviation of SLR residuals by up to 13% for Galileo FOC satellites. This method is based on bias aligning the offsets to single-photon SLR stations that were treated as a reference. The proposed approach increased the internal consistency of the SLR dataset, facilitating the detection of orbit modeling issues. The standard deviation of SLR residuals of the best individual solution versus the combined solution equals 13/13, 15/17, 17/17, 18/19 mm for Galileo-FOC, -IOV, GLONASS-K1B, -M, respectively. Therefore, the combined solution can be considered equal in quality compared to the best individual AC solutions. Searching for patterns in SLR residuals for different satellite-Sun-Earth geometries revealed that some orbit modeling issues are not fully diminished for individual ACs. Eventually, our findings suggest that the delivered combined orbit product may be considered the best solution overall, as it benefits from the best individual solutions for each satellite type.

Estimating inter-satellite link Ka-band antenna phase center offsets to improve BDS-3 satellite orbit determination

Satellite antenna phase center offset (PCO) correction is essential not noly for L-band ground-tracking GNSS high-precision data processing, but also for Ka-band inter-satellite link (ISL) observations. The Ka-band ISL and L-band satellite-ground observations of 90 days from 25 November 2020 to 23 February 2021 are used for the Ka-band PCO estimation and precise orbit determination (POD) validation for BDS-3 satellites. The results show that the standard deviations of the estimated Ka-band PCOs for 29 BDS-3 satellites are 1.73, 2.63 and 3.58 cm in the X, Y and Z components, respectively. The accuracy of the Ka-band PCO estimates for medium Earth orbit satellites (MEO) is better than those of inclined geosynchronous orbit (IGSO) and geostationary Earth orbit (GEO) satellites due to a better observation geometry. The POD validation results indicate that using the estimated Ka-band PCOs can improve the ISL-based orbit determination performance. Using Ka-band ISL and L-band satellite-ground observations with 40 global stations for combined POD, the orbit overlap comparison shows that the average three-dimensional (3D) root mean square (RMS) error for 29 BDS-3 satellites is decreased from 12.6 to 6.6 cm with a reduction of 47.6%, and satellite laser ranging (SLR) validation for 4 BDS-3 satellites shows an accuracy improvement from 6 to 5 cm in RMS. When using ISL and only 7 regional ground stations are used, the average 3D RMS of orbit overlap errors is decreased from 15.9 to 10.9 cm with a reduction of 31.6%, and the RMS of SLR residuals for 4 BDS-3 satellites is decreased from 6.7 to 5.6 cm. For ISL-only POD, the average 3D RMS of orbit overlap errors is decreased from 18.0 to 12.4 cm with a reduction of 31.1%, and the RMS of SLR residuals for 4 BDS-3 satellites is decreased from 7.1 to 6.9 cm.

Assessment of Swarm Kinematic Orbit Determination Using Two Different Double-Difference Methods

The Swarm mission aims to study the principle and change regularities of the Earth’s magnetic field. Precise orbit determination is essential to the successful implementation of the mission and relevant scientific research. This article focuses on using two different double-difference methods to improve the accuracy of Swarm kinematic orbit determination. The accuracy of the kinematic orbit determination relies entirely on the space-borne observation data, independent of any dynamic parameters. The article analyzes the data quality of the Swarm space-borne global positioning system (GPS) receiver and presents a detailed introduction to the data pre-processing algorithms. The double-difference observation gathering and the applied orbit determination strategy using two different double-difference methods are discussed. The results of the kinematic orbits under different solar cycle conditions are presented, along with an evaluation based on analysis of GPS carrier phase residuals, subtracting from the post-processed orbits, and assessment with satellite laser ranging (SLR) measurements. The results show that the accuracy of the kinematic orbit determination is at the centimeter level for the three Swarm satellites’ orbit solutions. The daily root mean square (RMS) values of the three satellites’ phase residuals remain at around the 6 mm level. The RMS values of the position residuals between the kinematic orbits and the reduced dynamic orbits released by the European Space Agency (ESA) are at about the 2–3 cm level. The external evaluation with SLR measurements shows a good agreement with the ESA level, with the RMS values of the SLR residuals for kinematic orbits around 2 cm.

Improving Precise Orbit Determination of LEO Satellites Using Enhanced Solar Radiation Pressure Modeling

AbstractPrecise orbit knowledge is a fundamental requirement for low Earth orbit (LEO) satellites. High‐precision non‐gravitational force modeling directly improves the overall quality of LEO precise orbit determination (POD). To address the potential systematic errors in solar radiation pressure (SRP), we introduce observed radiation data and modeled physical effects to describe the real in‐flight environment of satellites. Time‐dependent solar irradiance data and a highly physical shadow model are considered for SRP modeling. We develop an advanced thermal reradiation model for satellite solar panels. A set of improved non‐gravitational force models is performed for LEO POD, and we discuss the benefits of the enhanced dynamic models on orbit quality and dependence on empirical parameters. The Gravity Recovery and Climate Experiment Follow‐On (GRACE‐FO), Jason‐3, and Haiyang‐2B missions are selected for the POD process. Estimated empirical acceleration and scale parameters and independent satellite laser ranging (SLR) are used to validate the final orbit solutions. The magnitude of empirical acceleration estimated in POD is reduced by 19% with the enhanced dynamic modeling, and the estimated scale factor for the SRP converges to stable and reasonable level. Furthermore, the steady‐state temperature model used in thermal reradiation can effectively reduce mismodeled effects in the SRP, and the systematic linear dependency revealed by the SLR residuals is significantly reduced for the GRACE‐C and Jason‐3 satellites, with improvements of approximately 61% and 49%, respectively. Overall, advances are made in the explicit modeling of non‐gravitational forces to pursue superior satellite orbits, suggesting a more dynamic orbit solution.

Precise orbit determination and baseline consistency assessment for Swarm constellation

Precise orbit determination (POD) and precise baseline determination (PBD) of Swarm satellites with 4 years of data are investigated. Ambiguity resolution (AR) plays a crucial role in achieving the best orbit accuracy. Swarm POD and PBD based on single difference (SD) AR and traditional double difference (DD) AR methods are explored separately. Swarm antenna phase center variation (PCV) corrections are developed to further improve the orbit determination accuracy. The code multipath of C1C, C1W and C2W observations is first evaluated and clear variations in code noise related to different receiver settings are observed. Carrier phase residuals of different time periods and different loop tracking settings of receiver are studied to explore the effect of ionospheric scintillation on POD. The reduction of residuals in the polar and geomagnetic equator regions confirms the positive impact of the updated carrier tracking loops (TLs) on POD performance. The SD AR orbits and orbits with float ambiguity (FA) are compared with the Swarm precise science orbits (PSOs). An average improvement of 27 %, 4 % and 16 % is gained in along-track, cross-track and radial directions by fixing the ambiguity to integer. For Swarm-A/B and Swarm-B/C formations, specific days are selected to perform the DD AR-based POD during which the average distance of the formation satellites is less than 5000 km. Satellite laser ranging (SLR) observations are employed to validate the performance of FA, SD AR and DD AR orbits. The consistency between the SD AR orbits and SLR data is at a level of 10 mm which shows an improvement of 25 % when comparing with the FA results. An SLR residuals reduction of 15 % is also achieved by the DD AR solution for the selected days. Precise relative navigation is also an essential aspect for spacecraft formation flying missions. The closure error method is proposed to evaluate the baseline precision in three dimensions. A baseline precision of 1–3 mm for Swarm-A/C formation and 3–5 mm for Swarm-A/B and Swarm-B/C satellite pairs is verified by both the consistency check and closure error method.

Enhanced precise orbit determination of BDS-3 MEO satellites based on ambiguity resolution with B1C/B2a dual-frequency combination

Broadcasting the multi-frequency signals is one of the outstanding advantages of BDS-3 (BeiDou Navigation Satellite System). In addition to the legacy signals B1I and B3I, we are interested in the performance of two new signals B1C and B2a. For the purpose of improving the orbit accuracy of BDS-3 backbone MEOs (Medium Earth Orbits), this study focuses on Precise Orbit Determination (POD) using B1C/B2a dual-frequency measurements. First, the performance of POD will be benefit from an increasing observation quantity and qualities. Observations from the public MGEX (Multi-GNSS Experiment) network are statistically collected, and the number of stations able to track B1C and B2a signals is analyzed. In recent years, the sustainable growth number of stations is up to more than one hundred, which brings the new chance to improve the POD accuracy based on B1C/B2a linear combination. Then, the multipath and post-fit residuals of observations are compared between B1C/B2a and B1I/B3I to investigate the observation qualities. Second, as a key strategy, the ambiguities are successfully fixed during POD processing of BDS-3 MEOs. The ambiguity fixing percentage is as high as 93.1% on average using B1C/B2a combination. Compared with float solutions, the improvements of ambiguity fixed orbits reach 36.2%, 19.9% and 36.6% on the radial, along-track and cross-track components, respectively. For the fixed solutions of B1I/B3I and B1C/B2a, it is noticed that the orbits achieve further benefits from B1C/B2a combination for 8.4%, 6.9% and 8.0% on the radial, along-track and cross-track components, respectively. Finally, the performance of B1C/B2a combination POD is validated by orbital overlapping, the orbit signal-in-space range error (SISRE) analysis, independent technology Satellite Laser Ranging (SLR) and Precise Point Positioning (PPP), respectively. SLR residuals show the radial orbital accuracy has an improvement of 8.0% for C20, C21, C29 and C30 on average when employing B1C/B2a measurements. Compared with B1I/B3I, both static and kinematic PPP based on B1C/B2a orbits indicate the improvement on the vertical component. And a faster and more stable convergence performance is also observed using B1C/B2a combination.

Station-dependent satellite laser ranging measurement corrections for TOPEX/Poseidon

The TOPEX/Poseidon (T/P) altimetry mission with its main objectives to monitor variations of the global and regional sea level as well as ocean circulation is a milestone in Earth observation. The spacecraft was launched in 1992 and it is the predecessor mission of the Jason series and Sentinel-6A. For laser ranging measurements from the Earth to the spacecraft, T/P was equipped with a non-ideally designed annular retroreflector array. Its large dimensions of over 160 cm in diameter caused huge optical phase centre variations which limit the orbit accuracy. In this study, we developed a continuous, analytical correction function that counteracts both, phase centre variations at the spacecraft and station-related laser ranging measurement errors such as range biases. For most Satellite Laser Ranging (SLR) ground stations that tracked the T/P spacecraft, an individual set of six correction parameters is estimated. The parameters are valid for the entire mission or, in some cases, for a defined period. The developed function uses the observation’s viewing angles to determine a correction value which is added to the range measurement. Applying the measurement correction reduces the overall T/P mission root mean square fit of SLR residuals from 33.78 cm to 1.97 cm (1.59 cm for SLR core stations). External orbits based on the joint analysis of two different space-geodetic techniques are used to validate the quality of our improved orbit solution. The comparisons show good agreement. The mean values of the radial, transverse, and normal components differ by 0.0,-0.1, and -0.2 cm, respectively. To investigate the impact of the corrected orbit on sea level computations, a single-satellite crossover analysis is performed. When using the corrected measurements, the standard deviation of crossover differences reduces from 16 cm to 6 cm. The computed coefficients of the measurement correction function are provided publicly and can be used in any software package to obtain T/P orbits.

Estimating surface optical properties and thermal thrust for Galileo satellite body and solar panels

Precise orbit determination of GNSS (Global Navigation Satellite System) satellites requires accurate models of perturbing forces acting on the spacecraft, i.e., solar radiation pressure (SRP) and thermal radiation forces. With the officially published satellite metadata, the analytical box-wing model is usually used to describe most of the SRP accelerations and the rest is assumed to be compensated by estimating ECOM/ECOM2 (Empirical CODE Orbit Model) parameters. However, we find that the precision of Galileo satellite orbits shows notable degradation inside eclipse seasons for 3-day-arc solutions and 24-h predictions. For instance, the RMS (root-mean-square) of orbit misclosures increases by about a factor of two in the eclipse season when using the box-wing model as the a priori and the 5-parameter ECOM model on top. The reason is proven to be mostly due to ignoring imbalanced thermal radiation forces (i.e., radiator emission and thermal radiation of solar panels) as satellite thermal properties are unknown. These imbalanced thermal effects cannot be fully absorbed by the ECOM/ECOM2 parameters inside eclipse seasons because the earth’s shadowing of a satellite in orbit causes periodic changes of the thermal environment. To cope with this problem, we first estimate satellite optical and thermal parameters as part of orbit determination based on Galileo tracking data covering 1 year. Then, we add physical thermal radiation models for radiators and solar panels as part of the a priori model and evaluate the performance of different ECOM models in Galileo satellite orbit determination. As shown by orbit misclosures, 24-h orbit predictions and SLR (Satellite Laser Ranging) residuals, the 7-parameter ECOM2 model performs better than the 5-parameter ECOM and the 9-parameter ECOM2 model for Galileo satellites. When using the 7-parameter ECOM2 model on top, the impact of the radiator emission and the thermal radiation of solar panels on Galileo satellite orbits is about 1 and 2 cm, respectively, inside eclipse seasons for 3-day-arc solutions.

BeiDou-3 orbit and clock quality of the IGS Multi-GNSS Pilot Project

Within the Multi-GNSS Pilot Project (MGEX) of the International GNSS Service (IGS), precise orbit and clock products for the BeiDou-3 global navigation satellite system (BDS-3) are routinely generated by a total of five analysis centers. The processing standards and specific properties of the individual products are reviewed and the BDS-3 orbit and clock product performance is assessed through direct inter-comparison, satellite laser ranging (SLR) residuals, clock stability analysis, and precise point positioning solutions. The orbit consistency evaluated by the signal-in-space range error is on the level of 4–8 cm for the medium Earth orbit satellites whereas SLR residuals have RMS values between 3 and 9 cm. The clock analysis reveals sytematic effects related to the elevation of the Sun above the orbital plane for all ACs pointing to deficiencies in solar radiation pressure modeling. Nevertheless, precise point positioning with the BDS-3 MGEX orbit and clock products results in 3D RMS values between 7 and 8 mm.

Satellite laser ranging to GNSS-based Swarm orbits with handling of systematic errors

Satellite laser ranging (SLR) retroreflectors along with GNSS receivers are installed onboard numerous active low earth orbiters (LEOs) for the independent validation of GNSS-based precise orbit determination (POD) products. SLR validation results still contain many systematic errors that require special handling of various biases. For this purpose, we derive methods of reducing systematic effects affecting the SLR residuals to LEO Swarm satellites. We test solutions incorporating the estimation of range biases, station coordinate corrections, tropospheric biases, and horizontal gradients of the troposphere delays. When estimating range biases once per day, the standard deviation (STD) of Swarm-B SLR residuals is reduced from 10 to 8 mm for the group of high-performing SLR stations. The tropospheric biases estimated once per day, instead of range biases, further reduce the STD of residuals to the level of 6 mm. The systematic errors that manifest as dependencies of SLR residuals under different measurement conditions, e.g., elevation angle, are remarkably diminished. Furthermore, introducing troposphere biases allows for the comparison of the orbit quality between kinematic and reduced-dynamic orbits as the GPS-based orbit errors become more pronounced when SLR observations are freed from elevation-dependent errors. Applying tropospheric biases in SLR allows obtaining the consistency between the POD solution and SLR observations that are two times better than when neglecting to model of systematic effects and by 29% better when compared with solutions considering present methods of range bias handling.

Validation and evaluation on B1IB3I-based and B1CB2a-based BDS-3 precise orbits from iGMAS

The quality analysis of precise orbit products are important for the BeiDou global navigation satellite system (BDS-3) to provide precise positioning, navigation and timing (PNT) services to users. The International GNSS Monitoring and Assessment System (iGMAS) is simultaneously been provided with B1IB3I-based and B1CB2a-based BDS-3 precise orbits since September 29, 2019. However, a comprehensive comparison and analysis of these two types of orbit products is still missing. In this contribution, we study the performance of these precise orbits by analyzing the day boundary discontinuity and Satellite Laser Ranging (SLR). Moreover, we compare these precise orbits with the BDS-3 orbit products from the multi-GNSS Experiment (MGEX) analysis center. The test period is more than two years from September 29, 2019 to October 21, 2021. We found that: First, compare them with the MGEX BDS-3 precise orbits provided by GFZ and WHU, the average root mean square (RMS) value of the position deviation in the 3D directions orbit component is 0.081 ∼ 0.154 m for MEO satellites and 0.202 ∼ 0.255 m for IGSO satellites. The ranks of different BDS-3 orbit consistency are MEO SECM > MEO CAST > IGSO and we observe that the orbit consistency in eclipse period is better than the ones in non-eclipse period expect for rare cases; Secondly, the RMS of SLR residuals is 5.0 ∼ 6.0 cm, and compared to the orbit in non-eclipse period, the decrease of orbit in eclipse period in percentage is 1.72%∼8.00%. In addition, the radial performance of B1IB3I-based precise orbits for MEO CAST satellites is better than that of B1CB2a-based precise orbits, but it is opposite for MEO SECM satellites; Thirdly, combine with the CODE’s new solar radiation pressure model, namely ECOM2, by using the day boundary discontinuity, the preliminary results illustrate that the average RMS value of the day boundary discontinuity for BDS-3 MEO satellites in the 3D directions orbit component is 0.546 ∼ 0.639 m, and it is 0.881 ∼ 0.897 m for BDS-3 IGSO satellites. The performance ranking of day boundary discontinuity of BDS-3 precise orbits in the 3D directions orbit component is B1IB3I-based MEO CAST > B1IB3I-based MEO SECM > B1CB2a-based MEO CAST > B1CB2a-based MEO SECM > B1CB2a-based IGSO > B1IB3I-based IGSO satellites.

Evaluation of strategies for the ultra-rapid orbit prediction of BDS GEO satellites

ABSTRACT The quality of BeiDou Navigation Satellite System (BDS) Geostationary Earth Orbit (GEO) ultra-rapid products is unsatisfactory because GEO satellites are nearly stationary relative to ground stations. To optimize the quality of these ultra-rapid orbit products, we investigated the effects of the fitting arc length, an a priori Solar-Radiation Pressure (SRP) model, and the along-track empirical acceleration on the prediction of BDS GEO satellite orbits. The predicted orbit arcs of 24-h were evaluated through comparisons with the corresponding observed orbit arc and Satellite Laser Ranging (SLR) observations. In both eclipse and non-eclipse seasons, accuracy of the orbit predictions obtained using a 48-h fitting arc length were better than those obtained using 24-h and 72-h fitting arc lengths. Although the overlapping precision of predicted orbits exhibited no obvious improvement when an a priori SRP model was employed, the systematic bias in the SLR residuals was significantly reduced. Specifically, the mean value of SLR residuals decreased from −0.248 m to −0.024 m during non-eclipse seasons and from −0.333 m to −0.041 m during eclipse seasons, respectively. In addition, when an empirical acceleration in the along-track direction was introduced, the three-Dimensional Root-Mean-Square (3D RMS) of overlapping orbits during eclipse seasons decreased from 2.964 to 1.080 m, which is comparable to that during non-eclipse seasons. Furthermore, the Standard Deviation (STD) of SLR residuals decreased from 0.419 to 0.221 m during eclipse seasons. The analysis of SRP estimates shows that the stability of SRP parameters was significantly enhanced after the introduction of along-track empirical acceleration in eclipse seasons. The optimal BDS GEO ultra-rapid orbit prediction products were yielded by using a 48-h fitting arc length, an a priori SRP model and an along-track empirical acceleration.

COST-G gravity field models for precise orbit determination of Low Earth Orbiting Satellites

Accurate time-variable Earth’s gravity field information is important for precise orbit determination (POD) of low Earth orbiters (LEO). In particular the POD of altimeter satellites benefits from accurate modelling of the time-variable gravity field. The Combination Service for Time-variable Gravity Fields (COST-G) provides monthly gravity fields based on a combination of GRACE/GRACE-FO derived monthly gravity field solutions from different analysis centers. These monthly solutions are available with a latency of 2–3 months. The monthly GRACE-FO gravity fields from June 2018 to June 2021 serve as base for a fitted signal model (FSM) of time-variable gravity that enables a few months prediction and therefore application in operational LEO POD. The combined monthly solutions themselves, as well as several variants of the FSM are used for GPS-based POD of Copernicus Sentinel-2 and −3 satellites flying at altitudes of 786km and 814,5km, respectively. The orbit results show a significant improvement when using both, the COST-G monthly and the new COST-G FSM gravity field models compared to results using existing long-term static gravity field models, including seasonal and trend estimates, as well. RMS values of the GPS carrier phase residuals are reduced up to 20%, orbit overlaps between 30–50%. Orbit validations performed by Satellite Laser Ranging (SLR) also show that SLR residuals are reduced by about 10%. The COST-G FSMs show a similar performance as the monthly gravity fields. The LEO POD in particular benefits from the more realistic mass trend estimates in river basins with strong non-periodic inter-annual variations compared to the predicted trends of outdated long-term gravity fields. The FSM thus complements the monthly gravity fields as a new COST-G product. It is foreseen to update the COST-G FSM quarterly.

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.