Inclined Geosynchronous Orbit Research Articles

Orbit Determination of High Orbit Spacecraft based on BDS-3 using Extended Kalman Filter

The positioning of highly elliptical orbit (HEO) spacecraft is very important in deep space exploration and other missions. The Beidou Navigation Satellite System (BDS-3), which was completed in 2020, has put forward new ideas for using GNSS technology to make up for the lack of earth-based TT&C systems. At the same time, the selection of the gravity field order and the perturbation forces must be considered when using dynamical methods for HEO spacecraft positioning. In this paper, simulation experiments are carried out based on the Two-Line Orbital Element (TLE) data of BDS. Two HEO satellites with different semimajor-axis (67509.6 km, 212857.3 km) are selected as the target spacecraft orbits. The single point positioning (SPP) method and Extended Kalman Filter (EKF) method are used for HEO spacecrafts positioning respectively. The experiments show that a 20×20 order gravity field model can be chosen for the orbital integral of HEO spacecraft orbits, and the integration error caused by solar gravity and lunar gravity should be taken into account with the change of orbital altitude. Adding Geostationary Orbit (GEO) and Inclined GeoSynchronous Orbit (IGSO) navigation satellites improved the visibility of these two orbits by 19.45% and 16.64%, respectively, and improved the GDOP value by 14.1% and 3.8%, respectively, compared to observing MEO satellites only. The use of the EKF method can effectively improve the positioning accuracy of the HEO spacecraft compared to the SPP method. For the spacecraft in HEO-1, the RMS value of the positioning errors of the EKF method was 67.82 m, and the accuracy of each direction was improved by 44.28%, 38.85%, and 38.62%, respectively. For the spacecraft in HEO-2, the positioning errors in the x direction of EKF method was within 12 km and within 2 km in the y and z directions. The accuracy of these three directions was improved by 84.91%, 79.24%, and 58.64%, respectively.

Effects of Equatorial Plasma Bubbles on Multi-GNSS Signals: A Case Study over South China

Equatorial plasma bubbles (EPBs) occur frequently in low-latitude areas and have a non-negligible impact on navigation satellite signals. To systematically analyze the effects of a single EPB event on multi-frequency signals of GPS, Galileo, GLONASS, and BDS, all-sky airglow images over South China are jointly used to visually determine the EPB structure and depletion degree. The results reveal that scintillations, or GNSS signal fluctuations, are directly linked to EPBs and that the intensity of scintillation is positively correlated with the airglow depletion intensity. The center of the airglow depletion often corresponds to stronger GNSS scintillation, while the edge of the bubble, which is considered to have the largest density gradient, corresponds to relatively smaller scintillation instead. This work also systematically analyzes the responses of multi-constellation and multi-frequency signals to EPBs. The results show that the L2 and L5 frequencies are more susceptible than the L1 frequency is. For different constellations, Galileo’s signal has the best tracking stability during an EPB event compared with GPS, GLONASS, and BDS. The results provide a reference for dual-frequency signal selection in precise positioning or TEC calculation, that is, L1C and L2L for GPS, L1C and L5Q for Galileo, L1P and L2C for GLONASS, and L1P and L5P for BDS. Notably, BDS-2 is significantly weaker than BDS-3. And inclined geosynchronous orbit (IGSO) satellites have abnormal data error rates, which should be related to the special signal path trajectory of the IGSO satellite.

Orbit determination of Earth-Moon libration point navigation constellation based on Inter-satellite links

