Next Article in Journal
Practical Limitations of Using the Tilt Compensation Function of the GNSS/IMU Receiver
Previous Article in Journal
Compound Jamming Recognition Based on a Dual-Channel Neural Network and Feature Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time LEO Satellite Clocks Based on Near-Real-Time Clock Determination with Ultra-Short-Term Prediction

1
National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China
2
University of Chinese Academy of Sciences, Beijing 100094, China
3
Key Laboratory of Time Reference and Applications, Chinese Academy of Sciences, Xi’an 710600, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(8), 1326; https://doi.org/10.3390/rs16081326
Submission received: 5 March 2024 / Revised: 8 April 2024 / Accepted: 8 April 2024 / Published: 10 April 2024
(This article belongs to the Section Satellite Missions for Earth and Planetary Exploration)

Abstract

:
The utilization of Low Earth Orbit (LEO) satellites is anticipated to augment various aspects of traditional GNSS-based Positioning, Navigation, and Timing (PNT) services. While the LEO satellite orbital products can nowadays be produced with rather high accuracy in real-time of a few centimeters, the precision of the LEO satellite clock products that can be achieved in real-time is less studied. The latter, however, plays an essential role in the LEO-augmented positioning and timing performances. In real-time, the users eventually use the predicted LEO satellite clocks, with their precision determined by both the near-real-time clock precision and the prediction time needed to match the time window for real-time applications, i.e., the precision loss during the prediction phase. In this study, a real-time LEO satellite clock determination method, consisting of near-real-time clock determination with ultra-short-term clock prediction is proposed and implemented. The principles and strategies of this method are discussed in detail. The proposed method utilized Kalman-filter-based processing, but supports restarts at pre-defined times, thus hampering continuous bias propagation and accumulation from ancient epochs. Based on the method, using Sentinel-3B GNSS observations and the real-time GNSS products from the National Center for Space Studies (CNES) in France, the near-real-time LEO satellite clocks can reach a precision of 0.2 to 0.3 ns, and the precision loss during the prediction phase is within 0.07 ns for a prediction time window from 30 to 90 s. This results in a total error budget in the real-time LEO satellite clocks of about 0.3 ns.

1. Introduction

With the Global Positioning System (GPS), the Global Navigation Satellite System (GLONASS) in Russia, the BeiDou Navigation Satellite System (BDS), and the Galileo now providing global services, the Global Navigation Satellite System (GNSS) development has entered a new era. However, due to the inherent characteristics of GNSS, it still faces the challenge of being unable to address all the service requirements using a specific navigation technology. As a result, the navigation systems of countries are targeting higher service accuracy, more diverse functionalities, and more reliable services, initiating the planning and construction of the next generation of navigation systems [1].
In recent years, due to the unique advantages of orbital characteristics and signal strength, augmentation with Low Earth Orbit (LEO) satellites has gained attention and favor in the global satellite navigation field, promising to be an incremental development in the next generation of satellite navigation systems. Worldwide, there is the active implementation and deployment of LEO satellite constellation plans, such as Musk’s Starlink [2], the OneWeb [3], the Iridium [4], the Xona’s Pulsar of the United States [5], the Kepler system of Germany [6], and the Centispace [7], Hongyan and Hongyun [8] of China, involving tens of thousands of LEO satellites. Some of these constellations have the capability to actively broadcast navigation signals, such as Xona’s Pulsar and the Centispace [7], aiming to augment the performance of GNSS, providing users with high-precision, fast convergence, low-cost, and highly reliable satellite Positioning, Navigation and Timing (PNT) services [9]. The LEO-augmented positioning results using real LEO satellite navigation signals have, e.g., been achieved with the CentiSpace satellites [7].
For cm-level high-accuracy LEO-augmented PNT services, high-precision real-time LEO satellite clock products are indispensable. Therefore, in addition to the LEO satellite real-time orbital products that have been intensively studied over the past decade [10,11,12,13], high-precision LEO satellite clock products are a prerequisite for integrating the LEO navigation augmentation system with GNSS to provide high-performance PNT services. Their contributions to the LEO satellite Signal-In-Space-Ranging-Error (SISRE) to ground users are unignorable.
In contrast to GNSS satellites that are predominantly positioned in Medium Earth Orbits (MEOs), LEO satellites have an altitude ranging from a few hundred kilometers to approximately 1500 km. This results in significantly smaller footprints compared to that of GNSS satellites [14]. When employing the same clock determination strategy as that for the GNSS satellite clocks using a network of ground stations, frequent reinitializations and discontinuous estimation cannot be avoided due to the difficulties in building dense network stations in remote areas and oceans [11,14]. Luckily, LEO satellites serve as users of GNSS satellite signals and can perform continuous clock determination, treating the clock bias as a time-independent white noise estimation parameter [15]. When processing LEO satellite clocks using batch least-squares adjustment in reduced-dynamic mode, a high precision of about 0.15 ns can be achieved in near-real-time when using high-quality GNSS real-time products [16]. This, however, implies a processing time of about 5 min for 24 h processing with a sampling interval of 30 s, even with a powerful processing unit [12]. Considering a time interval between subsequent processing rounds of, e.g., another 5 min, a clock prediction of 10 min is needed to catch up with the real-time applications.
The mid-term LEO satellite clock prediction turns out to be a more challenging task for LEO satellite clocks compared to that for the GNSS satellite clocks. The estimable LEO satellite clocks exhibit diverse systematic effects, which contain complicated relativistic effects, temperature-related hardware biases, influences caused by the South Atlantic Anomaly, and some other reasons that are not yet fully clarified [17,18]. Wang et al. utilized real data from GRACE-FO and Sentinel-3B to identify the factors impacting the stability of LEO satellite clocks and introduced a prediction model that takes into account the systematic effects [18]. A substantial reduction was demonstrated in the prediction errors compared to polynomial fitting, however, with a precision loss of approximately 0.4 to 0.5 ns at a 10 min prediction time. Wu et al. attempted to remove the complicated relativistic effects on LEO satellite clocks for a better clock prediction, leading to a precision loss of about 0.3 to 0.5 ns at a prediction time point of 10 min [19]. Ge et al. proposed to perform the LEO satellite clock prediction with least-squares harmonic estimation and obtained a precision loss of about 0.3 ns within prediction arcs of 15 min [20]. In short, LEO satellite clock prediction over more than 10 min already leads to significant precision loss. This urges the employment of more efficient LEO satellite clock determination methods with ultra-short-term prediction for real-time applications.
Processing the LEO satellite clocks in a filter-based manner delivers more efficient but less precise clock solutions. Lots of previous studies have discussed the filter-based LEO satellite Precise Orbit Determination (POD), with the clocks automatically produced as byproducts [11,21]. The focus, however, was mostly put on the POD accuracy due to the requirements of other applications instead of navigation. Using simulation data, a real-time estimation approach for the LEO satellite clock based on ground tracking stations was proposed for the first time [22]. Yang et al. introduced a method to estimate the LEO satellite clocks by integrating ground observations from specific regions with onboard observations from LEO satellites. Without using a ground network, the kinematic clock estimation achieved a precision of over 0.3 ns using GNSS products, while incorporating a simulated ground network improved the clock precision to 0.15 ns [23].
In this contribution, using solely the GNSS observations tracked onboard LEO satellite, a real-time clock estimation strategy is introduced containing filter-based near-real-time clock estimation and ultra-short-term clock prediction. The filter-based processing is proposed to be restarted at pre-defined times to hamper bias propagation and accumulation from ancient epochs, allowing also for convergence through processing arcs that are long enough. Using real data from Sentinel-3B, the effects of using different processing arc lengths and time intervals between adjacent processing sessions are discussed, and the total error budget consisting of the near-real-time clock estimation, ultra-short-term clock prediction and the total real-time clock precision achieved under different cases are analyzed.
The study begins with an overview of the LEO satellite real-time clock determination method, containing the near-real-time clock determination and the ultra-short-term clock prediction. The near-real-time clock determination is elaborately discussed, particularly emphasizing the selection and utilization of the GNSS real-time products. Afterward, the processing strategy for real-time clock determination is introduced with near-real-time clock precision and efficiency achieved with processing arc lengths discussed. Subsequently, the total error budget of the real-time LEO satellite clock precision is assessed for the near-real-time estimation part, the ultra-short-term prediction part, and then combined. The conclusions and discussions are given at the end.

