Satellite Clock Error and Orbital Solution Error Estimation for Precise Navigation Applications ()
In this paper, the satellite signal time of transmission is precisely estimated by considering the clock correction parameters like satellite clock bias (a0), satellite clock drift (a1), satellite clock drift rate (a2), transmitted as part of the navigation message, signal emission time at the antenna and the pseudorange measured between the satellite and the receiver.
The orbital solution of each satellite is estimated using the broadcast ephemerides, which differ from their true satellite positions as shown in Figure 1. This deviation is due to uncertainty in the gravitational model, limited accuracy of the orbit representations and inadequately modelled surface forces on the satellites like solar radiation pressure and atmospheric drag. In this paper, the orbital solution errors are estimated by comparing the broadcast ephemerides with the precise ephemerides obtained by JPL, USA. Here for error calculation, the JPL estimations are taken as the benchmark assuming that these estimates
Figure 1. Satellite ephemeris errors.
are of high precision. The errors are estimated and analysed for ephemeris data collected on 11th March 2011 at Department of Electronics and Communication Engineering, Andhra University College of Engineering, Vishakhapatnam.
2. Estimation of Satellite Clock Error and Orbital Solution Error
The GPS receiver uses the same PRN codes which are transmitted by the 32 satellites for determining the distance between each satellite and receiver. This distance is also known as pseudorange. For precise navigation solution the computed pseudorange needs to be corrected for the errors like satellite clock error, tropospheric error, multipath errors, ionospheric errors etc. [1]. Among these the satellite clock error has a major impact on the pseudorange and is given by the following Equation (1)
(1)
where, 
= Measured range (meters); 
= True range (meters); 
= Satellite clock error (sec); 
= Velocity of light = 3 ´ 108 meters/sec; from Equation (1) it is evident that satellite clock error of 1 microsecond will lead to 300 meters error in the pseudorange [2].
2.1. Satellite Clock Error
All satellites contain atomic clocks that control all onboard timing operations, including broadcast signal generation. Although, these clocks are highly stable still they lack of perfect synchronization between the timing of the satellite broadcast signals and GPS system time. Satellite clock correction terms (
,
and
) and time of clock
, account for this lack of synchronization. All these parameters are obtained from navigation file. The on board clock stability is about 1 to 2 parts in 1013 over a period of one day [3]. These errors are common to all users observing the same satellite. The satellite clock error is caused by the satellite oscillator not being synchronised to true time (GPS time).
The deviation of a particular clock from GPS system time is modeled as a quadratic function of time. The parameters of this model are estimated, uploaded to the satellite, and are broadcasted. The coefficients
, 
and 
in the below equation are called the bias, drift, and aging parameters of the satellite clock [4]. The estimation of satellite clock errors are as given below.
The satellite time of transmission is computed as
(2)
where, 
= satellite time of transmission (sec); 
= Signal reception time at receiver (sec); P = Pseudorange (meters); c = speed of light (meters/sec).
The time at which the signal received from the satellite and the pseudorange are read from the broadcast ephemerides and observation data files respectively. The term 
in Equation (2) gives the signal propagation time. This propagation time is subtracted from the signal reception time computed at the receiver and the residual time will be the instant at which the signal started from the satellite.
The satellite clock offset 
is calculated as,
(3)
where, 
= Satellite time of transmission (sec); 
= Time of clock (sec).
The satellite time of transmission is calculated using Equation (2) and time of clock is available in navigation data.
The satellite clock error 
is modelled using the second order polynomial as given below. The correction parameters bias, drift and aging are available in navigation data.
(4)
where, 
= clock bias (sec); 
= clock drift (sec/sec); 
= frequency drift (i.e. aging) (sec/sec2); 
= satellite clock offset (sec).
The GPS time of transmission corresponds to the signal emission time from the satellite and is computed as below,
(5)
where, 
= Corrected GPS time of transmission (sec); 
= satellite time of transmission (sec); 
= Satellite clock error (sec).
Finally the satellite position is estimated using the corrected time of transmission 
and the broadcast orbital parameters.
2.2. Broadcast Orbital Solution Error
The GPS navigation message file which contains the broadcast ephemerides gives the Keplerian parameters needed to compute the coordinates and clock correction for each satellite. This ephemeris data also contains the time information needed to correct that satellite’s clock to match with the GPS time. Error in satellite clock and orbital parameters results in estimation of incorrect orbital solution.
2.3. JPL Orbital Solution Interpolation
The Lagrange interpolation algorithm is used to find the intermediate values for the points (known as the data nodes) at a given interval. The precise satellite positions estimated by JPL for an epoch interval of 15 minutes where as broadcast satellite positions are estimated for an epoch interval of 15 seconds. The precise 3D satellite positions, Sk at different regular intervals of time are used to interpolate the satellite positions in between the known intervals by using the Lagrange interpolation algorithm [5],
(6)
where, 
are Lagrange interpolating polynomials defined as,
(7)
where, x = Lagrange multipliers; n = Number of intervals.
The above Lagrange interpolating polynomial for time interval t is written as the nth polynomial
These calculated polynomials are substituted in equation (6) to get the interpolated precise satellite positions.
3. Results and Discussion
The satellite clock errors and orbital solution errors are estimated and the analysis of the error which are supported by the relevant graphs and the tables are presented in this paper.
The broadcast orbital solution error and satellite clock error are estimated using the data collected by the receiver from 02:00 hours to 02:39 hours with an epoch interval of 15 seconds. During this observation period of 39 minutes (total of 157 epochs) out of 31 satellites, 9 (02, 05, 12, 15, 18, 21, 25, 26, 29) satellites were visible. Though the errors are computed and analyzed for all the visible satellites, the results corresponding to a particular Satellite Vehicle Pseudo Random Noise (SVPRN 26) are presented in this paper. To compute the positional error of the satellite, the positions estimated by the receiver are compared with the one estimated by JPL, assuming the JPL estimates are of high precision.The precise ephemerides estimated by JPL are interpolated using the Lagrange interpolation algorithm to generate satellite position for SVPRN 26 for 157 epochs with an epoch interval of 15 seconds. Tables given in this paper details 10 epochs of data and figures are plotted for 157 epochs.
The broadcast ephemerides are collected on 11th March 2011 from the dual frequency GPS receivers located at Department of Electronics and Communication, Andhra University College of Engineering, Visakhapatnam (17.73˚N/ 83.319˚E). The precise ephemerides are collected from Jet Propulsion Laboratory (JPL), USA for the same day.
3.1. Satellite Clock Error
Table 1 details the satellite clock offset and satellite clock error of SVPRN26 for 10 epochs. Figure 2 and figure 3 shows these variations over 157 epoch and it is observed that the satellite clock offset is varied from −0.07262 to 2340 seconds and the corresponding clock error is varied from −0.00014192 to −0.00014194 seconds.
3.2. Orbital Solution Error
Table 2 details the error in each of the coordinates of the broadcast satellite position with that of the precise satellite position of SVPRN26 for 10 epochs. Figure 4 shows this variation over 157 epochs and it is observed that the satellite positions in the x-coordinate is deviated by 17.02 to −5.939 meters. Similarly, the error is varied from 2.973 to 63.89 meters and −235.8 to −218.1 meters for yand z-coordinates respectively.
Table 1. satellite clock offset and error of SVPRN 26.
Figure 2. Satellite clock offset for SVPRN 26.
Figure 3. Satellite clock error for SVPRN26.
Table 2. Orbital solution error for SVPRN 26.
4. Conclusions
This paper reveals the importance of the satellite clock error and the orbital solution error.
Figure 4. Orbital solution error for SVPRN 26.
Satellite clock error. Over 157 epochs of data analysed, it is observed that the satellite clock offset varied from −0.07262 to 2340 seconds and the corresponding clock error varied from −0.00014192 to −0.00014194 seconds. Hence the clock error estimation needs to be modeled for precision applications, (e.g. CAT I/II aircraft landings, missile navigation).
Broadcast Orbital Solution Errors.
The error in satellite positions inturn affects the accuracy of navigation solution. Broadcast satellite position varied from precise positions and the deviation in x-coordinate is 17.02 to −5.939 meters over 157 epochs. Similarly the error varied from 2.973 to 63.89 meters and −235.8 to −218.1 meters for yand z-coordinates respectively. These positional errors have significant effect on the critical applications, (e.g. studies of the crustal dynamics of the earth).
Acknowledgements
The work undertaken in this paper is supported by Ministry of Science and Technology, Department of Science and Technology (DST), Government of India, New Delhi, under Woman Scientist Scheme(WOS-A), Vide Saction letter No. SR/WOS-A/ET-04/2013.