Earth-Moon libration point navigation constellation has a merit of fully covering the cislunar space and providing navigation services with merely a small number of satellites. Furthermore, this type of constellation itself can achieve orbit determination by using only its Inter-satellite links (ISL) owing to the well-known gravitational asymmetry near the libration points. This paper adopts a representative four-satellite navigation constellation including three satellites located near the libration point L3, L4, L5, and an extra one in Distant Retrograde Orbit (DRO), and studies the orbit determination accuracy of the constellation under two conditions: establishing ISLs within itself (called Solely libration point satellites) and establishing additional ISLs with BeiDou Navigation Satellite System (BDS) (called Libration + BeiDou satellites). The simulation results indicate that, for the Solely libration point satellites scenario with the 1 m range error and solar radiation pressure (SRP) error with the level of 10% deviation between ball model and macro model, the final orbit accuracies for L3, L4, L5, and DRO satellites can respectively achieve 11.3 m, 12.7 m, 12.6 m, and 7.3 m. Reducing the link interval can significantly shorten the arc length to achieve final accuracy, whileas it has less impact on the final orbital accuracy. When the link interval is set as 240, 60, 10, and 2 min, the arc lengths at the final accuracy of better than 15 m are about 35, 29, 23, and 20 days, respectively. To balance the communication burden and the orbit determination arc length, the recommended ISL building interval is 10 min. Errors coming from range measurement model and dynamic model, in addition, are key factors in orbit determination, which can result in meters to tens of meters of orbital accuracy. For the Libration + BeiDou satellites scenario, the augmented ISLs to BDS can substantially improve constellation orbital accuracy and shorten the orbit determination arc length. When the DRO satellite attached to the Geostationary Earth Orbit (GEO), Inclined Geosynchronous Orbit (IGSO), and Medium Earth Orbit (MEO) satellites of BDS respectively, the improved accuracy are 37.1%, 61.2%, and 42.2%, while the shortened arc length are 20.1%, 28.3%, and 28.1%, respectively. The research findings presented in this paper can serve as a reference for the construction and assessment of navigation constellations in the Earth-Moon system.

Gravity field recovery of inter-satellite links between Beidou navigation satellite system (BDS) and LEO based on geodesy and time reference in space (GETRIS)

The goal of this contribution is to construct the inter-satellite links (ISLs) between Global Navigation Satellite System (GNSS) and Low Earth Orbit (LEO) satellite to recover the Earth’s gravity field based on the concept of Geodesy and Time Reference in Space (GETRIS). Taking the constellation of the Beidou Navigation Satellite System (BDS) as an example, several types of high-low satellite-to-satellite tracking (H-L SST) modes combining BDS and GETRIS are used to research the tracking mode of high orbit satellites and the distribution of LEO satellite constellations to explore the optimal mode for gravity field recovery. The effects of multiple H-L SST modes on the accuracy of gravity field recovery are studied. In addition, the contribution of each error source is also analyzed. Finally, this study varies LEO satellite constellations to show reduced errors in the temporal gravity field solutions for the improved ISLs. The results showed that the Inclined GeoSynchronous Orbit (IGSO) satellite greatly improved the Geosynchronous orbit (GEO) satellite tracking LEO satellite mode in the high latitude regions. Compared with the Medium Earth Orbit (MEO) satellite tracking LEO satellite, the error of gravity field recovered by combining with GEO, IGSO and MEO can be reduced by 16 %. The largest error is the non-tidal signal, followed by the ocean tidal model difference, and the smallest is the instrument error. For a single LEO satellite, the combination model is more sensitive to the satellite orbital altitude, the cumulative error was decreased from 5.08 mm (550 km) to 0.33 mm (250 km). For LEO constellations, adding inclined orbit satellites can improve the accuracy of the solution to some extent. But the optimal solution can be obtained by the constellation composed of polar orbit satellites with different Right Ascension of Ascending Node (RAAN), the cumulative errors are 0.27 mm (two satellites), 0.19 mm (three satellites), and 0.19 mm (four satellites).

Investigating the Global Performance of the BDS-2 and BDS-3 Joint Real-Time Undifferenced and Uncombined Precise Point Positioning Using RTS Products from Different Analysis Centers

