Study of Fast and Reliable Time Transfer Methods Using Low Earth Orbit Enhancement

: The Global Navigation Satellite System (GNSS) can be utilized for long-distance and high-precision time transmission. With the ongoing development of low Earth orbit (LEO) satellites and the rapidly changing geometric relationships between them, the convergence rate of ambiguity parameters in Precise Point Positioning (PPP) algorithms has increased, enabling fast and reliable time transfer. In this paper, GPS is used as an experimental case, the LEO satellite constellation is designed, and simulated LEO observation data are generated. Then, using the GPS observation data provided by IGS, a LEO-enhanced PPP model is established. The LEO-augmented PPP model is employed to facilitate faster and more reliable high-precision time transfer. The application of the LEO-augmented PPP model to time transfer is examined and discussed through experimental examples. These examples show multiple types of time transfer links, and the experimental outcomes are uniform. GPS + LEO is compared with exclusive GPS time transfer schemes. The clock offset of the time transfer link for the GPS + LEO scheme converges more swiftly, meaning that the time required for the clock offset to reach a stable level is the briefest. In this paper, standard deviation is employed to assess stability, and Allan deviation is utilized to assess frequency stability. The re-sults show that the clock offset stability and frequency stability achieved by the GPS + LEO scheme are superior within the convergence time range. Controlled experiments with different numbers of satellites for LEO enhancement indicate that time transfer performance can be improved by increasing the number of satellites. As a result, augmenting GPS tracking data with LEO observations enhances the time transfer service compared to GPS alone.


Introduction
The technology of time transfer has resolved the problem of synchronizing clocks between two locations.Currently, time transfer can be executed through a variety of methods, including fiber optics, broadcast radio signals, telephones, networks, and more.The Global Navigation Satellite System (GNSS) time transfer method, which employs radio signals for time transfer, offers extensive coverage, high accuracy, and ease of operation [1].In recent years, GNSS time transfer based on Precise Point Positioning (PPP) solutions has experienced rapid development [2].The PPP method combines pseudorange and carrier-phase observation data, with precision products, such as precise orbits and clock deviations, to achieve highprecision time transfer services [3].Zhang et al., achieved PPP time transfer with an accuracy of 0.1-0.2ns and frequency stability of 1 × 10 −15 -2 × 10 −15 in half a day using International GNSS Service (IGS) precise satellite clock differential products [4].However, attaining PPP time transfer accuracy at the sub-nanosecond level typically necessitates a substantial amount of convergence time, spanning several tens of minutes [5].In practical applications of GNSS time transfer, the method's real-time reliability is affected by environmental interference on the receiving terminal when receiving satellite signals [6,7].
Numerous scholars have engaged in research aimed at enhancing the performance of PPP time transfer.Zhao et al. [8] proposed a method to enhance PPP time transfer stability by constraining receiver clock offset, resulting in a significant short-term stability improvement.Xu et al. [9] created a multi-frequency ionosphere elimination PPP time transfer model for the BeiDou Satellite Navigation System (BDS).The tri-frequency PPP time transfer achieved an accuracy of 0.61 ns and the dual-frequency accuracy was 0.75 ns, demonstrating the advantages and feasibility of BDS tri-frequency PPP in time transfer applications.Lv et al. [10] performed PPP time transfer by incorporating multi-GNSS observations and ambiguity fixing techniques.The results demonstrated that integrating multi-system observations into PPP time transfer reduced the standard deviation by 38.1% compared to the traditional GPS-only approach.
In addition to the aforementioned methods that enhance PPP time transfer performance through parameter constraints, multiple frequencies, and multi-system integration, the advent of low Earth orbit (LEO) satellites has introduced a novel avenue for performance enhancement.Since the 1990s, numerous organizations have launched constellations of low Earth orbit (LEO) communications satellites, with the LEO Iridium U.S. communications satellite constellation being the most prevalent example.[11].Compared with GNSS satellites, LEO satellites offer an increased signal strength and a faster geometry change due to their lower orbital altitudes.This confers robust anti-spoofing and anti-jamming capabilities, ensuring normal continuous signal reception in signal-obscuring environments [6,12].The ongoing evolution of LEO satellites has led to new advancements in Global Navigation Satellite Systems.Reid et al. [13,14] established the viability of developing integrated positioning, navigation, and timing (PNT) systems with LEO satellites by analyzing signal strength, propagation times, and the geometric dynamics between LEO satellites.This confirms that LEO-augmented GNSS systems are achievable.Ke et al. [15] formulated a PPP model for ionospheric removal using GPS and LEO satellites.They confirmed that GPS combined with LEO satellites allows for a shorter PPP convergence time than just processing sole GPS data.Ge et al. [16] created simulation software to simulate GNSS and LEO observations from ground stations.Their simulation outcomes indicated that LEO-augmented GNSS could reduce the PPP convergence time to 5 min.Su et al. [17] assessed the performance of BDS-3 enhanced by an LEO satellite constellation through data simulation, and the findings demonstrated that the PPP convergence time was decreased from approximately 30 min to 1 min.Zhao et al. [2] appraised the efficacy of the integrated LEO/GNSS PPP under adverse environmental conditions and executed comparative experiments with observational simulations, revealing that the convergence velocity and precision of both fixed and floating solutions derived from the LEO enhancement of GNSS were enhanced.Li et al. [18] developed three types of LEO constellations with 60,192 and 288 satellites and simulated the observation data with GNSS observations over the ground station.Concurrently, the LEO constellation was designed to transmit navigation signals at three frequencies: L1, L2 and L5.The results demonstrated that the introduction of LEO constellations can markedly enhance the performance of both single-system and multi-GNSS PPP AR.Li et al. [19] studied the feasibility of providing real-time PPP services with LEO-augmented GNSS.The results confirmed that LEO-augmented GNSS can accomplish spatial signal ranging with an error at the centimeter level.
Based on the aforementioned studies, the integration GNSS with LEO satellites can enhance the performance of PPP.However, a scarcity of studies have applied LEO enhancement to time transfer.Consequently, this study concentrates on the application of LEO-augmented GNSS to time transfer.This paper employs GPS data as a case study.By designing a constellation of LEO satellites and simulating satellite observations, we conducted an analysis of the performance of GPS time transfer enhanced by LEO, incorporating actual GPS observation data.

