Performance Assessment of Multi-GNSS PPP Ambiguity Resolution with LEO-Augmentation

: The fast motion of low Earth orbit (LEO) satellites provides rapid geometric changes in a short time, which can accelerate the initialization of precise point positioning (PPP). The rapid convergence of ambiguity parameters is conducive to the rapid success of ambiguity fixing. This paper presents the performance of single-and four-system combined PPP Ambiguity Resolution (AR), enhanced with an ambiguity-float solution LEO. Two LEO constellations were designed: L was a typical polar orbit constellation, with a higher number of visible satellites at high latitudes than at low and middle latitudes; and M was designed to compensate for the lack of visible satellites at low and middle latitudes. The ground observation data of the LEO satellites at the MGEX stations were simulated. Because the global navigation satellite systems (GNSSs) were fully operational, the GNSS data were real observation data from the MGEX stations. Based on the daily observation datasets collected at 258 stations in the global MGEX observation network over three days (from January 1 to January 3 2022), in addition to the LEO simulation data, we evaluated the positioning performance of LEO ambiguity-float solution-enhanced PPP ambiguity resolution and compared it with LEO-enhanced PPP. The L+M mixed constellation was able to reduce the time to first fix (TTFF) of the four-system combined PPP-AR to 5 min, and four LEO satellites were sufficient to achieve this. L+M mixed constellation was able to redu ce the convergence time of the four-system combined PPP to 2 min. Unlike PPP-AR, PPP required more LEO satellites for augmentation to saturate.