Compared to the traditional ionospheric-free (IF) precise point positioning (PPP) model, the undifferenced and uncombined (UU) PPP has the advantages of lower observation noise and the ability to obtain ionospheric information. Thanks to the IGS (International GNSS Service), real-time service (RTS) can provide RT vertical total electron content (VTEC) products, and an enhanced RT UU-PPP based on the RT-VTEC constraints can be achieved. The global performance of the BeiDou Navigation Satellite System-2 (BDS-2) and BDS-3 joint RT UU-PPP using different RTS products was investigated. There is not much difference in the RTS orbit accuracy of medium earth orbit (MEO) satellites among all analysis centers (ACs), and the optimal orbit accuracy is better than 5, 9, and 7 cm in the radial, along-track, and cross-track directions, respectively. The orbit accuracy of inclined geosynchronous orbit (IGSO) satellites is worse than that of MEO satellites. Except for CAS of 0.46 ns, the RTS clock accuracy of MEO satellites for other ACs achieves 0.2–0.27 ns, and the corresponding accuracy is about 0.4 ns for IGSO satellites. In static positioning, due to the limited accuracy of RT-VTEC, the convergence time of the enhanced RT UU-PPP is longer than that of RT IF-PPP for most ACs and can be better than 25 and 20 min in the horizontal and vertical components, respectively. After convergence, the 3D positioning accuracy of the static RT UU-PPP is improved by no more than 8.7%, and the optimal horizontal and vertical positioning accuracy reaches 3.5 and 7.0 cm, respectively. As for the kinematic mode with poor convergence performance, with the introduction of RT-VTEC constraints, the convergence time of RT UU-PPP can be slightly shorter and reaches about 55 and 60 min in the horizontal and vertical components, respectively. Both the horizontal and vertical positioning accuracies of the kinematic RT UU-PPP can be improved and achieve around 7.5 and 10 cm, respectively.

Analysis of BDS inter-satellite link ranging performance

The Ka-band inter-satellite link (ISL) is a key technology of BeiDou Global Navigation Satellite System (BDS-3) that can improve the system performance. Based on 30 days of ISL observations, we analyze the BDS-3 ISL ranging performance and the related influencing factors, such as satellite manufacturer, satellite orbit type and link type. The fitting residuals of the ISL geometry-free observations are used to characterize the measurement noise, and the post-fit ISL residuals from precise orbit determination (POD) and precise clock estimation (PCE) are used to further validate the ranging performance. The fitting residuals follow a normal distribution, and the average root mean square (RMS) value for all 421 links is approximately 2.5 cm. The satellite manufacturers, China Academy of Space Technology (CAST) and Shanghai Engineering Center for Microsatellites (SECM), have an obvious impact on the ISL measurement noise. The RMS of the ISL fitting residuals for the links established by the two CAST satellites is only 1.45 cm, while it can reach 3.34 cm for the links established by the two SECM satellites. This correlation between ISL ranging performance and the satellite manufacturer is also found in the POD and PCE post-fit ISL residuals. The RMS of PCE residuals is about 3.1 cm, while the POD residuals are larger with the RMS of 5.5 cm due to the effect of orbital errors, especially for the Inclined Geosynchronous Orbit (IGSO) and Geostationary Orbit (GEO) satellites. Overall, the ranging accuracy of the BDS-3 ISL is better than 5 cm according to the above analyses. Furthermore, the BDS-3 satellite orbit and clock results are evaluated by internal consistency check and external comparison to analyze the application of ISL observations.

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).

A Study of Solar Radiation Pressure Model in BDS-3 Precise Orbit Determination

For the BDS (BeiDou Navigation Satellite System) satellites, the solar radiation pressure is the major non-gravitational perturbation the satellites suffer. Due to the influence of multiple factors, accurate modeling the solar radiation pressure acting on these satellites is difficult to accomplish. The solar radiation pressure perturbation is an essential source of error in the precise orbit determination (POD) and orbit prediction processes of the BDS satellites. As the ECOMC (Empirical CODE Orbit Model 1 and 2 Combined) model adopts the advantageous features of the ECOM1 (Empirical CODE Orbit Model 1) model and the ECOM2 (Empirical CODE Orbit Model 2) model, a greater number of parameters have been introduced into the model, leading to strong correlation between the ECOMC model parameters. In response to the deficiency of the ECOMC model, the authors collected the BDS-3 satellite precise ephemeris data that were provided by Wuhan University (WHU) from January 2019 to April 2022. And by these ephemeris data, the ECOMC model’s 13 parameters were obtained through the dynamic orbit fitting method. By conducting these time sequence diagrams analysis of the ECOMC model parameters, this thesis affords the selection strategies of the solar radiation pressure parameters of the ECOMC model for BDS-3 inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellite. And the rationality of the selection strategies of the solar radiation pressure parameters for the ECOMC model was verified by orbit fitting and orbit prediction experiments. The results show that using the improved ECOMC model gains the best orbit fitting effects for BDS-3 IGSO and MEO satellite, and can also enhance the accuracy of medium and long-term orbit prediction of BDS-3 IGSO and MEO satellite.