LEO Constellation Design and Observation Simulation
To investigate the time transfer performance of the GNSS enhanced by LEO satellites, due to the lack of existing observation data of low orbit satellites, this study simulates the observation data of low orbit satellites.The simulation of observation data is the inverse process of positioning calculation [20].The simulated observation data of low orbit used in this study are obtained by introducing various types of errors into the geometric distance between the satellites and the receivers.Since GNSS observations use GPS, for ease of interoperability, the simulated observation frequency for the LEO satellites is the same as that of GPS, 1575.42MHz (L1) and 1227.60 MHz (L2).Before the simulation, fixed station coordinates, receiver clock bias, and zenith tropospheric delay are obtained through a static PPP solution in preparation for subsequent simulations.
First, the LEO satellite constellation is designed with simulated orbits and clock discrepancies.This study simulates 160 low Earth orbit satellites, with 70 distributed across six polar orbits and 90 across ten inclined orbits.The inclination of the polar orbits is 90 • , while the inclination of the inclined orbits is 60 • [21,22].The low Earth orbit satellite group adopts a Walker constellation, which consists of multiple circular orbit satellites with the same orbital altitude, with the satellite orbital height being 1000 km.For the simulation of the low Earth orbit satellite orbits, the initial positions and initial velocities of each satellite are calculated based on the aforementioned LEO constellation parameters, and then, the orbital dynamics are integrated to obtain full-day orbit data.The orbits are then transformed from the geocentric inertial coordinate system to the Earth-centered, Earth-fixed coordinate system (ECEF), with the orbital data interval designed to be 30 s for interoperability with IGS products.The simulation of orbit errors is achieved by combining cosine trend terms, constant trend terms, and Gaussian noise.The clock bias of the low Earth orbit satellites is simulated using the GPS satellite clock biases provided by IGS.Without considering the system time bias between LEO and GPS, the clock bias of GPS satellites is randomly substituted for the low Earth orbit satellites.The errors published by IGS are calculated clock biases from the IF combination (a linear combination of L1/L2), and the simulated clock errors absorb the pseudorange hardware delays of dual-frequency ionosphere-free combination satellites.It should be noted that the simulation results for orbit and clock errors may not be absolutely accurate, so further consideration should be given to systematic and random errors after combining satellite orbits and satellite clock errors.To simplify the process, noise of different magnitudes is directly added to the orbital components and satellite clock errors.
Geometric distances are calculated using known receiver coordinates, satellite orbits, and satellite clock errors.Then, various error terms are calculated using error models.Specifically, the ionospheric delay is calculated using the projection function and total electron content from the global ionospheric maps published by the Center for Orbit Determination in Europe (CODE).For the simulation of tropospheric delay, the wet component of the troposphere delay is calculated using static PPP estimation and the GMF (Global Mapping Function) projection function, while the dry component of the troposphere delay is calculated using the Saastamonien empirical model and GMF projection function [11,23].The initial measurement noise for code and phase is simulated as Gaussian noise with a standard deviation of 0.5 m and 3 mm, respectively, and is inversely correlated with elevation angle.Earth rotation errors, relativistic effects, and tides are corrected using existing corresponding models.Finally, all types of errors are added to the geometric distance between the satellite and ground station to obtain the final observed data at the ground station.