2. Method

In this study, the real-time LEO satellite clocks are proposed to be determined and provided to users with the following steps:
(1)
Set a timer and loop through the steps (2)–(4);
(2)
Use real-time GNSS orbits, clocks and LEO satellite on-board GNSS observation to determine the near-real-time LEO satellite clocks based on the kinematic model, including GNSS product selection and the clock re-referencing;
(3)
Predict the LEO satellite clocks in the ultra-short term, ensuring that users obtain LEO satellite clock products in real time;
(4)
Broadcast the low-order polynomial parameters to users at regular intervals.
The flowchart of real-time clock determination is illustrated in Figure 1. The principles of the steps are described briefly in the following subsections.

2.1. LEO Satellite Clock Determination in Near-Real-Time

LEO satellite clock determination can be realized through a reduced-dynamic POD model or a kinematic POD model. The reduced-dynamic model provides higher accuracy in both the POD and the clock determination [12], while the kinematic model combined with, e.g., Kalman-filtering, achieves higher computational efficiency and better flexibility during, e.g., orbital maneuvers. The foundation for both models is the GNSS raw observation equations, with the differences of combining and not combining dynamic models. Therefore, this sub-section starts with GNSS raw observation equations, which can be formulated as:
P i s = ρ i s + c d t r d t s + I i s + d i s d r , i + ε i , p
L i s = ρ i s + c d t r d t s + λ i ( N i s + b r , i b i s ) I i s + ε i , L
where P i s and L i s represent the pseudorange and phase observations, respectively. i indicates the frequency, and s indicates the GNSS satellite. c is the speed of light. d t r and d t s are the LEO satellite clock bias and the GNSS satellite clock bias, respectively. I i s is the slant ionospheric delay. d i s and d r , i represent the pseudorange hardware delays of the GNSS satellite and the LEO satellite, respectively. λ i is the wavelength of the frequency i , and N i s is the carrier phase ambiguity. b r , i and b i s denote the carrier phase hardware delays of the LEO satellite and the GNSS satellite, respectively. ε i , p and ε i , L are the measurement noise and multipath effects of the pseudorange and carrier phase observations, respectively. Other correctable error terms such as tidal effects, phase center offsets and variations, etc., are not specifically described in Equations (1) and (2).
The kinematic model determines the orbit and clocks of LEO satellites using only the GNSS observational data and related real-time GNSS products, without introducing any dynamic models. The kinematic model has the advantages of simplicity, low computational requirements, short processing time, and the ability to maintain good precision. Determining the LEO satellite clocks and orbit using the kinematic model is similar to the PPP (Precise Point Positioning) principle. However, compared to PPP, the GNSS observations received by LEO satellites are not affected by the tropospheric delays.
In this study, the kinematic model is implemented using the Ionosphere-Free (IF) pseudorange and carrier phase observations. The Observed-Minus-Computed (O-C) terms can be expressed as:
P I F = A x Δ r K + c × d t ̄ r + ε P I F
L I F = A x Δ r K + c × d t ̄ r + λ N ¯ I F + ε L I F
with
d t ̄ r = d t r + d r , I F c + d r e l
Δ r K = Δ r x , Δ r y , Δ r z
λ I F = c f 1 + f 2
N ̄ I F = f 1 2 f 1 2 f 2 2 N 1 f 2 2 f 1 2 f 2 2 N 2 + b r , I F b I F s + d r , I F λ I F
where the pseudorange and phase observations O-C terms for the IF combinations are represented by P I F and L I F , respectively. f 1 and f 2 correspond to the first and the second frequencies used, respectively. A x is the design matrix containing partial derivatives of the observations with respect to the orbital coordinates. Δ r K is the orbital correction vector based on the initial orbits, with Δ r x ,   Δ r y and Δ r z representing the increments of the LEO satellite position in the X , Y and Z directions of the Earth-Centered Earth-Fixed (ECEF) system, respectively. d t ̄ r represents the estimable receiver clocks with the inclusion of the IF receiver code bias d r , I F and the relativistic effects d r e l . Recall that d t r denotes the LEO satellite clock bias. It is assumed here that the IF code hardware delay of the GNSS satellite is already included in the GNSS satellite clock products. Depending on the specific type of pseudorange dual-frequency observations used, Differential Code Bias (DCB) corrections may be applied. For multi-GNSS processing, Inter-System Biases (ISBs) need to be additionally considered [13,17].
The Kalman filter is a commonly used parameter estimation method in the field of satellite navigation. It first requires the establishment of a linearized kinematic observation model and a state-transition model for filtering. Subsequently, by recursion, the parameters to be estimated are determined. For detailed information, see [24]. In this study, ambiguity parameters are constrained as constants if there is no cycle slip.
This study aims to determine LEO satellite clocks timely, which involves firstly the near-real-time clock determination as introduced above. Therefore, the choice of real-time GNSS products needs to be carefully considered. Currently, many analysis centers of the International GNSS Service (IGS) and other analysis systems like the International GNSS Monitoring and Assessment System (iGMAS) provide rich real-time GNSS orbit and clocks products [25,26,27,28,29], such as those offered by the National Centre for Space Studies (CNES) in France [30], The Federal Agency for Cartography and Geodes (BKG) in Germany [31], the GMV Aerospace and Defense (GMV) in Spain [32], the Wuhan University in China [33] and the GeoForschungsZentrum (GFZ) in Germany [34]. In this study, the selection of real-time GNSS products is based on the following steps:
(1)
Analyze the availability and continuity of each real-time product over a past period.
(2)
Analyze the accuracy of GNSS satellite orbits, the precision of the GNSS satellite clocks, and their Signal-In-Space Range Error (SISRE) over a past period through comparison with the final GNSS products of, e.g., the Center for Orbit Determination in Europe (CODE) [35,36].
(3)
Based on the results of the first two steps, identify a set of optional real-time GNSS products and establish an initial ranking. The ranking list is updated on a daily basis.
(4)
Upon each processing round, conduct an assessment of the availability/continuity of the real-time products of GNSS satellites in the corresponding processing period. For those passing the pre-defined thresholds for availability and continuity, select the real-time GNSS products according to the ranking list, so that real-time GNSS products with both good completeness and precision can be used for the real-time determination of LEO satellite clocks.
Figure 2 shows, e.g., the precision of real-time GPS satellite clocks from the BKG, CNE, GMV, GFZ, and WHU on Day of Year (DOY) 358 to 364, 2023 with a sampling interval of 10 s. It can be seen that the average clock precision from all the analysis centers is better than 0.2 ns, with BKG and WHU exhibiting the best performances during the test period, approximately at 0.12 ns.
Using the above-mentioned selection strategy for real-time GNSS products, the estimated LEO satellite clocks determined may face two issues:
(1)
Inconsistent time references for LEO satellite clocks are determined in each session. This is caused by the different time references of different real-time GNSS satellite clocks provided by different analysis centers.
(2)
Poor stability of the real-time LEO satellite clock time reference due to the poor stability of the time reference of the real-time GNSS satellite clocks, as shown in Figure 3 (red). In Figure 3, the real-time time reference is calculated as the epoch mean difference between all the usable real-time GPS satellite clocks from the CNES and the CODE final GPS satellite clocks [18].
With inconsistent and unstable time references for LEO satellite clocks determined during different sessions, jumps between sessions and degraded clock prediction results can occur. To solve these issues, the LEO satellite clock estimates need to be re-referenced to a stable time reference directly after the near-real-time clock determination in each session. In real-time, the time reference is proposed to be aligned to a selected satellite clock in the corresponding GNSS with good stability and completeness, avoiding at least jumps (red) in Figure 3. After aligning to a stable GPS satellite clock, as shown in the blue line in Figure 3, the standard deviation (STD) of the time reference is reduced from above 10 ns to 1.5 ns, with all the large jumps removed (blue).
Figure 4 shows the Sentinel-3B satellite clocks determined using the CNES real-time clock products. It can be observed that when using the original real-time CNES clock products, the LEO satellite clocks (red) exhibit significant fluctuations. This is attributed to the unstable reference of real-time CNES clocks as shown in Figure 3. Aligning the CNES real-time clocks to a stable GPS satellite clock and then using it to determine the Sentinel-3B satellite clocks (blue), results in a substantial improvement in stability. The improvement in the stability of the time reference does not affect the clock estimation precision but will influence the LEO satellite clock prediction [18].

