A set of Global Navigation Satellite Systems (GNSS) satellite orbit and clock offset are an essential prerequisite for precise application. However, abrupt changes in accuracy at the boundaries are prevalent in products provided by international GNSS services, resulting in decreased orbit interpolation precision near the daily boundary. In addition, the effect of this phenomenon is reflected in the deterioration of accuracy and the fluctuations in subsequent applications. In this study, time-weighted and equal-weighted calibrated methods were utilized for adjacent Global Positioning System (GPS) satellite orbits and the orbit variations were then corrected for the clock offset to ensure their consistency. The calibration method is evaluated based on the accuracy and smoothness of post-processing kinematic precise point positioning (PPP) and low earth orbit (LEO) precise orbit determination (POD) near the day boundary. In a variety of scientific applications, the results indicate that the proposed calibration method can effectively reduce the excessive differences near the day boundary between adjacent days. Near the boundary, maximum improvements for post-processing kinematic PPP, dynamic LEO precision orbit, kinematic LEO precision orbit are 41.5%, 9.4%, and 20.5%, respectively.

Multi Global Navigation Satellite System (GNSS) Precise Point Positioning (PPP) has become the mainstream of PPP technology. Due to the differences in the coordinates and time references of each GNSS, multi-GNSS PPP must include additional Inter-System Bias (ISB) parameters to ensure compatibility between different GNSSs. Therefore, research on the characteristics of ISB is also essential. To analyze the short- and long-term time characteristics of multi-GNSS ISBs, as well as their relationship with receiver type and receiver antenna type, the Undifferenced and Uncombined (UDUC) PPP model of Global Positioning System (GPS), BeiDou navigation satellite system (BDS), and Galileo satellite navigation system (Galileo) was rigorously derived, and the physical of ISBs was elaborated in depth. ISB parameters were estimated and analyzed using 31 days of data from the 31 Multi-GNSS Experimental stations (MGEX). The results indicate that: (1) the ISB value is dependent on the station receiver type, receiver antenna type, analysis center product utilized, and GNSS system. (2) The short-term time characteristics of ISB-COM, ISB-WUM, and ISB-GBM are similar for the same station but not for the long term. In addition, ISBs are more stable in the short term. (3) There is little correlation between the ISB time characteristics, the receiver type, and the receiver antenna type, and the day-boundary discontinuity(DBD) on the ISB can be ignored for the concussive days’ process.