LEO-Augmented PPP Model
In this study, we utilized both GNSS observation data and simulated LEO satellite observation data.The GNSS observation data encompassed pseudorange observations as well as carrier-phase observations.The PPP observation equation is constructed as follows: where the superscript G indicates the GNSS system type; the superscript S indicates a specific satellite; the subscript f indicates the carrier frequency of the satellite signal; the subscript R indicates the receiver; P G,S f ,R is the pseudorange observation; ρ G,S R is the geometric distance between the station and the satellite; c is the speed of light under vacuum conditions; δT R is the receiver clock offset, which refers to the time deviation between the system time and the receiver clock; δT G,S is the satellite clock offset; τ G,S ion, f ,R is the ionospheric delay at the f carrier frequency; τ trop is the tropospheric delay; τ G delay_P, f ,R is the hardware delay related to the receiver code, τ G,S delay_P, f is the code hardware delay on the satellite side, with the code calibration delay in meters; ε G,S P, f ,R is the pseudo-distance observation noise; Φ G,S f ,R is the carrier-phase observation; λ G,S f is the wavelength of the satellite carrier signal; N G,S f ,R is integer ambiguity; τ G delay_Φ, f ,R is the hardware delay for the receiver phase, τ G,S delay_Φ, f is the hardware delay for the phase on the satellite end, with phase calibration delays in cycles; and ε G,S Φ, f ,R is the carrier-phase observation noise.In this study, to eliminate the first-order ionospheric delay of the pseudorange and carrier-phase observables, the ionosphere-free (IF) PPP model was employed as follows: where f G,S 1 and f G,S 2 are the frequencies of different carrier signals.Therefore, the IF-PPP model observation equation is as follows: The ionosphere-free combination code pseudorange hardware delay at the satellite end and the receiver end is absorbed by the ambiguity parameters.So, Formulas (5) and ( 6) are rewritten as Formulas (7) and (8), respectively.
The reparameterized receiver clock error and ambiguity are shown in Formula (9): where τ G UPD,IF,R is the uncalibrated phase delay at the receiver end; τ G,S UPD,IF is the uncalibrated phase delay at the satellite end; and N G,S IF,R is the IF combined ambiguity parameter, which contains the hardware delay deviation at the satellite and receiver ends.
By combining the simulated observational data from the LEO satellites, the PPP model for LEO enhancement was constructed as follows: Among them, the L superscript represents the LEO satellite.

Data Source
To investigate the impact of the GNSS PPP model with LEO-augmented GNSS time transfer performance, we used the GPS observation data provided by the IGS for day 1 in 2022 and the related precise products.A cohort of 13 stations was curated based on equipment type and geographic location.The station coordinates were fixed during PPP solving, and the station information is shown in Table 1 below.Utilizing the aforementioned stations, ten time transfer links were established.The time transfer link information is shown in Table 2, and the spatial distribution of the ten time transfer links is illustrated in Figure 1.
With regard to the LEO satellite observation data, the ground station of the MGEX observation network was used to simulate the LEO satellite ground observation data based on the above LEO satellite observation simulation method.The daily average number of visible satellites in the LEO constellation on the first day of 2022 is shown in Figure 2. The relevant data used in the study are presented in Table 3. Utilizing the aforementioned stations, ten time transfer links were established.The time transfer link information is shown in Table 2, and the spatial distribution of the ten time transfer links is illustrated in Figure 1.With regard to the LEO satellite observation data, the ground station of the MGEX observation network was used to simulate the LEO satellite ground observation data based on the above LEO satellite observation simulation method.The daily average number of visible satellites in the LEO constellation on the first day of 2022 is shown in Figure 2. The relevant data used in the study are presented in Table 3.

Experimental Solution
To analyze the effects of LEO enhancement on GPS time transfer performance, a sole GPS time transfer experiment was carried out.We used the following two experimental