An improved satellite-induced code bias correction model for BDS un-GEO satellites in consideration of elevation and azimuth angles

The satellite-induced code bias (SICB) resulting in code-phase divergences of nearly 1.0 m has been observed in the code measurements of BeiDou navigation satellite system (BDS). Since the SICB is proved to be elevation- and frequency-dependent, several SICB correction models have been proposed to remove the systematic bias, but most of them only consider elevation angles, overlooking the influence of azimuth angles, which would limit the accuracy of the correction models, especially for Medium Earth Orbit (MEO) satellites. To deal with this problem, we proposed an improved SICB piecewise correction model as a function of elevation and azimuth angles and estimated parameters for each BDS un-GEO (Geostationary Earth Orbit) satellite and each available frequency with node-separation of 1°. The epoch-difference method is employed to remove the effect of stable terms, such as phase-ambiguities, receiver- and satellite-dependent hardware delays. Moreover, to reduce the effect of code noise and multipath at low-elevation conditions, the threshold of elevation angle for the improved model is determined based on the statistics of Multipath (MP) combinations. Besides, weight-average algorithm is used to obtain the optimal MP deviations and improve the accuracy of proposed correction model in consideration of signal-to-noise (SNR). The experiment results with real GNSS data show that the SICB models are more suitable for BDS-2 MEO satellites than Inclined Geosynchronous Orbit (IGSO) satellites. The code measurements on B1 signal are more susceptible to the SICB than that on B3 signal. After applying the improved SICB correction models, the RMS errors of MP series are significantly reduced for BDS-2 un-GEO satellites. However, for BDS-3 satellites, only centimeter-level SICB deviations can be detected for all frequency-bands, and the RMS of MP variations have slight changes after applying the SICB corrections, indicating that the systematic bias in MP combinations is not obvious, especially for BDS-3 MEO satellites.

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.

The Extraction and Characterization of Pseudorange Multipath Based on BDS-3 Multi-Frequency Observations.

Global Navigation Satellite System (GNSS) observations are subject to various errors during their propagation process. A reasonable correction of these errors can improve the positioning, navigation, and timing (PNT) service capability. The impact of multipaths on pseudorange observations can reach a decimeters or even meters level. However, their mechanism is complex and there is currently no universally accepted high-precision correction model. The correlation between the pseudorange multipaths (MP) of BDS-2 satellites and satellite elevation has been confirmed, while there have been fewer analyses of the MP characteristics for different frequencies of BDS-3 satellites. The broadcasting of multi-frequency observations in BDS-3 should theoretically make the extracted MP more accurate compared to traditional methods. Based on this, in this contribution, a multi-frequency MP extraction algorithm based on the least squares principle is proposed, which can simultaneously eliminate the influence of higher-order ionospheric delay. The analytical expression for only eliminating first-order ionospheric delay is successfully derived. Subsequently, the characteristics of the MPs extracted from different frequency combinations and the impact of combination noise on the extraction accuracy are discussed. The influence of second-order ionospheric delay on the MPs for each frequency under different combination noises, as well as the periodic behavior exhibited in long-term observations of the BDS-3 medium earth orbit (MEO) and inclined geosynchronous orbit (IGSO) satellites, are also analyzed. Finally, the correlations between the MPs of each frequency of BDS satellite and elevation are quantitatively analyzed based on observations from 35 stations. Overall, this work has positive implications for the study of the MP characteristics of BDS-3 and subsequent modeling efforts.

Low-Cost BDS Reflectometry for Real-Time Water Surface Retrieval

