1. Introduction
The precise orbit and clock products of global navigation satellite system (GNSS) are of great importance for GNSS applications, such as precise point positioning (PPP), GNSS meteorology etc. Currently, the precise GNSS orbit and clock products provided by the international GNSS service (IGS) are generated from the global ground stations. Using ground observations, GPS and GLONASS orbits present a 1D accuracy which is better than 3 cm (
http://www.igs.org/article/official-igs-citation-updated). However, compared with GPS and GLONASS, the orbit accuracy of the new GNSS constellation (such as BDS) is relatively poor [
1]. In the case of BDS constellation, this poor quality mainly results from: (1) special attitude control mode [
2], (2) inaccurate solar radiation pressure (SRP) modeling [
3,
4], (3) relatively few, unevenly distributed station networks, and (4) unfavorable geometry conditions [
5]. In particular, the worst orbit accuracy of 3~4 m can be found for BDS GEO satellites because of their stationary characteristics [
3,
6], which leads to a weak tracking geometry.
Ideally, a globally well-distributed tracking station network is needed to ensure adequate observation strengths for all satellites at every epoch in the GNSS orbit determination. Meanwhile, relative motion between the satellite and ground stations is also necessary to yield a changeable tracking geometry. However, these conditions are not fully available for the BDS constellation due to the limited ground stations and special constellation design of the BDS. Fortunately, the weakness in observability and tracking geometry can be overcome by introducing low earth orbiters (LEO) as moving stations. LEO satellites have become a very important platform in earth observing and communication over the past two decades. Along with the proposing of many large LEO constellation projects, such as Iridium, SpaceX, OneWeb, and Hongyan etc., LEO satellites have been recently gaining widespread attention for the prospect of a LEO-augmented global navigation satellite systems (GNSSs). The constellation, with a large number of LEO satellites, is expected to enhance the positioning service of GNSS, notably reducing the long convergence time of the precise point positioning technique [
7]. Many studies on the availability, design, and performance of LEO-augmented GNSS positioning have been conducted by many researchers [
8,
9,
10,
11,
12,
13]. The aforementioned LEO constellations (such as Hongyan), which aim to provide augmented navigation services, usually carry a state-of-the-art GNSS receiver to obtain the precise knowledge of LEO orbits. These abundant onboard GNSS observations from hundreds (or even thousands) of LEO satellites can be used to improve the orbit accuracy of GNSS satellites through integrated processing with ground observations. The fast motion of LEO can greatly improve the observation geometry of GNSS satellites. Furthermore, the strong algebraic correlations between LEO orbit, GNSS orbits, and the ground reference frame will be introduced in the process of integrated POD.
The integrated precise orbit determination (also called the one-step method) refers to simultaneously estimating the GNSS and LEO orbit parameters as well as other parameters (for example, earth rotation parameters, etc.) using the onboard and ground observations in one parameter estimation process. This method was first validated and investigated by Zhu et al. [
14] using CHAMP and GRACE onboard GPS data. The integrated POD of GPS and LEO satellites have also been reported and analyzed in many other papers. Bookmkamp and Dow [
15] developed a double-differenced, integrated algorithm based on CHAMP and Jason satellites. König et al. [
16] integrated GPS ground observations with CHAMP, GRACE, and SAC-C data to estimate the orbits and the Earth system parameters simultaneously. The results show that the integrated solution achieves better accuracy for both orbits and Earth parameters. Geng et al. [
17] presented an initial comparison between the performances of integrated POD with different ground networks. Zoulida et al. [
18] reported an increase of about 0.7% in the ambiguities fixed rate and better GPS orbits when Jason-2 data were introduced into the orbit determination. Männel and Rothacher [
19] jointly processed a network of 53 global ground stations and four LEO satellites to assess the improvement on geocenter coordinates. Zeng et al. [
20] used a few LEO satellites to discuss the impact of LEO with different orbital altitudes and inclinations, as well as ground station networks, on the GPS orbit determination. In terms of the GPS and BDS dual-system integrated POD, Li et al. [
21] performed an integrated POD of FengYun-3C (FY-3C), GPS, and BDS satellites using FY-3C onboard GPS and BDS observations as well as global and regional ground measurements. The results evidence the contribution of FY-3C onboard observations for BDS orbits, especially for BDS GEO. Based on the Iridium satellite constellation and global network, Li et al. [
10] tested the performance of integrated POD with the GPS, BDS, and LEO constellations using simulated data. Their initial results demonstrated the significant accuracy improvements of GPS and BDS orbits as a result of the LEO constellation.
The aforementioned works evidence the contribution of onboard GNSS observations to the GNSS POD, and shed light on the performance of the integrated POD. However, these studies mainly focus on the integrated POD with a few LEO satellites, since no onboard GNSS observations from large LEO constellations are available at present. Additionally, the impact of different LEO configurations, such as LEO number, orbital altitude, and orbit type, on integrated POD has not yet been well studied. Motivated by the aforementioned studies, we propose a LEO constellation-augmented multi-GNSS precise orbit determination method and assess the LEO-augmented POD performance with different configurations of ground stations and LEO constellations. Within this paper, we mainly focus on the geometric factors that limit the accuracy of GNSS orbit determination. Six LEO constellations with different orbital altitudes, satellite numbers, and orbit types are designed to evaluate the influence of LEO constellations on the GNSS POD. Several station networks with different station number and distributions are employed to discuss the dependency of the integrated POD on ground stations. The paper is structured as follows.
Section 2 starts by introducing GNSS and LEO constellations we used, and gives a detailed description of the sophisticated simulation processing methods used for GNSS observations. Then, the LEO constellation-augmented four-system precise orbit determination algorithm and strategy are introduced in
Section 3.
Section 4 analyzes LEO-augmented GNSS POD results under different tracking conditions in detail. A discussion of our result is given in
Section 5. Finally, conclusions are provided in
Section 6.
4. Result Analysis and Discussion
In this section, the performance of the multi-GNSS integrated POD with different LEO constellations and different station network schemes is analyzed. Several impact factors including LEO numbers, orbital types, and altitude, as well as the numbers and distribution of ground stations, are discussed in detail. The previously simulated orbits are regarded as the true values, and the differences between our estimated orbits and true orbits provide a method by which to assess the accuracy of the integrated POD results.
4.1. Integrated POD with Full LEO Constellation
4.1.1. Integrated POD with Different Numbers of LEO
We firstly performed the integrated POD using global network from MGEX. The ground observations from the 120-station MGEX network are processed as a reference. In order to make a fair comparison, only a 65-station network with a good global distribution from the MGEX we selected is used to perform the integrated POD (shown as red dots with a black circle edge in
Figure 2). The onboard observations from 60, 66, and 96 polar-orbiting LEO satellites with altitudes of 1000 km (see
Figure 1a,b,d) are processed to investigate the influence of LEO numbers on the integrated POD.
Figure 3,
Figure 4,
Figure 5 and
Figure 6 present the average 3D Root Mean Square (RMS) values of orbit differences between integrated POD solution and true orbits for all satellites from GPS, GLONASS, Galileo, and BDS constellations respectively. The result shows that by using only the ground multi-GNSS observations, the majority of GNSS satellites can achieve a 3D orbit accuracy of better than 4 cm, which is comparable to the accuracy of the 24-h orbit recovered from the real ground tracking measurements. This demonstrates that the error models and simulation strategy are close to the real situation.
Table 4 computes the average 3D RMS of orbit differences for four systems. The worst orbit quality can be found in BDS GEO, and the average 3D RMS values of orbit differences even exceeds 2 m. This is reasonable, because the stationary and regional coverage feature of GEO satellites lead to poor tracking geometry, greatly hampering the orbit determination for GEO satellites.
As shown in
Figure 3,
Figure 4,
Figure 5 and
Figure 6, all satellites providing global coverage can achieve a sub-centimeter level orbit with the support of LEO onboard observations, while the orbit accuracy of regional coverage satellites (BDS GEO and IGSO) is at the centimeter level. It can be seen that LEO satellites can make more contributions to GNSS orbits than the same numbers of ground stations, since they can not only provide high-quality observations, but also improve the geometry diversity. The orbit of GNSS satellites obtains a remarkable accuracy improvement after introducing the LEO observations into the POD processing. The accuracy improvement with respect to the ground-based POD result can reach over 70% for all of 60 LEO-, 66 LEO-, and 96 LEO-augmented POD schemes (as shown in
Table 4). This indicates that the onboard data from the large LEO constellations can be an important supplement to the GNSS POD, even under the condition of a global network. We find that the BDS GEO satellites present the largest accuracy improvement when the integrated POD is implemented; the improvement percentage even reaches up to about 98%. This inspiring result evidences the contribution of LEO to improve the tracking geometry of GEO satellites.
Furthermore, it can be seen that smaller orbit differences can be achieved when more LEO satellites are introduced into the integrated POD. The 96-LEO solution achieves a slightly better orbit accuracy than the other two LEO-augmented POD solutions. However, compared with the slight accuracy improvement of the 96-LEO solution, introducing additional 30 LEO satellites into the processing of the integrated POD leads to a huge burden on the computational efficiency because of the larger number of unknown parameters to be solved. The average computation time of integrated POD for the 96-LEO scheme is about 38 h, which is almost three times more than that of 60-LEO and 66-LEO schemes (about 13 h).
4.1.2. Integrated POD with Different LEO Orbital Type
To assess the impact of LEO orbital type, two typical LEO orbits, nearly-polar and sun-synchronous orbits, are adopted in this study. Both of LEO constellations consist of 60 LEO satellites flying at an altitude of 1000 km (see
Figure 1a,b). The ground observations from a global network with 65 stations are used. The results are shown in
Figure 7. It can be seen that the sun-synchronous-orbiting constellation presents a slightly stronger enhancement to the GNSS orbits than the polar-orbiting constellation. With the sun-synchronous-orbiting data of 60 satellites, the orbit of GNSS satellites can achieve an accuracy similar to that of the 96-polar-orbiting scheme, and the orbit accuracy of BDS GEO for the sun-synchronous scheme is even better than that for the 96-polar-orbiting scheme.
As reported by Li et al. [
9], the polar-orbiting constellation outperforms the sun-synchronous-orbiting constellation in terms of the convergence time of PPP. However, the situation is reversed when the LEO constellation is employed to enhance the GNSS precise orbit determination. In order to determine the reason for this, we computed the actual number of LEO satellites used for the GNSS POD and calculated an orbit dilution of precision (ODOP) value for every GNSS satellite using a similar algorithm like the position dilution of precision (PDOP) for ground stations. Assuming one GNSS satellite can be tracked by
ground stations simultaneously, the ODOP of this GNSS satellite can be calculated as follows:
where
and
are the positions of ground station and GNSS satellite respectively.
is the distance vector with the direction from GNSS satellite to ground station. The ODOP value describes the geometric strength of observations for GNSS satellites. The better the geometric strength, the smaller the ODOP value.
Figure 8 shows the used LEO numbers and ODOP values for C02 (GEO) and G13 (MEO) respectively. The two 60-LEO schemes contribute an almost comparable number of LEO satellites to the orbit determination for BDS GEO C02. The average number of the used sun-synchronous-orbiting LEO satellite is 21.5, which is slightly larger than polar-orbiting LEO number, i.e., 20.6. In terms of ODOP, the ODOP values of the sun-synchronous-orbiting constellation present a stronger fluctuation with a larger amplitude than those of polar-orbiting constellation. For BDS GEO satellites, the changeless tracking geometry is the key limitation of their orbit determination. The main contribution of LEO satellites is to bring a fast variation to the tracking geometry for BDS GEO. The faster fluctuation of ODOP indicates that the sun-synchronous-orbiting constellation can bring more variations to the GEO tracking geometry, which can yield more pronounced orbit accuracy improvements for GEO satellites. For the 96-LEO solution, more available LEO satellites are introduced compared with the two 60-LEO solutions. Though the 96-LEO solution presents smaller ODOP values, a more moderate variation can be observed in GEO ODOP values for the 96-LEO solution compared with the 60 sun-synchronous-orbiting solution. The slow variation of ODOP values for the 96-LEO scheme may be the reason why the 96-LEO solution achieves a relatively worse orbit quality of GEO than the 60 sun-synchronous-orbiting solution. In the case of satellite G13, the differences between the ODOP variation of the two 60-LEO solutions is small, which leads to a relatively small discrepancy in the corresponding orbit accuracy, i.e., 1 mm.
4.1.3. Integrated POD with Different LEO Altitude
Table 5 lists the result of the integrated POD with LEO satellites at different altitudes. The 60 polar-orbiting LEO satellites at altitudes of 600, 1000, and 1400 km are considered (see
Figure 1a). It can be seen that there are almost no orbit accuracy differences among the solutions. The integrated POD with higher LEO satellites presents a slightly higher orbit quality than that with lower LEO satellites, except for the BDS GEO satellites. For BDS GEO, the best orbit accuracy is achieved by the 600 km solution. Given the motionless characteristic of GEO satellites, the LEO satellites orbiting at an altitude of 600 km can contribute more improvements to the tracking geometry for BDS GEO due to their faster motion, which may be the reason for the better performance of the 600-km solution in GEO orbit accuracy. In addition, it should be noted that the LEO satellites at low altitude usually suffer from complex SRP modelling. Also, the trajectory of low-altitude LEO is perturbed by a stronger drag force, which will increase the difficulty of drag force modeling and lead to accelerated orbital decay for LEO. The difficulty of air drag and SRP modelling may negatively affect the result of the integrated POD. Fortunately, with the proper modeling of these forces, such as estimating the drag parameter with a shorter interval, low-altitude LEO satellites can obtain a high precision orbit result, which has been evidenced in many low-altitude LEO missions, such as Swarm [
38].
4.1.4. Integrated POD with Different Ground Network
In order to evaluate the impact of station quantity and distribution on the LEO-augmented POD, fewer global stations are selected to perform the integrated POD with different LEO constellations. We firstly adopted about 22 global stations from iGMAS (see
Figure 2). Typically, the three LEO schemes, i.e., 60 polar-orbiting LEO, 60 sun-synchronous-orbiting LEO, and 96 polar-orbiting LEO constellations at the altitude of 1000 km, are chosen.
Figure 9 illustrates the integrated POD results for GPS, GLONASS, Galileo, and BDS using iGMAS ground observations. For comparison, we also plot the GNSS orbit results determined from the 65 MGEX stations and 60 sun-synchronous-orbiting LEO satellites. The corresponding statistical values are provided in
Table 6. With regard to the numbers provided in
Table 4, a clear degradation in the orbit accuracy can be recognized for all GNSS satellites compared with the POD result using MGEX network when only ground iGMAS observations are processed. This is reasonable, because fewer ground stations are employed to recover GNSS satellites orbit. Similar to the MGEX solution, the significant orbit accuracy improvement for all satellites can be recognized after the incorporation of LEO onboard observations into the GNSS orbit determination. Certainly, BDS GEO satellites exhibit the largest orbit accuracy improvement. The different LEO constellations present a similar performance to that in the MGEX solution. As shown in
Table 6, the sun-synchronous-orbiting constellation yields slightly better performance in LEO-augmented POD compared to the two polar-orbiting constellations.
Meanwhile, it should be noted that although a sparse iGMAS network is employed to perform the integrated POD, the orbit accuracy of all three iGMAS integrated POD schemes is comparable to that of the corresponding MGEX solutions. The performance of the iGMAS and MGEX solutions mainly differs in their estimated BDS GEO orbits. Compared with the MGEX solution, the integrated POD with the iGMAS network provides a relatively worse orbit for the BDS GEO. The result demonstrates that the dependence of GNSS orbit determination on the ground stations can be largely reduced when introducing a large number of LEO satellites into the POD processing. Serving as moving stations, LEO satellites can not only provide a large amount of onboard observations, but also bring evident improvements to the tracking geometry for GNSS satellites.
The above results indicate that the onboard observations play an important role in reducing the contribution of ground observations to GNSS orbit determination when introducing a large number of LEO satellites. To further investigate the dependency of the integrated POD on the ground stations, we design three integrated POD schemes with only 8 or 4 ground stations involved. The distribution of the small number of stations we selected is shown in
Figure 10. The sun-synchronous-orbiting constellation with 60 LEO satellites is adopted due to its best performance in the previous study. It should be noted that the earth rotation parameters are fixed when we perform the integrated POD with a few stations.
As shown in
Figure 11, the estimated orbit of the integrated POD with 8 regional stations can achieve a centimeter-level accuracy, which is better than the POD result only using ground MGEX observations (see
Table 4). This indicates that with the assistance of 60 LEO satellites, only using observations from regional stations can already obtain a relatively high-quality satellite orbit. A significant reduction in orbit difference can be recognized for all the GNSS satellites when the 8 regional stations are replaced by 8 global stations. We can find that the orbit accuracy of the integrated POD with 8 global stations is almost in line with the accuracy of the corresponding MGEX and iGMAS integrated POD schemes. The accuracy improvement can be attributed to the better performance of the globally distributed stations in anchoring the whole constellations compared to the regional stations. Slight orbit degradation for GNSS satellites is observed when the number of global stations is reduced to 4. The result shows that the POD of GNSS satellites has less dependence on the number of ground stations and is more sensitive to the distribution of ground stations after the inclusion of the large LEO constellation.
4.2. Integrated POD with Partial LEO Constellation
The previous results demonstrate the strong enhancement of the entire LEO constellation to the GNSS orbit estimation. However, it can be seen that the integrated POD suffers from a long computation time when a large number of LEO satellites are introduced. The computation time of the 60-LEO scheme is about 13 h, while for 96-LEO scheme, it reaches about 38 h. Only moderate accuracy differences are obtained for the 60-, 66- and 96-LEO constellations. In order to balance the conflict between orbit accuracy and computational efficiency, it is worthwhile to investigate the integrated POD with a partial LEO constellation, and to explore how many LEO satellites can be used as an auxiliary to improve the orbit accuracy of GNSS satellites. In this section, we process the observation from 65 MGEX stations and a partial LEO constellation. The 60 polar-orbiting LEO constellation and 60 sun-synchronous-orbiting LEO constellation at an altitude of 1000 km are chosen.
Figure 12 shows the average 3D RMS of GNSS orbit differences for the integrated POD as a function of the number of LEOs. We find that with only adding 10 LEO satellites, BDS GEO satellite orbits can already achieve an accuracy of better than 10 cm, while the orbit differences of other satellites are below 2 cm. The corresponding computation time is less than 2 h (shown in
Figure 13). This result indicates that the introduction of a small number of LEO satellites can not only obtain a relatively high accuracy GNSS orbit, but that it also limits the computation time, which is very beneficial to the GNSS application with high timeliness requirements. For both, the orbit differences of GNSS satellites decrease monotonically as the number of introduced LEO satellites increases. However, with more LEO satellites included, the improvement of orbit accuracy gradually becomes smaller. A reduction of about 0.5–4 cm can be observed in the orbit differences when the LEO increases from 10 to 20, whereas the additional 10 LEO satellites can only contribute to an accuracy improvement of less than 0.2 cm as the LEO increases from 40 to 50. On the other hand, the increase of LEO satellites introduces a large number of parameters to be estimated, resulting in a steep rise in computational time (shown in
Figure 13). For example, the computation time of the integrated POD greatly rise from 9.46 h to 13.42 h when the LEO number increases from 50 to 60. But the addition of 10 LEO satellites only yields an orbit accuracy improvement of about 1 mm for the majority of GNSS satellites.
In addition, it can be seen that the sun-synchronous LEO satellites exhibit a stronger enhancement to the GNSS orbits than those in nearly-polar orbit, no matter how many LEO satellites are introduced. The accuracy differences between the sun-synchronous LEO scheme and nearly-polar LEO scheme is very small for the GPS, GLONASS, Galileo, BDS IGSO, and MEO satellites. This is because the relative motion between these satellites and ground stations already yields good geometric diversity, so that the small contribution differences in tracking geometry between sun-synchronous LEO and nearly-polar LEO satellites have little impact on the orbit determination of these satellites. However, differences of more than 1 cm are visible in the BDS GEO orbits between the two solutions. Evidently, the combination of global stations and sun-synchronous LEO satellites result in better quality GEO orbits. This is because LEO satellites in sun-synchronous orbit can give rise to faster variations of observation geometry for BDS GEO satellites (see
Figure 8).
5. Discussion
Serving as moving stations, the LEO satellites can significantly strengthen the tracking geometric diversity of GNSS satellites and improve the accuracy of the GNSS precise orbit determination, particularly in the case of regional network or sparse global network. Although many studies have reported the contribution of LEO onboard observations to the GNSS POD, few studies have focused on the performance of the integrated POD with a large LEO constellation. Meanwhile, the impact factors of the integrated POD, such as LEO number, orbital height, and orbit type, is barely analyzed. In this study, we presented the GNSS orbit solution derived in an integrated processing of the ground network and the large LEO constellation using simulated data. Six LEO constellations are designed to study the performance of the integrated POD with different tracking configurations. The results demonstrate that the orbit accuracy of GNSS satellites can be dramatically improved when a large number of LEO satellites is introduced; the more LEO satellites, the better the orbit accuracy. Notably, BDS GEO orbits present an accuracy of a few centimeters with the recruitment of LEO satellites. Moreover, the LEO satellite in sun-synchronous orbit can contribute more to GNSS orbits than nearly polar-orbiting satellites, since they give rise to a stronger observation geometry. Furthermore, the orbit altitude of LEO satellites is found to have no evident impact on the GNSS satellites except BDS GEO satellites. In addition, it can be seen that the integrated POD with the entire LEO constellation suffers from a long computational time because of the abundant GNSS observations to be processed and the large number of parameters to be estimated. The long computation time of the integrated POD cannot meet the timeliness requirement of the real-time precise positioning. This weakness in terms of computational efficiency can be overcome by just introducing a certain number of LEO satellites.
Although our study presents an optimistic integrated POD result due to the simulated data, the impact of introducing LEO constellations to improve GNSS orbit determination is clearly demonstrated. Indeed, not only tracking geometry, but many other factors, such as attitude model, solar radiation pressure and drag models, and antenna calibrations, constitute key limitations in real-world LEO and GNSS POD. Fortunately, with the efforts of many scientists, the force and observation models have been rigorously refined in recent years, which can reduce the impact of model errors on the GNSS POD as much as possible. Using the current models, GNSS satellites (except BDS GEO) can achieve centimeter-level orbits [
1], though more efforts are needed to further refine the force and observation models. Meanwhile, the force model errors do not represent a weakness in the observability and tracking geometry, which is a separate issue from the tracking geometry that needs to be further investigated. In addition, except for the factors mentioned above, there are still many issues which need to be taken into consideration when implementing the integrated POD using real GNSS observations. For example, the incorporation of a large number of real multi-GNSS onboard observations may introduce unexpected unmodeled errors into the orbit recovery process. Meanwhile, different from ground receivers, the performance of space-borne receivers can be more easily affected by the space environment due to the high altitude, i.e., 600 km~1000 km, where LEO satellites are flying. As discussed in Xiong et al. [
39], large equatorial plasma irregularities can degrade the tracking capability of the Swarm onboard receivers, leading to severe signal loss at low latitudes. This loss of signal for space-borne receivers undoubtedly has negative effects on the integrated POD. Moreover, except for GNSS data, many other types of satellite tracking data are not considered in our study. Once the construction of the navigation augmentation LEO constellation is complete, the new ranging observations from LEO satellites to ground stations can be employed. These new types of satellite tracking data are expected to remove the potential systemic bias in the estimated orbits and further improve the orbit accuracy of both GNSS and LEO satellites.
6. Conclusions
This paper is devoted to investigating the integrated precise orbit determination with large LEO constellations. We performed LEO-augmented multi-GNSS precise orbit determination in this study. Based on simulated GNSS observations, several integrated POD schemes are implemented to investigate the performance of integrated POD under different tracking conditions. The potential influence factors of the integrated POD, including LEO satellites number, altitude, and orbit type as well as ground station number and distribution, are analyzed and discussed in detail.
The result shows that joint processing ground observations from the global network and the large LEO constellation onboard observations can significantly improve the orbit accuracy for GNSS satellites, especially for BDS GEO satellites. The accuracy improvement percentage with respect to the ground-based POD results can reach over 70% for all the integrated POD schemes with 60-, 66-, 96-LEO satellites. The largest orbit accuracy improvement of over 98% can be recognized for BDS GEO, since the fast motion of LEO satellites brings tremendous variations to the tracking geometry of GEO satellites. By incorporating a large number of LEO satellites, BDS GEO satellites can obtain a centimeter-level orbit. Compared with the 60- and 66-LEO schemes, a slightly better orbit quality is observed in the 96 LEO scheme due to the introduction of more LEO satellites. However, the increase of the involved LEO satellites results in a sharp rise in the computational time of the integrated POD because of more unknown parameters to be solved. We also present the impact of the LEO orbit type on the integrated POD. With the same number of LEO satellites, the sun-synchronous-orbiting constellation presents a stronger enhancement to GNSS orbits than the polar-orbiting constellation. The improvement can be attributed to the more rapid variation of GNSS tracking geometry provided by a sun-synchronous-orbiting constellation. In terms of the LEO altitude, the orbit altitude of LEO satellites is found to have little influence on the enhancement to GNSS orbits except for with BDS GEO satellites. Benefiting from the faster motion of lower LEO satellites, the 600 km scheme achieves a better orbit accuracy for BDS GEO than the 1000 km and 1400 km schemes.
In order to assess the impact of station number and distribution, global and regional networks with different numbers of stations are employed. With the inclusion of large LEO constellation, the integrated POD appears to be more sensitive to the distribution of ground stations, rather than the station number. This is because in the integrated POD, the dependence of GNSS orbit determination on ground stations is largely reduced by the onboard GNSS observations when a large LEO constellation is considered. Based on 8 regional stations, the orbit accuracy of the integrated POD with a sun-synchronous-orbiting constellation is at the centimeter level. However, a subcentimeter accuracy can be recognized for the orbits of GNSS satellites when 8 globally distributed stations are adopted, which is comparable to the result of the integrated POD with a denser global network of 65 stations.
Although GNSS orbit estimations can greatly benefit from the incorporation of a full LEO constellation, the integrated POD suffers from long computation time due to the addition of a large number of LEO satellites. To balance the conflict between orbit accuracy and computational efficiency, the integrated POD with a partial LEO constellation is investigated. The orbit accuracy of GNSS satellites improves gradually as the LEO satellites increase. The more LEO satellites, the better the orbit accuracy. However, with an increasing number of LEOs, the accuracy improvement of GNSS orbits becomes smaller. Meanwhile, the increase of LEO satellites results in a steep rise in computational time. Considering both orbit accuracy and computational efficiency, we prefer 40 LEO satellites in a sun-synchronous orbit.
With the rapid development of the LEO constellations, we expect an improvement of GNSS POD when jointly processing GNSS observation data from reference sites and LEOs. Continued efforts for force modeling, algorithm optimization, and scheme design are required to achieve better performance of the integrated POD.