Experimental Solution
To analyze the effects of LEO enhancement on GPS time transfer performance, a sole GPS time transfer experiment was carried out.We used the following two experimental schemes for the comparative analysis of time transfer performance: Solution A: Time transfer using a sole GPS system.Solution B: Time transfer using GPS combined with LEO satellites.
To verify the applicability of LEO-augmented GPS time transfer at different latitudes and its performance with different numbers of satellites, we designed time transfer links and simulated different numbers of satellites.

Example Analysis 4.1. Analysis of LEO-Augmented GPS Time Transfer Convergence Performance
Three representative stations, namely areg, harb, and ons1, were selected to analyze the receiver clock offset calculated by LEO-augmented PPP.The three stations-areg, harb, and ons1-used rubidium, cesium, and hydrogen, respectively, as their clock frequencies.The first 30 min receiver clock offsets of the three stations are shown in Figure 3a-c, and the first-order difference results of the clock offsets are shown in Figure 4a-c.In the calculation, the three stations involved an average of eight low-orbit satellites.According to the clock offset variation diagram of the three stations, receiver clock offsets converge faster under LEO enhancement than under the sole GPS.To quantify the convergence time, by calculating the first-order difference of the receiver clock offset, the mean and standard deviation were calculated using the converged first-order difference data, and the convergence confidence interval was obtained by adding and subtracting the mean from two times the standard deviation.The convergence criterion was that the first-order difference result of the receiver clock offset should be within the range of the convergence confidence interval for 10 min.The red lines in Figure 4a-c represent the boundaries of the convergence confidence interval.The results of the statistical convergence time calculation are presented in Table 4 and Figure 5.According to the statistics, the convergence period was significantly shorter compared to processing sole GPS data and moreover significantly different at the three stations.The convergence speed increased the most (65.0%) at station harb, with its cesium clock; second-most (50.0%) at station ons1, with its hydrogen clock; and least (18.8%) at station areg, with its rubidium clock.gence time calculation are presented in Table 4 and Figure 5.According to the statistics, the convergence period was significantly shorter compared to processing sole GPS data and moreover significantly different at the three stations.The convergence speed increased the most (65.0%) at station harb, with its cesium clock; second-most (50.0%) at station ons1, with its hydrogen clock; and least (18.8%) at station areg, with its rubidium clock.boundaries of the convergence confidence interval.The results of the statistical convergence time calculation are presented in Table 4 and Figure 5.According to the statistics, the convergence period was significantly shorter compared to processing sole GPS data and moreover significantly different at the three stations.The convergence speed increased the most (65.0%) at station harb, with its cesium clock; second-most (50.0%) at station ons1, with its hydrogen clock; and least (18.8%) at station areg, with its rubidium clock.gence time calculation are presented in Table 4 and Figure 5.According to the statistics, the convergence period was significantly shorter compared to processing sole GPS data and moreover significantly different at the three stations.The convergence speed increased the most (65.0%) at station harb, with its cesium clock; second-most (50.0%) at station ons1, with its hydrogen clock; and least (18.8%) at station areg, with its rubidium clock.In light of the aforementioned conclusions, it can be posited that the receiver clock offset obtained by LEO-augmented PPP can effectively enhance convergence speed.

Receiver
The clock offset of a time transfer link calculated for the established 10 links are shown in Figure 5.The time transfer link clock offset had the same performance as the receiver clock offset.The time transfer link achieved using LEO-augmented PPP exhibited a faster performance compared to that attained by the GPS system (Figure 6).Concurrently, the 10 links had different latitude distributions, baseline lengths, and clock sources.Although the different time transfer links had different characteristics, they exhibited similar performances.A first-order differential calculation was applied to the clock offset sequence of the time transfer link, and the convergence confidence interval was determined.The results are presented in Figure 7.The first-order difference of the clock offset clearly showed faster convergence after LEO enhancement.The results of the convergence time calculations are presented in Table 5 and Figure 8.