2.2. LEO Satellite Clock Prediction in Ultra-Short-Term

As described in [19,37], different LEO satellite clock prediction algorithms do not lead to significant differences in ultra-short-term prediction, i.e., a few minutes or shorter. In this contribution, the LEO satellite clocks are predicted with periodic terms and polynomials, expressed as [18,19]:
C ^ l k t p t 0 = a ^ 0 + a ^ 1 t p t 0 + + a ^ m ( t p t 0 ) m + j = 1 k A ^ j s i n 2 π T ^ j t p t 0 + φ ^ j
where C ^ l k denotes the predicted clocks, t p denotes prediction epoch, t 0 is the initial epoch, a ^ i ( i = 1 , , m ) is the polynomial fitting coefficient, k denotes the number of periodic terms, T ^ j , A ^ j and φ ^ j denotes the period, the amplitude and the phase of the periodic terms, respectively. The determination of the coefficients can be referred to [19].

2.3. LEO Satellite Clock Broadcasting for Real-Time Applications

The GNSS satellite clocks are typically broadcast as low-order polynomial coefficients to users either in navigation messages with a low frequency via satellite links, or in real-time streams with high frequency via Internet links [12]. Like the clocks of GNSS satellites, predicted LEO satellite clocks (see Section 2.2) can be fitted with low-order polynomials and provided to users. Due to the bad LEO satellite clock prediction behaviors over the mid- to long-term [38,39], broadcasting clock coefficients over low-frequency navigation messages like GNSS is not a good option for LEO satellites. It is suggested to either significantly increase the update rate of the broadcast clock messages, or transfer the clocks via real-time streams with the Internet. An update interval of a few seconds can be considered.

3. Processing Strategies