Introduction
Precise point positioning (PPP) [1] technology has been widely used in many fields, such as precision agriculture, precise mapping, deformation monitoring, water vapor inversion, precise navigation, and positioning [2], owing to its high positioning accuracy, simple structure, and flexible operation (without any reference station). However, it usually takes tens of minutes to achieve centimeter-level positioning accuracy. Reliable ambiguity resolution (AR) [3,4] also requires a longer observation time, due to the obstacle of atmospheric errors, especially ionospheric delays. Regional augmentation precise point positioning, also known as precise point and real time kinematic positioning (PPP-RTK) [5], utilizes precise atmospheric corrections (generated using the reference network) to accelerate ambiguity resolution [6]. Unfortunately, there are places where reference stations cannot be established or available, such as in oceans, deserts, or developing areas.
Many efforts have been made to improve PPP performance in terms of observations, such as combining multiple global navigation satellite systems (GNSSs) and using multifrequency signals. GPS/GLONASS combined PPP can improve the convergence speed [7], but the positioning accuracy is not as good as the GPS-only PPP when the orbit and clock offset quality of the GLONASS satellite is poor [8]. At the same time, GPS/GLONASS combined PPP can also accelerate the convergence of ambiguity parameters, thereby shortening the time to first fix (TTFF) of ambiguity [9,10]. For GPS/BDS/GLONASS/Galileo four-system combined PPP, compared with single-system PPP, multisystem combined positioning not only converges faster but also improves the availability of PPP technology in special environments, as well as the reliability of positioning results [11][12][13]. In addition to multisystem fusion, multifrequency observations can improve the performance of PPP. Based on simulated GPS three-frequency data, Geng et al. [14] conducted PPP fixed-solution experiments and pointed out that three-frequency PPP could achieve ambiguity-fixed solutions in a few minutes. Elsobeiey et al. [15] conducted PPP float solution experiments based on real GPS three-frequency data, and the results showed that the convergence time of three-frequency PPP could be shortened by 10%. Li et al. [16] found that the convergence speed and positioning accuracy of three-frequency PPP ambiguity resolution was better than that of three-frequency PPP and dual-frequency PPP-AR in the initialization stage, through a BDS three-frequency uncombined PPP-AR experiment. Li et al. [17] evaluated the performance of a BDS/Galileo combination with three-frequency PPP-AR, and the results showed that multi-GNSS fusion could further improve the position estimation capability of a single system. Not only were the observations themselves found to be beneficial for reducing the convergence time in PPP, but so too was the incorporation of external ionospheric corrections from multiscale network configurations [18,19]. However, due to the slow geometric change in the satellites' position relative to the ground station, PPP initialization still took approximately 10 min.
In addition to multi-GNSS and multifrequency satellites, low Earth orbit (LEO) satellites are another way to improve PPP performance. In recent years, commercial companies have proposed their own LEO constellation, consisting of hundreds to thousands of LEO satellites that would provide broadband Internet services worldwide. Reid et al. [20] demonstrated the potential of the LEO constellation as a navigation system, which was also confirmed in the first generation of the Iridium constellation. This constellation was the first low Earth orbit satellite constellation to integrate communication and navigation functions. It comprises 66 satellites distributed in six orbital planes, with an orbital altitude of 781 km and an orbital inclination of 86.4°. Compared with the ordinary GNSS medium-Earth orbit (MEO) satellite, the LEO constellation has the following advantages: First, the LEO satellite has a lower orbital height, being 20 times closer to Earth than the MEO satellite, in addition to a shorter signal transmission path, and a smaller degree of power loss. Therefore, the signal power received on the ground is stronger, and the antiinterference and anti-deception performance is better. For example, the signal power of the Iridium satellite is 1000 times that of GPS [21]. Therefore, even in highly shielded environments (such as indoors, in forests, under bridges, and inside canyons), users can still receive continuous LEO navigation signals. Second, the LEO satellite moves faster than the ground station, and the continuous observation duration is only 10-20 min. The Doppler frequency shift phenomenon is evident, which is favorable for the carrier phase cycle slip detection [22]. At the same time, the fast motion of LEO satellites provides rapid geometric changes in a short period of time, which can accelerate the initialization of PPP, thus shortening the convergence time of precise point positioning.
To prove the contribution of LEO to PPP, many researchers have conducted numerous LEO-enhanced precision point positioning experiments based on simulated observations. Ke et al. [23] tested the performance of LEO-enhanced GPS PPP and compared it with that of GPS or GPS/GLONASS combination schemes. The results showed that the convergence time of the LEO-enhanced GPS and GPS/GLONASS combined scheme was reduced by 51.31% and 3.93%, respectively, compared with the GPS-only scheme, indicating that the contribution of the LEO satellite to the convergence speed of the PPP was greater than that of a MEO satellite. Ge et al. [24] found that LEO-enhanced GNSSs (GPS/BDS/Galileo) could reduce the convergence time of PPP to 5 min. Li et al. [25] studied the performance of LEO-enhanced GNSSs (GPS/GLONASS/BDS/Galileo) under LEO constellations with different satellite numbers (60,96,192,288), orbital heights (1000 km, 600 km), and orbital types (polar orbit, sun-synchronous orbit), and proposed a LEO constellation scheme consisting of 192 satellites with an orbital altitude of 1000 km, which could significantly shorten the convergence time of GRCE PPP by more than 70%. Li et al. [26] analyzed the performance of PPP ambiguity resolution with LEO-augmented multi-GNSS (GPS/GLONASS/BDS/Galileo), and the results showed that with the enhancement of 60, 192, and 288 LEO constellations, the TTFF of GREC PPP-AR was shortened from 7.1 min to 4.8 min, 1.1 min, and 0.7 min, respectively, and the positioning accuracy was improved by approximately 60%, 80%, and 90%, respectively. The performance of the threefrequency LEO-enhanced PPP-AR was better than that of dual-frequency. Zhao et al. [27] found that in harsh environments, the float solution and fixed solution performance of triple-frequency uncombined BDS/GPS/LEO PPP were improved to different degrees compared with that of only BDS, which they ascertained through six sets of shielding experiments using different elevation and azimuth angles.
The above studies show that the LEO constellation can significantly improve the performance of PPP and PPP-AR; however, few studies have examined the contribution of the LEO float solution to PPP-AR. As mentioned above, the fast motion of LEO satellites provides rapid geometric changes in a short time, which can accelerate the initialization of PPP, and the fast convergence of ambiguity parameters is conducive to the rapid success of ambiguity fixing. Therefore, this study focused on the performance of the LEO float solution to enhance PPP-AR. The remainder of this article is structured as follows: The two LEO constellations we designed and the LEO observation simulation on the Multi-GNSS Experiment (MGEX) stations are introduced in Section 2. Because GNSSs (GPS/GLONASS/BDS-3/Galileo) were fully operational, the GNSS data are real observation data from MGEX stations. The functional model of LEO-enhanced GNSSs PPP and PPP-AR is introduced in Section 3. In Section 4, we introduce the positioning performance of single-system and the four-system combined PPP-AR enhanced using the LEO ambiguity-float solution, and then compare it with LEO-enhanced PPP. A summary of our results is presented in Section 5.

