When satellites enter into the noon maneuver or the shadow crossing regimes, the actual attitudes will depart from their nominal values. If improper attitude models are used, the induced-errors due to the wind-up effect and satellite antenna PCO (Phase Center Offset) will deteriorate the positioning accuracy. Because different generations of satellites adopt different attitude control models, the influences on the positioning performances deserve further study. Consequently, the impact of three eclipsing strategies on the single-system and multi-GNSS (Global Navigation Satellite System) Precise Point Positioning (PPP) are analyzed. According to the results of the eclipsing monitor, 65 globally distributed MGEX (Multi-GNSS EXperiment) stations for 31-day period in July 2017 are selected to perform G/R/E/C/GR/GREC PPP in both static and kinematic modes. The results show that the influences of non-nominal attitudes are related to the magnitude of the PCO values, maximum yaw angle differences, the duration of maneuver, the value of the sun angle and the satellite geometric strength. For single-system, using modeled attitudes rather than the nominal ones will greatly improve the positioning accuracy of GLONASS-only and BDS-only PPP while slightly contributions to the GPS-only and GALILEO-only PPP. Deleting the eclipsing satellites may sometimes induce a longer convergence time and a worse solution due to the poor satellite geometry, especially for GLONASS kinematic PPP when stations are located in the low latitude and BDS kinematic PPP. When multi-GNSS data are available, especially four navigation systems, the accuracy improvements of using the modeled attitudes or deleting eclipsing satellites are non-significant.

In this paper, we propose a new reduced-dynamic (RD) method by introducing the second-order time-difference position (STP) as additional pseudo-observations (named the RD_STP method) for the precise orbit determination (POD) of low Earth orbiters (LEOs) from GPS observations. Theoretical and numerical analyses show that the accuracies of integrating the STPs of LEOs at 30 s intervals are better than 0.01 m when the forces (<10⁻⁵ ms⁻²) acting on the LEOs are ignored. Therefore, only using the Earth’s gravity model is good enough for the proposed RD_STP method. All unmodeled dynamic models (e.g., luni-solar gravitation, tide forces) are treated as the error sources of the STP pseudo-observation. In addition, there are no pseudo-stochastic orbit parameters to be estimated in the RD_STP method. Finally, we use the RD_STP method to process 15 days of GPS data from the GOCE mission. The results show that the accuracy of the RD_STP solution is more accurate and smoother than the kinematic solution in nearly polar and equatorial regions, and consistent with the RD solution. The 3D RMS of the differences between the RD_STP and RD solutions is 1.93 cm for 1 s sampling. This indicates that the proposed method has a performance comparable to the RD method, and could be an alternative for the POD of LEOs.

GPS-based, satellite-to-satellite tracking observations have been extensively used to elaborate the long-scale features of the Earth’s gravity field from dedicated satellite gravity missions. We proposed compiling a satellite gravity field model from Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite accelerations directly estimated from the onboard GPS data using the point-wise acceleration approach, known as the carrier phase differentiation method. First, we composed the phase accelerations from the onboard carrier phase observations based on the sixth-order seven-point differentiator, which can eliminate the carrier phase ambiguity for Low Earth Orbiter (LEO). Next, the three-dimensional (3D) accelerations of the GOCE satellite were estimated from the derived phase accelerations as well as GPS satellite ephemeris and precise clock products. Finally, a global gravity field model up to the degree and order (d/o) 130 was compiled from the 71 days and nearly 2.5 years of 3D satellite accelerations. We also recovered three gravity field models up to d/o 130 from the accelerations derived by differentiating the kinematic orbits of European Space Agency (ESA), Graz, and School of Geodesy and Geomatics (SGG), which was the orbit differentiation method. We analyzed the accuracies of the derived accelerations and the recovered gravity field models based on the carrier phase differentiation method and orbit differentiation method in time, frequency, and spatial domain. The results showed that the carrier phase derived acceleration observations had better accuracy than those derived from kinematic orbits. The accuracy of the recovered gravity field model based on the carrier phase differentiation method using 2.5 years observations was higher than that of the orbit differentiation solutions for degrees greater than 70, and worse than Graz-orbit solution for degrees less than 70. The cumulative geoid height errors of carrier phase, ESA-orbit, and Graz-orbit solutions up to degree and order 130 were 17.70cm, 21.43 cm, and 22.11 cm, respectively.