The official launch of the Chinese BeiDou Navigation Satellite System with global coverage (BDS-3) presents significant opportunities for various applications, including precision agriculture and autonomous driving, among others. With its global spatial coverage and hybrid space constellation comprising geosynchronous Earth orbit (GEO), inclined geosynchronous orbit (IGSO), and medium Earth orbit (MEO) satellites, BDS can significantly contribute to various GNSS remote sensing applications that require real-time, precise water surface height measurements with high temporal and spatial resolution, such as in tidal monitoring. In this paper, we propose a carrier-phase-based method for BDS Reflectometry (BDS-R) to precisely retrieve real-time water surface height. Firstly, the BDS-R altimetry method is introduced, along with a detailed explanation of the data processing procedures. Secondly, a quality control method tailored to the characteristics of low-cost BDS devices is developed. Thirdly, a land altimetry experiment is conducted to evaluate the precision of BDS-R and analyze the specific contribution of the BDS hybrid constellation. Finally, a water surface altimetry experiment validates the real-time monitoring capabilities for low-cost BDS-R. The results indicate that low-cost BDS-R can achieve real-time centimeter-level water level monitoring with a temporal resolution of 1 s in lakefront environments. The performance of BDS-R can be significantly improved by the BDS hybrid constellation, particularly IGSOs. It is concluded that low-cost BDS-R has great potential for promoting ground-based GNSS remote sensing applications.

BDS Orbit Maneuver Detection Based on Epoch-Updated Orbits Estimated by SRIF

The orbit maneuver detection is crucial in Global Navigation Satellite System (GNSS) precise orbit determination, which is necessary for adjusting data processing strategies. The frequency of orbit maneuvers for the BeiDou Navigation System is significantly higher than that of other navigation systems, especially for geosynchronous orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites. We propose a novel real-time and postprocessing method for detecting orbit maneuvers for BeiDou satellites based on the orbit differences between the epoch-updated orbit estimated using square root information (SRIF) and the predicted orbit according to the precise orbit estimated during non-maneuver period, as well as the orbital state difference during maneuver and non-maneuver periods. This method has significant advantages over using observation residuals and it is not affected by observation outliers, thus improving the accuracy and timeliness of orbit maneuver detection. We demonstrated that 32 orbit maneuver events of BeiDou satellites were successfully detected in 2022, of which 1 was for medium Earth orbit (MEO), 7 were for IGSO with an average detected maneuvering time of 7–8 min, and 24 were for GEO satellites with an average detected time of 4–5 min. Moreover, our method can be easily integrated into current real-time filter-based precise orbit determination (POD) processing without any extra task line, which simplifies the overall data processing. The data used in this method can be accessed easily, including GNSS observation data, broadcast ephemeris, and other open-source information files.

Determination of Earth rotation parameters by Beidou navigation satellite system

Earth Rotation Parameters (ERP) play an important role in astrometry, geophysics, and space geodesy. Beidou Navigation Satellite System (BDS) has three orbital types of satellites, i.e., Geostationary orbit (GEO), Inclined Geosynchronous Orbit (IGSO), and Medium Earth Orbit (MEO) satellites. The orbital diversity of BDS constellation provides the opportunity to study the influence of system-specific satellite constellation with different orbital characteristic on ERP. In this article, different strategies of ERP estimation are conducted in terms of different solar radiation pressure (SRP) models and different orbital types of satellites. For ERP derived from BDS MEO observations, compared with C04 products from International Earth Rotation and Reference System Service, the results show that the accuracy of ERP with extended Empirical CODE Orbit Model (ECOM) with 9 parameters (ECOM2) outperforms other SRP models (e.g., ECOM with 5 parameters (ECOM1) and ECOM1 with an empirical prior model (EMP + ECOM1)). In EMP + ECOM1 solutions, the application of a prior EMP model would intensify the influence of draconitic year in pole coordinates spectrum, while the visible signal of a period close to 7.4 days with ECOM1 solution can be mitigated. Moreover, the bias of B0 parameter in EMP + ECOM1 would cause a mean offset in Length Of Day (LOD). The accuracy of ERP is also affected by the BDS constellation characteristic. For pole coordinate, MEO-only solution with ECOM2 model shows the best accuracy, while the LOD achieves the best solution when 3 orbital types are employed. Additionally, using the combined constellation of MEO, IGSO, and GEO satellites can decorrelate some of the systematic signals, e.g., the spurious signals of LOD at 2.35-day and 3.5-day in MEO-only solutions are reduced. Likewise, the increase of orbital planes can reduce the amplitude of 3 cycles per draconitic year (close to 121 days) in the pole coordinate spectrum, which is the main error source in pole coordinates estimated from BDS MEO constellation. However, some signals around 10 ∼ 24, 33– and 50-day are found in the Y-pole coordinate and its rate.