LEO Constellation Design and Ground Observation Simulation
In this study, two LEO constellations were designed: L and M, both of which were Walker constellations [28]. L consisted of 70 satellites distributed in six orbital planes, with an orbital altitude of 1000 km and an orbital inclination of 90°. M consisted of 90 satellites distributed in 10 orbital planes, with an orbital altitude of 1000 km and an orbital inclination of 60°. The main difference between the two LEO constellations was their orbital inclinations. Figure 1 shows the daily average number of visible satellites of the LEO constellation on 1 January 2022, at a cutoff elevation angle of 7°, with L, M, and L+M from top to bottom. The number of visible satellites from the L constellation in high latitudes was much greater than that in middle and low latitudes, because the L constellation is a polar orbit constellation with an inclination of 90°. The M constellation compensated for the lack of visible satellites in the low and middle latitudes of the L constellation, but it had no coverage in the polar regions, which was also caused by orbital inclination. This study conducted LEO ground-observation simulations for all stations in the MGEX observation network. The simulation data included two frequencies, L1 (1575.42 MHz) and L2 (1227.60 MHz). Errors, such as satellite clock offsets, receiver clock offsets, relativistic effect, phase wind-up, tides, satellite and receiver antenna phase center offsets and variations, tropospheric delay, Earth rotation, and measurement noise, were considered. The dry and wet components of the tropospheric delay in the zenith direction were calculated using the Saastamoinen model [29], and the GMF projection function [30] was used to calculate the total tropospheric delay in the signal propagation direction. The measurement noise of the pseudorange and the carrier phase were 0.5 and 0.003 meters, respectively.

LEO-Enhanced GNSSs PPP
where the subscript r is the station number, j represents the frequency point, the superscript s is the satellite number, and K represents GPS, GLONASS, BDS-3, or Galileo.
, sK r  indicates the geometric distance between the station and satellite. c is the speed of light in a vacuum, r dt is the receiver clock offset, and are the measurement noise errors of the code pseudorange and carrier phase, respectively. There were also some potential error items, such as Earth rotation, satellite and receiver antenna phase center variations and offsets, antenna wind up, solid tide, ocean tide, atmospheric tide, and relativistic effect, which were corrected according to existing models [31]; therefore, they have not been listed in the above equation.
To eliminate the first-order ionospheric delay of the code pseudorange and carrier phase observation, the ionospheric-free (IF) combination of  2  2  2  2  2  , , , 2  2  2  2  2  , , , where , 1 sK f and , 2 sK f denote the different frequency points. In this study, GPS used L1 and L2, GLONASS used G1 and G2, BDS-3 used B1 and B3, and Galileo used E1 and E5a, which was consistent with the International GNSS Service (IGS) analysis center. The observation equation after the combination is ,IF   ,  ,IF   ,  ,  ,  ,  ,  ,IF  ,IF  IF   ,  ,  , To reduce satellite orbit and clock errors, precision ephemeris and precision satellite clock products are generally used in PPP. Because the precision satellite clock error is generated from the IF combined pseudorange observation data of the GNSS ground tracking network, the c of the IF combination at the satellite end. The code hardware delay at the receiver end varied with different satellite systems (such a difference is called intersystem bias), and it will be absorbed by the receiver clock error; therefore, the receiver-clock offsets of different satellite systems are different. The ambiguity parameter is closely related to the phase hardware delay. They all have constant characteristics, which vary between different satellites, making it difficult to separate them from each other. Simultaneously, owing to the introduction of the satellite clock error based on pseudorange observation, the code hardware delay at the satellite end was introduced into the phase observation equation, being then absorbed by the ambiguity parameter. Because the pseudorange observation equation and phase observation equation work together to estimate the receiver clock error, the code hardware delay at the receiver end was also introduced into the phase observation equation and absorbed by the ambiguity parameter. Therefore, because the satellite clock error was corrected by the precise product, Equation (5) where the reparameterized receiver clock error and ambiguity are shown in Equations (12)- (15). In addition, the three-dimensional position of the receiver and the tropospheric wet delay in the zenith direction have been estimated: , , , Comprehensively, considering the error source of the LEO simulation data, the observation equation of LEO-enhanced GNSSs PPP is as follows: ,  ,IF   ,  ,IF   ,  ,IF   ,  ,IF   ,  ,  ,  ,IF   ,  ,  ,  ,  ,IF  IF  ,IF ,, In the formula, the superscript L represents the LEO system. The state vector of unknowns is r , which will be estimated using a Kalman filter approach. Because the tropospheric dry and wet delays of the LEO simulation data are calculated using the Saastamoinen model (see Section 2), the LEO system does not need to estimate the tropospheric delay parameters, but instead directly uses the model to correct them. A stochastic model based on elevation angle [32] has been adopted, and the cutoff elevation angle was 7 degrees.