Figure 5 depicts a processing timeline of the real-time LEO satellite clock determination based on near-real-time clock determination and ultra-short-term clock prediction. The entire processing flow is initiated on a scheduled basis. The time interval between the launching of two subsequent processing sessions is represented as T .
In Figure 5, the black line corresponds to the near-real-time clock determination part of LEO satellite clocks, while the red line represents the ultra-short-term clock prediction part. The initial time of the entire process is T 0 , the processing arc length of the GNSS observations is T o b s , and the prediction period is T p r e .   B is the end of T o b s . C represents the time point when the near-real-time clock determination is completed, D represents the time point when the clock prediction is completed, and E represents the end epoch of the predicted clocks. The total processing time for the entire workflow is denoted as T p r o c e s s ( B D ), while T p r o c e s s is mainly the time taken to determine the near-real-time LEO satellite clocks due to the very short time needed for clock prediction, i.e., typically less than 1 s. According to Figure 5, users can use real-time clock products within the interval [ D , E ], and when the products from the next session are available, [ D + T , E + T ].
In general, T p r e ( B E ) must be longer than the sum of T p r o c e s s ( B D ) and T to provide users with real-time LEO satellite clocks; T p r o c e s s is typically shorter than T to avoid overlapped processing and to save computational resources; T p r e needs to be determined according to T and T p r o c e s s . When updating the clock products in each session regularly, as shown in Figure 5, the clock products within the interval [ D , D + T ] are used by real-time users, with the precision representing the real-time LEO satellite clock precision available to users. To achieve a high real-time clock precision, one needs to, on the one side, obtain high-precision near-real-time clocks as the basis for good prediction and on the other side, attempt to suppress T p r o c e s s ( B D ) to shorten the prediction time window, here mainly shorten the BC part for the near-real-time clock determination.
As a representative example, the LEO satellite Sentinel-3B flying at about 810 km [40] is used to show the relationships between the length of the processing arcs and the processing time for the near-real-time clock estimation (the B C part) within T p r o c e s s . Using its 30 s GPS L1/L2 observations, the processing is performed in 288 sessions with the starting time on DOY 228, 2018, shifted by 5 min in each round. The averaged processing time required for the near-real-time clock determination ( B C in Figure 5) is shown in Table 1. For a conservative assessment, a relatively poor processing unit with an operating frequency of the server of 2.1 GHz was used for this test.
From Table 1 it can be observed that reducing the processing arc ( T o b s ) from 24 to 6 h, the processing time is shortened from 24 to 17 s, i.e., with a shortening of almost 30%. Comparing the Kalman-filter-based kinematic solutions employing real-time products from the CNES with the batch least-squares reduced-dynamic clock solutions using the CODE final products (denoted as R), reducing the T o b s from 24 to 4 h leads to almost no changes in the near-real-time clock determination. Further reducing T o b s can slightly reduce the processing time, but could lead to a degradation in the clock precision. Considering the similar processing time and resulting clock precision between T o b s of 4 h and 6 h, in this study, T o b s is set to the longer processing arc of 6 h for safety reasons, with T p r o c e s s (including the determination and the prediction) equal to 18 to 20 s using a processing unit described above. When calculating the STDs, a 4.42-sigma outlier exclusion was performed.
In addition to T o b s and T p r o c e s s , other time points in Figure 5 are also important. T , for example, can be set to 60 to 300 s depending on the real-time precision needed. A short T leads to a short prediction time needed to match the real-time window [ D , D + T ], and delivers a higher precision of the real-time clocks of the LEO satellite.
In summary, the real-time determination process of LEO satellite clocks in this study is as follows:
(1)
At the current time point T , acquire observational data in [ T T o b s , T ], and obtain real-time GNSS products selected according to Section 2.1;
(2)
Based on the strategies outlined in Table 2, determine the LEO satellite clocks in near-real-time;
(3)
According to the prediction method detailed in Section 2.2, predict the LEO satellite clocks over T p r e .
(4)
Perform a second-order polynomial fitting using the predicted clocks, and broadcast the fitted polynomial coefficients to users;

4. Test Results

Real-time LEO satellite clocks determined based on near-real-time clock determination and ultra-short-term clock prediction are tested and analyzed in this section. Firstly, near-real-time LEO satellite clocks are assessed. Based on the clock predictions tested for different periods, the real-time clocks are then analyzed with their total error budget. GPS L1/L2 observations with a sampling interval of 30 s from Sentinel-3B on DOY 227–229, 2018 are used for the processing.

4.1. Near-Real-Time LEO Satellite Clocks

Using the CNES real-time GNSS satellite orbital and clock products, with the processing arc length T o b s set to 6 h, and the time interval T between subsequent processing sessions tested for 60, 120, 180, 240, and 300 s as explained in Section 3, near-real-time Sentinel-3B satellite clocks are determined.
The blue lines in Figure 6 connect the near-real-time Sentinel-3B satellite clocks during the last T period over each processing arc, i.e., the most recently updated near-real-time satellite clocks to be used for prediction. The Kalman-filter-based daily clocks based on the kinematic model utilizing the real-time products provided by the CNES (red line), and the daily clocks based on the batch least-squares reduced-dynamic model using the CODE final products (yellow line) are given for comparison. The LEO satellite clocks are re-referenced to the same time reference, i.e., that of the CODE final clocks, for a fair comparison with:
Δ t p r e f = s = 1 m p ( t p s 1 t p s 2 ) m p
Δ t p s , l e o = t p s , l e o t 0 p s , l e o t p r e f
where Δ t p r e f is the epoch mean difference between the real-time and the CODE final GPS clocks at epoch p , and m p is the number of the compared GPS satellites at epoch p . t p s 1 and t p s 2 are satellite clock biases from the real-time and the CODE final products at epoch p , respectively, and t p s , l e o and t 0 p s , l e o are the LEO satellite clocks determined based on real-time and the CODE final products at epoch p , respectively. Δ t p s represents the re-referenced clock difference at epoch p .
As depicted in Figure 6, the trends of the connected near-real-time Sentinel-3B satellite clocks determined during the last T period of different arcs are generally consistent with the two types of daily solutions. Compared to the daily solutions based on the Kalman-filter-based kinematic model using the CNES products (yellow), the connected near-real-time Sentinel-3B clocks (blue) have fewer outliers. When using the 6 h arcs, the near-real-time clocks determined during the last T period are fully converged. The connected clocks determined from the different sessions do not exhibit obvious jumps at boundaries either, since T is much shorter than the processing arc length of 6 h, and shifting the observation arc by T each time does not lead to significant solution differences at the ending period of the processing arc.
Figure 7 shows the differences between the connected near-real-time Sentinel-3B clocks (blue lines in Figure 6) and the daily clock solutions using the reduced-dynamic model with the CODE final products (yellow lines in Figure 6) over the test period. The results are shifted by N = T 60 1 × 2 ns for each shifted T for a better demonstration. As presented in Figure 7, it is noticeable that the trend of the clock errors is highly consistent for different T . As T increases, the STDs of the clock errors are slightly improved.
Table 3 lists the daily precision of the connected near-real-time clocks for different T . From Table 3, it can first be observed that the connected near-real-time clocks over T based on 6 h Kalman-filter-based processing exhibit very similar precision to the daily kinematic clocks using the same CNES real-time products, i.e., with STD of the differences at 0.01 to 0.03 ns. Secondly, compared with the clocks determined by reduced-dynamic with CODE final products serving as the “true” clocks in this study, the STD of connected near-real-time clocks is about 0.2 to 0.3 ns. Having a look at the last column of Table 3, increasing the T can slightly improve the connected near-real-time clock precision, i.e., with an improvement within 0.02 ns when increasing T from 30 to 300 s. To reduce the required prediction period for real-time applications, a short T of 1 to 2 min should be a good choice.
Figure 8 shows the Modified Allan Deviations (MDEVs) [41] of the connected near-real-time clock errors using different T . As shown in Figure 8, when the averaging time τ is between the range of [100, 300], MDEV decreases with increasing T . This corresponds to the conclusions from Table 3, i.e., a large T is beneficial for the continuity of the clock solutions, especially over the short term, due to fewer connections of the results. For a T of 300 s, the MDEV is slightly increased due to the disturbance of gross errors, which is consistent with its slightly larger STD as shown in Figure 7.