Low-thrust station-keeping towards exploiting the inclined geosynchronous dynamics

Small satellite constellations in Inclined Geosynchronous Orbits (IGSO) provide an alternative to the increasingly crowded Geostationary Orbit belt. Satellite platforms equipped with full-electric propulsion systems are of great capacity for the advantages of significant propellant saving. This work focuses on the IGSO constellation design and its low-thrust station-keeping strategies. An analytical dynamical model based on a triple-averaged Hamiltonian formulation with respect to the ecliptic plane is first applied for a preliminary selection of the frozen orbits, which can serve as the baseline orbit for the whole constellation. Afterwards, the impact of the inclination on the manifolds in the eccentricity vector space is analyzed and the station-keeping slot for IGSO satellites is defined. Then, the low-thrust arc is designed to make the drifting arc closed, which maximize the long-term behaviors due to the resonances. Moreover, by feeding a three-impulse solution to the initial costate guess, an efficient solver to guarantee the convergence of the indirect optimization is proposed. Finally, the efficiency of the constellation design approach and the low-thrust SK strategy are validated via numerical simulations.

A two-step method to detect broadcast ephemeris unavailable periods caused by BeiDou satellite orbital maneuvers

The space constellation of the BeiDou navigation satellite system (BDS) is a hybrid constellation containing medium earth orbit (MEO) satellites, geostationary earth orbit (GEO) satellites, and inclined geosynchronous orbit (IGSO) satellites. Due to the geosynchronous characteristics of GEO and IGSO, GEO satellites and IGSO satellites often need to perform orbital maneuvers, which can affect the signal-in-space (SIS) availability performance of BeiDou satellites. A two-step detection method for BeiDou satellite orbital maneuvers has been proposed in this paper. The first step is to identify orbital maneuvers based on time series analysis of broadcast ephemeris, and the second step is to verify orbital maneuvers based on bidirectional orbit prediction. The two-step detection method was used to detect the orbital maneuvers of BeiDou satellites in 2019. Through the double guarantees of identification and verification, the detection accuracy of BeiDou satellite orbital maneuvers has been effectively improved. And the orbital maneuver detection results are continued to be used to assess the SIS availability of BeiDou satellites. The results show that the availability loss of GEO satellite orbital maneuvers is about 0.45%–1.07%, and the availability loss of IGSO satellite orbital maneuvers is about 0.12%–0.19%.

Signal-in-Space Range Error of the Global BeiDou Navigation Satellite System and Comparison with GPS, GLONASS, Galileo, and QZSS

