Original Research ARTICLE
Relativistic Time Transfer for Inter-satellite Links
- 1Department of Astronomy, School of Astronomy and Space Sciences, Nanjing University, Nanjing, China
- 2Shanghai Key Laboratory of Space Navigation and Position Techniques, Shanghai, China
- 3Key Laboratory of Modern Astronomy and Astrophysics, Nanjing University, Ministry of Education, Nanjing, China
Inter-Satellite links (ISLs) will be an important technique for a global navigation satellite system (GNSS) in the future. Based on the principles of general relativity, the time transfer in an ISL is modeled and the algorithm for onboard computation is described. It is found, in general, satellites with circular orbits and identical semi-major axes can benefit inter-satellite time transfer by canceling out terms associated with the transformations between the proper times and the Geocentric Coordinate Time. For a GPS-like GNSS, the Shapiro delay is as large as 0.1 ns when the ISL passes at the limb of the Earth. However, in more realistic cases, this value will decrease to about 50 ps.
Inter-Satellite links (ISLs), also called crosslinks, will be a promising technique for enhancing the reliability and integrity of a global navigation satellite system (GNSS). Since the early stage of the Global Positioning System (GPS), ISLs have been a part of it. ISLs are applied to transfer the Nuclear Detonation Detection System data and the autonomous navigation and ranging data between space vehicles (Sonntag, 1997). For the Block IIR satellites of GPS, the pseudorange measurements with ISLs are used by the onboard computers to update the stored navigation messages, including the ephemeris and clock states, which endows the satellites with the ability to maintain accuracy for a specified period of time after loss of contact with the control segment (Abusali et al., 1998). The network system of ISLs and its uses for the next-generation GPS was discussed (Maine et al., 2003). The European Space Agency (ESA) also assessed the potential improvements on a GNSS for navigation and dissemination performance by introducing ISLs for ranging and communication (Fernández, 2011). It was officially reported1 that China is developing and testing the new-generation Beidou Navigation Satellite System, which will have a new system of navigation signaling and ISLs.
In order to fulfill these practical purposes, time signals of onboard clocks together with other information need to be transfered among the satellites through ISLs for ranging and computing clock offsets. In the procedure of time transfer, Einstein's general relativity (GR) has “long since passed from the realm of theoretical physics to the realm of engineering design” (Nelson, 2011). In fact, GR is an inevitable part of a GNSS (Ashby, 2003; Han et al., 2011). In the light of the principles of GR, it is necessary to abandon the Newton's absolute time for different kinds of times: proper times and coordinate times (Misner et al., 1973; Landau and Lifshitz, 1975).
Theoretically, the readings of an ideal clock give the proper time τ, which is an observable and only belongs to the clock itself. However, there exists no such thing as an ideal clock. An atomic clock, which is widely used on the ground and in space, drifts from the ideal clock with some known and unknown factors. The coordinate times cannot be measured directly. However, they might be used as variables in the equations of motion of celestial and artificial bodies and light rays. The coordinate times are connected with the proper time through the four-dimensional spacetime interval, which changes the way of clock synchronization and time transfer significantly (Petit and Wolf, 2005; Nelson, 2007, 2011). Experiments involving time/frequency transfer might also be used for testing theories of gravity (Samain, 2002; Cacciapuoti and Salomon, 2009; Christophe et al., 2009; Wolf et al., 2009; Christophe et al., 2012; Deng and Xie, 2013a,b, 2014; Hees et al., 2014; Zhang et al., 2014; Delva et al., 2015; Xie and Huang, 2015).
The relativistic time transfers in various contexts, such as in the vicinity of the Earth (Klioner, 1992; Petit and Wolf, 1994; Wolf and Petit, 1995; Petit and Wolf, 1997; Kouba, 2002, 2004; Petit and Wolf, 2005; Nelson, 2007, 2011) and in the Solar System (Nelson, 2007, 2011; Deng, 2012; Pan and Xie, 2013, 2014, 2015), have been intensively studied and discussed. In such contexts, one clock is on-board a space vehicle and the other clock or a user is usually on the ground. However, for ISLs, each end of a link is a satellite with clocks in orbital motion, which makes the description of the time of light propagation in the ISL different from previously mentioned cases. In the present work, we will model the time transfer through ISLs and discuss its algorithm for onboard computation, where the relativistic effects are fully taken into account.
The paper is organized as follows. Section 2 is devoted to modeling the time transfer through ISLs. Some mathematical details for this purpose can be found in the Appendix. The algorithm for onboard computation is also discussed. Finally, in Section 3, we summarize our results and discuss their implication.
2. General Relativistic Description of Time Transfer in an ISL
We consider an ISL between two satellites 1 and 2 in the vicinity of the Earth. It is also assumed that these two satellites belong to a GNSS like GPS. The proper times of their onboard clocks are τ1 and τ2. Satellite 2 emits electromagnetic signals encoded with its proper time of this emission of τ2 and other necessary information, such as its ephemeris and clock offset. Satellite 1 receives these signals at its proper time τ1, decodes out τ2 and other messages and computes τ1 − τ2. In the framework of the International Astronomical Union (IAU) 2000 Resolutions for general relativistic reference systems (Soffel et al., 2003), we describe this problem in the Geocentric Celestial Reference System (GCRS), which is physically adequate to describe processes occurring in the vicinity of Earth. For the time transfer in the ISL, we can have
where T1 and T2 are the time coordinates of GCRS, called the Geocentric Coordinate Time (TCG).
The first and third terms on the r.h.s of Equation (1) account for the relativistic four-dimensional transformations between the proper time and TCG. Considering a post-Newtonian development and keeping all terms larger than one part in 1018, the first term on the r.h.s of Equation (1) reads (Wolf and Petit, 1995; Petit and Wolf, 1997)
Here, ϵ = c−1 and c is the speed of light; UE(X1) is the Newtonian gravitational potential of the Earth evaluated at the position of satellite 1, X1; V1 is the velocity of satellite 1 in GCRS; Ū(x) is the Newtonian gravitational potential of external masses (except the Earth) evaluated at the coordinate of x; xE is the coordinate of the Earth's center of mass in the Solar System Barycentric Celestial Reference System (BCRS); the Q term is related to the four-acceleration of the Earth's center of mass in the external gravitational field due to its mass quadrupole and it is numerically at the order of ~ 10−11 m s−2. It was found in Kouba (2004) that, at the pico-second (ps) precision level (or the inaccuracy of frequency of 10−16), the effects of external masses can be neglected in time transfer for GPS. Therefore, we can simplify Equation (2) as
and we get a similar equation for the third term in Equation (1) as
It was also pointed out in Kouba (2004) that only the first oblateness term J2 in UE is necessary for the orbital altitude of GPS satellites. In the GPS relativistic transformations between the proper time and TCG, both the orbit perturbations and the Earth gravity field oblateness J2, which is a main cause of the perturbations, need to be taken into account together (Kouba, 2004).
The second term in Equation (1) describes the propagation of the electromagnetic signal from satellite 2 to satellite 1 in TCG. For a one-way signal transmission, we can have the coordinate time interval as (Blanchet et al., 2001)
Here, ϵF1 is the Euclidean geometric effect and is the Shapiro time delay (Shapiro, 1964; Weinberg, 1972) where we only consider the monopole component of the Earth's potential. R, X, and are quantities associated with the propagation of the signal so that
Substituting Equations (3–5) into Equation (1), we can have
Leaving the details of integration in the Appendix, we can finally have
where the terms due to Euclidean geometry, ΔτE, time transformations for satellite 1 and 2, ΔτS1 and ΔτS2, and the Shapiro delay, ΔτSh, are
In order to calculate all of the above terms, models of clock offsets and ephemerides for both satellites are required. In principle, according to the model of clock offset used in GNSS data processing (Ashby and Spilker, 1996), the difference between the proper time τ and T in TCG can be determined as
where c0, c1, and c2 are the coefficients of the clock offset model, a is the orbital semi-major axis of a satellite, e is the orbital eccentricity and E is the eccentric anomaly. All these quantities are available in the navigation messages.
For the Euclidean term ΔτE, T1, and T2 can be obtained by the clock offset models by Equation (17), and X1(T1) and X2(T2) can be calculated by the ephemerides, which are basically represented in the form of orbital elements and are also available in the navigation messages. If we consider a GPS-like GNSS that a1 ≈ a2 ≈ a and e1 ≈ e2 ≈ 0, where a1 and a2 are the orbital semi-major axes of satellite 1 and 2, e1 and e2 are the orbital eccentricities and a ≈ 26559.7 km, we can find that |X1(T1)| ≈ |X2(T2)| ≈ a, and ΔτE ≲ 0.172 s, where R⊕ is the radius of the Earth. In order to maximize |X1(T1)−X2(T2)|, the ISL needs to pass through the Earth, which never happens in practice. In a more realistic situation, when X1(T1) is perpendicular to X2(T2), ΔτE is about 0.125 s.
For the time transformation terms, when the oblateness of the Earth J2 can be neglected, it can be found that (Nelson, 2011)
where Δτ = τ1 − τ2 and E1 and E2 are the eccentric anomalies of the two satellites. When satellite 1 and 2 are set in the orbits with almost the same size and close to circles, ΔτS1 + ΔτS2 will nearly cancel out each other. In the case that J2 has to be taken into account, additional corrections will be needed and their mathematical expressions can be found in Ashby (2003); Kouba (2004); Nelson (2011). For a GPS satellite with an inclination of i = 55.0°, the drift rate of the secular term due to J2 is at the level of −3.4 × 10−17 and the amplitude of the periodic effect is 24 ps. For a GPS-like system but with different inclinations, the maximum value of the drift rate is about 5.2 × 10−15 for i = 0° and the maximum value of the amplitude of the periodic effect is 36 ps for i = 90°. Because of its dependence of the orbital inclination, the argument of perigee and the true anomaly, the J2-associated component in ΔτS1 + ΔτS2 will generally not cancel out each other.
For the Shapiro time delay, if we consider that a1 ≈ a2 ≈ a and e1 ≈ e2 ≈ 0, we can find that the maximum value of ΔτSh = 0.12 ns when the ISL is passing the limb of the Earth. However, in more probable cases, this value will be much smaller. For example, when X1(T1) is perpendicular to X2(T2), ΔτSh is about 52 ps.
It was reported in Hadas and Bosy (2015) that, for GPS, real-time clocks accuracy is 0.28 ns and orbits accuracy is 5 cm (0.17 ns). Based on the above estimation, the J2 effects and the Shapiro time delay can be neglected for ISLs in a GPS-like system.
3. Conclusions and Discussion
In the framework of GR, we investigate the time transfer for ISLs, which will be a promising technique for a GNSS in the future. We give a model describing the difference between two onboard clocks of an ISL and discuss its algorithm for onboard computation. It is found that the difference between the proper times of two satellites contains the contributions from Euclidean geometry, time transformations between the proper time and TCG for these satellites and the Shapiro delay. For a GPS-like system, the net effect of the time transformations is close to 0 because they cancel out each other and the J2 effects in these transformations are less than several tens of ps. The Shapiro delay is about 52 ps in such a system. For current clocks accuracy of GPS, the effects of J2 and the Shapiro delay can be neglected. After all of the differences among the clocks in a GNSS are obtained, the ephemeris and the clock states can be evaluated and updated to maintain the system.
In practice, some additional factors need to be taken into account as well, such as the time delay caused by the Earth's atmosphere and instruments. Due to the limited capability of on-board computation, the question about how to optimize the constellation of a GNSS and the measurements through ISLs is another important issue. More detailed and specific case studies are definitely required.
In the future, when optical clocks and laser links (Chou et al., 2010a; Chou et al., 2010b; Predehl et al., 2012; Bloom et al., 2014) are deployed and applied in space, the external gravitational effLects can no longer be neglected and will have to be carefully included (Blanchet et al., 2001; Soffel et al., 2003; Petit and Luzum, 2010). Moreover, an accurate processing of ISL data might also offer a new opportunity for the application of GNSS in fundamental physics and astronomy (Wolf and Petit, 1997; Bertolami and Páramos, 2011; Kentosh and Mohageg, 2012; Aoyama et al., 2014).
The author confirms being the sole contributor of this work and approved it for publication.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Editorial Office is aware of co-publications from the Associate Editor and Reviewer SB, and confirms the review was carried out to the highest degree of professionalism, and is free of bias.
I acknowledge the helpful comments and suggestions from my referees. The work is funded by the National Natural Science Foundation of China under Grant No. 11573015 and the Opening Project of Shanghai Key Laboratory of Space Navigation and Position Techniques under Grant No. 14DZ2276100.
Abusali, P. A. M., Tapley, B. D., and Schutz, B. E. (1998). Autonomous navigation of global positioning system satellites using cross-link measurements. J. Guid. Control Dynam. 21, 321–327. doi: 10.2514/2.4238
Aoyama, S., Tazai, R., and Ichiki, K. (2014). Upper limit on the amplitude of gravitational waves around 0.1 Hz from the Global Positioning System. Phys. Rev. D 89:067101. doi: 10.1103/PhysRevD.89.067101
Ashby, N. (2003). Relativity in the global positioning system. Living Rev. Relativ. 6, 1. doi: 10.12942/lrr-2003-1. Available online at: http://relativity.livingreviews.org/Articles/lrr-2003-1/
Ashby, N., and Spilker, J. (1996). “Introduction to relativistic effects on the global positioning system,” in Global Positioning System: Theory and Applications, Vol. 1, eds B. Parkinson, J. Spilker Jr., P. Axelrad, and P. Enge (Reston, VA: American Institute of Aeronautics and Astronautics), 623–697.
Bloom, B. J., Nicholson, T. L., Williams, J. R., Campbell, S. L., Bishof, M., Zhang, X., et al. (2014). An optical lattice clock with accuracy and stability at the 10−18 level. Nature 506, 71–75. doi: 10.1038/nature12941
Chou, C. W., Hume, D. B., Koelemeij, J. C. J., Wineland, D. J., and Rosenband, T. (2010a). Frequency comparison of two high-accuracy al+ optical clocks. Phys. Rev. Lett. 104:070802. doi: 10.1103/PhysRevLett.104.070802
Christophe, B., Spilker, L. J., Anderson, J. D., André, N., Asmar, S. W., Aurnou, J., et al. (2012). OSS (Outer Solar System): a fundamental and planetary physics mission to Neptune, Triton and the Kuiper Belt. Exp. Astron. 34, 203–242. doi: 10.1007/s10686-012-9309-y
Delva, P., Hees, A., Bertone, S., Richard, E., and Wolf, P. (2015). Test of the gravitational redshift with stable clocks in eccentric orbits: application to Galileo satellites 5 and 6. Class. Quant. Grav. 32:232003. doi: 10.1088/0264-9381/32/23/232003
Fernández, F. A. (2011). Inter-satellite ranging and inter-satellite communication links for enhancing GNSS satellite broadcast navigation data. Adv. Space Res. 47, 786–801. doi: 10.1016/j.asr.2010.10.002
Hees, A., Folkner, W. M., Jacobson, R. A., and Park, R. S. (2014). Constraints on modified Newtonian dynamics theories from radio tracking data of the Cassini spacecraft. Phys. Rev. D 89:102002. doi: 10.1103/PhysRevD.89.102002
Nelson, R. A. (2007). “Relativistic time transfer in the solar system,” in Frequency Control Symposium, 2007 Joint with the 21st European Frequency and Time Forum. IEEE International (Geneva), 1278–1283.
Pan, J.-Y., and Xie, Y. (2013). Relativistic transformation between τ and TCB for Mars missions: Fourier analysis on its accessibility with clock offset. Res. Astron. Astrophys. 13, 1358–1362. doi: 10.1088/1674-4527/13/11/006
Pan, J.-Y., and Xie, Y. (2015). Relativistic algorithm for time transfer in Mars missions under IAU Resolutions: an analytic approach. Res. Astron. Astrophys. 15, 281–292. doi: 10.1088/1674-4527/15/2/011
Predehl, K., Grosche, G., Raupach, S. M. F., Droste, S., Terra, O., Alnis, J., et al. (2012). A 920-kilometer optical fiber link for frequency metrology at the 19th decimal place. Science 336:441. doi: 10.1126/science.1218442
Soffel, M., Klioner, S. A., Petit, G., Wolf, P., Kopeikin, S. M., Bretagnon, P., et al. (2003). The IAU 2000 resolutions for astrometry, celestial mechanics, and metrology in the relativistic framework: explanatory supplement. Astron. J. 126, 2687–2706. doi: 10.1086/378162
Wolf, P., Bordé, C. J., Clairon, A., Duchayne, L., Landragin, A., Lemonde, P., et al. (2009). Quantum physics exploring gravity in the outer solar system: the SAGAS project. Exp. Astron. 23, 651–687. doi: 10.1007/s10686-008-9118-5
Xie, Y., and Huang, Y. (2015). Spacecraft Doppler tracking with possible violations of LLI and LPI: upper bounds from one-way measurements on MEX. Res. Astron. Astrophys. 15:1751. doi: 10.1088/1674-4527/15/10/012
Zhang, Y.-F., Zhang, X.-Z., Liu, J.-H., Huang, Y., and Xie, Y. (2014). Spacecraft Doppler tracking with possible violations of LLI and LPI: preliminary bounds on LLI from Mars Express. Res. Astron. Astrophys. 14, 1201–1206. doi: 10.1088/1674-4527/14/9/010
Appendix: Derivation of Equation (12)
Integrating Equation (11), we can have
Since Equation (4) can lead to
and Equation (5) can give
we can combine some terms in Equation (A1) as
where different forms for Shapiro time delay can also be found in Nelson (2011). By putting all of the three equations above into Equation (A1), we can finally obtain that
In fact, the difference between the proper time and coordinate time can be neglected in the terms with the orders of ϵ2 and ϵ3 in the above equation.
Keywords: general relativity, time transfer, inter-satellite link, time comparison, global navigation satellite system
Citation: Xie Y (2016) Relativistic Time Transfer for Inter-satellite Links. Front. Astron. Space Sci. 3:15. doi: 10.3389/fspas.2016.00015
Received: 03 February 2016; Accepted: 12 April 2016;
Published: 26 April 2016.
Edited by:Christophe Le Poncin-Lafitte, Observatoire de Paris, France
Reviewed by:Sven Zschocke, Technical University Dresden, Germany
Stefano Bertone, University of Bern, Switzerland
Copyright © 2016 Xie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yi Xie, email@example.com