Metrology in the Galileo Navigation System The Experience of the Italian National Metrology Institute

I.Sesia, G.Signorile, G.Cerretto, E.Cantoni, P.Tavella A.Cernigliaro, A.Samperi Optics Division, INRiM, Turin, Italy DASS Division, aizoOn, Turin, Italy [email protected] [email protected], [email protected]

Abstract—Timekeeping is crucial in Global Navigation system, the research on new technologies, calibrations, Satellite Systems (GNSS), being the positioning accuracy directly evaluation of uncertainty, and dissemination of accurate time. related to a time measurement. As a consequence, the typical expertise of time metrology laboratories is necessary in many This paper presents the experience of the Italian National different aspects of a navigation system. This paper presents the Institute of Metrological Research (INRiM) time laboratory in experience of INRIM in the development of the Galileo the Galileo experimental and validation phases. navigation system from the earlier studies till the very recent In Orbit Validation phase showing how the time metrology practice II. THE GALILEO EXPERIENCE has been useful in the understanding of time aspects of INRiM has been involved in the Galileo system since 1999 navigation and showing as well how the navigation service and participated to different phases of the project. perspective has stimulated new ideas and better understanding of time measures. The first experimental phase was the Galileo System Test Bed Version 1 (GSTB V1) in 2002 in which INRiM participated, Keywords—time metrology; atomic ; steering; time together with the British NPL and the German PTB scales; GNSS timing; timekeeping; space clocks; system noise laboratories, to generate the experimental Galileo , the reference time scale of the system, obtained from the I. INTRODUCTION Experimental Precise Timing Station located in INRIM. GSTB The European Union (EU) and the European Space Agency V1 supported the validation of the on-ground algorithms for (ESA) are building the European satellite navigation system, orbit determination and time synchronization. In that phase the Galileo. In addition to the positioning service, Galileo will be main learned lesson was the difficulty of generating a reference able to disseminate globally accurate time. time scale, synchronized to UTC, in a complete automatic and unmanned way, as Galileo requested. It appeared clear that not Currently there are 4 Galileo satellites already in orbit, only optimal algorithms and devices are necessary, but that launched in 2011 and 2012 for the In Orbit Validation (IOV) robustness was the issue and different monitoring and control phase. Two experimental satellites were launched in 2005 and checks had to be devised [2]. This project constituted an 2008 and operated in the GIOVE Mission phase. The essential contribution to the understanding of the Galileo remaining satellites of the system will be launched between mission segment development. 2014 and 2018, in order to complete the constellation for the Initial Operational Capability (IOC) and the Full Operational The next phase, from 2005 to 2011, was the Galileo System Capability (FOC) phases [1]. Test Bed Version 2 (GSTB-V2), then named GIOVE Mission [3-4], in which INRiM role consisted in the realization of the During these , the role of European metrological institutions, which cooperated and worked together with space experimental Galileo System Time, driven by the free running industries and the ESA, has been significant and fruitful. Active Hydrogen Maser connected to an experimental Galileo receiver, named GIEN, located at INRiM, and the National Metrology Institutes (NMI) have the instruments, characterization of the clocks of two experimental satellites expertise, and competencies required to ensure the high GIOVE-A and GIOVE-B, in addition to the ground station performances of time measurements as necessary in navigation clocks. systems. In fact, NMI time laboratories realize their national Ground and space clocks performances were evaluated time scales, called UTC(k), and are hence used to face the through systematic monitoring and campaign activities. Both issues related to the accuracy, reliability, and robustness of the short and long term clock behaviors were characterized timing systems. The international reference time is the [5-7]. Coordinated (UTC), and UTC(k) are local real Particular care was given to the detection and identification of time approximation as realized by the k laboratory. feared events in satellite clocks, as for example the frequency Time metrology is implicated in several aspects of a navigation jumps that sometimes could be observed over long data system: from the reference time scale definition to the design periods, as well as to the correlation of clocks’ behavior and of specific algorithms for clock characterization and anomaly clock internal telemetries [8-9]. Moreover, specific analysis detection, as well as in the evaluation of the measurement were performed on clock prediction strategies [10]. The results of this phase have been used for the risk mitigation after the publication of the BIPM Circular T containing the for the IOV subsequent phase [5]. final UTC data. Fig 1 shows the |UTC-GST|(mod 1s) offset in the period May-Oct 2013, which remained less than 10 ns. From 2011 to 2013, within the In Orbit Validation Phase A product of the TVF was the daily prediction of the (IOV), INRiM role consisted in the support to the In Orbit Test offset UTC −GST to be transmitted in the Galileo navigation (IOT) after the launch of the first 4 IOV Galileo satellites and message to disseminate UTC to the user. The true offset in the design and development of the Time Validation Facility, UTC −GST is known only a posteriori, after BIPM Circular T to assess and characterize all the timing aspects of the Galileo publication, one a posteriori. system [11]. The In Orbit Validation satellites were launched During the test period, as showed in Fig 2, the prediction error on October 2011 and October 2012 from the Guiana Space remained below ±8 ns [11, 19]. Centre, and their design was already representative of that of the final Galileo constellation. The fundamental time metrology aspects in the support to the Galileo experimental phases will be discussed in the next sections.