Based on the standard of precise multisystem orbit and clock offset, the evaluation methods of satellite orbit, satellite clock offset and the signal-in-space range error (SISRE) from broadcast ephemeris were given in this study. The SISRE was used as a key performance index to compare the global BeiDou navigation satellite system (BDS-3) with the regional navigation satellite system (BDS-2), Global Positioning System (GPS), Galileo, Globalnaya Navigazionnaya Sputnikovaya Sistema (GLONASS), and Quasi-Zenith Satellite System (QZSS). Through the combination of domestic tracking stations and inter satellite links, the orbit determination accuracy of BDS-3 has been significantly improved. The radial (R), tangential (T) and normal (N) direction orbit accuracies of BDS-3 medium Earth orbit (MEO) satellites were 0.08, 0.37, and 0.34 m, respectively, which were considerably better than those of other satellite navigation systems. The R-, T-, and N-direction orbit accuracies of BDS-3 inclined geosynchronous orbit (IGSO) satellites were 0.17, 0.59, and 0.56 m, respectively, which were better than those of QZSS IGSO satellites. The R-, T-, and N-direction orbit accuracies of BDS-3 geostationary earth orbit (GEO) satellites were 0.27, 3.42, and 1.29 m, respectively, revealing considerable accuracy advantages in the R and N directions compared with QZSS GEO, although performance in the T direction was inferior. For the clock offset, the BDS-3 MEO satellites had an accuracy of 1.83 ns, which was noticeably superior, similar, slightly inferior, and superior to those of BDS-2, GPS, Galileo, and GLONASS. The clock offset accuracies of BDS-3 IGSO and GEO satellites were 3.13 and 3.17 ns, respectively, which were slightly better than those of QZSS IGSO and GEO satellites but did not have an obvious advantage over BDS-2. The primary reason could be that the launch time of the two types of BDS-3 satellites was close to the analysis period, and therefore, the clocks may be still in the aging process. Considering only orbit error, the average SISRE (orb) of BDS-3 MEO, Galileo, GPS, and GLONASS was 0.10, 0.28, 0.57, and 1.03 m, respectively. The SISRE (orb) of BDS-3 IGSO was 0.19 m, which was marginally greater than that of BDS-3 MEO but smaller than those of the other systems and considerably smaller than that (0.90 m) of QZSS IGSO. The SISRE (orb) of BDS-3 GEO was 0.60 m, which was comparable to that of GPS (0.59 m) and considerably smaller than that of QZSS GEO (0.91 m) satellites. Considering both orbit and clock offset error, the SISRE of BDS-3 MEO satellites was 0.52 m, which was slightly greater than that of Galileo (0.40 m), marginally smaller than that of GPS (0.59 m), and considerably smaller than that of GLONASS (2.33 m). The BDS-3 IGSO satellites had a SISRE of 0.90 m, which was greater than those of GPS and Galileo but comparable to that of QZSS IGSO (0.92 m). The average SISRE of BDS-3 GEO satellites was 1.15 m, which was greater than those of GPS and Galileo and slightly greater than that of QZSS0.95 m GEO satellites, but smaller than that of GLONASS.

Anomaly Detection Method of BDS Signal-in-Space Based on Autoregressive Distributed Lag Model

The signal-in-space (SIS) anomaly of BeiDou navigation satellite system (BDS) is an important factor affecting its high accuracy SIS quality assessment. Detecting and eliminating SIS anomaly is not only an important method to build SIS fault model of BDS, but also helps to guarantee the integrity of BDS navigation and positioning. Based on the problem that the traditional empirical threshold method cannot accurately identify the start and end times of anomalies in anomaly detection, which leads to anomaly detection leakage, a combined detection method based on autoregressive distributed lag model and empirical threshold is proposed in this paper. Before the calculation, the spurious anomalies of SIS are removed by data cleaning. The high-precision SIS ranging error (SISRE) is recovered by Space State Representation (SSR) correction number, and then projected to the user’s line of sight direction, and the anomaly detection threshold is determined by using the combined threshold of empirical threshold and autoregressive distributed lag (ARDL) model. The feasibility and effectiveness of the method were analyzed by using the data collected in 2021. The test results show that, compared with the traditional threshold method, the proposed method can more accurately detect the start and end points of SIS anomalies caused by clock anomalies, which is also the detection accuracy is improved. In addition, the anomaly detection method proposed in this paper is used to count the anomalies throughout the year, and the results show that the highest frequency of anomalies is found in geostationary orbit (GEO) and inclined geosynchronous orbit (IGSO), and these anomalies are mainly caused by satellite clocks.

Studies and Analysis of Combining BDS-3/GNSS Ultra-Rapid Orbit Products from Different IGS Analysis Centers