PPP-AR
Owing to the complexity of the frequency division multiple access system, only the code division multiple access system was considered for the ambiguity-fixed solution. In Equation (13), the reparametrized real-value ambiguity of the IF combination can be expressed as a combination of wide-lane ambiguity and narrow-lane ambiguity [4]: The wide-lane ambiguity In this formula, In the above formula, the narrow-lane FCB absorbs the wide-lane FCB. Narrow-lane FCBs can be obtained using an estimation process similar to that of wide-lane FCBs.
After accurate estimation of WL and NL FCBs, PPP-AR can be realized. According to Equation (20), the IF combination ambiguity is reversed, and the intersatellite single difference is used to eliminate the narrow-lane FCB at the receiver end to obtain Equation The LAMBDA search method was used to fix the integer ambiguity of the narrowlane , NL st f , and the ratio value was set to 2. Taking the above equation as the constraint of the IF combined observation equation, the ambiguity-fixed solution of the PPP is calculated.

Analysis
To analyze the performance of the LEO float solution to enhance PPP-AR, observation data of 258 stations in the MGEX observation network from January 1 to 3, 2022 were selected, as shown in Figure 2. The experiments of single-system and four-system combined PPP-AR enhanced using LEO ambiguity-float solution were carried out for these stations and compared with LEO-enhanced PPP. The detailed data processing strategy is shown in Table 1. In Figure 2, the red dots indicate that we conducted LEO-enhanced GPS PPP and PPP-AR experiments on this site, the green dots indicate that we conducted LEOenhanced Galileo PPP and PPP-AR experiments on this site, the blue dots indicate that we conducted LEO-enhanced GLONASS PPP experiments on this site, and the yellow dots indicate that we conducted LEO-enhanced BDS-3 PPP experiments on the site. The presence of four colors in one dot indicates that we carried out LEO-enhanced four-system combination PPP and PPP-AR experiments on the site. In all experiments, the reference coordinates of the experimental results were the coordinate values in the SNX weekly solution file provided by the IGS. Supplementary information, such as the precise satellite orbit and the clocks (post-processed products), were provided by the IGS data center of the German Research Centre for Geosciences.