Stability Analysis of LEO-Augmented GPS Time Transfer
For the sustained and dependable application of GNSS time transfer, the reliability of the clock offset in the time transfer link is of paramount importance.Considering that LEO-augmented PPP can improve the convergence speed of time transfer links, the standard deviation of the clock offset obtained by the two methods in the convergence interval was examined.As evident from the results presented in Table 6 and Figure 9, the standard deviation of the clock offset for the LEO-augmented time transfer link exhibited a marked improvement within the convergence interval, when compared to that of GPS.Therefore, LEO-augmented PPP provides a higher stability in the initial convergence period of time transfer and can better ensure the stability and reliability of real-time fast convergence time transfer.The convergence speed of the time transfer link constructed by the rubidium clock link and the hydrogen-cesium hybrid clock group was the slowest.The performance of the rubidium clock link was consistent with the convergence speed of the receiver clock offset, which reflected the atomic clock performance.The kiru-dif1 link had the greatest increase in convergence speed among the cesium atomic clock links and was the link with the shortest baseline length among the four cesium atomic clock links.In the hydrogen atomic clock link, convergence time was consistently within four minutes.Therefore, LEO-augmented PPP exhibited different convergence speed improvements in time transfer links with different atomic clock and baseline lengths.There were no significant differences among time transfer links at different latitudes because the LEO constellations had wide latitude coverage.Thus, there was no significant latitude difference in the convergence speed of clock offsets in the LEO-augmented PPP time transfer links.

Stability Analysis of LEO-Augmented GPS Time Transfer
For the sustained and dependable application of GNSS time transfer, the reliability of the clock offset in the time transfer link is of paramount importance.Considering that LEOaugmented PPP can improve the convergence speed of time transfer links, the standard deviation of the clock offset obtained by the two methods in the convergence interval was examined.As evident from the results presented in Table 6 and Figure 9, the standard deviation of the clock offset for the LEO-augmented time transfer link exhibited a marked improvement within the convergence interval, when compared to that of GPS.Therefore, LEO-augmented PPP provides a higher stability in the initial convergence period of time transfer and can better ensure the stability and reliability of real-time fast convergence time transfer.The time and frequency signals of the time transfer link are continuous.Time and frequency signals are not always fixed due to the influence of frequency source performance variation characteristics.However, this alteration has the same impact on the two data processing schemes [1], whether long-term or short-term.The Allan deviation is employed to assess the frequency stability of the time transfer clock offset sequence in the convergence time span.The Allan deviation results for each time transfer link are shown in Figure 10.The blue curve represents the Allan deviation obtained by the sole GPS time transfer scheme, and the orange curve represents the Allan deviation obtained by the GPS + LEO time transfer scheme.The Allan deviation obtained by the GPS + LEO time transfer scheme obviously shows a smaller variance.Thus, adding low-orbit satellites to the sole GPS data improves the frequency stability of the time transfer clock error sequence within the convergence time range.The Allan deviation results for the 10 links were calculated and averaged.As depicted in Table 7, the frequency stability was improved after LEO enhancement over different sampling time intervals, with a minimum improvement of 25.7%.Consequently, within the time transfer convergence interval, LEO-enhanced PPP The time and frequency signals of the time transfer link are continuous.Time and frequency signals are not always fixed due to the influence of frequency source performance variation characteristics.However, this alteration has the same impact on the two data processing schemes [1], whether long-term or short-term.The Allan deviation is employed to assess the frequency stability of the time transfer clock offset sequence in the convergence time span.The Allan deviation results for each time transfer link are shown in Figure 10.The blue curve represents the Allan deviation obtained by the sole GPS time transfer scheme, and the orange curve represents the Allan deviation obtained by the GPS + LEO time transfer scheme.The Allan deviation obtained by the GPS + LEO time transfer scheme obviously shows a smaller variance.Thus, adding low-orbit satellites to the sole GPS data improves the frequency stability of the time transfer clock error sequence within the convergence time range.The Allan deviation results for the 10 links were calculated and averaged.As depicted in Table 7, the frequency stability was improved after LEO enhancement over different sampling time intervals, with a minimum improvement of 25.7%.Consequently, within the time transfer convergence interval, LEO-enhanced PPP time transfer can effectively improve the frequency stability of the time transfer link's clock error sequence.The principal rationale is the fast geometry change in the low-orbit satellites, which allows for a faster and more stable estimation of the receivers' position and clock offset.
ADEV Tau(s) Tau(s)   To further investigate the influence of the enhancement level of LEO satellites on time transfer performance, experiments were designed with different numbers of LEO satellites.Four scenarios were considered for the LEO enhancement level, corresponding to 1, 2, 3, and 4 LEO satellites.Regarding the LEO satellite number control method involved in this paper, the selection strategy of three LEO satellites is taken as an example to illustrate the following: When the number of LEO satellites is more than 3, the first three satellites are retained according to the elevation angle.If a retained LEO satellite has an elevation angle lower than the cut-off elevation angle, the satellite with the higher elevation angle is selected for replacement.
Simultaneously, three time transfer links, namely pie1-ons1, harb-sydn, and ohi2-syog, were selected for conducting experiments on the LEO enhancement level.The diagrams depicting clock offsets for the pie1-ons1 link across the four scenarios are detailed in Figure 11, with the first-order differences shown in Figure 12.Similarly, the clock offset diagrams for the harb-sydn link under the four scenarios are shown in Figure 13, with the corresponding first-order difference results displayed in Figure 14.Lastly, the clock offset diagrams for the ohi2-syog link in the four scenarios are depicted in Figure 15, and the first-order difference results are exhibited in Figure 16.