4.2. Real-Time LEO Satellite Clocks

As shown in Figure 5, the interval [ D , D + T ] represents the actual time window when the predicted LEO satellite clocks need to be used by users in real time. In other words, the real-time clocks of interest to users are those within the prediction period [ T p r o c e s s , T p r o c e s s + T ]. This implies that the prediction time T p r e needs to be at least longer than T p r o c e s s + T . As mentioned before, T o b s in Figure 5 is set to 6 h, and T p r o c e s s amounts to about 18–20 s. For possible delays in the GNSS real-time products and transmission of the LEO onboard GNSS observations, T p r o c e s s is set to 30 s in this study. In this sub-section, T is tested for 60, 120, 180, 240 and 300 s, and T p r e is consequently set to 90, 150, 210, 270, and 330 s.
Using the connected near-real-time clocks determined based on different T in Section 4.1, and according to the corresponding T p r e , predictions are performed based on Equation (9) to obtain the real-time Sentinel-3B satellite clocks. The fitting time is tested from 90 to 1800 s with a step of 30 s. The most suitable one, as shown in Table 4, is the one delivering the smallest STD in prediction errors. The clocks determined based on the batch least-squares reduced-dynamic model with the final CODE products are taken as the reference.
The daily precision of the predicted and then connected real-time clocks within [ T p r o c e s s , T p r o c e s s + T ] is listed in Table 4. The precision of the near-real-time clocks is about 0.2 to 0.3 ns. As the T p r e increases from 90 to 330 s, the precision loss due to prediction is from about 0.06 to 0.17 ns, with the total error budget of the real-time clocks amounting to 0.2 to 0.35 ns. As T p r e increases, the precision of real-time clocks gradually decreases due to the slightly increasing prediction errors. As shown in Figure 9, the STD of the near-real time clocks decreases slightly with the increasing T (blue bars), while the precision loss due to prediction increases substantially when increasing the T (red bars). Therefore, it can be inferred that the total error budget of real-time LEO satellite clocks (yellow bars) based on near-real-time determination and ultra-short-term prediction mainly comes from the prediction, and a small T is preferred if possible.
In summary, there are several key factors affecting the final error budget of the real-time LEO satellite clocks:
(1)
The processing time T p r o c e s s is mainly dependent on the processing arc length T o b s . For the very similar T p r o c e s s , it is preferable to have a longer T o b s to enhance the short-term stability of the near-real-time clocks. The time interval between subsequent processing sessions ( T ) is independent of the precision of the near-real-time clocks but directly determines the prediction time T p r e and the usable prediction window [ T p r o c e s s , T p r o c e s s + T ] for real-time applications. The latter is directly related to the predicted clock precision, i.e., the precision of real-time LEO satellite clocks.
(2)
The precision of the connected near-real-time LEO satellite clocks of the last T period within the 6 h processing arcs is similar to the kinematic daily solutions. The major differences in the real-time clocks for different T come from the clock prediction, more concretely, from the different prediction times. Further research is needed to improve the ultra-short-term prediction of LEO satellite clocks, i.e., within 5 min.
(3)
When setting T to 60 s, a real-time clock precision around or lower than 0.3 ns can be achieved.

5. Discussion and Conclusions

The utilization of LEO satellites is anticipated to enhance various aspects of traditional GNSS-based PNT services. These enhancements include improvements in coverage, signal strength, reliability, initialization time, and landing power. The improved performance of LEO satellites positions them as valuable assets for future real-time high-precision PNT users, particularly in applications like autonomous driving and unmanned aerial vehicles. However, in practical applications, LEO satellites also face some challenges. For instance, the determination of real-time LEO satellite clocks still awaits further investigation due to challenges in mid- to long-term LEO satellite clock prediction.
To address this need, this study introduces real-time LEO satellite clock determination based on efficient high-precision near-real-time clock determination and ultra-short-term clock prediction. The method connects the most updated clock estimates from different rounds of Kalman-filter-based solutions and performs ultra-short-term prediction using the clock estimates of the last few minutes, hampering unnecessary bias propagation and accumulations of ancient epochs.
Using Sentinel-3B GPS L1/L2 observations and the CNES real-time GNSS products, near-real-time satellite clocks were determined and assessed. The results showed that the precision of the near-real-time LEO satellite clocks is generally at 0.2 to 0.3 ns. A processing arc of 6 h can be selected considering the efficiency and the precision. The time interval T between subsequent processing sessions affects the prediction time windows needed for real-time applications. When setting the time interval to 60 s, it is possible to obtain real-time LEO satellites with a precision of about 0.3 ns, with a precision loss within 0.07 ns during a prediction time of 90 s. A longer T induces larger prediction errors due to the longer prediction time needed.
Although notable precision is achieved in estimating LEO satellite clocks in this study, there are still challenges in real-world applications. For instance, the timely acquisition of observation data and GNSS real-time products is essential, as the corresponding latencies would directly add to the prediction time required to catch up with the real-time applications, which requires further research in the clock prediction of a longer time. Furthermore, different rounds of the processing could lead to the selection of different clock references, and eventually result in jumps in the time reference between subsequent rounds. This might require re-referencing the clock estimates to a stable time reference before broadcasting to users.

Author Contributions