The third generation of China’s BeiDou Navigation Satellite System (BDS-3) began to provide global service at the end of 2018, and the completion of BDS was announced in July 2020, which includes GEO (Geostationary Earth Orbit), IGSO (Inclined Geosynchronous orbit), and MEO (Medium Earth Orbit) satellites. The resulting BDS orbits vary due to the inconsistent alignment strategies used by different analysis centers. Therefore, it is necessary to study the method of determining the BDS-3 orbit products combined from different analysis centers. In this research, the accuracy of the combined orbits for BDS-3 and other GNSS systems is evaluated and analyzed. To verify the reliability of the orbit combination method proposed in this paper, the GPS orbit is first selected for verification. Compared to the analysis centers, the mean Signal in Space User Ranging Error (SISURE) for GPS combined orbits is significantly reduced, and the mean SISURE of combined orbits for Block IIF, Block IIR, and Block III is 4.15 mm, 5.43 mm, and 5.63 mm, respectively. This demonstrates the effectiveness of the orbit combination method in this research. Besides, the accuracy of the combined orbits is improved by the ERP correction, and the mean RMS of the without Earth rotation Parameters (ERP) correction orbits and with ERP correction orbits is 4.78 mm and 4.53 mm, respectively. This demonstrates that orbit consistency corrections should be considered when performing orbit combinations. Compared to the GFZ orbits, the accuracy of the combined orbits has improved for GPS, GLONASS, GALILEO, and BDS is 8.2%, 9.9%, 9.9%, and 5.5%, respectively. It shows that the orbit combination method improves the orbital accuracy compared to the individual analysis center orbits. The mean RMS of the combined orbits for GPS, GLONASS, GALILEO, BDS MEO, and BDS IGSO is 1.7 cm, 2.61 cm, 2.52 cm, 2.59 cm, and 4.90 cm, respectively. The results demonstrate that the accuracy of the combined orbit for the BDS-3 MEO satellite is already similar to other systems; an orbit combination also available for the BDS-3 satellite.

Analysis of Precise Orbit Determination of BDS-3 MEO and IGSO Satellites Based on Several Dual-Frequency Measurement Combinations

The Chinese BeiDou-3 navigation satellite system (BDS-3) is capable of transmitting both old B1I, B3I signals and new B1C, B2a, B2b signals. Current BDS-3 precise orbits are generally calculated using a B1I/B3I combination considering overlap with the BeiDou-2 navigation satellite system (BDS-2). In this contribution, the observation quality of BDS-3 medium earth orbit (MEO) satellites and inclined geosynchronous orbit (IGSO) satellites are analyzed based on three aspects, i.e., carrier to noise ratio (C/N0), pseudo-range noise and pseudo-range multipath (MP). The C/N0 of the MEO satellite is 2~3 dB higher than that of the IGSO satellite at the same elevation angle. Meanwhile, the order of the Root Mean Square (RMS) values of both pseudo-range noise and MP is B1I < B1C < B3I < B2a ≈ B2b. Three kinds of combinations, i.e., B1CB2a, B1CB2b and B1IB3I, are selected for the BDS-3 precise orbit determination (POD) experiment. Orbits are assessed by the orbit-only signal-in-space range error (SISRE) computed between pairs of the three kinds of combinations in this contribution, CODE and GFZ final orbits. Orbit-only SISRE assessment shows that B1CB2a/CODE, B1CB2b/CODE, B1CB2a/GFZ and B1CB2b/GFZ are at the same level with CODE/GFZ, and the orbit-only SISRE is at the level of 5 cm for MEOs and 9 cm for IGSOs, respectively. Meanwhile, B1IB3I/CODE and B1IB3I/GFZ are about 1–2 cm worse. Inter-solution comparison between B1CB2a, B1CB2b and B1IB3I also indicate that B1CB2a and B1CB2b have good consistency, while B1IB3I shows poor performance. Satellite laser ranging (SLR) residuals indicate that the mean RMS is 3–4 cm for the four BDS-3 MEOs for CODE final orbit, GFZ final orbit, B1CB2a and B1CB2b combinations, while the mean RMS value for B1IB3I combination is a few millimeters worse, at approximately 4–5 cm.