Analysis of the Influence of LEO Enhancement Degree on Time Transfer Performa
To further investigate the influence of the enhancement level of LEO satelli transfer performance, experiments were designed with different numbers of lites.Four scenarios were considered for the LEO enhancement level, correspo 2, 3, and 4 LEO satellites.Regarding the LEO satellite number control method i this paper, the selection strategy of three LEO satellites is taken as an example t the following: When the number of LEO satellites is more than 3, the first thre are retained according to the elevation angle.If a retained LEO satellite has an angle lower than the cut-off elevation angle, the satellite with the higher eleva is selected for replacement.
Simultaneously, three time transfer links, namely pie1-ons1, harb-sydn, syog, were selected for conducting experiments on the LEO enhancement leve grams depicting clock offsets for the pie1-ons1 link across the four scenarios a in Figure 11, with the first-order differences shown in Figure 12.Similarly, the c diagrams for the harb-sydn link under the four scenarios are shown in Figure 1 corresponding first-order difference results displayed in Figure 14.Lastly, the c diagrams for the ohi2-syog link in the four scenarios are depicted in Figure 1 first-order difference results are exhibited in Figure 16.As the LEO satellite count increases, the effect of LEO enhancement was further strengthened, resulting in an enhanced convergence speed for the time transfer link.To examine the convergence speed, the convergence times of the three links in the four scenarios were tabulated, as shown in Table 8 and Figure 17.Under the four different LEO enhancement conditions, the trend in convergence of the three time transfer links was consistent, and convergence speed increased with satellite number.Combined with the data in Table 5, the convergence speed of the time transfer link clock offset under different LEO enhancement conditions was calculated and compared with that of the GPS.The results are summarized in Table 9.The largest improvement in convergence speed was observed in the harb-sydn link, which uses cesium as its clock frequency source.A sole LEO satellite increased convergence speed by 13.6%, whereas four satellites improved it by up to 77.3%.The ohi2-syog link with the mixed clock group showed little change in  As the LEO satellite count increases, the effect of LEO enhancement was further strengthened, resulting in an enhanced convergence speed for the time transfer link.To examine the convergence speed, the convergence times of the three links in the four scenarios were tabulated, as shown in Table 8 and Figure 17.Under the four different LEO enhancement conditions, the trend in convergence of the three time transfer links was consistent, and convergence speed increased with satellite number.Combined with the data in Table 5, the convergence speed of the time transfer link clock offset under different LEO enhancement conditions was calculated and compared with that of the GPS.The results are summarized in Table 9.The largest improvement in convergence speed was observed in the harb-sydn link, which uses cesium as its clock frequency source.A sole LEO satellite increased convergence speed by 13.6%, whereas four satellites improved it by up to 77.3%.The ohi2-syog link with the mixed clock group showed little change in convergence speed according to the number of satellite.With one satellite, the convergence speed was almost the same as that of GPS.With three and four, the convergence speed increased by 33.3%.The aforementioned results further corroborate that LEO-augmented PPP has different degrees of applicability to time transfer links depending on the clock frequency source.fer link in the convergence interval is analyzed under the condition of different numbers of LEO satellites, and the results are shown in Table 10.The findings demonstrate that the stability of the time transfer link in the convergence interval increased with the number of satellites.Accordingly, with the increase in the number of LEO satellites, the combination of GPS + LEO has a better geometric shape, which is more conducive to improving the performance of PPP time transfer.