Formal analysis, M.W., K.W. and X.Y., Methodology, M.W. and K.W., Writing—original draft preparation, M.W., Writing—review and editing, M.W., K.W. and J.W., Software V1.0 M.W., K.W., J.W., J.L., B.C., Z.Z. and W.X., Funding acquisition M.W. and K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation of China (No. 12003041), the National Time Service Center, the Chinese Academy of Sciences (No. E167SC14), the International Partnership Program of the Chinese Academy of Sciences (Grant No. 021GJHZ2023010FN), and the Special Research Assistant Funding Project, Chinese Academy of Sciences (No. 110400T0XW).

Data Availability Statement

The data of Sentinel-3B were obtained from the ESA via https://scihub.copernicus.eu/gnss/#/home, accessed on 1 February 2024. The final products of the Center for Orbit Determination in Europe (CODE) were obtained from ftp://ftp.aiub.unibe.ch/CODE/, 1 February 2024. The CNES real-time products were obtained from http://www.ppp-wizard.net/products/REAL_TIME/, 1 February 2024.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wang, L.; Li, D.; Chen, R.; Fu, W.; Shen, X.; Jiang, H.; Hao, J. Low earth orbiter (LEO) navigation augmentation: Opportunities and challenges. Strateg. Study CAE 2020, 22, 144–152. [Google Scholar]
  2. Starlink. Available online: https://www.starlink.com/ (accessed on 26 February 2024).
  3. Eutelsat OneWeb. Available online: https://oneweb.net/ (accessed on 26 February 2024).
  4. Iridium Satellite Communications. Available online: https://www.iridium.com/ (accessed on 26 February 2024).
  5. Xona Space Systems. Available online: https://www.xonaspace.com/ (accessed on 26 February 2024).
  6. Kepler. Available online: https://www.kepler.global/conf/system/ (accessed on 26 February 2024).
  7. Li, W.; Yang, Q.; Du, X.; Li, M.; Zhao, Q.; Yang, L.; Qin, Y.; Chang, C.; Wang, Y.; Qin, G. LEO augmented precise point positioning using real observations from two CENTISPACE™ experimental satellites. GPS Solut. 2024, 28, 44. [Google Scholar] [CrossRef]
  8. Zhang, X.; Ma, F. Review of the development of LEO navigation-augmented GNSS. Acta Geod. Cartogr. Sin. 2019, 48, 1073. [Google Scholar]
  9. Yi, H.; Lei, W.; Wenju, F.; Haitao, Z.; Tao, L.; Beizhen, X.; Ruizhi, C. LEO navigation augmentation constellation design with the multi-objective optimization approaches. Chin. J. Aeronaut. 2021, 34, 265–278. [Google Scholar]
  10. Hauschild, A.; Tegedor, J.; Montenbruck, O.; Visser, H.; Markgraf, M. Precise onboard orbit determination for LEO satellites with real-time orbit and clock corrections. In Proceedings of the 29th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2016), Portland, OR, USA, 12–16 September 2016; pp. 3715–3723. [Google Scholar]
  11. Li, K.; Zhou, X.; Guo, N.; Zhao, G.; Xu, K.; Lei, W. Comparison of precise orbit determination methods of zero-difference kinematic, dynamic and reduced-dynamic of GRACE-A satellite using SHORDE software. J. Appl. Geod. 2017, 11, 157–165. [Google Scholar] [CrossRef]
  12. Wang, K.; Liu, J.; Su, H.; El-Mowafy, A.; Yang, X. Real-time LEO satellite orbits based on batch least-squares orbit determination with short-term orbit prediction. Remote Sens. 2022, 15, 133. [Google Scholar] [CrossRef]
  13. Hauschild, A.; Montenbruck, O.; Steigenberger, P.; Martini, I.; Fernandez-Hernandez, I. Orbit determination of Sentinel-6A using the Galileo high accuracy service test signal. GPS Solut. 2022, 26, 120. [Google Scholar] [CrossRef]
  14. Wang, K.; Yang, X.; El-Mowafy, A. Visibility of LEO Satellites under Different Ground Network Distributions. In Proceedings of the 35th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2022), Denver, CO, USA, 19–23 September 2022; pp. 2478–2491. [Google Scholar]
  15. Montenbruck, O.; Ramos-Bosch, P. Precision real-time navigation of LEO satellites using global positioning system measurements. GPS Solut. 2008, 12, 187–198. [Google Scholar] [CrossRef]
  16. El-Mowafy, A.; Wang, K.; Li, Y.; Allahvirdi-Zadeh, A. The impact of orbital and clock errors on positioning from LEO constellations and proposed orbital solutions. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2023, 48, 1111–1117. [Google Scholar] [CrossRef]
  17. Montenbruck, O.; Hackel, S.; Wermuth, M.; Zangerl, F. Sentinel-6A precise orbit determination using a combined GPS/Galileo receiver. J. Geod. 2021, 95, 109. [Google Scholar] [CrossRef]
  18. Wang, K.; El-Mowafy, A. LEO satellite clock analysis and prediction for positioning applications. Geo Spat. Inf. Sci. 2022, 25, 14–33. [Google Scholar] [CrossRef]
  19. Wu, M.; Wang, K.; Liu, J.; Zhu, Y. Relativistic effects of LEO satellite and its impact on clock prediction. Meas. Sci. Technol. 2023, 34, 095005. [Google Scholar] [CrossRef]
  20. Ge, H.; Wu, T.; Li, B. Characteristics analysis and prediction of Low Earth Orbit (LEO) satellite clock corrections by using least-squares harmonic estimation. GPS Solut. 2023, 27, 38. [Google Scholar] [CrossRef]
  21. Suesser-Rechberger, B.; Krauss, S.; Strasser, S.; Mayer-Guerr, T. Improved precise kinematic LEO orbits based on the raw observation approach. Adv. Space Res. 2022, 69, 10. [Google Scholar] [CrossRef]
  22. Yang, Z.; Liu, H.; Qian, C.; Shu, B.; Zhang, L.; Xu, X.; Zhang, Y.; Lou, Y. Real-time estimation of low Earth orbit (LEO) satellite clock based on ground tracking stations. Remote Sens. 2020, 12, 2050. [Google Scholar] [CrossRef]
  23. Yang, Z.; Liu, H.; Wang, P.; Xu, X.; Qian, C.; Shu, B.; Zhang, Y. Integrated kinematic precise orbit determination and clock estimation for low Earth orbit satellites with onboard and regional ground observations. Meas. Sci. Technol. 2022, 33, 125002. [Google Scholar] [CrossRef]
  24. Yang, Y. PPP Based Sequential Orbit Determination for Satellite in Low Earth Orbit. Ph.D. Thesis, Northwestern Polytechnical University, Xi’An, China, 2015. [Google Scholar]
  25. IGS Real-Time Service. Available online: https://igs.org/rts/ (accessed on 26 February 2024).
  26. International GNSS Monitoring & Assessment System. Available online: http://www.igmas.org/Product/Info/ (accessed on 26 February 2024).
  27. Bock, H.; Dach, R.; Jäggi, A.; Beutler, G. High-rate GPS clock corrections from CODE: Support of 1 Hz applications. J. Geod. 2009, 83, 1083–1094. [Google Scholar] [CrossRef]
  28. Rülke, A.; Agrotis, L.; Caissy, M.; Habrich, H.; Neumaier, P.; Söhne, W.; Weber, G. IGS Real-Time Service-Status And Future Developments. In Proceedings of the EGU General Assembly Conference, Vienna, Austria, 27 April–2 May 2014; p. 7429. [Google Scholar]
  29. Hadas, T.; Bosy, J. IGS RTS precise orbits and clocks verification and quality degradation over time. GPS Solut. 2015, 19, 93–105. [Google Scholar] [CrossRef]
  30. Kazmierski, K.; Sośnica, K.; Hadas, T. Quality assessment of multi-GNSS orbits and clocks for real-time precise point positioning. GPS Solut. 2018, 22, 11. [Google Scholar] [CrossRef]
  31. BKG—Bundesamt für Kartographie und Geodäsie. Available online: https://www.bkg.bund.de/ (accessed on 26 February 2024).
  32. Tobías, G.; Calle, J.D.; Navarro, P.; Rodríguez, I.; Rodríguez, D. magicGNSS’ Real-Time POD and PPP Multi-GNSS Service. In Proceedings of the 27th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2014), Tampa, FL, USA, 8–12 September 2014. [Google Scholar]
  33. Guo, L.; Zhao, Q.; Xu, X.; Tao, J.; Zhang, Q.; Qu, Z.; Chen, G.; Wang, C. Real-time orbit and clock products at Wuhan University to support Multi-GNSS applications. In Proceedings of the IGS Workshop, Wuhan, China, 29 October–2 November 2018. [Google Scholar]
  34. Männel, B.; Brandt, A.; Nischan, T.; Brack, A.; Sakic, P.; Bradke, M. GFZ Rapid Product Series for the IGS; GFZ Data Services; GFZ: Potsdam, Germany, 2020. [Google Scholar]
  35. Chen, B.; Wang, K.; El-Mowafy, A.; Yang, X. Real-Time GNSS satellite SISRE and its integrity for LEO satellite Precise Orbit Determination. In Proceedings of the 2024 International Technical Meeting of The Institute of Navigation (ION ITM 2024), Long Beach, CA, USA, 22–25 January 2024. [Google Scholar]
  36. Prange, L.; Orliac, E.; Dach, R.; Arnold, D.; Beutler, G.; Schaer, S.; Jäggi, A. CODE’s five-system orbit and clock solution—The challenges of multi-GNSS data analysis. J. Geod. 2017, 91, 345–360. [Google Scholar] [CrossRef]
  37. Larson, K.M.; Ashby, N.; Hackman, C.; Bertiger, W. An assessment of relativistic effects for low Earth orbiters: The GRACE satellites. Metrologia 2007, 44, 484. [Google Scholar] [CrossRef]
  38. Ko, U.-D.; Tapley, B.D. Computing the USO frequency instability of GRACE satellites. In Proceedings of the 2010 IEEE Aerospace Conference, Big Sky, MT, USA, 6–13 March 2010; pp. 1–8. [Google Scholar]
  39. Hudson, A.I.; Iyanu, G.H.; Wang, H.; Bloch, M.B.; McClelland, T.; Terracciano, L. Long-term Performance of a Space-system OCXO. In Proceedings of the 2019 Joint Conference of the IEEE International Frequency Control Symposium and European Frequency and Time Forum (EFTF/IFC), Orlando, FL, USA, 14–18 April 2019; IEEE: Piscataway, NJ, USA, 2019. [Google Scholar]
  40. Fletcher, K. Sentinel-3: ESA’s Global Land and Ocean Mission for GMES Operational Services; ESA Communications: Oakville, ON, Canada, 2012. [Google Scholar]
  41. Riley, W. Handbook of Frequency Stability Analysis; NIST special publication 1065; NIST: Gaithersburg, MD, USA, 2008.
