Abstract
Limiting factors for the precise orbit determination (POD) of low-earth orbit (LEO) satellite using dual-frequency GPS are nowadays mainly encountered with the in-flight phase error modeling. The phase error is modeled as a systematic and a random component each depending on the direction of GPS signal reception. The systematic part and standard deviation of random part in phase error model are, respectively, estimated by bin-wise mean and standard deviation values of phase postfit residuals computed by orbit determination. By removing the systematic component and adjusting the weight of phase observation data according to standard deviation of random component, the orbit can be further improved by POD approach. The GRACE data of 1–31 January 2006 are processed, and three types of orbit solutions, POD without phase error model correction, POD with mean value correction of phase error model, and POD with phase error model correction, are obtained. The three-dimensional (3D) orbit improvements derived from phase error model correction are 0.0153 m for GRACE A and 0.0131 m for GRACE B, and the 3D influences arisen from random part of phase error model are 0.0068 m and 0.0075 m for GRACE A and GRACE B, respectively. Thus the random part of phase error model cannot be neglected for POD. It is also demonstrated by phase postfit residual analysis, orbit comparison with JPL precise science orbit, and orbit validation with KBR data that the results derived from POD with phase error model correction are better than another two types of orbit solutions generated in this paper.
1. Introduction
Since spaceborne dual-frequency Global Positioning System (GPS) receiver was successfully applied in CHAllenging Minisatellite Payload (CHAMP) mission [1, 2], more and more low-earth orbit (LEO) satellite missions have been equipped with dual-frequency GPS receiver for precise navigation, such as Gravity Recovery And Climate Experiment (GRACE) mission [3], Gravity field and steady-state Ocean Circulation Explorer (GOCE) mission [4], and TanDEM-X mission [5]. Due to the characteristic of good continuity, high precision, and low cost, the spaceborne dual-frequency GPS receiver has become a primary instrument for precise orbit determination (POD) of LEO satellite. By the processing of GPS observation data, 3dimensional (3D) in-flight position information of satellite can be obtained with centimeter-level precision [6–8].
At present, the reduced dynamic orbit determination method [8–10] is widely used for LEO POD based on dual-frequency GPS, which combines the geometric strength of the GPS observations with the orbit dynamical model constraints. In this way, the available accuracy of the GPS measurements may be fully exploited without sacrificing the robustness offered by dynamic orbit determination techniques. Thanks to numerous improvements in the final GPS orbit and clock products provided by the International GNSS Service (IGS) [11], in the orbit dynamical models of LEO satellite, and in the approaches of data processing and orbit determination, the accuracy of LEO POD based on dual-frequency GPS can be steadily improved.
Note that the noise level of GPS carrier phase measurement of GPS receiver can reach 1-2 mm, which is much more accurate than GPS pseudo code observation. Therefore, the quality of phase observation will directly determine the accuracy of LEO POD. However, there exist different kinds of errors in the phase observation, such as phase measurement error, receiver antenna phase center location error, and near-field multipath. The phase error belongs to mixed error formed by random error and systematic error, and the error character is complex. Therefore, the researches on phase error modeling and correction provide a valid way to further improve the precision of orbit determination. The residual approach [7, 12–15] is a widely used method for in-flight calibrations of phase errors. In this approach, calibrations are derived as bin-wise mean values from GPS carrier phase postfit residuals obtained from orbit determination. The JASON-1 orbits [12] have already been successfully improved by this approach. Meanwhile, the Jet Propulsion Laboratory (JPL) [13] applies this approach to the GRACE satellites, which serves as primary data source to derive phase center variations (PCVs) for the GPS transmitter antennas. In addition, Montenbruck et al. [14], Jäggi et al. [7], and Bock et al. [15] use this approach to calibrate the PCVs of GPS receiver antennas onboard GRACE, TerraSAR-X, and GOCE, respectively. In this approach, only bin-wise mean values of phase postfit residuals are considered for the in-flight calibrations, and the systematic part in phase error can be partly calibrated. However, the random part is less considered, that is, the bin-wise standard deviation information of phase postfit residuals is neglected.
Subsequently, we mainly focus on the whole modeling of phase error and assess its impact on GRACE POD. In this paper, the phase error is modeled as a systematic and a random component each depending on the direction of signal reception. The systematic part and the standard deviation of random part in the direction of signal reception are, respectively, estimated. The final LEO orbit can be determined by removing the systematic part and adjusting the weight of phase data. For this investigation, the GRACE satellites are selected. The GRACE mission [3], launched on March 17, 2002, is a joint partnership between the National Aeronautics and Space Administration (NASA) in the United States and Deutsches Zentrum für Luft-und Raumfahrt (DLR) in Germany. It consists of two identical formation flying spacecraft (GRACE A and GRACE B) in a near polar, near circular orbit with an initial altitude of about 500 km. The spacecraft has a nominal separation of 220 km. The key science instruments onboard both spacecraft include a BlackJack GPS receiver, a SuperSTAR accelerometer, a star tracker, a K-band ranging (KBR) system, and a satellite laser ranging (SLR) retroreflector. The BlackJack GPS receiver exhibits a representative noise level of 1 mm for L1 and L2 carrier phase measurements.
Section 2 focuses on the in-flight phase error modeling, parameter estimation of error model, and POD with phase error correction. In Section 3, the phase error model is applied to the real GPS observation data of GRACE and its impact on POD is discussed. Section 4 shows the conclusions.
2. In-Flight GPS Phase Error Modeling
2.1. Observation Equation
In order to eliminate the first order ionosphere delay, dual-frequency ionosphere-free combination observations are always adopted. For the pseudocode and carrier phase observations, the ionosphere-free combination yields where subscript “IF” denotes ionosphere-free combination, subscripts 1 and 2 denote different frequencies, is pseudo code ionosphere-free combination, is phase ionosphere-free combination, is the carrier frequency, is the real signal traveling time from GPS satellite to LEO satellite which can be obtained by iterative calculation, is geometric distance between the mass center position of GPS satellite and LEO satellite at signal transmission epoch and signal reception epoch , respectively, and is light velocity. δt is the clock offset of LEO satellite, is the ambiguity of phase ionosphere-free combination, and is a series of corrections which can be denoted as where is the clock correction of GPS satellite at epoch , is the relativity correction of GPS satellite j, is the phase center offset correction of GPS satellite j, and is the ionosphere-free phase center correction of LEO satellite. These corrections have been studied by other scholars [16] and are not discussed in this paper.
For convenience, has to be linearized as where where is the mass center position of GPS satellite in conventional inertial reference frame (CIRF) at epoch , is the approximate mass center position of LEO satellite in CIRF, and is a line of sight (LOS) vector.
2.2. Phase Error Modeling
For the convenience of description, the definition of Antenna-Fixed Coordinate System (AFCS) [14] is given at first. The origin O is mechanical antenna reference point (ARP). The positive -axis coincides with the mechanical symmetry axis and points along the boresight direction. The -axis and -axis point from the mechanical ARP into the respective directions, which depend on the specific mounting of the antennas. Take AFCS of GRACE satellites for instance. The -axis of AFCS coincides with -axis of Satellite Body Coordinate System (SBCS). The -axis of AFCS is in opposite direction with -axis of SBCS and -axis of AFCS completes a right-handed coordinate system. In AFCS, the azimuth angle of a vector is defined as an angle between the projection of in XOY-plane and positive -axis and is counted in a counter-clockwise sense from the -axis to the -axis. The elevation angle of a vector is defined as an angle between and XOY-plane.
As the in-flight GPS phase error consists of GPS receiver antenna phase center location error, phase observation error, near-field multipath, and all other nonmodeled phase errors, it can be composed of systematic part and random part, which both depend on the direction of GPS signal reception. Assuming that the LOS vector of GPS reception signal in CIRF is , the phase error on Li band along this direction can be modeled as where and are the systematic part and random part of phase, respectively, is a random variable from normal distribution with mean value of 0 and standard deviation of , and and are independent if .
From (2.6), the ionosphere-free phase error model can be expressed as
2.3. Parameter Estimation of Phase Error Model
Assuming that the GPS ionosphere-free carrier phase observation value at epoch is , LOS vector is , and the ionosphere-free phase observation model value from (2.2) is , then the ionosphere-free phase error value can be expressed as
From (2.8), the ionosphere-free phase error can be directly estimated by the postfit residuals of orbit determination. The LEO orbits can be directly obtained by reduced dynamic orbit determination approach. As the phase error depends on the direction of GPS reception signal, that is, the azimuth angle and elevation angle of the LOS vector in AFCS, phase error can be estimated by azimuth/elevation bins. In this study, the phase postfit residuals are sorted in azimuth/elevation bins of . The mean value and standard deviation of all the phase postfit residuals fallen into the region of are obtained. The values of and are both selected as 5° in this paper. If the azimuth angle and elevation angle of the LOS vector are fallen into the region of , is selected as the estimation of phase systematic part and is the estimation of standard deviation of phase random part . They are denoted as follows:
2.4. POD with Phase Error Correction
From (2.8) and (2.9), we can get
According to (2.10), the orbit parameters can be improved by reduced dynamic orbit determination approach in a second step after removing the systematic part and adjusting the weight of phase observation data according to .
Assuming that the phase outliers have been removed, the elevation cutoff angle is , the azimuth angle and elevation angle of phase observation data are and , respectively, the standard deviation of all the phase postfit residuals is , and the weights of phase observation data are set as follows.(1)If , the weight of phase observation data is 0.(2)When , and are fallen into the region of , if > 0, the weight of phase observation data is , if = 0, the weight of phase observation data is 0.
The flow of data processing for LEO POD is shown in Figure 1.
3. Numerical Analysis for GRACE Satellites
3.1. Data Sets and Processing Strategies
The data sets used here include GPS observation data (GPS1B), spacecraft attitude data (SCA1B), KBR data (KBR1B), JPL precise science orbit data (GNV1B) of GRACE A and GRACE B [17] from GeoForschungsZentrum (GFZ), and the final GPS orbits and the 30 s high-rate satellite clock corrections from the Center for Orbit Determination in Europe (CODE) [18]. The data cover a period of 31 days from January 1 to January 31 of 2006. The phase center offsets of GPS receiver antennas onboard GRACE satellites [19] in respective SBCS are listed in Table 1.
The LEO orbit determination is implemented in the separate software tools as part of the NUDT Orbit Determination Software 1.0. The GPS observation data processing consists of GPS observation data preprocessing [16], reduced dynamic orbit determination with medium precision [16, 19], GPS observation data editing [16, 19], and reduced dynamic orbit determination with high precision. Both reduced dynamic orbit determinations with medium and high precision make use of zero-difference (ZD) reduced dynamic batch Least-SQuares (LSQ) approach. In this approach, the individual spacecraft positions at each measurement epoch are replaced by the spacecraft trajectory model in (2.1) and (2.2). We make use of known physical models of the spacecraft motion to constrain the resulting position estimates and introduce three empirical acceleration components to absorb all other nonmodeled perturbation. As the atmospheric density and solar activity are difficult to model accurately, atmospheric drag coefficient and solar radiation pressure coefficient in the models of atmospheric drag and solar radiation pressure are also estimated. The orbit dynamical models and reference frames used for reduced dynamic orbit determination are listed in Table 2.
Three empirical acceleration components in the radial (R), transverse (T) and normal (N) direction are denoted as [25] where is satellite latitude, , , ,, , and are empirical acceleration coefficients, and , , and are unit vectors in directions of R, T, and N, which are three axes of RTN coordinate system.
All the parameters are estimated by batch LSQ approach [16, 19]. The parameters estimated in orbit determination are listed in Table 3. Orbit determination process is typically conducted in 24-hour data batches, and Adams-Cowell multistep integration method [26] is used for orbit integration.
3.2. Parameter Estimation Results Adopting Phase Error Model
At first, all the GPS observation data of GRACE A and GRACE B are processed to obtain the orbits using POD without phase error model correction. The mean values and standard deviations of phase postfit residuals are stored with a resolution of . The mean value patterns and standard deviation patterns of GRACE satellites are shown in Figure 2. It is shown by Figure 2 that the mean value patterns and standard deviation patterns of both satellites exhibit the similar distributions. In addition, it is shown by mean value patterns and standard deviation patterns from different periods that the distributions of mean value patterns and standard deviation patterns keep relatively steady and can be assumed to keep constant with time.
(a) Mean value pattern of phase postfit residuals of GRACE A
(b) Mean value pattern of phase postfit residuals of GRACE B
(c) Standard deviation pattern of phase postfit residuals of GRACE A
(d) Standard deviation pattern of phase postfit residuals of GRACE B
3.3. Impact of Phase Error Modeling on GRACE POD
In this section, besides the GRACE orbit solutions computed by POD without phase error model correction as mentioned in Section 3.2, another two types of GRACE orbit solutions derived from POD with mean value correction of phase error model and POD with phase error model correction are also analyzed. The mean value correction of phase error model is just the residual approach. Phase postfit residual analysis, orbit comparison, and orbit validation with KBR data are used to analyze the impact of phase error modeling on GRACE POD.
3.3.1. Phase Postfit Residual Analysis
Phase postfit residual analysis can be used to measure the consistency of the applied models with the GPS observation data. The Root Mean Square (RMS) errors of ionosphere-free phase postfit residuals of POD are shown in Table 4.
It is found by statistic results that the mean values of phase postfit residuals are close to 0 m (see Figure 3). From Table 4, the RMS values of phase postfit residuals derived from POD with phase error model correction are 9.96 mm for GRACE A and 9.28 mm for GRACE B, which are less than those computed by POD without phase error model correction and POD with only mean value correction of phase error model. It is shown by these results that the applied models are in better consistency with GPS observation data using POD with phase error model correction.
(a) Ionosphere-free phase postfit residuals of POD without phase error model correction for GRACE A on January 15, 2006
(b) Ionosphere-free phase postfit residuals of POD without phase error model correction for GRACE B on January 15, 2006
(c) Ionosphere-free phase postfit residuals of POD with mean value correction of phase error model for GRACE A on January 15, 2006
(d) Ionosphere-free phase postfit residuals of POD with mean value correction of phase error model for GRACE B on January 15, 2006
(e) Ionosphere-free phase postfit residuals of POD with phase error model correction for GRACE A on January 15, 2006
(f) Ionosphere-free phase postfit residuals of POD with phase error model correction for GRACE B on January 15, 2006
3.3.2. Orbit Comparison
In an effort to obtain some information on the orbit improvement arisen from phase error model correction and the quality of GRACE orbits computed in this paper, internal and external quality metrics are used for orbit comparison. The sampling intervals of orbits are 30 s. In internal quality metric, the orbits generated by POD with phase error model correction (identifier “POD + PMC”) are compared with the orbits obtained by POD without phase error model correction (identifier “POD + NOPMC”) and the orbits obtained by POD with mean value correction of phase error model (identifier “POD + MVC”), respectively, which can be used to obtain the improvement of POD with phase model correction. In external quality metric, the orbits computed in this paper are compared with the JPL precise science orbits (identifier “JPL”), which shows the consistency with JPL precise science orbits. JPL precise science orbits represent one of highest precision of GRACE orbit solutions and are created by processing zero-difference ionosphere-free pseudo code and carrier phase data with the GIPSY-OASIS software package, which are distributed along with the GRACE GPS data as a part of GRACE Level 1B product. The RMS values of orbit comparisons computed in R, T, and N direction and 3D (see Figure 4) are listed in Table 5.
(a) POD + NOPMC – POD + PMC (GRACE A)
(b) POD + NOPMC – POD + PMC (GRACE B)
(c) POD + MVC – POD + PMC (GRACE A)
(d) POD + MVC – POD + PMC (GRACE B)
(e) POD + NOPMC – JPL (GRACE A)
(f) POD + NOPMC – JPL (GRACE B)
(g) POD + MVC – JPL (GRACE A)
(h) POD + MVC – JPL (GRACE B)
(i) POD + PMC – JPL (GRACE A)
(j) POD + PMC – JPL (GRACE B)
It is shown by the results of internal quality metric that the 3D orbit improvements arisen from phase error model correction for GRACE A and GRACE B are 0.0153 m and 0.0131 m, respectively. The most improvements are all in T direction, 0.0118 m and 0.0097 m, which is because T direction corresponds to flight direction, the phase error perpendicular to the flight direction are to some extent absorbed by the carrier phase ambiguities and clock offset, but the phase error along the flight direction will almost be absorbed by orbit parameters in orbit determination. In addition, the comparisons between the orbits derived from POD with mean value correction of phase error model and POD with phase error model correction reflect the influence arisen from random part of phase error model, and the 3D RMS values for GRACE A and GRACE B are 0.0068 m and 0.0075 m, respectively. Therefore, the random part of phase error model cannot be neglected for POD.
In external quality metric, the 3D RMS values of orbit comparison results are reduced from 0.0553 m to 0.0527 m for GRACE A and from 0.0515 m to 0.0502 m for GRACE B by POD with phase error model correction. These are slightly better than the comparison results between the orbits derived from POD with mean value correction of phase error model and JPL precise science orbits. It demonstrates that orbits computed by POD with phase error model correction have the higher consistency with JPL precise science orbits.
3.3.3. Orbit Validation with KBR Data
KBR system is one of the key scientific instruments onboard the GRACE satellites, which measures the one-way range change between the twin GRACE satellites with a precision of about 10 μm for KBR range at a 5-second data interval. Due to the high precision of KBR data, the relative position accuracy of the GRACE satellites can be validated. The relative positions computed by POD without phase error model correction, POD with mean value correction of phase error model, POD with phase error model correction, and JPL precise scientific orbits are validated by KBR data, respectively, (see Figure 5), and the average standard deviations of KBR comparison residuals are shown in Table 6.
From Figure 5 and Table 6, we can see that the average K-band standard deviation of the relative positions computed by POD results with phase error model correction is the least of all three types of GRACE orbit solutions generated in this paper. In addition, it is also less than the average K-band standard deviation of the relative positions computed by JPL precise science orbits. These results show that phase error model correction could remove the phase errors of both satellites and relative position accuracy could be obviously improved.
4. Conclusions
The phase error model has been set up and its model parameters have been estimated. The impact of phase error modeling on GRACE POD has been analyzed. The following conclusions are drawn from above results and discussions.
The whole phase error model is set up. The phase error is modeled as a systematic and a random component each depending on azimuth and elevation of signal reception in AFCS. The systematic part and standard deviation of random part in phase error model can both be estimated by the bin-wise mean and standard deviation values of phase postfit residuals obtained by direct orbit determination. By removing the systematic part and adjusting the weight of phase observation data according to standard deviation of random part, it leads to improved orbit parameters in a POD procedure.
The 31-day GRACE data of January 2006 are processed, and three types of GRACE orbit solutions by POD without phase error model correction, POD with mean value correction of phase error model, and POD with phase error model correction are generated. The impact of phase error modeling on LEO POD is analyzed by a number of tests, including phase postfit residual analysis, orbit comparison, and orbit validation with KBR data, respectively. At first, the phase postfit residual analysis shows that the RMS values of ionosphere-free phase postfit residuals of GRACE satellites derived from POD with phase error model correction are the least in the three types GRACE orbit solutions, and the applied models are in the best consistency with GPS observation data. Secondly, two quality metrics, the internal and the external, are used for orbit comparisons. In the internal metric, the 3D orbit improvements derived from phase error model correction for GRACE A and GRACE B are 0.0153 m and 0.0131 m, and the 3D influences arisen from random part of phase error model are 0.0068 m and 0.0075 m for GRACE A and GRACE B, respectively. As a result, the random part of phase error model cannot be neglected for POD. From the external quality metric, the 3D RMS values of orbit comparison residuals between the orbits computed by POD with phase error model correction and JPL precise science orbits are 0.0527 m for GRACE A and 0.0502 m for GRACE B, which are also better than another two types of orbit solutions. It shows that the orbit solutions of GRACE satellites with phase error model correction are in higher consistency with JPL precise science orbits. At last, the average K-band standard deviation of the relative positions computed by POD results with phase error model correction is the least of all three types of GRACE orbit solutions generated in this paper, and it is also better than that derived from JPL precise science orbits.
As the aforementioned, the phase error modeling is helpful to improve precision of orbit determination, and the phase error in GPS phase observation data can be reduced. This research will have some theoretical and engineering significance on phase error correction of LEO POD.
Acknowledgments
The authors are grateful to GeoForschungsZentrum (GFZ) for providing GPS observation data, attitude data, KBR data, and precise science orbit data of GRACE mission. The authors are grateful to the Center for Orbit Determination in Europe (CODE) for providing GPS orbit solutions and clock corrections. This study is supported by the National Natural Science Foundation of China (Grant no. 61002033 and no. 60902089) and Open Research Fund of State Key Laboratory of Astronautic Dynamics of China (Grant no. 2011ADL-DW0103).