Quality of FCB Products
We used observation data from 152 globally distributed MGEX stations to estimate GCE FCBs, in addition to evaluating the daily stability of NL FCBs, as shown in Figure 3. NL FCBs in Galileo were the most stable, with all NL FCBs remaining within 0.1 cycles in a day. Among GPS satellites, the G23 satellite showed a range of 0.3 cycles, while other satellites remained stable within 0.1 cycles. For BDS-3, the stability was relatively poor. Satellites C34 and C46 had the largest changes, up to 0.6 cycles, while other satellites had changes within 0.4 cycles. Overall, BDS-3 exhibited the worst NL FCB stability among the three systems. The average STD of the NL FCBs for G, C, and E were 0.016, 0.054, and 0.014 cycles, respectively.

PPP and PPP-AR Enhanced with Different LEO Constellations
The positioning performance of both PPP and PPP-AR (enhanced using the L, M, and L+M constellations) was evaluated using ground observation simulation data of different LEO constellations. MGEX station ASCG (7.92°S, 14.33°W) was selected as a typical example. The results of the single-system and four-system combined float solution, in addition to the fixed-solution PPP enhanced using different LEO constellations, for this station on 3 January 2022, are shown in Figure 4a,b. It can be seen that the addition of the LEO satellite could significantly accelerate the convergence of PPP, and effectively accelerated the first fix of PPP-AR. The convergence speeds in the east (E), north (N), and up (U) directions were clearly improved, and the fixed solutions exhibited similar results. Compared with single-system PPP, four-system combined PPP converges faster, leaving less room for improvement of LEO satellites. Similarly, the improvement of the fixed solution was not as significant as that of the float solution. The L+M mixed constellation was the best, followed by the M constellation, and then the L constellation.   Figure 5 shows the average convergence time of PPP enhanced using different LEO constellations. In this study, the convergence time was set as the time when the positioning error in the E, N, and U directions was less than 0.1 m and lasted for 10 min. By introducing the L, M, and L+M constellations, the convergence time of single-system PPP was shortened from 29.1 to 9.7, 5.7, and 2.3 min, and the convergence time of the four-system combined PPP could be shortened from 11.4 to 6.1, 3.0, and 1.8 min. The L+M constellation showed the shortest times, followed by the M constellation, and then the L constellation. This was because the L constellation had few visible satellites in almost the entire midlow latitudes, and the M constellation had only a blank coverage in the polar region. The hybrid constellation combined the advantages of the L constellation in high latitudes and the M constellation in middle and low latitudes, with an abundance of visible satellites on a global scale, as shown in Figure 1. For the positioning scheme with only the BDS-3 system, the convergence speed of PPP improved significantly. After the introduction of the L, M, and L+M constellations, the convergence time was reduced from 45 min to 2 min, which is a 96% reduction.  Figure 6 shows the average TTFF of the PPP-AR enhanced with different LEO constellations. In this study, TTFF was defined as the time required for at least five epochs to successfully achieve a fixed ambiguity solution. By introducing the L, M, and L+M constellations, the TTFF of single-system PPP could be shortened from 14.7 to 10.4, 9.4, and 7.9 min, respectively, and the time to first fix of the four-system combination PPP could be shortened from 8.6 to 6.8, 4.8, and 4.5 min, respectively. Similar to the float solution, the L+M constellation was the best, followed by the M constellation, and then the L constellation. The fixed solution showed only a 46% improvement, compared with a 92% improvement for the single-system float solution.  Table 2 shows the average positioning accuracy of different LEO constellation-augmented PPP and PPP-AR in the east, north, and up directions, respectively. Again, the reference coordinate of this study was the coordinate value in the SNX weekly solution file provided by the IGS, and the positioning accuracy was represented by the positioning error of the last epoch. The positioning accuracy of ambiguity-float solution PPP in E, N, and U directions was 0.40, 0.34, and 1.20 cm, respectively, and that of ambiguity-fixed solution PPP was 0.28, 0.27, and 1.00 cm, respectively. With the help of the LEO constellation, the positioning accuracy of the float and fixed solutions was improved to varying degrees. The L+M constellation in the E and U directions was superior to the M and L constellations. However, the overall improvement was very limited. Even in the U-direction with the poorest positioning accuracy, the improvement was less than half a centimeter.

