The failure to detect anomalies and maneuvering of the orbits of navigation satellite sensors will deteriorate the performance of positioning and orbit determination. Motivated by the influence of the frequent maneuvering of BDS GEO and IGSO satellites, this paper analyzes the limitations of existing methods, where BDS orbit maneuvering and anomalies can be detected, and develops a method to solve this problem based on the RMS model of orbit mutual differences proposed in this paper. The performance of this method was assessed by comparison with the health flag of broadcast ephemeris, precise orbit products of GFZ, the O-C values of a GNSS station and a conventional method. The results show that the performance of the method developed in this paper is better than that of the conventional method when the periodicity and trend items are obvious. Meanwhile, three additional verification results show that the method developed in this paper can find error information in the merged broadcast ephemeris provided by iGMAS. Furthermore, from the testing results, it can be seen that the detection of anomaly and maneuvering items do not affect each other based on the robust thresholds constructed in this paper. In addition, the precise orbit of the maneuvering satellites can be determined under the circumstances that the maneuver information detected in this paper is used, and the root mean square (RMS) of orbit overlap comparison for GEO-03/IGSO-03 in Radial, Along, Cross, 1D-RMS are 0.7614/0.4460 m, 1.8901/0.3687 m, 0.3392/0.2069 m, 2.0657/0.6145 m, respectively.

By August 2016, 5 new-generation BeiDou satellites (BeiDou-3) have successfully been launched. The observations of a very limited number of 9 International GNSS (Global Navigation Satellite System) Monitoring and Assessment Service (iGMAS) stations and 52 Multi-GNSS Experiment (MGEX) stations from 16 July to 14 August 2016 are processed to determine the orbits of BeiDou-3 and BeiDou-2 satellites, respectively. The internal consistency and satellite laser ranging (SLR) validations are conducted for the orbit validation. BeiDou-3 MEO (Medium Earth Orbit) (C33 and C34) have larger root mean square (RMS) values than those BeiDou-3 IGSO (C31 and C32), whereas BeiDou-2 MEO satellites have smaller RMS values than the BeiDou-2 IGSO satellites. Furthermore, BeiDou-3 IGSO and BeiDou-2 satellites have RMS values at identical levels, whereas BeiDou-3 MEO satellites have larger RMS values than the BeiDou-2 MEO satellites. The RMS residuals are approximately 10 cm in the radial component and approximately 25 cm in the along component for BeiDou-3 IGSO satellites. For BeiDou-3 MEO satellites, the RMS residuals are approximately 40 cm in the radial component and approximately 60 cm in the along component. The SLR validation reports that the orbit radial component can reach an accuracy on the level of 1 decimeter and 4 decimeters for BeiDou-3 IGSO and MEO, respectively.

In this study, the characteristics of the total electron content (TEC) fluctuations and their regional differences over China were analyzed by utilizing the rate of the TEC index (ROTI) based on GPS data from 21 reference stations across China during a solar cycle. The results show that there were significant regional differences at different latitudes. Strong ionospheric TEC fluctuations were usually observed at lower latitudes in southern China, where the occurrence of TEC fluctuations demonstrated typical nighttime- and season-dependent (equinox months) features. This phenomenon was consistent with the ionospheric scintillation characteristics of this region. Additionally, compared to low-latitude China, the intensity of TEC fluctuations over mid-latitude China was significantly weaker, and the occurrence of TEC fluctuations was not a nighttime-dependent phenomenon. Moreover, the intensity of TEC fluctuations was much stronger during high solar activity than during low solar activity. Furthermore, the summer-dependent characteristics of TEC fluctuations gradually emerged over lower mid-latitude areas as equinox characteristics weakened. Similar to the equinox characteristics, the summer-dependent characteristics gradually weakened or even disappeared with the increasing latitude. Relevant discussions of this phenomenon are still relatively rare, and it requires further study and analysis.

An analytical solar radiation pressure (SRP) model, IGGBSPM (an abbreviation for Institute of Geodesy and Geophysics BeiDou Solar Pressure Model), has been developed for three BeiDou satellite types, namely, geostationary orbit (GEO), inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO), based on a ray-tracing method. The performance of IGGBSPM was assessed based on numerical integration, SLR residuals and analyses of empirical SRP parameters (except overlap computations). The numerical results show that the integrated orbit resulting from IGGBSPM differs from the precise ephemerides by approximately 5 m and 2 m for GEO and non-GEO satellites, respectively. Moreover, when IGGBSPM is used as an a priori model to enhance the ECOM (5-parameter) model with stochastic pulses, named ECOM + APR, for precise orbit determination, the SLR RMS residual improves by approximately 20–25 percent over the ECOM-only solution during the yaw-steering period and by approximately 40 percent during the yaw-fixed period. For the BeiDou GEO01 satellite, improvements of 18 and 32 percent can be achieved during the out-of-eclipse season and during the eclipse season, respectively. An investigation of the estimated ECOM D0 parameters indicated that the β-angle dependence that is evident in the ECOM-only solution is no longer present in the ECOM + APR solution.

Offering real-time precise point positioning (PPP) services for global and large areas based on global navigation satellite systems (GNSS) has drawn more and more attention from institutions and companies. A precise and reliable satellite orbit is a core premise for multi-GNSS real-time services, especially for the GPS and GLONASS, which are undergoing modernization, whereas the Galileo, BDS and QZSS have just fulfilled the construction stage. In this contribution, a real-time precise orbit determination (POD) strategy for the five operational constellations based on the hourly updated ultrarapid orbit prediction method is presented. After combination of 72 h arc through three adjacent 24 h arc normal equations, the predicted orbits are finally generated (hourly updated). The POD results indicate that the mean one-dimensional (1-D) root mean square (RMS) values compared with the Deutsches GeoForschungsZentrum (GFZ) final multi-GNSS orbits are approximately 3.7 cm, 10.2 cm, 5.8 cm, 5.7 cm, 4.1 cm and 25.1 cm for GPS, BDS IGSOs, BDS MEOs, GLONASS, Galileo and QZSS NONE GEOs, respectively. The mean 1-D RMS values of the hourly updated ultrarapid orbit boundary overlapping comparison are approximately 1.6 cm, 6.9 cm, 3.2 cm, 2.7 cm, 1.8 cm and 22.2 cm for GPS, BDS IGSOs, BDS MEOs, GLONASS, Galileo and QZSS NONE GEOs, respectively. The satellite laser ranging (SLR) validation illuminates that the mean RMS values are approximately 4.53 cm and 4.73 cm for the four MEOs of BDS-3 and four BDS-2 satellites, respectively.

The predicted parts of ultra-rapid orbits are important for (near) real-time Global Navigation Satellite System (GNSS) precise applications; and there is little research on GPS/GLONASS/BDS/Galileo/QZSS five-system ultra-rapid precise orbit determination; based on the one-step method and double-difference observation model. However; the successful development of a software platform for solving five-system ultra-rapid orbits is the basis of determining and analyzing these orbits. Besides this; the different observation models and processing strategies facilitate to validate the reliability of the various ultra-rapid orbits. In this contribution; this paper derives the double-difference observation model of five-system ultra-rapid precise orbit determination; based on a one-step method; and embeds this method and model into Bernese v5.2; thereby forming a new prototype software platform. For validation purposes; 31 days of real tracking data; collected from 130 globally-distributed International GNSS Service (IGS) multi-GNSS Experiment (MGEX) stations; are used to determine a five-system ultra-rapid precise orbit. The performance of the software platform is evaluated by analysis of the orbit discontinuities at day boundaries and by comparing the consistency with the MGEX orbits from the Deutsches GeoForschungsZentrum (GFZ); between the results of this new prototype software platform and the ultra-rapid orbit provided by the International GNSS Monitoring and Assessment System (iGMAS) analysis center (AC) at the Institute of Geodesy and Geophysics (IGG). The test results show that the average standard deviations of orbit discontinuities in the three-dimension direction are 0.022; 0.031; 0.139; 0.064; 0.028; and 0.465 m for GPS; GLONASS; BDS Inclined Geosynchronous Orbit (IGSO); BDS Mid-Earth Orbit (MEO); Galileo; and QZSS satellites; respectively. In addition; the preliminary results of the new prototype software platform show that the consistency of this platform has been significantly improved compared to the software package of the IGGAC.