Conclusions
This study introduces a novel method for enhancing GNSS time transfer through the use of LEO satellites.The LEO satellite constellation was designed following a Walker constellation.The observational data from the LEO satellites were simulated based on data from all stations in the MGEX observation network.The GNSS observation data incorporated GPS data and precision products offered by the IGS.Utilizing these data, an enhanced PPP with LEO was developed.The receiver clock offset and time transfer link clock offset were additionally calculated.This study investigated the convergence performance and reliability of LEO-augmented GPS time transfer through experimental examples.Four scenarios were designed to assess the impact of LEO enhancement on time transfer reliability.Based on the comprehensive analysis and discussion presented above, the following conclusions were formulated:  Simultaneously, the standard deviation of the clock offset sequence of the time transfer link in the convergence interval is analyzed under the condition of different numbers of LEO satellites, and the results are shown in Table 10.The findings demonstrate that the stability of the time transfer link in the convergence interval increased with the number of satellites.Accordingly, with the increase in the number of LEO satellites, the combination of GPS + LEO has a better geometric shape, which is more conducive to improving the performance of PPP time transfer.

Conclusions
This study introduces a novel method for enhancing GNSS time transfer through the use of LEO satellites.The LEO satellite constellation was designed following a Walker constellation.The observational data from the LEO satellites were simulated based on data from all stations in the MGEX observation network.The GNSS observation data incorporated GPS data and precision products offered by the IGS.Utilizing these data, an enhanced PPP with LEO was developed.The receiver clock offset and time transfer link clock offset were additionally calculated.This study investigated the convergence performance and reliability of LEO-augmented GPS time transfer through experimental examples.Four scenarios were designed to assess the impact of LEO enhancement on time transfer reliability.Based on the comprehensive analysis and discussion presented above, the following conclusions were formulated: (1) The receiver clock offset calculated using the LEO-augmented PPP model demonstrated accelerated convergence.The receiver clock offsets for three representative stations-areg, harb, and ons1-were examined.The convergence velocity of the receiver clock offsets was contrasted with those derived solely from GPS data.The convergence speed post-LEO augmentation was enhanced, with varying degrees of enhancement observed across different links.A convergence confidence interval was established to determine convergence time.The results showed that harb's cesium clock had the greatest increase in convergence speed at 65.0%.ons1's hydrogen clock had an increase of 50.0%, while areg's rubidium clock had the slowest increase at The findings demonstrate that the frequency stability obtained by the PPP model based on LEO satellite fusion in the convergence interval was significantly better than that of the sole GPS.Consequently, the PPP model enriched with LEO satellite data positively influences the reliability of time transfer.(4) Experiments involving varying quantities of LEO satellites demonstrated that an increase in LEO satellite data correspondingly improved the enhancement degree.This was mainly manifested as faster convergence speeds and higher stability in the convergence interval.
In summary, the time transfer method proposed in this study can effectively enhance the performance of time transfer.This will be beneficial for further advancing the development of time-frequency application services and industries, especially for mobile time synchronization terminals; it can improve their application efficiency and product performance.

Figure 1 .
Figure 1.Global distribution of stations and time transfer links.

Figure 1 .Figure 2 .
Figure 1.Global distribution of stations and time transfer links.

Figure 2 .
Figure 2. Average number of daily visible satellites of the LEO constellation on the first day of 2022.

