Modelling GPS Observables for Time Transfer Marek Ziebart Department of Geomatic Engineering University College London
Presentation structure Overview of GPS Time frames in GPS Introduction to GPS observables Modelling the observables Applications, accuracies, issues
GPS segments I: Space Segment ~30 satellites in orbit, six orbital planes, 55 inclination Orbit near circular, 26,000km semi-major axis Satellites move ~4kms -1, orbital period ~12 hours On board rubidium and caesium clocks Thermally stabilised bus Several generations of satellite on orbit
GPS segments II: Control Segment Limited network (5) of tracking and uplink stations (about to grow) Ensemble of atomic clocks at Master Control Station, Colorado Springs -> defines system time System time steered to remain an integer number of seconds offset from UTC Control segment determines predicted satellite trajectories and clock behaviour
GPS Segments III: User Segment
Satellite clocks Freely running drift steadily from system time, not steered Satellite oscillators detuned before launch to offset first order relativistic time dilation Satellite clock offset and drift from system time estimated by control segment and made available to user via satellite navigation message
General context of problem In positioning getting accurate absolute position is generally difficult Where is my receiver?
Relative positioning is far easier Measure multiple ranges to satellites Rover at unknown point Compute relative position vector Reference receiver placed at known point New point for which We need to calculate position Point with known co-ordinates
Determining absolute time from GPS is complex What is the UTC time by my clock?
Determining relative time from GPS is more straightforward My clock is 4 seconds slower than yours.
Time frames in GPS GPST (Colorado) UT1/UTC link between orbital reference frame (inertial space) and Earth-centred, Earth-fixed terrestrial frames TDB, TDT, GAST, TAI User time frame (timing laboratories in this case)
GPS signals 19 cm and 24 cm (L-band) wavelength carriers Pseudo-random noise (PRN) codes (30 m and 300 m chip spacing) Navigation Message predicted orbit, predicted clock behaviour, offset from UTC
GPS observables L1 and L2 carrier phase L1 and L2 Doppler count C/A, P1 and P2 pseudoranges L2C as of yesterday
Basic mathematical model: Pseudoranges P i = x x + r s r s P δ t c δt c + I + T + M + ε satellite clock offset from system time receiver clock offset from system time satellite antenna phase centre position ground antenna phase centre position measured pseudorange
Basic mathematical model: signal errors ionospheric delay multipath P i = x x + t c t c I T M r s δ + δ + + + + ε r s dry and wet tropospheric delay Target parameter noise
Basic mathematical model: Phase Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ phase multipath phase ambiguity measured phase
Modelling I: Receiver position ε δ δ + + + + + + = Φ Φ M T I N c t c t x x s r s r i Position of antenna phase centre at the instant of measurement Ideally treated as a known quantity
Instantaneous antenna phase centre position Elevation dependent phase centre variations (L1, L2) Mean L1/L2 phase centre offset L2 phase centre offset L1 phase centre offset antenna reference point (ARP) antenna pre-amplifier vector from ground mark to ARP ground mark Earth and ocean tide correction Published mean coordinates of station
Kinematic Earth tide variation, Lhasa: 1 day, perihelion 0.2 0.1 Earth tide correction (m) 0-0.1-0.2-0.3-0.4-0.5-0.6 0.8m range (on this day) -0.7-0.8 time (MJD)
Modelling II: Transmitter position Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ Position of satellite transmit antenna phase centre at the instant of signal transmission Known to varying degrees of accuracy Requires knowledge of orbit and satellite attitude (Block IIA only)
Orbit accuracies and latency Broadcast ephemeris: 2-5m accuracy, real-time access IGS predicted : 10 20 cm accuracy, quasi-real-time IGS rapid : ~5 cm radial accuracy, on the day IGS final : <5 cm radial accuracy, two weeks latency Note: 30 cm orbit error => 1 ns time error
Modelling III: Satellite Clock Error Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ Crude estimate of satellite clock error broadcast in navigation message Precise applications require estimation of residual offset Parameterised as offset, drift and ageing Can be eliminated by common view differencing
7.0 6.0 5.0 Difference between OSPF broadcast satellite clock offset and final IGS solution satellite clock error (m) 4.0 3.0 2.0 1.0 0.0 Circa 3 metres (10 ns) -1.0 1 2 3 5 6 7 8 9 10 11 13 14 16 18 20 21 25 26 27 28 29 30 31 PRN number
Common view differencing Φ 1 - Φ 2 eliminates satellite clock errors Φ 1 Φ 2
Modelling IV: Phase ambiguity Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ Integer number of cycles biasing the phase range observable Constant parameter for a phase connected arc Can be determined from a single epoch of data
Atmospheric models Ionosphere: Troposphere: from 9km 100km depth to at 1000km poles, 16km Driving at parameter: equator driving electron parameters: content, Neutral solar gas UV density, flux, solar water cycle, vapour coupling distribution with solar magnetic field
Modelling IV: Ionospheric refraction Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ Dispersive (frequency dependent) effect for microwaves Highly variable, chaotic delay (2m -> 80m) Can be eliminated (to first order) using dual frequency ionosphere-free combination
Modelling IV: Tropospheric delays Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ Dry and wet components Total delay ~2-3m Dry component modelled successfully using surface meteorological data Wet component driven by water vapour density most approaches estimate the delay
Wet component Tropospheric delay estimation Most up to date approaches estimate Zenith delay and two horizontal gradient parameters Slowly varying effect can be modelled using a random walk in time Still a major source of uncertainty ~5-10 cm
Modelling V: Multipath Φ i = x r x s + δ t c δt c + N + I + T + M + ε r s Φ Signal degradation due to the physical environment of the receiver and the transmitter High and low frequency noise/bias Hot area of research
Approaches to dealing with Multipath Hardware mitigation: ground planes, choke rings, antenna gain pattern Correlator windows Smart choice of antenna siting Modelling: sidereal filter (GPS constellation geometry repeats daily) Ray tracing Statistical approaches
Data pre-processing: pseudoranges Read RINEX Gross error check Ionosphere-free combination Cycle slip detection Tropospheric delay Multipath detection Phase Smoothing Test Data
Pseudo-range time estimation: single station -4.0E-08-5.0E-08 Receiver Clock Offset (Seconds) -6.0E-08-7.0E-08-8.0E-08-9.0E-08-1.0E-07 averaged clock offset (OSPF) averaged clock offset (SP3) IGS clock (truth model) ~3m -1.1E-07-1.2E-07 0 500 1000 1500 2000 2500 3000 3500 Seconds of Hour
Project tracking network = atomic clocks (H-masers, caesium, rubidium) = conventional quartz clocks
Pseudo-range time estimation: network approach -8.80E-08-9.00E-08-9.20E-08 estimated clock offset (seconds) -9.40E-08-9.60E-08-9.80E-08-1.00E-07-1.02E-07-1.04E-07-1.06E-07-1.08E-07 averaged clock offset IGS receiver clock (truth model) Global network solution, GLPS, clocks only, weighted RMS scatter: ~0.2m or ~1ns -1.10E-07 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 Time (minutes)
Receiver clock offset (Goldstone H-maser, USA): Carrier phase solution 0.0E+00-1.0E-08 Estimated accuracy: ~0.1ns -2.0E-08-3.0E-08-4.0E-08-5.0E-08-6.0E-08-7.0E-08-8.0E-08-9.0E-08 0 2400 4800 12300 14700 17100 19500 21900 24300 26700 29100 31500 33900 36300 38700 41100 43500 45900 48300 50700 53100 55500 57900 60300 62700 65100 67500 69900 79800 82200 84600 receiver clock offset (seconds) Time (seconds of day)
Conclusions and Outlook Many ways of using GPS observables in time transfer 1 ns real-time accuracy with pseudo-ranges and networks is feasible 0.1 ns post-processed results are routine (two week latency) 0.1 ns daily solution may be feasible now Careful modelling of observables underpins all applications many ways to fail, many ways to improve performance