The fusion of low earth orbit (LEO) constellation and Global Navigation Satellite Systems (GNSS) can increase the number of visible satellites and optimize spatial geometry, which is expected to improve the performance of precise point positioning (PPP) ambiguity resolution (AR). In addition, the multi-frequency signals of LEO satellites can bring a variety of observation combinations, which is potential to further improve the efficiency of PPP AR. In this contribution, multi-frequency PPP AR was achieved with the augmentation of different LEO constellations. Three types of LEO constellations were designed with 60, 192, and 288 satellites. Moreover, the corresponding observation data were simulated with the GNSS observations over the ground stations. The LEO constellations were designed to transmit navigation signals on three frequencies: L1, L2, and L5 at 1575.42, 1227.6, and 1176.45 MHz, respectively, which are consistent with the GPS signals. For PPP AR, the uncalibrated phase delay (UPD) products of GNSS and LEO were estimated first. Furthermore, the quality of UPD products was also analyzed. The research findings show that the performance of estimated LEO UPD is comparable to that of GNSS UPD. Based on the UPD products, LEO-augmented multi-GNSS PPP AR can be achieved. Numerous results show that the performance of single-system and multi-GNSS PPP AR can be significantly improved by introducing the LEO constellations. The augmentation performance is more remarkable in the case of increasing LEO satellites. The time to first fix (TTFF) of the GREC fixed solution can be shortened from 7.1 to 4.8, 1.1, and 0.7 min, by introducing observations of 60-, 192-, and 288-LEO constellations, respectively. The positioning accuracy of multi-GNSS fixed solutions is also improved by about 60%, 80%, and 90% with the augmentation of 60-, 192-, and 288-LEO constellations, respectively. Compared to the dual-frequency solutions, the triple-frequency LEO-augmented PPP fixed solution presents a better performance. The TTFF of GREC fixed solutions is shortened to 33 s with the augmentation of 288-LEO constellation under the triple-frequency environment. It is worth indicating that the 288-satellite LEO-only PPP AR was conducted in dual-frequency and triple-frequency modes, respectively. The averaged TTFFs of both modes are 71.8 s and 55.2 s, respectively. It indicates that LEO constellation with 288 satellites is capable of achieving high-precision positioning independently and shows an even better performance than GNSS-only solutions.

Global navigation satellite system (GNSS) orbits are traditionally determined by observation data of ground stations, which usually need even global distribution to ensure adequate observation geometry strength. However, good tracking geometry cannot be achieved for all GNSS satellites due to many factors, such as limited ground stations and special stationary characteristics for the geostationary Earth orbit (GEO) satellites in the BeiDou constellation. Fortunately, the onboard observations from low earth orbiters (LEO) can be an important supplement to overcome the weakness in tracking geometry. In this contribution, we perform large LEO constellation-augmented multi-GNSS precise orbit determination (POD) based on simulated GNSS observations. Six LEO constellations with different satellites numbers, orbit types, and altitudes, as well as global and regional ground networks, are designed to assess the influence of different tracking configurations on the integrated POD. Then, onboard and ground-based GNSS observations are simulated, without regard to the observations between LEO satellites and ground stations. The results show that compared with ground-based POD, a remarkable accuracy improvement of over 70% can be observed for all GNSS satellites when the entire LEO constellation is introduced. Particularly, BDS GEO satellites can obtain centimeter-level orbits, with the largest accuracy improvement being 98%. Compared with the 60-LEO and 66-LEO schemes, the 96-LEO scheme yields an improvement in orbit accuracy of about 1 cm for GEO satellites and 1 mm for other satellites because of the increase of LEO satellites, but leading to a steep rise in the computational time. In terms of the orbital types, the sun-synchronous-orbiting constellation can yield a better tracking geometry for GNSS satellites and a stronger augmentation than the polar-orbiting constellation. As for the LEO altitude, there are almost no large-orbit accuracy differences among the 600, 1000, and 1400 km schemes except for BDS GEO satellites. Furthermore, the GNSS orbit is found to have less dependence on ground stations when incorporating a large number of LEO. The orbit accuracy of the integrated POD with 8 global stations is almost comparable to the result of integrated POD with a denser global network of 65 stations. In addition, we also present an analysis concerning the integrated POD with a partial LEO constellation. The result demonstrates that introducing part of a LEO constellation can be an effective way to balance the conflict between the orbit accuracy and computational efficiency.

The rapid growing number of earth observation missions and commercial low-earth-orbit (LEO) constellation plans have provided a strong motivation to get accurate LEO satellite position and velocity information in real time. This paper is devoted to improve the real-time kinematic LEO orbits through fixing the zero-differenced (ZD) ambiguities of onboard Global Navigation Satellite System (GNSS) phase observations. In the proposed method, the real-time uncalibrated phase delays (UPDs) are estimated epoch-by-epoch via a global-distributed network to support the ZD ambiguity resolution (AR) for LEO satellites. By separating the UPDs, the ambiguities of onboard ZD GPS phase measurements recover their integer nature. Then, wide-lane (WL) and narrow-lane (NL) AR are performed epoch-by-epoch and the real-time ambiguity–fixed orbits are thus obtained. To validate the proposed method, a real-time kinematic precise orbit determination (POD), for both Sentinel-3A and Swarm-A satellites, was carried out with ambiguity–fixed and ambiguity–float solutions, respectively. The ambiguity fixing results indicate that, for both Sentinel-3A and Swarm-A, over 90% ZD ambiguities could be properly fixed with the time to first fix (TTFF) around 25–30 min. For the assessment of LEO orbits, the differences with post-processed reduced dynamic orbits and satellite laser ranging (SLR) residuals are investigated. Compared with the ambiguity–float solution, the 3D orbit difference root mean square (RMS) values reduce from 7.15 to 5.23 cm for Sentinel-3A, and from 5.29 to 4.01 cm for Swarm-A with the help of ZD AR. The SLR residuals also show notable improvements for an ambiguity–fixed solution; the standard deviation values of Sentinel-3A and Swarm-A are 4.01 and 2.78 cm, with improvements of over 20% compared with the ambiguity–float solution. In addition, the phase residuals of ambiguity–fixed solution are 0.5–1.0 mm larger than those of the ambiguity–float solution; the possible reason is that the ambiguity fixing separate integer ambiguities from unmodeled errors used to be absorbed in float ambiguities.