Figure 3 .
Figure 3. Receiver clock offset time series of stations (black denotes GPS alone; blue, GPS + LEOs.

FirstFigure 4 .
Figure 4.The first-order difference of the receiver clock offset time series of the stations (black denotes GPS alone; blue, GPS + LEOs).(a) The first-order difference results of the clock offset sequence for the areg station, (b) The first-order difference results of the clock offset sequence for the harb station, (c) The first-order difference results of the clock offset sequence for the ons1 station.

Figure 6 .
Figure 6.Time transfer link clock offset time series (black denotes GPS alone; blue GPS + LEOs).(a) The clock offset sequence of the areg-mcil time transfer link, (b) The clock offset sequence of the gold-pie1 time transfer link, (c) The clock offset sequence of the pie1-ons1 time transfer link, (d) The clock offset sequence of the ons1-bor1 time transfer link, (e) The clock offset sequence of the bor1-usud time transfer link, (f) The clock offset sequence of the kiru-dlf1 time transfer link, (g) The clock offset sequence of the dlf1-harb time transfer link, (h) The clock offset sequence of the harb-sydn time transfer link, (i) The clock offset sequence of the kiru-syog time transfer link, (j) The clock offset sequence of the ohi2-syog time transfer link.

Figure 6 .
Figure 6.Time transfer link clock offset time series (black denotes GPS alone; blue GPS + LEOs).(a) The clock offset sequence of the areg-mcil time transfer link, (b) The clock offset sequence of the gold-pie1 time transfer link, (c) The clock offset sequence of the pie1-ons1 time transfer link, (d) The clock offset sequence of the ons1-bor1 time transfer link, (e) The clock offset sequence of the bor1-usud time transfer link, (f) The clock offset sequence of the kiru-dlf1 time transfer link, (g) The clock offset sequence of the dlf1-harb time transfer link, (h) The clock offset sequence of the harb-sydn time transfer link, (i) The clock offset sequence of the kiru-syog time transfer link, (j) The clock offset sequence of the ohi2-syog time transfer link.

Figure 8 .
Figure 8.Time transfer link convergence time histogram.

Figure 8 .
Figure 8.Time transfer link convergence time histogram.

Figure 9 .
Figure 9. Histogram of standard deviation of time transfer link clock offsets in the convergence interval.

Figure 9 .
Figure 9. Histogram of standard deviation of time transfer link clock offsets in the convergence interval.

Figure 10 .
Figure 10.The Allan deviation results of the time transfer link in the convergence interval.

Figure 10 .
Figure 10.The Allan deviation results of the time transfer link in the convergence interval. 044

Figure 12 .
Figure 12.First-order difference diagram of pie1-ons1 time transfer link clock offset in the four enhancement scenarios.

Figure 13 .
Figure 13.Clock offset diagram of harb-sydn time transfer link in the four enhancement scenarios.

Figure 13 .Figure 13 .FirstFigure 14 .
Figure 13.Clock offset diagram of harb-sydn time transfer link in the four enhancement scenarios.

Figure 14 .
Figure 14.First-order difference diagram of harb-sydn time transfer link clock offset in the four enhancement scenarios.

Figure 17 .
Figure 17.Convergence time histograms of the three links across the four enhancement scenarios.

Figure 17 .
Figure 17.Convergence time histograms of the three links across the four enhancement scenarios.

18 .
8%.Therefore, LEO-augmented PPP affected clock offset convergence speeds differently on stations with different atomic clock frequency sources, reflecting the performance of the atomic clocks.The performance of atomic clocks directly determines the accuracy of time measurement, which in turn affects the measurement results of satellite signals and is indirectly reflected in the later data processing.(2) The analysis focused on the convergence performance of the clock offsets of the ten time transfer links.The convergence speed of all the time transfer link clock offsets exhibited uniform performance.The convergence speed achieved with the LEOaugmented PPP was superior to that of the GPS-only model.PPP model based on LEO satellites fusion has different improvements in the convergence speed of time transfer links of different atomic clock types, and the improvement in the convergence speed of time transfer links of different baseline lengths is also different.There were no obvious differences in time transfer links at different latitudes.The enhanced convergence speed has been demonstrated to be more advantageous for securing the speed and real-time performance of GNSS time transfer.(3) The standard deviation of the clock offsets for the time transfer link, as derived from both satellite constellation schemes during the convergence time range, was calculated.The PPP model based on LEO satellites fusion did significantly improve stability in the convergence interval.Additionally, the Allan deviation was employed to assess the frequency stability of the time transfer link clock offset in the convergence interval.

Table 1 .
List of selected measurement station information.

Table 2 .
Time transfer link information.

Table 2 .
Time transfer link information.

Table 4 .
Receiver clock offset convergence time improvement (min).

Table 5 .
Convergence time and percentage of improvement in time transfer links (min).
sequence for the pie1-ons1 time transfer link, (d) The first-order difference results of the clock offset sequence for the ons1-bor1 time transfer link, (e) The first-order difference results of the clock offset sequence for the bor1-usud time transfer link, (f) The first-order difference results of the clock offset sequence for the kiru-dlf1 time transfer link, (g) The first-order difference results of the clock offset sequence for the dlf1-harb time transfer link, (h) The first-order difference results of the clock offset sequence for the harb-sydn time transfer link, (i)The first-order difference results of the clock offset sequence for the kiru-syog time transfer link, (j) The first-order difference results of the clock offset sequence for the ohi2-syog time transfer link.

Table 6 .
Standard deviation of time transfer link clock offset in the convergence interval (ns).

Table 6 .
Standard deviation of time transfer link clock offset in the convergence interval (ns).

Table 7 .
Mean Allan deviation results.

Table 8 .
Convergence time statistics of the three links across the four enhancement scenarios (min).

Table 9 .
Convergence time percentage of improvement in the three links across the four enhancement scenarios.

Table 10 .
Standard deviations of convergence interval clock offsets of the three links across the four enhancement scenarios (ns).

Table 9 .
Convergence time percentage of improvement in the three links across the four enhancement scenarios.

Table 10 .
Standard deviations of convergence interval clock offsets of the three links across the four enhancement scenarios (ns).