LEO-Enhanced PPP and PPP-AR in Different Latitudes
Considering that the enhancement of performance at different locations may vary, the 258 stations were divided into low-, middle-, and high-latitude zones for analysis. Figure 7a shows the average convergence time of PPP enhanced using the same LEO constellation at different latitudes. The average convergence time of the L constellation-enhanced PPP in low, middle, and high latitudes was 14.1, 7.8, and 3.3 min, respectively. The higher the latitude, the faster the convergence speed. This was because the polar orbital constellation was mostly at higher latitudes, so it provided more visible satellites, as shown in Figure 1a. The average convergence time of the M constellation in the middle-latitude region was 4.6 min, which was slightly faster than the 5.4 min of the low-latitude region and the 7.0 min of the high-latitude region. This was because the M constellation provided the largest number of visible satellites in midlatitudes, as shown in Figure 1b. The average convergence times of the L+M constellation in low, middle, and high latitudes were 2.8, 2.0, and 1.9 min, respectively. The convergence speed became faster with increasing latitude. Compared with the L and M constellations, the L+M constellation presented the least significant latitude difference in the acceleration of convergence time. Figure 7b shows the average TTFF of the PPP-AR enhanced using the same LEO constellation at different latitudes. Unlike ambiguity-float solutions, the TTFF of ambiguity-fixed solution PPP enhanced using the three LEO constellations showed the same latitude difference as that of the non-LEO constellation-enhanced PPP-AR. This was because the degree of enhancement of fixed-solution PPP enhanced using LEO was weaker than that of the float solution. After the introduction of LEO satellites, it did not obviously follow that the more LEO satellites there were, the shorter the TTFF would be, as shown in Section 4.2.1.

PPP and PPP-AR Enhanced by Different LEO Satellite Numbers
To further study the relationship between PPP and PPP-AR enhanced with the LEO constellation and the number of LEO satellites, we designed four experiments to control the number of LEO satellites. It should be noted that the LEO constellation used in this experiment was the L+M mixed constellation. Two LEO satellites have been taken as an example, to illustrate the method of controlling the number of LEO satellites in this study: when there were more than two LEO satellites, the first two satellites were retained according to the elevation angle. If one of the retained LEO satellites fell (i.e., the elevation angle became lower than the cutoff angle), the satellite with the higher elevation angle was selected for replacement. Figure 9 shows the convergence times of one, two, three, four, and 'all' LEO satellites that enhanced the single-system PPP at different latitudes. One LEO satellite hardly accelerated the convergence, whereas two, three, four, and 'all' LEO satellites could accelerate the convergence time of the single-system PPP from approximately half an hour to 8.7, 4.4, 3.1, and 2.3 min, respectively. The convergence speed increased nonlinearly with the increase in the number of LEO satellites, and gradually reached saturation. Meanwhile, the latitude difference in the convergence time almost disappeared with the control of the number of LEO satellites. This also explained why the latitude difference in convergence time in Figure 7a was caused by the difference in the number of visible LEO satellites.  Figure 10 shows the TTFF of one, two, three, four, and 'all' LEO satellites with enhanced single-system PPP-AR at different latitudes. As with the float solution, one LEO satellite hardly accelerated the first fixation of ambiguity, while two, three, four, and 'all' LEO satellites accelerated the TTFF of the single-system PPP-AR from approximately 15 min to 9.0, 7.5, 6.8, and 6.4 min, respectively. However, unlike the float solution, the TTFF for the fixed-solution PPP did not decrease significantly with an increase in the number of LEO satellites, and was close to saturation when there were four LEO satellites. This also indicated that the latitude difference of the TTFF in Figure 7b was not caused by the latitude difference in the number of visible LEO satellites, but by the latitude difference of the first fixed time itself.  (4), and all (-) LEO satellites enhanced single-system PPP-AR in different latitudes.

Conclusions and Discussion
This paper focused on the positioning performance of single-system and four-system combined PPP-AR enhanced with the LEO ambiguity-float solution, to determine if the float solution could realize PPP-AR in a shorter TTFF. We designed two LEO constellations, L and M. L was a typical polar orbit constellation, with a larger number of visible satellites at high latitudes than at low and middle latitudes. M compensated for the shortage of visible satellites in L at low and middle latitudes, but it had no coverage in the polar regions. The ground observation data of the LEO satellites at the MGEX stations were simulated. Because GNSSs were fully operational, the GNSS data were real observation data from the MGEX stations. This study evaluated the positioning performance of LEO ambiguity-float solution-enhanced PPP-AR and compared it with LEO-enhanced PPP, based on the daily observation dataset collected at 258 stations in the global MGEX observation network over three days, from 1 January to 3 January 2022, in addition to LEO simulation data.
In our estimated global GCE FCBs products, the NL FCBs of most GPS and Galileo satellites changed within 0.1 cycles in a day. For BDS-3, the stability was relatively poor. Satellites C34 and C46 had the largest changes, up to 0.6 cycles, while other satellites had changes within 0.4 cycles. The average STD of the NL FCBs for G, C, and E were 0.016, 0.054, and 0.014 cycles, respectively.
By introducing L, M, and L+M constellations, the convergence time of single-system PPP was shortened from 29.1 to 9.7, 5.7, and 2.3 min, and the TTFF was shortened from 14.7 to 10.4, 9.4, and 7.9 min, respectively. The convergence time of the four-system combination PPP was shortened from 11.4 to 6.1, 3.0, and 1.8 min, and the TTFF was shortened from 8.6 to 6.8, 4.8, and 4.5 min. The L+M mixed constellation was the best solution for shortening convergence time of PPP and TTFF, followed by the M, and then the L constellation.
Considering that the enhancement performance of the stations at different locations may vary, the 258 stations were divided into low-, middle-, and high-latitude zones for analysis. For ambiguity-float solutions, the convergence time of the L constellation-enhanced PPP in low, middle, and high latitudes was 14.1, 7.8, and 3.3 min, respectively. The higher the latitude, the faster the convergence speed. The convergence time of the M constellation-enhanced PPP in low, middle, and high latitudes was 5.4, 4.6, and 7.0 min, respectively. The convergence time of the L+M mixed constellations in low, middle, and high latitudes was within 3 min, and the latitude difference was not significant. Unlike ambiguity-float solutions, for the TTFF, the ambiguity-fixed solution PPP enhanced using the three LEO constellations showed the same latitude difference as that of the non-LEO constellation-enhanced PPP-AR.
To further study the relationship between both PPP and PPP-AR enhanced with the LEO constellation and the number of LEO satellites, we controlled the number of LEO satellites in the L+M mixed constellation. For ambiguity-float solutions, the convergence speed increased nonlinearly with an increase in the number of LEO satellites, and gradually reached saturation. Meanwhile, the latitude difference in the convergence time almost disappeared with the control of the number of LEO satellites. For fixed-solution PPP, the TTFF decreased with an increase in the number of LEO satellites (which was not significant) and was close to saturation for four LEO satellites.
We recommend the use of a mixed low-orbit constellation that integrates both inclined orbit satellites and polar orbit satellites, which have good enhancement performance not only in high latitudes but, more importantly, in low and middle latitudes as well. Since the world population is mostly distributed in the middle and low latitudes, the shorter convergence time in this region is more meaningful than that it would be in the high latitudes. In addition, the LEO float solution was employed to enhance PPP-AR. This method was undoubtedly a better choice when the ground tracking network of the LEO system was not perfect enough to solve the LEO FCB to fix the LEO satellite ambiguity, which can reduce the time to first fix of the four-system combined PPP-AR to 5 min on the global scale (thus reducing the contingency caused by selecting typical examples). In addition, it does not require too many LEO satellites; four is enough. Similarly, Hong et al. [36] also used this method to study the rapid precise point positioning enhanced with low Earth orbit satellites, pointing out that (based on the small-scale European stations) the time to first fix using the G single system as well as the G/C/E combined system with LEO satellites added could be reduced by 86.4% and 82.8%, respectively, with respect to systems without additional LEO satellites. Data Availability Statement: The raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.