III. TIME METROLOGY ACTIVITIES A. Timing Infrastructure The most recent project in which INRiM has been involved from 2010 to 2013 is the design and development of the Time Validation Facility (TVF), in the frame of the Time and

Geodetic Validation Facility (TGVF) project. The TVF had two main tasks [11]: Fig. 1. GST validation: UTTC-GST( modulo 1s) time offset • the operation as a preliminary Galileo Time Service (May-October 2013) Provider (TSP) during IOV phase, computing the steering corrections needed to align the Galileo System Time to the international reference time UTC • the validation of the whole timing infrastructure of the IOV phase, including space and ground clocks and Galileo System Time (GST) performances. To accomplish these objectives, the TVF has worked in cooperation with some European laboratories: NPL (UK), OP (F), ORB (Belgium), PTB (Germany), ROA (Spain), and, recently, SP (Sweden), contributing with their UTC(k) national time scales. B. Reference Time Scale Validation In order to validate the Galileo reference time scale, the TVF collected the data from the different European time laboratories involved in the project, as well as from external Fig. 2. UTC-GST prediction error (May-October 2013) entities as the BIPM. The difference between each laboratory time scale and GST, namely UTC(k)-GST, was measured C. System Clocks Characterization using time transfer data obtained at the time laboratories and An additional task of the TVF was the characterization and at the Galileo Precise Time Facility (PTF) realizing GST, monitoring of the satellite and ground station clocks of the using the Two Way Satellite Time and Frequency Transfer system. Similar task had already been performed during the (TWSTFT) stations and the GPS timing receiver techniques. GIOVE Mission phase of the Galileo project [5]. Starting from these offsets, a real time approximation of System clocks were analyzed by means of proper UTC −GST, named UTCapprox−GST , was obtained. Finally, mathematical algorithms, exploiting the methods and software such offset was used to compute the steering corrections regularly used by INRiM for institutional purposes [12-19, required to align the GST to UTC. These corrections were 20]. These activities were particularly important during the In regularly sent to the Galileo PTF since February 2013, in order Orbit Test phases performed after the launch of the Galileo to maintain the time offset |UTC−GST|(mod 1s) within the IOV satellites. By evaluating the frequency offset of the clock specified limits, as asked by Galileo requirements. The TVF on board with respect to ground clock, in particular for the verified the fulfilment of such requirement on a monthly basis, high performance space Passive Hydrogen Maser, it has been possible to verify the effect of general relativity shifting the From the difference between the two previous quantities, it is frequency of the on board clocks [5]. possible to compute UTC(SIS) −UTC(ORB) which can be Additionally the clock on board showed sometimes behavior compared with the reference value of UTC−UTC(ORB) not typically observed on ground, for example the frequency published in the BIPM Circular T, as well as with the BIPM values could suddenly “jump” to another value. This atypical rapid data for UTCr-UTC(ORB), available with a shorter behavior called for the development of additional clock latency. Fig 4 shows the preliminary results of UTC analysis tools, also studied for INRiM institutional activities, dissemination, where the accuracy level of Galileo broadcast that could be integrated in a real time monitoring of GNSS values was 10 ns. satellite clocks [16-18] Galileo Navigation message is broadcasting also the GPS to D. Time Dissemination Galileo time offset (GGTO), which is fundamental for GPS- Galileo started in April 2013 to broadcast the conversion Galileo interoperability [23]. The TVF validated the broadcast parameters required to obtain the offset between GST and GGTO comparing it with the same quantity evaluated at TVF UTC, as well as the GST to GPS Time Offset (GGTO) [21]. using the measures of the PTF GPS receiver directly The Time Validation Facility verified the Galileo time connected to GST, and with the GGTO computed by dissemination service, analyzing the broadcast navigation processing the measures of the ORB receiver pseudo-range message. To this purpose, the UTC-GST (mod 1s) offset measurements, as illustrated in Fig 5. predicted and commanded by TVF was routinely compared with the value reconstructed at user level through the conversion parameters broadcast in the navigation message. Fig 3 shows the comparison between the broadcast and reference value computed by TVF during the period June - August 2013. The differences between the two values are due to the rounding of the parameters reported in the navigation message, and the delay in the upload of the computed parameters.

Fig. 4. Signal in Space Validation: comparison among UTC(SIS)- UTC(ORB), UTC-UTC(ORB) and UTCr-UTC(ORB)