Figure 1. The flowchart of real-time LEO satellite clock determination.
Figure 1. The flowchart of real-time LEO satellite clock determination.
Remotesensing 16 01326 g001
Figure 2. Precision of the real-time GPS satellite clocks from different analysis centers.
Figure 2. Precision of the real-time GPS satellite clocks from different analysis centers.
Remotesensing 16 01326 g002
Figure 3. The time reference of the CNES real-time GPS satellite clocks on DOY 227, 2018 before and after re-referencing to the satellite clock of G05.
Figure 3. The time reference of the CNES real-time GPS satellite clocks on DOY 227, 2018 before and after re-referencing to the satellite clock of G05.
Remotesensing 16 01326 g003
Figure 4. Sentinel-3B satellite clocks estimated using the CNES real-time clocks on DOY 227, 2018 without and with alignment to the satellite clocks of G05.
Figure 4. Sentinel-3B satellite clocks estimated using the CNES real-time clocks on DOY 227, 2018 without and with alignment to the satellite clocks of G05.
Remotesensing 16 01326 g004
Figure 5. Processing timeline of the real-time LEO satellite clock determination based on near-real-time clock determination and ultra-short-term clock prediction.
Figure 5. Processing timeline of the real-time LEO satellite clock determination based on near-real-time clock determination and ultra-short-term clock prediction.
Remotesensing 16 01326 g005
Figure 6. Comparison between the connected near-real-time Sentinel-3B clocks at the last T period over each processing arc, and different daily solutions, (a): T is 60 s; (b): T is 120 s; (c): T is 180 s; (d): T is 240 s; (e): T is 300 s. (K denotes the daily clock solutions using the Kalman-filter-based kinematic model with the CNES real-time products, RD denotes the daily clock solutions using the reduced-dynamic model with the CODE final products).
Figure 6. Comparison between the connected near-real-time Sentinel-3B clocks at the last T period over each processing arc, and different daily solutions, (a): T is 60 s; (b): T is 120 s; (c): T is 180 s; (d): T is 240 s; (e): T is 300 s. (K denotes the daily clock solutions using the Kalman-filter-based kinematic model with the CNES real-time products, RD denotes the daily clock solutions using the reduced-dynamic model with the CODE final products).
Remotesensing 16 01326 g006
Figure 7. The errors between the connected near-real-time Sentinel-3B clocks at the last T period for each processing arc, and the daily clock solutions using the reduced-dynamic model with the CODE final products to determine the LEO satellite clocks.
Figure 7. The errors between the connected near-real-time Sentinel-3B clocks at the last T period for each processing arc, and the daily clock solutions using the reduced-dynamic model with the CODE final products to determine the LEO satellite clocks.
Remotesensing 16 01326 g007
Figure 8. MDEVs of the connected near-real-time Sentinel-3B clocks at the last T period over each processing arc.
Figure 8. MDEVs of the connected near-real-time Sentinel-3B clocks at the last T period over each processing arc.
Remotesensing 16 01326 g008
Figure 9. Error budget of the real-time LEO satellite clocks for Sentinel-3B from DOY 227 to 229, 2018.
Figure 9. Error budget of the real-time LEO satellite clocks for Sentinel-3B from DOY 227 to 229, 2018.
Remotesensing 16 01326 g009
Table 1. The average processing time ( T p r o c e s s ) and the precision (STD) of the near-real-time clock estimates using different processing arc lengths ( T o b s ). (K denotes Kalman-filter-based daily solution using the kinematic model and the CNES real-time products, RD denotes the batch least-squares reduced-dynamic solutions using the CODE final products).
Table 1. The average processing time ( T p r o c e s s ) and the precision (STD) of the near-real-time clock estimates using different processing arc lengths ( T o b s ). (K denotes Kalman-filter-based daily solution using the kinematic model and the CNES real-time products, RD denotes the batch least-squares reduced-dynamic solutions using the CODE final products).
T o b s (h) T p r o c e s s (s)STD (vs. K) (ns)STD (vs. R) (ns)
2424------
12210.0150.19
8190.0150.19
6180.0150.20
4170.0180.20
2170.0480.22
Table 2. Processing strategy.
Table 2. Processing strategy.
TypeParametersProcessing Strategies
Near-real-time clock determinationObservationsUndifferenced IF code and carrier phase combination
(GPS: L1/L2)
Sampling30 s
Elevation cut-off angle
GNSS orbits and clocksSelected (Section 2.1)
GNSS satellite antenna PCO/PCVsigs20.atx
LEO satellite antenna PCO/PCVsOperator supplied
Phase wind-upCorrected
EstimatiorKalman filter, Kinematic
T o b s 6 h
Clock Prediction T p r e Longer than T p r o c e s s + T
T 60/120/180/240/300 s
Broadcast intervalOptional
Table 3. Daily STDs between the connected near-real-time clocks for different T and different daily clocks (K denotes Kalman-filter-based daily kinematic clocks using the CNES real-time products, R denotes the batch least-squares reduced-dynamic daily clocks using the CODE final products).
Table 3. Daily STDs between the connected near-real-time clocks for different T and different daily clocks (K denotes Kalman-filter-based daily kinematic clocks using the CNES real-time products, R denotes the batch least-squares reduced-dynamic daily clocks using the CODE final products).
DOYSTD (K vs. R) (ns) T (s)STD (vs. K) (ns)STD (vs. R) (ns)
2270.211600.0250.218
1200.0210.218
1800.0180.218
2400.0150.215
3000.0150.201
2280.297600.0250.304
1200.0210.300
1800.0180.300
2400.0140.299
3000.0250.301
2290.237600.0280.246
1200.0240.245
1800.0200.245
2400.0170.241
3000.0280.230
Table 4. The daily STDs of the connected real-time clocks according to different T (R denotes the daily reduced-dynamic clocks based on batch least-squares adjustment using the CODE final products).
Table 4. The daily STDs of the connected real-time clocks according to different T (R denotes the daily reduced-dynamic clocks based on batch least-squares adjustment using the CODE final products).
DOY T (s) T p r e (s)Fitting Time (s)STD of Near-Real-Time Clocks Used for Fitting (ns)Precision Loss (ns)STD (vs. R) (ns)
22760901200.2190.0640.229
1201501800.2190.0840.246
1802101800.2220.1060.247
2402702700.2110.1270.264
3003303000.2230.1560.297
22860901200.3040.0690.311
1201501800.3010.0920.326
1802101800.3010.1150.329
2402702700.2990.1380.341
3003303000.3030.1690.342
22960901200.2470.0680.266
1201501800.2450.0900.279
1802101800.2460.1270.295
2402702700.2430.1360.312
3003303000.2280.1730.327
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, M.; Wang, K.; Wang, J.; Liu, J.; Chen, B.; Xie, W.; Zhang, Z.; Yang, X. Real-Time LEO Satellite Clocks Based on Near-Real-Time Clock Determination with Ultra-Short-Term Prediction. Remote Sens. 2024, 16, 1326. https://doi.org/10.3390/rs16081326

AMA Style

Wu M, Wang K, Wang J, Liu J, Chen B, Xie W, Zhang Z, Yang X. Real-Time LEO Satellite Clocks Based on Near-Real-Time Clock Determination with Ultra-Short-Term Prediction. Remote Sensing. 2024; 16(8):1326. https://doi.org/10.3390/rs16081326

Chicago/Turabian Style

Wu, Meifang, Kan Wang, Jinqian Wang, Jiawei Liu, Beixi Chen, Wei Xie, Zhe Zhang, and Xuhai Yang. 2024. "Real-Time LEO Satellite Clocks Based on Near-Real-Time Clock Determination with Ultra-Short-Term Prediction" Remote Sensing 16, no. 8: 1326. https://doi.org/10.3390/rs16081326

APA Style

Wu, M., Wang, K., Wang, J., Liu, J., Chen, B., Xie, W., Zhang, Z., & Yang, X. (2024). Real-Time LEO Satellite Clocks Based on Near-Real-Time Clock Determination with Ultra-Short-Term Prediction. Remote Sensing, 16(8), 1326. https://doi.org/10.3390/rs16081326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop