The emerging potential of low-Earth-orbit (LEO) satellite-based Positioning, Navigation, and Timing services has increased the need for real-time, stable, and accurate orbit determination techniques. Here, we propose a method for estimating sub-meter-level LEO satellite orbits using Global Navigation Satellite System (GNSS) code pseudorange observations. To mitigate ionospheric delay, a dual-frequency ionosphere-free combination was applied, while code-carrier smoothing was employed to reduce code observation noise. A satellite weighting model based on Signal-in-Space Range Error was developed to reflect the orbit and clock error characteristics of different GNSS, and a robust weighting scheme was applied to alleviate the impact of observation outliers. Further, Galileo High Accuracy Service corrections compensated for orbit, clock and code bias errors. The algorithm was validated using the GNSS observation data collected from the Sentinel-6A satellite on 10 August 2023. Each successively applied technique gradually improved orbit determination accuracy, achieving up to a 51% reduction in 3D root mean square error (RMSE). The final RMSE values in the radial, along-track, cross-track, and 3D components were 39.4, 18.8, 23.5, and 49.6 cm, respectively. Temporal analysis showed no distinct periodicity in orbit errors and no significant correlation with satellite visibility or ground track.
Lee et al. (Tue,) studied this question.