Article Prediction of BeiDou Satellite Orbit Maneuvers to Improve the Reliability of Real-Time Navigation Products
Zhiwei Qin 1 , Le Wang 1,*, Guanwen Huang 1, Qin Zhang 1, Xingyuan Yan 2, Shichao Xie 1, Haonan She 1, Fan Yue 1 and Xiaolei Wang 3
1 School of Geology Engineering and Geomatics, Chang’an University, 126 Yanta Road, Xi’an 710054, China; [email protected] (Z.Q.); [email protected] (G.H.); [email protected] (Q.Z.); [email protected] (S.X.); [email protected] (H.S.); [email protected] (F.Y.) 2 School of Geospatial Engineering and Science, Sun Yat-sen University, No. 135, Xingang Xi Road, Guangzhou 510275, China; [email protected] 3 School of Earth Sciences and Engineering, Hohai University, 1 Xikang Road, Nanjing 210098, China; [email protected] * Correspondence: [email protected]; Tel.: +86-029-8233-9043
Abstract: The positioning, navigation, and timing (PNT) service of the Global Navigation Satellite System (GNSS) is developing in the direction of real time and high precision. However, there are some problems that restrict the development of real-time and high-precision PNT technology. Satellite orbit maneuvering is one of the factors that reduce the reliability of real-time navigation products, especially the high-frequency orbit maneuvering of geostationary earth orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites. The BeiDou Navigation Satellite System (BDS) constellation is designed to contain GEO, IGSO, and medium earth orbit (MEO). These orbit maneuvers bring certain difficulties for data processing, especially for BeiDou satellites, such as decreased real-time service performance, which results in real-time navigation products including unusable maneuvered Citation: Qin, Z.; Wang, L.; Huang, satellites. Additionally, the performance of real-time navigation products will decrease because G.; Zhang, Q.; Yan, X.; Xie, S.; She, H.; the orbit maneuvers could not be known in advance, which diminishes the real-time PNT service Yue, F.; Wang, X. Prediction of BeiDou performance of BDS for users. Common users cannot obtain maneuvering times and strategies owing Satellite Orbit Maneuvers to Improve the Reliability of Real-Time to confidentiality, which can lead to a decline in the BDS real-time service performance. Thus, we Navigation Products. Remote Sens. propose a method to predict orbit maneuvers. BDS data from the broadcast ephemeris were analyzed 2021, 13, 629. https://doi.org/ to verify the availability of the proposed method. In addition, the results of real-time positioning 10.3390/rs13040629 were analyzed by using ultra-rapid orbit products, demonstrating that their reliability is improved by removing maneuvered satellites in advance. This is vital to improve the reliability of real-time Academic Editor: Ali Khenchaf navigation products and BDS service performance. Received: 20 January 2021 Accepted: 5 February 2021 Keywords: BeiDou satellites; orbit maneuvers; prediction; real-time navigation products Published: 9 February 2021
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in 1. Introduction published maps and institutional affil- iations. The positioning, navigation, and timing (PNT) technology of the Global Navigation Satellite System (GNSS) grew out of nothing through development. The PNT service of GNSS is developing in the direction of real time and high precision. However, there are some problems restricting the development of real-time and high-precision PNT technology that decrease the reliability of real-time navigation products. Satellite orbit maneuvering is Copyright: © 2021 by the authors. one of the factors that reduce the reliability of real-time navigation products, especially the Licensee MDPI, Basel, Switzerland. high-frequency orbit maneuvers of geostationary earth orbit (GEO) and inclined geosyn- This article is an open access article chronous orbit (IGSO) satellites. It is necessary to improve the reliability of the PNT service distributed under the terms and of GNSS in order to satisfy the demands of high-reliability and high-precision navigation conditions of the Creative Commons Attribution (CC BY) license (https:// products for users. Thus, it is vital to enhance the reliability of real-time navigation prod- creativecommons.org/licenses/by/ ucts (such as orbit, clock, etc.). Since 27 December 2012, BeiDou-2 has provided continuous 4.0/). PNT service for users in the entire Asia-Pacific region, with a constellation of five GEO,
Remote Sens. 2021, 13, 629. https://doi.org/10.3390/rs13040629 https://www.mdpi.com/journal/remotesensing Remote Sens. 2021, 13, 629 2 of 18
five IGSO, and four MEO satellites. The second phase, BeiDou-3, will provide PNT service for global users at the end of 2020 [1–3]. The successful launch of the last satellite is a milestone of the construction completed for the BeiDou-3 constellation. At present, there are 30 BeiDou-3 satellites shining in the sky, consisting of three GEO, three IGSO, and 24 MEO satellites. BDS is the only GNSS that contains GEO and IGSO satellites in the nominal satellite constellation, which is the highlight of this GNSS. The MEO satellites offer complete global coverage, sharing numerous common features with GPS and GLONASS. Despite the higher orbit of the GEO and IGSO satellites as compared with the legacy MEO constellations, BDS achieves fully adequate signal strength and competitive code and car- rier tracking performance [4,5]. Together with GEO, IGSO and MEO satellites can provide fully global PNT service for worldwide users [6–10]. The GEO and IGSO satellites are considered relevant and valuable complements to GNSS services in the Asia-Pacific region, playing an important role [11–14]. They not only enhance satellite visibility and PNT avail- ability for users in China and surrounding areas, but also improve PNT precision for global users. The estimation accuracy of the earth orientation parameters can be improved based on the different orbital altitudes of the BDS satellites [6]. However, due to gravitational perturbations from the Earth, Sun, and Moon, as well as other various perturbations to the satellites, orbit maneuvers are required at regular intervals to maintain predesignated positions. GEO and IGSO satellites are maneuvered more frequently than MEO satellites because they are geosynchronous and geostationary. When orbit maneuvering, the orbital elements are changed, which causes the status variable to change after maneuvering. The broadcast ephemeris includes the health status of GNSS satellites; their health identifier will be marked as 0 for healthy, and otherwise the satellite is unhealthy. The satellite health identifier from broadcast ephemeris is also marked as unhealthy for nonmaneu- vered anomalies. Where information on the maneuvered satellites are unavailable owing to confidentiality, it brings several problems for GNSS data processing. With the development of Multi-GNSS, the position users can remove the unhealthy satellites directly, and the po- sition accuracy could not be decreased. However, the orbit maneuvers and unmaneuvered abnormalities of satellites cannot be distinguished by using the SV health from broadcast ephemeris. They both would be marked as unhealthy. In addition, the unmaneuvered could not affect the precise orbit determination. However, the precise orbit of maneuvered satellites cannot be determined and predicted because of the extra maneuvered thrust. As the orbit maneuvers and unmaneuvered abnormalities cannot be distinguished, the ultra-rapid precise orbit products published by IGS analysis centers include the unusable maneuvered satellites. These problems decrease the reliability of high-precision real-time navigation products of BDS, such as the ultra-rapid precise orbit products published by the analysis center of the international GNSS Monitoring and Assessment System (iGMAS) and the Multi-GNSS Experiment (MGEX). The time information of maneuvered satellites is unavailable in advance, so the real-time navigation products contain the unusable orbits of maneuvered satellites. Unusable real-time navigation products that include maneuvered satellites will decrease the real-time PNT performance of BDS. In order to improve the reliability of the real-time navigation products, the time information of orbit maneuvers must be obtained in advance. The predicted information for orbit maneuvering can be provided for position strategy changes and precise orbit determination. Thus, it is vital to provide orbit maneuver information in order to remove unusable maneuvered satellites, which is crucial for real-time navigation product solutions and real-time positions. Previous studies have investigated orbit maneuvers, yielding certain useful results. For orbit maneuvering of an inertial device to rendezvous and dock with another space- craft, Lemmens et al. [15] and Li et al. [16] presented a maneuver detection method for LEO satellites based on historical two-line element (TLE) data. LEO satellite maneuvers can be relatively easily handled with precise Global Positioning System (GPS) or TLE data. Compared to LEO satellites, the maneuvering of navigation satellites is significantly more difficult to investigate due to their high altitudes and the missing positions of ma- neuvered satellites. These methods cannot be used to analyze GNSS orbit maneuvers. Remote Sens. 2021, 13, 629 3 of 18
Huang et al. [17] and Cao et al. [18] used observations from the Chinese Area Positioning System (CAPS) to estimate the orbits of maneuvered GEO satellites. However, IGSO and MEO satellites cannot be continuously monitored via CAPS, so common users cannot obtain relevant observations. These methods bring certain difficulties when analyzing the characteristics of maneuvered orbits due to data limitations, especially for IGSO and MEO satellites. Huang et al. detected the start time of orbital maneuvers by using single point positioning (SPP) technology and residual values of a pseudo-range [19,20]. Qin et al. detected the orbit maneuvering time period by backward prediction of orbit and pseudo-range and determined the orbit for maneuvered satellites [21,22]. Numerous stud- ies have investigated orbit maneuver detection and determination; however, there is a lack of studies specifically focusing on predicting orbit maneuvers. This is why we propose a prediction method for orbit maneuvers. In this study, we analyze BeiDou satellite orbit maneuvers to characterize orbital elements, including the relationship with the orbital semimajor axis. First, we analyzed long-term variations in the orbital semimajor axis and proposed a prediction method for east–west maneuvering. The experiment results based on broadcast ephemeris were analyzed to demonstrate the validity of the prediction method for orbit maneuvers. In addition, we analyzed the results of real-time positioning based on ultra-rapid precise orbit products from the German Research Centre for Geosciences (GFZ), and we present our conclusions and implications.
2. Orbital Maneuver Analysis and Prediction Method Orbital elements have changed continuously due to perturbations from multiple planets and spacecraft, which cause long-term satellite drift. Satellites tend to gradually deviate from their designed orbits due to various perturbations, including perturbations from solar radiation pressure, Earth’s nonspherical gravity, and the gravity of the Sun, Moon, Venus, Mars, and other planets. External forces are solar, and other forces come from places such as nonspherical bodies. The total forces are calculated by:
= + + Ft ∑ Fgi Fs ∑ Fni (1)
where Ft is total forces affecting satellites, Fgi is gravity from other planets, Fs is the force
from solar radiation pressure, and Fni is the force of nonspherical body [23]. There are two modes of orbit maneuvers: east–west (in plane) maneuvers can retain the shape of the orbit plane, and north–south (out of plane) maneuvers can retain the inclination of the orbit plane. Orbit maneuvers in the north–south direction are much less frequent than maneuvers in the east–west direction. Thus, orbit maneuvers in the north–south direction are not discussed in this study. A schematic diagram of the gravitational force from the Sun and other planets affecting satellites is shown in Figure1. In Figure1, the orange circle and triangle denote the Sun, the yellow circle denotes Venus, and the dotted yellow line is its orbit around the Sun. The red circle denotes Mars, and its orbit is illustrated by a dotted red line. Earth is denoted by a green circle, and the dotted green line denotes its orbit around the Sun. The blue line denotes the designed orbit for satellites around the Earth, and the blue object is a satellite. Other planets (Mercury, Jupiter, Saturn, Uranus, Neptune, and Moon) are not shown in Figure1. The effect of instantaneous gravity from planets on satellites is denoted by arrows. The dotted orange arrow denotes the joint forces of gravity and solar pressure from the Sun. The dotted yellow, red, and green arrows denote the gravity from Venus, Mars, and Earth, respectively. Supposing the force affecting satellites is only from the Earth and the density of the Earth is uniform, satellites could run perennially around the Earth in their designed orbits; the blue arrow denotes the velocity of satellites. In fact, there are various forces acting on satellites from the gravity of multiple planets, solar pressure, and nonspherical gravity from the Earth. Thus, the actual instantaneous joint forces acting on satellites are different with the gravity from the Earth, which is denoted by the dotted black arrow, and the actual instantaneous velocity is expressed by the solid black arrow. Satellite orbits are Remote Sens. 2021, 13, 629 4 of 18
Remote Sens. 2021, 13, x FOR PEER REVIEW 4 of 18 continuously changed due to perturbations from multiple planets and other factors, as shown in Figure2.
Remote Sens. 2021, 13, x FOR PEER REVIEW 5 of 18
Figure 1.Figure Schematic 1. Schematic diagram diagramof gravitational of gravitational force from force planets from affecting planets affecting satellites. satellites.
In Figure 1, the orange circle and triangle denote the Sun, the yellow circle denotes Venus, and the dotted yellow line is its orbit around the Sun. The red circle denotes Mars, and its orbit is illustrated by a dotted red line. Earth is denoted by a green circle, and the dotted green line denotes its orbit around the Sun. The blue line denotes the designed orbit for satellites around the Earth, and the blue object is a satellite. Other planets (Mer- cury, Jupiter, Saturn, Uranus, Neptune, and Moon) are not shown in Figure 1. The effect of instantaneous gravity from planets on satellites is denoted by arrows. The dotted or- ange arrow denotes the joint forces of gravity and solar pressure from the Sun. The dotted yellow, red, and green arrows denote the gravity from Venus, Mars, and Earth, respec- tively. Supposing the force affecting satellites is only from the Earth and the density of the Earth is uniform, satellites could run perennially around the Earth in their designed orbits; the blue arrow denotes the velocity of satellites. In fact, there are various forces acting on satellites from the gravity of multiple planets, solar pressure, and nonspherical gravity from the Earth. Thus, the actual instantaneous joint forces acting on satellites are different with the gravity from the Earth, which is denoted by the dotted black arrow, and the ac- tual instantaneous velocity is expressed by the solid black arrow. Satellite orbits are con- tinuously changed due to perturbations from multiple planets and other factors, as shown in Figure 2.
Figure 2. FigureSchematic 2. Schematic diagram ofdiagram satellite of orbit satellite drift orbit due todrift various due to perturbations various perturbations from multiple from planetsmultiple and spacecraft. planets and spacecraft. In Figure2, the green circle denotes the Earth, which is the focal point of satellite In Figure orbits.2, the green The initial circle orbit denotes of satellites the Earth, is shown which theis the elliptical focal point blue curve,of satellite line a is the orbital orbits. The initialsemimajor orbit of axis,satellites line isb isshown the orbital the elliptical semiminor blue axis,curve, and line line c isis the the orbital half focal length of semimajor axis,the line satellite is the orbit. orbital The semiminor light blue partaxis, denotesand line the is normal the half range focal of length satellite of orbits. When the satellite orbit.satellites The light deviate blue from part the denote normals the orbit normal range, range they of can satellite be adjusted orbits. by When orbit maneuvering, satellites deviateusing from their the propulsionnormal orbit system. range, they The greencan be arrow adjusted denotes by orbit the maneuvering, instantaneous gravity force using their propulsionacting on system. satellites The from green the arrow Earth. denotes The black the arrow instantaneous denotes the gravity instantaneous force joint forces acting on satellitesfrom from multiple the Earth. planets The acting black onarrow satellites, denotes which the instantaneous leads to satellites joint tendingforces to gradually from multiple deviateplanets fromacting their on normalsatellites, orbital which range. leads The to satellites shape of thetending satellite to gradually orbits also changes with deviate from their normal orbital range. The shape of the satellite orbits also changes with time. This can be shown by the pure curve, which denotes the satellite orbit that needs to be adjusted though orbit maneuvering; line ′ is the orbital semimajor axis, line ′ is the orbital semiminor axis, and line ′ is the half focal length of the satellite orbit. The varia- tion in the orbital semimajor axis ∆ is the difference value of the norm of and ′, which is indicated by the dotted red line in Figure 2. Satellites periodically traverse the sky around the center of the Earth, which allows the perturbations from multiple planets to yield periodicity. The comprehensive perturbations from planets on satellites can be counteracted on short-term time scales [24]. The short-term periodical variation of the semimajor axis is not discussed for orbit maneuver prediction in this study. The compre- hensive perturbations also cause long-term satellite drift, which leads to variation of the orbital semimajor axis. The orbital semimajor axis can be expressed as follows: