Basic performance and future developments of BeiDou global navigation satellite system

Abstract

The core performance elements of global navigation satellite system include availability, continuity, integrity and accuracy, all of which are particularly important for the developing BeiDou global navigation satellite system (BDS-3). This paper describes the basic performance of BDS-3 and suggests some methods to improve the positioning, navigation and timing (PNT) service. The precision of the BDS-3 post-processing orbit can reach centimeter level, the average satellite clock offset uncertainty of 18 medium circular orbit satellites is 1.55 ns and the average signal-in-space ranging error is approximately 0.474 m. The future possible improvements for the BeiDou navigation system are also discussed. It is suggested to increase the orbital inclination of the inclined geostationary orbit (IGSO) satellites to improve the PNT service in the Arctic region. The IGSO satellite can perform part of the geostationary orbit (GEO) satellite’s functions to solve the southern occlusion problem of the GEO satellite service in the northern hemisphere (namely the “south wall effect”). The space-borne inertial navigation system could be used to realize continuous orbit determination during satellite maneuver. In addition, high-accuracy space-borne hydrogen clock or cesium clock can be used to maintain the time system in the autonomous navigation mode, and stability of spatial datum. Furthermore, the ionospheric delay correction model of BDS-3 for all signals should be unified to avoid user confusion and improve positioning accuracy. Finally, to overcome the vulnerability of satellite navigation system, the comprehensive and resilient PNT infrastructures are proposed for the future seamless PNT services.

Introduction

The BeiDou navigation satellite system (BDS) provides more services compared with other global navigation satellite systems (GNSSs). The BeiDou global navigation satellite system, (BDS-3), is the third step of China’s satellite navigation system construction [1, 2]. In addition to the positioning, navigation and timing (PNT) service provided by all GNSSs, BDS-3 also provides regional message communication (1000 Chinese characters per time) and global short message communication (40 Chinese characters per time), global search and rescue service (SAR), regional precise point positioning (PPP) service, embedded satellite-based augmentation service (BDSBAS), and space environment monitoring function [3,4,5]. By the end of 2019, 28 BDS-3 satellites had been successfully launched, including 24 medium circular orbit (MEO) satellites, 3 inclined geostationary orbit (IGSO) satellites and 1 geostationary orbit (GEO) satellite. However, only 18 MEO satellites can provide services thus far, and other satellites are in test phase.

The BDS-3 baseline system with 18 MEO satellites has begun providing initial services to global users on December 27, 2018, and at least five satellites are visible for global users. Many scholars have described the BDS-3 status and the main service function indexes [2,3,4,5,6], including the constellation design, service type, navigation signal system, space–time datum and orbit determination method of BDS-3 system, etc. The space–time datum, signal-in-space quality, accuracy of satellite broadcast ephemeris, the accuracy of the post-processing precise orbit, time synchronization accuracy, satellite clock offset accuracy and the basic PNT service performance of the BDS-3 baseline system have be calculated and analyzed [3,4,5, 7,8,9,10]. Because BDS-3 satellites are equipped with inter-satellite links (ISL), more research achievements have been made in the fields of ISL supported determination of satellite orbit and clock offset [3,4,5, 11,12,13]. According to the preliminary calculation results, the ranging accuracy of the ISLs of BDS-3 is about 4 cm; if the satellite orbits are determined using only regional station observations, the three-dimensional orbit accuracy of the overlapping arc is about 60 cm; with ISL measurements, the orbit accuracy is about 30 cm, the 24-h orbit prediction accuracy is also raised from 140 to 51 cm [3,4,5], and the radial accuracy can reach 10 cm, as evaluated through laser observations.

The performance of BDS-3 presented by different scholars, using different data sources, are highly similar to each other. From the aspect of signal, the ratio bias of signal component effective power is better than 0.25 dB, the S curve bias is less than 0.3, the signal-in-space accuracy calculated with post-processing ephemeris and broadcast ephemeris is approximately 0.5 m (root mean square, RMS), and the signal-in-space continuity and availability are approximately 99.99% and 99.78%, respectively [7]. The timing accuracy is better than 19.1 ns (95%), which satisfies the system service performance specifications. Compared with BDS-2, BDS-3 exhibits significant improvement in system coverage, spatial signal accuracy, availability and continuity [2,3,4,5]. At the user’s end, service improvements such as the BDS-3 signal design and optimal receiving mode [14], precision orbit determination [3,4,5, 9, 15], and precise point positioning technology [16] have been reported in relevant literature and will not be repeated herein.

The main performance evaluation methods and formulae are given in the second section. The current performance of BDS-3 by the end of August 2019 is presented and analyzed in the third section. Under normal circumstances, the main performance indicators provided by BDS-3 can satisfy or surpass the design indicators, which are comparable to the service performance of other GNSSs. The possible improvement strategies for BDS-3 are presented in the fourth section from the aspects of system construction, satellite constellation, satellite payload, etc. Conclusions are presented in the final section of the paper.

Performance evaluation methods

The spatial signal accuracies of BDS-3 include broadcast orbit accuracy, broadcast clock offset accuracy, signal-in-space ranging error (SISRE), broadcast ionospheric delay model accuracy, etc. Specific evaluation methods are as follows:

Broadcast orbit accuracy: the post precise satellite positions are taken as references, and the accuracy of the positions calculated by the broadcast ephemeris is the evaluation object. The root mean square error (RMSE) of the broadcast orbit is evaluated using the differences (errors of the broadcast orbit) of the broadcast and precise satellite orbits in the same coordinate system, i.e. China Geodetic Coordinate System 2000 (CGCS2000) [17] and the same time scale, i.e. BDS time (BDT) [18]. The orbit error can be obtained using the following formula [19,20,21,22]

$$\Delta \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {R} = \left( {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {R}_{brd} + A \cdot PCO_{brd} } \right) - \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {R}_{pre}$$

(1)

where $\Delta \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {R}$ is the broadcast orbit error vector; $\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {R}_{brd}$ is the satellite position vector calculated by the broadcast ephemeris with the unit m; $A$ denotes the transformation matrix from the satellite-fixed coordinate system to the earth fixed system; $PCO_{brd}$ is the correction from the satellite antenna phase center to the satellite mass center with the unit m; $\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {R}_{pre}$ is the precise satellite position vector calculated by precise ephemeris with the unit m. And then the broadcast orbit accuracy in a specific period can be obtained from the statistics of differences of the satellite positions.

Broadcast clock offset accuracy: the post precise satellite clock offsets are taken as the references, and the accuracy of the satellite clock offsets calculated by the broadcast ephemeris is the evaluation object. The RMS of the broadcast clock offset is evaluated using the differences between the broadcast and precise clock offsets. The broadcast clock offset is calculated using the following formula

$$ \Delta \tilde{c}_{k}^{t} = \Delta c_{k}^{t} - \frac{{\mathop \sum \nolimits_{k = 1}^{N} \Delta c_{k}^{t} }}{N} $$

(2)

where $\Delta \tilde{c}_{k}^{t}$ is the broadcast clock offset of the satellite $k$ after subtracting the time datum difference at the epoch $t$; $\Delta c_{k}^{t}$ is the difference between the broadcast and precise clock offsets after time group delay (TGD) and antenna phase center corrections; $N$ denotes the number of satellites; $k$ denotes the satellite number; $t$ denotes the epoch time. $\Delta c_{k}^{t}$ is calculated using the following formula

$$c_{k}^{t} = B_{k}^{t} - \frac{{TGD_{{{\text{m}}1}} f_{{{\text{m}}1}}^{2} - TGD_{{{\text{m}}2}} f_{{{\text{m}}2}}^{2} }}{{f_{{{\text{m}}1}}^{2} - f_{{{\text{m}}2}}^{2} }} - \frac{{e_{k}^{pre} - e_{k}^{brd} }}{c} - \bar{C}_{k}^{t}$$

(3)

where $B_{k}^{t}$ is the clock offset of the $k$-th satellite calculated by the broadcast clock offset parameters at the epoch t; $\bar{C}_{k}^{t}$ is the precise clock offset of the $k$-th satellite at the epoch t with the unit $s$; $\frac{{TGD_{{{\text{m}}1}} f_{{{\text{m}}1}}^{2} - TGD_{{{\text{m}}2}} f_{{{\text{m}}2}}^{2} }}{{f_{{{\text{m}}1}}^{2} - f_{{{\text{m}}2}}^{2} }}$ is the TGD correction of the satellite clock offset obtained using the broadcast TGD parameters with the unit $s$; $TGD_{m1}$ is the time delay difference of on-board equipment corresponding to the $f_{m1}$ frequency signal with unit $s$; $f_{m1}$ and $f_{m2}$ denote the signal carrier frequencies adopted in the double differential ionosphere-free combination; $\frac{{e_{k}^{pre} - e_{k}^{brd} }}{c}$ is the difference between the antenna phase center corrections of the broadcast and precise clock offsets with the unit $s$; $e_{k}^{pre}$ and $e_{k}^{brd}$ denote the Z components of the satellite antenna phase center corrections adopted by the products of the precise clock offset and the broadcast clock offset respectively with the unit m. The accuracy of broadcast clock offset in the evaluation period can then be obtained through the statistical calculation of clock offset errors.

Signal-in-space ranging error (SISRE): the post orbit parameters are taken as the reference, and the evaluation object is set to be the projection accuracy of the combined error of the satellite orbit error and the clock offset on the sightline between the satellite and the user. It is usually calculated using the differences of the broadcast and post processed precise orbit parameters and the differences of the broadcast and post processed precise clock offsets [23,24,25]

$$SISRE = \sqrt {\left( {\alpha \delta_{R} - c\delta_{T} } \right)^{2} + \beta \left( {\delta_{A}^{2} + \delta_{C}^{2} } \right)}$$

(4)

where $\delta_{R}$, $\delta_{A}$ and $\delta_{C}$ denote the errors of broadcast orbit in the radial, tangential and normal components respectively with the unit m; $\delta_{T}$ denotes the clock offset error with the unit s; $\alpha$ denotes the contribution factor of the radial component error, and $\beta$ denotes that of the tangential and normal errors. The contribution factors of different GNSS satellites are listed in Table 1.

Table 1 Average contribution factors of satellites on SISRE (elevation mask angle: 5°)

Broadcast ionospheric model accuracy: the post-processing high-precision grid ionospheric delay model is taken as the reference model, and the accuracy of the broadcast ionospheric zenith delay value is the evaluation object. In this paper, the accuracy of BeiDou global broadcast ionospheric delay correction model (BDSGIM) is evaluated [26]. Four indexes are used for evaluation, including average bias, standard deviation, RMSE and correction percentage. The formulae are as follows [26, 27]:

$$bias = \frac{{\mathop \sum \nolimits_{i = 0}^{N} \left( {vTEC_{model}^{i} - vTEC_{ref}^{i} } \right)}}{N}$$

(5)

$$std = \sqrt {\frac{{\mathop \sum \nolimits_{i = 0}^{N} \left( {vTEC_{model}^{i} - vTEC_{ref}^{i} - bias} \right)^{2} }}{N - 1}}$$

(6)

$$RMS = \sqrt {\frac{{\mathop \sum \nolimits_{i = 0}^{N} \left( {vTEC_{model}^{i} - vTEC_{ref}^{i} } \right)^{2} }}{N}}$$

(7)

$$\left\{ {\begin{array}{*{20}l} {per_{i} = \left( {1 - \frac{{\left| {vTEC_{model}^{i} - vTEC_{ref}^{i} } \right|}}{{vTEC_{ref}^{i} }}} \right) \cdot 100\% } \hfill \\ {per = \left( {\frac{{\mathop \sum \nolimits_{i = 1}^{{N_{{\left( {per_{i} \ge 0} \right)}} }} per_{i} }}{{N_{{\left( {per_{i} \ge 0} \right)}} }}} \right) \cdot 100\% ,\left( {1 - \frac{{N_{{\left( {per_{i} \ge 0} \right)}} }}{N}} \right) \cdot 100\% } \hfill \\ \end{array} } \right.$$

(8)

where $i$ denotes the ith grid point; $N$ is the total number of grid points; $vTEC_{model}$ is the correction value calculated by the model at the grid point; $vTEC_{ref}$ is the reference value.

Latest performance of BDS-3

Broadcast satellite orbit accuracy

With the continuous improvement of the BDS-3 basic constellation, the stability and connectivity of the ISL continue to increase, and the performance of the basic constellation improves gradually. The basic service performance of BDS-3 first depends on the performance of the satellite orbit and satellite clock. The performance of the BDS-3 post-processing precision orbit was estimated from the observations of international GNSS monitoring assessment service (iGMAS) stations deployed all over the world until July 2019. It indicated that the accuracies of the radial, tangential and normal components of the satellite orbit were 1.5, 5.7 and 4.1 cm respectively (seen in Fig. 1a), and the accuracies were 8.0, 34 and 37 cm, respectively, when the system started offering global service in December 2018.

Fig. 1

a Precise orbit accuracy of the BDS-3 basic constellation. b Broadcast orbit accuracy of the BDS-3 basic constellation

As shown, the post-processing orbit accuracy of the BDS-3 satellites are increased by approximately five, six and nine times, respectively. The improvement in the normal accuracy is the most significant, as the ISL observation is more abundant and stable. The accuracy of the three components (radial, tangential and normal) of the broadcast orbit of the BDS-3 by referring to the post accurate orbits are 0.059, 0.323 and 0.343 m respectively (seen in Fig. 1b). Compared with the launch in December 2018, the broadcast orbit has also been significantly improved. Furthermore, with the significant improvement of post-processing orbit accuracy, the estimation reliability of broadcast ephemeris based on the post-processing precise orbit has also been increased accordingly.

Broadcast satellite clock offset accuracy

The uncertainty of the broadcast satellite clock offset is a main factor affecting the user PNT service. The accuracy of the satellite broadcast clock offset is evaluated using the differences of the broadcast satellite clock offsets and post-processing precise ones. The RMS are displayed in Fig. 2.

Fig. 2

Broadcast clock offset of the BDS-3 satellites

As shown in the figure, the average uncertainty of satellite clock offset of 18 MEO satellites is 1.55 ns, but the satellite clock offsets of M02, M06 and M15 satellites are 3.068, 2.555 and 2.832 ns, respectively, which are larger and unstable.

Estimation of SISRE

SISRE of the BDS-3 are estimated based on the post-processing ephemeris and broadcast ephemeris, which are shown in Fig. 3.

Fig. 3

SISRE of the BDS-3 satellites

Figure 3 shows that the average SISRE is approximately 0.474 m, and all the SISREs are in agreement with the satellite clock offsets. The SISRE of the M02, M06 and M15 satellites are significantly greater than those of other satellites; when BDS-3 started offering global service in December 2018, the SISRE was approximately 0.7 m.

Accuracy of ionospheric delay correction models

BDS-3 can broadcast four civil signals, including B1C, B1I, B2a/B2b and B3I. Among them, the B1I and B3I signals are broadcast by both BDS-2 and BDS-3. Therefore, the civil signals of BDS-3 adopt two types of ionospheric models, namely BDSK8 and BDGIM [28]. The former is used for B1I and B3I which is consistent with those of the BDS-2, while the latter is used for the new signals of BDS-3, i.e. B1a, B1C, B2a and B2b. To analyze the ionospheric correction accuracy of these two models, we used the global ionospheric maps model constructed by the Center for Orbit Determination in Europe (CODG model) as a reference, and compared the nine-month (from January to September) average accuracies of Nequick, GPSK8, BDGIM and BDSK8 models. The results are shown in Fig. 4, where “high latitude” ranges from 55° to 85°, “mid latitude” from 30° to 55°, and “low latitudes” from 0° to 30°.

Fig. 4

Evaluation of the ionospheric model accuracy in half a year based on CODG model

Figure 4 shows that the correction accuracy of BDGIM model is obviously better than that of all other models in the middle and low latitude zones; in the northern middle and low latitudes, the correction ratio of BDGIM is approximately 75%; in the southern middle latitude zones, the correction ratio is approximately 65%; in the high latitude zones, BDGIM is slightly worse than the Nequick model; the BDGIM performs better than the GPSK8 and BDSK8 model at any latitude.

The accuracy evaluation details of BDGIM are shown in Table 2, and the correction percentage and data rejection percentage are listed in the last row of low, middle and high latitude columns.

Table 2 Accuracy evaluation of BeiDou global ionospheric model

Table 2 shows that the BDGIM maintains a good correction accuracy, regardless of the ionospheric effect.

According to the performance evaluation of the BDS-3 orbit accuracy, satellite clock offset accuracy and SISRE, the basic functions and performance of BDS-3 satisfy the design requirements.

Future possible development

Although BDS-3 has developed many new designs and functions, such as satellite orbit design, service function innovation and payload improvements [3,4,5], there is still room for improvements in the future.

• 1. The PNT service performance in high latitude area is poor. The BDS-3 constellation is composed of 24 MEO, 3 GEO and 3 IGSO satellites. The inclination of MEO and IGSO satellites is 55°, which indicates that the service performance of BDS-3 will be significantly degraded in high latitude area. First, all visible BDS-3 satellites relative to the users in the polar regions are in low-elevation angle. Calculations show that the maximum altitude angle of the satellites is generally less than 50° when the users are in regions of latitude greater than 85° [29]. The low-elevation satellites will decrease the user ranging accuracy; next, the BDSK8 model is unavailable in the Arctic and Antarctic regions, and the accuracy of BDGIM ionospheric model is relatively low, which will degrade high-accuracy positioning. Furthermore, owing to the global climate change, sea ice in the polar regions of the Earth is melting rapidly. Especially after ice and snow in the Arctic has melted in the summer, the demand for PNT service in the Arctic region becomes more urgent considering the significant route value and rich resource reserve [29, 30].

To ensure the efficiency of all kinds of scientific research and the safety of transportation in the Arctic region, the satellite navigation service must be improved. For BDS-3, a reasonable and simplified method is to increase the inclination of IGSO satellites.

Using the actual constellation of BDS-3, we analyzed the user’s elevation angle variation according to changes in the orbit inclination of the IGSO satellites in high-latitude areas by adjusting the inclination of the IGSO constellation to 55°, 65° and 75°. Seven-day simulation data with a sampling interval of 10 min were used. To calculate the average dilution of precision (DOP), grids with the longitude interval of 2° and latitude interval of 1° were used, and the minimum altitude angle was set as 5°. The latitude of key test area ranged from 75° to 90°(N) whereas the longitude ranged from − 180° to 180°.

Table 3 shows that with a 10° increase of the inclination angle of the IGSO satellites in the high latitude region, the average altitude angle of the visible satellite increases by approximately 5°, and the maximum altitude angle increases by approximately 12°. With the increase in satellite elevation angle, the number of visible satellites will increase, and the effect of ionosphere will be weakened.

Table 3 Statistics of mean elevation angle with changes of the IGSO inclinations

• 2. An obvious “south wall effect” exists in BDSPPP and regional short message communication (RSMC) services based on GEO satellites. For users in the northern hemisphere, GEO satellites are always located south of the user, and for those at higher latitudes, the BDSPPP and short message communication service will likely be interrupted once obstacles appear in the south. The “south wall effect” is particularly serious for users of PPP, as 20–30 min are required to obtain the converged results after the PPP service recovers from the GEO satellites. Furthermore, the repeated convergence will result in a significant reduction in the application efficiency of BDSPPP in cities with high-rise buildings. Similarly, users in southern hemisphere will experience the “north wall effect”.

  To overcome the “south wall effect” of BDSPPP and RSMC provided by BeiDou GEO satellites, the most effective method is to use the BeiDou IGSO satellites to transmit rapid precise orbit and clock offset parameters, as well as the regional short message. IGSO and GEO satellites working together, can eliminate the “south wall effect” and effectively improve the featured service efficiency of BDS. Furthermore, if the inclination angle of IGSO satellites are increased and an appropriate ground reference station is built in the Arctic region, then the BDSPPP and RSMC services can be provided to users in the Arctic region.

• 3. Some satellites’ services are interrupted during the orbit maneuver of the satellites. Satellite orbit maneuver is inevitable. Almost every month, BeiDou GEO satellite must adjust its East-West orbit, and conduct a North-South orbit maneuver every half a year. IGSO and MEO satellites generally conduct an orbit maneuver every half a year. During the period of orbit adjustment, the satellite ephemeris cannot reflect the orbit maneuver quickly and is typically expressed as “unavailable” or “unhealthy”. Furthermore, the PNT service of the corresponding satellite will be suspended. The duration is approximately 7 h. For IGSO and MEO satellites, owing to their large number, low frequencies of adjustment and the sole functions, the PNT service will not be significantly affected in general. However, the BDSPPP and RSMC service provided by the three GEO satellites of the BDS-3 will significantly affect many authorized users owing to the unavailability caused by GEO satellites' orbit maneuver.

  As we know, the maneuvering satellites cannot provide normal service because the satellite ephemeris cannot accurately reflect the actual positions of the satellites. If inertial navigation system (INS) equipment was installed on the satellite, the orbit change in the satellite maneuver could be measured during the satellite maneuver, and thus the orbit parameters of the moving satellites could be obtained. However, one of the main problems is that during the non-maneuverable period of the satellite, the on-board INS must be calibrated using a precise satellite orbit or an error correction model must be established. Once the satellite enters the maneuvering state, INS observations and corresponding corrections can be used to provide the satellite position during the maneuvering period, which will improve the availability of the maneuvering satellites.

• 4. Datum drift and rotation of BeiDou satellite constellation will occur in autonomous navigation. If only ISL measurements are used in the autonomous navigation, then only the relative positions of constellation can be solidified. The overall drifts of the spatial datum and time datum cannot be determined or solidified, and neither can the overall rotation of the entire constellation. These types of systematic errors may not significantly affect high-accuracy differential positioning users; however, they will significantly affect the real-time navigation and timing service users.

  Regarding the drift and rotation of space and time reference of the BeiDou constellation in autonomous navigation, we may use high-accuracy and high stable hydrogen clocks or cesium clocks onboard as references, and then the autonomous time keeping of satellite clocks can be realized to reduce the time drift of the entire constellation. In addition, the onboard GNSS receiver can receive the signals of other GNSS satellites and determine the onboard satellite orbit; thus, the high-precision orbits of BeiDou satellites can be determined based on the joint adjustment with the autonomous orbit determination data. Therefore, the space and time references of the BeiDou satellite constellation can be maintained, and the datum drift and rotation of the autonomous orbit determination can be reduced.

• 5. Inconsistency of the BDS-3 ionospheric model. As previously analyzed, the B1I and B3I signals of BDS-2 are maintained in BDS-3, and the corresponding ionospheric model adopts the BDSK8 model. However, new signals such as B1C, B2a and B2b use the BDGIM ionospheric correction model. Different models result in the inconsistent of accuracies, which not only causes confusion among terminal manufacturers and users, but also reduces the ranging accuracy of B1I and B3I.

  Regarding the inconsistency of the BDS-3 ionospheric models, using the BDGIM as the unified model of all signals of the BDS-3 can reduce the confusion and inconvenience in users and terminal manufacturers. Since the accuracy of BDGIM model [28] is higher than that of BDSK8, the user equivalent range error can be improved slightly after the unification. Considering that BDS-2 is close to being obsolete, and to not affect the development of the BDS-3 terminal and the PNT service performance of BDS-3 users, we suggest unifying the ionospheric models for BDS-3 the soonest possible time for reducing loss to users and receiver manufacturers. In addition, the ionospheric delay correction model in high-latitude areas should be refined.

• 6. Vulnerability exists inherently in satellite navigation systems, such as the satellite constellation, the ground operational control system (OCS) and the signals. Once the core subsystem fails or malfunctions, e.g., power interruption, time system failure, and other core equipment failures, the PNT service may be interrupted.

To resolve the vulnerability of the BDS constellation, OCS and signals, we can expand PNT information sources by building a comprehensive PNT system [31,32,33], and then use the resilient PNT theoretical framework [34] to realize the resilient integration of multiple PNT sensors, and build the resilient function model, the resilient stochastic model, and resilient data fusion methodology, with the aim of realizing a seamless PNT service from deep space to deep sea, and from outdoor to indoor.

Conclusions

BDS-3 satisfied the requirements of design indexes in orbit determination accuracy, satellite clock accuracy, signal-in-space accuracy and PNT service performance. Additionally, the featured services such as BDSPPP, BDSBAS, regional message communication, and the global SAR function have wide application prospects.

In the future, the inclination angle of IGSO satellites may be increased to improve the service performance of BDS in polar regions. The IGSO satellites may be designed for providing the RSMC and BDSPPP services to overcome the “south wall effects”. The INS payloads are suggested to be added to various satellites, and thus the continuity and availability of satellite orbit parameters can be guaranteed in satellite maneuver. The high accuracy hydrogen clocks or cesium atomic clocks onboard might be used to control the time drift of satellite constellation during the autonomous navigation with the support of ISLs. The two commonly used BDS ionospheric models should be unified for all signals to reduce the users’ confusion and improve the range accuracy for the B1I and B3I signals. The comprehensive and resilient PNT infrastructure should be established for the seamless PNT services.