Fig. 3. Signal in Space Validation: Comparison between the broadcast GST-UTC(mod 1s) and the reference value computed at TVF level (from June 1st, 2013 to September 2nd, 2013)

The analysis of UTC Time and Frequency accuracy dissemination at user level was the result of the collaboration with the Royal Observatory of Belgium (ORB).

A GPS/Galileo receiver hosted, operated, and calibrated by Fig. 5. Signal in Space Validation: Comparison among the the broadcast ORB was connected to the Belgian time scale UTC(ORB), GGTO and the one estimated at user level, by TVF and allowed to obtain the [GST(SIS)−UTC(SIS)]p from the Signal In Space (SIS) that is broadcast navigation message, E. System Noise Evaluation and the difference between the local clock and the Galileo A typical metrologist job is to evaluate the noise of the System Time GST(SIS)−UTC(ORB) obtained by processing measurement system. Generally, in satellite navigation the pseudo range measurements of the ORB receiver [22]. systems, the orbit and clock estimates are generated by a least squares algorithms, named Orbit Determination and Time REFERENCES Synchronization [3,24], which process pseudorange measures to estimate satellite and ground station position and clock [1] A.Mudrak, "Galileo IOV Status and First Results", Joint UFFC, offsets. EFTF and PFM Symposium, July 2013, Prague The clock estimates represent the offset between each system [2] F. Cordara, R. Costa, L. Lorini, D. Orgiazzi, V. Pettiti, I. Sesia, P. clock and the ground reference clock. Such offset estimates Tavella, P. Elia, M. Mascarello, M. Falcone, J. Hahn, “Galileo System Test Bed V1: results on the Experimental Galileo System have to be considered “apparent” clock estimates since they Time”, in Proc. 19th European Frequency and Time Forum, include also the contribution of other instabilities besides the Besancon, 2005, on CD true clock behavior. The estimate of the system noise [3] Píriz, R., Fernández, V., Auz, A., Tavella, P., Sesia, I., Cerretto, introduced by this complex network algorithms to provide the G., Falcone, M., Navarro, D., Hahn, J., González, F., Tossaint, M., clock estimates was carried out [24] showing that the Gandara, M., "The Galileo System Test Bed V2 for Orbit and Clock Modeling," Proceedings of the 19th International Technical measurement system, based at the epoch on a limited number Meeting of the Satellite Division of The Institute of Navigation of ground stations, is barely capable to estimate the true (ION GNSS 2006), Fort Worth, TX, September 2006, pp. 549-562 quality of the on board passive H maser, being the noise of the [4] Sesia, I., Cerretto, G., Tavella, P., "Timing Activities at INRiM in measurement system at the same level of the maser instability. the Frame of the Galileo Project," Proceedings of the 40th Annual Precise Time and Time Interval Systems and Applications Meeting, In addition to the satellite clocks, we evaluated also the system Reston, VA, December 2008, pp. 641-656 noise level when assessing ground clock behavior, both in [5] P. Waller, F. Gonzalez, S.Binda, I.Sesia, I.Hidalgo, G.Tobias, Galileo and other GNSS systems. To this purpose, we P.Tavella, “The In Orbit Performances of GIOVE Clocks”, IEEE evaluated the system noise as the instability added to the Transactions on Ultrasonics, Ferroelectrics, and Frequency estimate of the offset between two Active Hydrogen Maser Control, VOL. 57, NO.3. March 2010 (AHM) clocks connected to ground stations of the network: [6] P.Waller, F.Gonzalez, S.Binda, I.Sesia, P.Tavella, I.Hidalgo, G.Tobias, “Update on the In-orbit Performance of GIOVE for example one in INRIM, Turin, Italy and one in USNO, Clocks”, in Proc. of EFTF-IFCS 2009, April 2009, Besançon, Washington, USA. Such H masers have an intrinsic stability France, pp. 388-392 which cannot be appreciated by the current state of the art time [7] Pierre Waller, Francisco Gonzalez, Stefano Binda, Daniel transfer systems, at least in the short term, as the measurement Rodriguez, Guillermo Tobias, Alice Cernigliaro, Ilaria Sesia, noise exceeds the H maser stability. The comparison of two Patrizia Tavella, "LONG-TERM PERFORMANCE ANALYSIS OF GIOVE CLOCKS", in Proc. of the 42th Precise Time and Time remote Active H masers allows setting a lower limit to the Interval (PTTI) Meeting, November 15th-18th 2010, Hyatt noise of the time transfer methods. Regency Reston Town Center, Reston, Virginia (USA) [8] Cerretto, Cueto, Delporte, Droz, Falcone, Fernandez, Gandara, Hahn, Hannes, Langley, Mandorlo, Navarro, Piriz, Rochat, Sesia, IV. CONCLUSIONS Tavella, Waller, “Time for GIOVE-A, the onboard Rubidium Clock Experiment”, GPS World, pp.64-69, May 2007 The INRIM time laboratory has been involved in Galileo [9] J. Hahn, F. Gonzalez, P. Waller, D. Navarro-Reyes, R. Píriz, A. activities since about 15 years. Time metrology algorithms and Mozo, V. Fernández, M. Cueto, P. Tavella, I. Sesia, "Giove-A expertise have been a support for the GNSS system Apparent Clock Assessment and Results", in Proceedings of 39th development, from the design to the validation phase. During Annual Precise Time and Time Interval (PTTI) Applications and Planning Meeting, November 2007,Long Beach, California, USA, the IOV phase of the Galileo project INRiM hosted and pp. 95-114 operated the Time Validation Facility. The fruitful cooperation [10] F. Gonzalez, A.Cernigliaro, P.Tavella, J.P Boyero, "Clock Strategy and support of different European laboratories have been key Experimentation with GIOVE Clocks", in Proc. of the 24th element for the success of the project. The paper aimed to European Frequency and Time Forum (EFTF), April 13th-16th show the timing results achieved in Galileo as well as the 2010, ESTEC; Noordwijk, The Nederland. cross fertilization that may arise between timekeeping and [11] I. Sesia, G.Signorile, G.Cerretto, E.Cantoni, P.Tavella, A.Cernigliaro, A.Samperi, “The first of the Galileo navigation systems. timekeeping and time dissemination: the role of the Time Validation Facility” ION PTTI, December 2013, Bellevue, Washington [12] G. Panfilo, P. Tavella, " Prediction based on ACKNOWLEDGMENT stochastic differential equations", Metrologia 45, 6, (2008) S108- The recent INRIM experience with the Time Validation S116 Facility is the result of the proactive and collaborative effort of [13] L. Galleani, P. Tavella, " Detection and identification of atomic many industrial and institutional partners. A special thanks to clock anomalies", Metrologia 45, 6, (2008) S127-S133 the European Space Agency, Thales Alenia Space France, the [14] P. Tavella, “Statistical and mathematical tools for atomic clocks”, Metrologia 45, 6, (2008) S183-S192 colleagues of the European UTC labs: NPL, OP, ORB, PTB, [15] I. Sesia, P. Tavella, "Estimating the Allan Variance in Presence of ROA, and SP; the PTF team and the Galileo Control Center Long Periods of Missing Data and Outliers”, Metrologia 45, 6, operators for their support. (2008) S134-S142 [16] L. Galleani, P. Tavella, “Dynamic Allan variance” , IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, UFFC, vol. 56, no. 3, March 2009, pp450-464 [17] I. Sesia, L. Galleani, P.Tavella “Application of the Dynamic Allan Variance for the Characterization of Space Clock Behavior” IEEE

Transactions on Aerospace and Electronic Systems, Vol. 47, N. 2 [22] P.DeDefraigne, W. Aerts, G. Cerretto, G. Signorile, E. Cantoni, I. April 2011 Sesia, P. Tavella, A. Cernigliaro, A. Samperi, J.M. Sleewaegen, [18] L.Galleani,P.Tavella, "Detection of Atomic Clock Frequency "Advances on the use of Galileo signals in time metrology: Jumps with the Kalman Filter," IEEE Transactions on Ultrasonics, calibrated time transfer and estimation of UTC and GGTO using a Ferroelectrics, and Frequency Control UFFC March 2012. vol. 59, combined commercial GPS-Galileo receiver", ION PTTI, no. 3, p. 504-509, March 2012. December 2013, Bellevue, Washington [19] J.Hahn, "Galileo IOV Timing Results", ICG-08 2013, November [23] R. Píriz, A. Mozo, G. Tobías, P. Tavella, I. Sesia, G. Cerretto, 2013, Dubai, UAE "GNSS interoperability: Offset Between Reference Time Scales and Timing Biases, Metrologia 45, 6, (2008) S87-S102 [20] A.Cernigliaro, I.Sesia, “Satellite Clocks Characterization and Monitoring for Global Navigation Satellite Systems”, Proc of the [24] G.Cerretto, N.Guyennon, I.Sesia, P.Tavella, F.Gonzalez, J. Hahn, V.Fernandez, V.Mozo, “Evaluation of the OD&TS System Noise XXX URSI General Assembly and Scientific Symposium of nd International Union of Radio Science, August 2011, Istanbul, in the Galileo GIOVE Mission”, Proceedings of the 22 European Turkey Frequency and Time Forum, Toulouse (France), April 2008 [21] http://www.esa.int/Our_Activities/Navigation/Galileo_starts_to_tel . l_UTC_the_world_s_time