Next Article in Journal
Impact of Winter Snowfall on Vegetation Greenness in Central Asia
Next Article in Special Issue
Potential Contributors to Common Mode Error in Array GPS Displacement Fields in Taiwan Island
Previous Article in Journal
Estimation of Maize Photosynthesis Traits Using Hyperspectral Lidar Backscattered Intensity
Previous Article in Special Issue
Reduced-Dynamic Precise Orbit Determination of Haiyang-2B Altimetry Satellite Using a Refined Empirical Acceleration Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GRACE-FO Antenna Phase Center Modeling and Precise Orbit Determination with Single Receiver Ambiguity Resolution

1
Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Space Star Technology Co., Ltd., Beijing 100194, China
4
GNSS Research Center, Wuhan University, 129 Luoyu Road, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(21), 4204; https://doi.org/10.3390/rs13214204
Submission received: 20 August 2021 / Revised: 10 October 2021 / Accepted: 14 October 2021 / Published: 20 October 2021
(This article belongs to the Special Issue BDS/GNSS for Earth Observation)

Abstract

:
Precise knowledge of the phase center location of the global navigation satellite system (GNSS) antenna is a prerequisite for precise orbit determination (POD) of the low Earth orbit (LEO) satellite. The phase center offset (PCO) and phase center variation (PCV) values for the LEO antenna obtained from ground calibration cannot reflect the error sources encountered in the actual spacecraft environment. PCV corrections are estimated by ionosphere free (IF) carrier phase post-fit residuals of reduced dynamic orbit determination. Ambiguity resolution (AR) plays a crucial role in achieving the best orbit accuracy. The single receiver AR concept is realized using wide-lane (WL) and narrow-lane (NL) bias products. Single difference (SD) observations between satellites are applied to remove the receiver dependent phase bias. SD AR and traditional double difference (DD) AR methods are applied to fix the ambiguities. The recovered SD and DD IF ambiguities are taken as pseudo-observations to constrain the undifferenced IF ambiguity parameters in the POD process. The LEO orbits based on float ambiguity (FA), SD, AR, and DD AR are investigated. One year’s data collected by the Gravity Recovery And Climate Experiment Follow-On (GRACE-FO) mission and GPS precise products provided by the Center for Orbit Determination in Europe (CODE) were analyzed. Precise orbit generated by the Jet Propulsion Laboratory (JPL), independent satellite laser ranging (SLR), and K-band ranging (KBR) measurements were utilized to assess the orbit accuracy. More than 98% of SD WL and 95% of SD NL ambiguities are fixed, which confirms the good quality of the bias products and correctness of the SD AR method. With PCV corrections, the average phase residuals of DD and SD AR solutions are 0.13 and 0.41 mm, which indicates improved consistency between applied models and observations. Compared with JPL’s orbit, the SD AR orbits achieve the accuracy of 6.0, 6.2, and 5.1 mm in along-track, cross-track, and radial directions. The SD AR solutions show an average improvement of 18.3% related to the FA orbits while 6.3% is gained by the DD AR approach. The root mean squares (RMSs) of SLR residuals for FA, DD AR, and SD AR solutions are 11.5, 10.2, and 9.6 mm, which validate the positive effect of AR on POD. Standard deviation (STD) of KBR residuals for SD AR orbits is 1.8 mm while 0.9 mm is achieved by the DD AR method. The explanation is that the phase bias products used for SD AR are not free of errors and the errors may degrade the KBR validation. In-flight PCV calibration and ambiguity resolution improve the LEO orbit accuracy effectively.

1. Introduction

Low Earth orbit (LEO) satellites are considered as key technologies for space missions due to their advantages of flexibility, redundancy, efficiency, and low cost. To fulfill scientific mission requirements, precise absolute or relative positions are required. The precise orbit determination (POD) capability, based on spaceborne Global Positioning System (GPS) observations, has been successfully verified on the TOPEX/Poseidon satellite [1,2]. The reduced dynamic POD technique [3,4] combined with the strength of GPS observations has shown its advantages in sampling rate and accuracy. Since then, more GPS receivers have been deployed on LEO satellites to meet the needs of various scientific missions for high precision orbit.
The quality of global navigation satellite system (GNSS) derived LEO orbits has steadily improved thanks to numerous improvements in the GNSS orbit and clock products provided by the International GNSS Service (IGS) [5], in the dynamic background models, and in modeling the carrier phase observations. LEO orbit accuracy of 1–3 cm can be fulfilled with float ambiguity (FA) [6,7]. With GPS precise orbit and clock products provided by IGS, CHAMP orbit calculated with ionosphere free (IF) observations achieves the accuracy of centimeter level [8]. GRACE orbit determined using GPS has been proven by independent satellite laser ranging (SLR) and K-band ranging (KBR) data, and 1 cm radial orbit accuracy has been achieved [9]. JASON’s GPS derived orbit is evaluated based on SLR residuals, dynamic orbit produced by SLR/DORIS data, and altimeter crossover tests [10,11], and the sub-centimeter radial accuracy is verified. Similar results are also obtained in other LEO satellite missions, such as GOCE [12,13], Sentinel [14,15], TerraSAR-X [7], HY [16], and FY-3 [17,18] satellites. Ambiguity resolution (AR) and in-flight antenna phase center variation (PCV) calibration are essential to fully exploit the precision of GPS observations for POD.
In order to fix the carrier phase ambiguities, a Kalman filter modeling the relative spacecraft dynamics has been developed for the GRACE mission and the double difference (DD) ambiguities are resolved to fully exploit the inherent measurement accuracy. Finally, the resulting GRACE orbit matches the KBR measurements with an accuracy of 1 mm [19]. A GRACE POD based on undifferenced and doubly differenced GPS data is studied [6]. Different baselines, including the space baseline between two GRACE satellites, the space-ground baselines consisting of LEO satellites and ground stations, and both types of baselines together, are processed. Results show that fixing of the GPS DD ambiguities has a significant impact on the space baseline [6]. With the orbit solution constrained by resolved DD ambiguities, a GRACE baseline accuracy of 2 mm is also achieved [20].
The double difference ambiguity resolution (DD AR) approach requires two satellites to construct the baseline. Several methods have been developed to resolve the integer ambiguity for the single receiver user. The idea of a single receiver integer ambiguity resolution forms the basis of PPPRTK. Taking the ionospheric delay as the unknown parameter, the common clock model [21,22] and distinct clocks model [23,24] are proposed to perform the single receiver ambiguity resolution. For IF formulation, the integer recovery clock (IRC) model [25,26], decoupled satellite clock (DSC) model [27], and uncalibrated phase delay (UPD) or fractional cycle bias (FCB) model [28,29] are presented. The S-system theory is applied by Teunissen to establish the linkage among different PPP-RTK methods [30,31]. Their differences are shown to lie (a) in the choice of S-basis; (b) in the choice of parameterization; (c) in the choice of whether or not to eliminate the ionospheric delay.
With wide-lane (WL) and phase biases provided by the Jet Propulsion Laboratory (JPL), the single receiver AR of GRACE and JASON-2 satellites is realized [32]. GRACE baseline accuracy is improved from 6 mm of float ambiguity solution to 2 mm of ambiguity fixed solution. Based on GPS orbit, clock, and WL bias products provided by Centre National d’Études Spatiales (CNES) [26], single receiver AR is also employed to generate the orbit of Sentinel-3A satellite and the root mean square (RMS) of SLR residuals decreases from 9 to 5 mm [15]. Using the observation specific bias (OSB) products provided by the Center for Orbit Determination in Europe (CODE), AR is applied to the POD of GRACE and Sentinel-3 satellites [33]. Ambiguity fixing improves the orbit quality with validation of KBR and SLR residuals, as well as the internal consistency between the reduced dynamic and kinematic orbits. With UPDs estimated via a global distributed network, kinematic orbits of Sentinel-3A and Swarm-A satellites are determined with AR [34]. SLR residuals show an improvement of 20% for AR solution when compared with the FA solution. Different AR solutions, including DD AR, single difference (SD) AR, and integrated SD and DD AR solutions are investigated with GRACE data to access their effects on orbit accuracy [35].
Another factor affecting the orbit accuracy of the LEO satellite is the phase center position of the spaceborne GNSS antenna. Nominal antenna models obtained from ground calibration have been made available for the antennas deployed on several space missions, such as GRACE [36], GOCE [37], Swarm [38], and TerraSAR-X [39]. Such nominal antenna models, however, do not reflect the influence of error sources, which are additionally encountered in the actual spacecraft environment, e.g., the influence of near-field multipath [36]. In-flight calibration of the LEO antenna is necessary for the stringent orbit accuracy requirement, especially for the altimeter mission. Two different approaches, the residual approach and the direct approach, can be used to determine the empirical PCV correction of the LEO antenna. Using the residual approach, the PCV map of JASON-1 satellite is created. The mean post-fit phase residual is decreased from 8 to 5 mm, and 1 cm radial accuracy is demonstrated [40]. PCV correction of the GOCE satellite is determined with 154 days of data and the consistency of reduced dynamic and kinematic orbits is improved when applying the PCV map [37]. An error in the given phase center offset (PCO) of Sentinel-1A antenna has been found by comparing different PCVs and different POD approaches [14]. The influence of relative PCV on precise baseline determination (PBD) of formation flying satellites is also studied [41]. With application of the generated relative PCVs of GRACE and GRACE-FO satellites, the consistency of KBR measurements is improved by 30%. The effects of antenna PCV on orbit determination and clock estimation for CentiSpace-1 satellite, which served as a navigation satellite, were analyzed [42].
The twin GRACE-FO satellites (named GRCC and GRCD hereafter), designed as a successor of GRACE, and jointly developed by the National Aeronautics and Space Administration (NASA) and the German Research Centre for Geosciences (GFZ), were launched from California, USA, on May 22, 2018 [43]. The main goals of GRACE-FO are to continuously provide high resolution monthly solutions of the Earth’s gravity field, surface mass change, and to measure the vertical temperature and humidity profiles of the Earth’s atmosphere [44]. GRACE and GRACE-FO satellites adopt the same satellite appearance design, and are equipped with the KBR system, GPS receiver, SLR retroreflector, star camera assembly (SCA), and other scientific instruments [44,45]. The appearance and instrument installation position of GRACE-FO satellites are illustrated in Figure 1 [46].
Effects of PCV correction on POD of GRACE-FO satellites have been explored [35,43]. However, impacts of different ambiguity resolution strategies and in-flight PCV calibration on POD and PBD still need more investigation. In this analysis, antenna PCV models of GRACE-FO satellites are developed to further exploit the POD accuracy. Single receiver and double difference integer ambiguity resolution methods are investigated and realized. With one year of data, GRACE-FO orbits based on different AR strategies, including FA, SD AR, and DD AR, are studied and evaluated with JPL’s orbit, SLR, and KBR measurements.
Section 2 provides details of the strategy and data adopted for GRACE-FO POD. The antenna PCV estimation approach and related results are presented in Section 3. Section 4 introduces the GNSS observation model and focuses on the SD and DD AR. POD and PBD results with PCV corrections and different AR methods were analyzed. The necessity of in-flight antenna calibration and ambiguity resolution is verified. Finally, discussions are made, followed by conclusions.

2. POD Strategy and Data Usage

2.1. POD Strategy

The reduced dynamic approach is employed to determine the orbits of GRACE-FO satellites. The Position And Navigation Data Analyst (PANDA) software [47], developed at the GNSS research center of Wuhan University, is modified and used for the POD process. The POD strategy is listed in Table 1. The macro model of GRACE-FO satellites [44] is applied to model the non-gravitational forces, which mainly result from atmospheric drag, solar, and earth radiation pressure. The atmospheric density values required for atmospheric drag modeling are obtained with the DTM94 model [48]. To account for deficiencies in DTM94 and the macro model, drag coefficient is estimated freely once per orbital revolution. Solar radiation pressure is calculated based on the satellite macro model and a scale factor is estimated per orbit determination arc. Additionally, one cycle per revolution (CPR) empirical accelerations in along-track and cross-track directions are estimated to compensate for deficiencies in the adopted force models. Other estimated parameters include receiver clock offsets and carrier phase ambiguities. The GRACE-FO antenna PCOs are applied according to the VGN1B product [44]. PCV corrections for IF carrier phase observations are calculated using the method described in Section 3.1.
According to the ambiguity resolution methods, three kinds of orbits are calculated. The first orbit is calculated with float ambiguities (donated as FA solution hereafter), the second solution is constrained by integer DD ambiguities (DD AR solution), and the third orbit is generated with the constraint of integer SD ambiguities (SD AR solution).

2.2. Data Usage

The twin GRACE-FO satellites fly a polar orbit with an altitude of 490 km. Over the mission lifetime, the two satellites will remain in co-planar orbits and the along-track separation is about 220 km. For POD, the GRACE-FO satellites are equipped with the new generation of the GNSS space science receiver, the TriG receiver. The receiver upgrades the capabilities offered by the BlackJack receiver, which was carried on the GRACE mission. GPS C/A, P1, and P2 pseudoranges, and associated carrier phase observations, namely LA, L1, and L2 can be provided by the receiver [54]. The GPS P1 and P2 pseudoranges and L1 and L2 carrier phase observations are used in this research.
One year of data in 2019 is processed. The data can be obtained from the GRACE-FO level 1B RL04 products, which are available at ftp://isdcftp.gfz-potsdam.de/grace-fo/ (accessed on 17 October 2021). The products also include the satellite attitude (SCA1B), precise orbits provided by JPL (GNV1B), biased inter-satellite ranges measured by KBR system (KBR1B), and the positions of the GPS antenna phase center (VGN1B). The GNV1B and KBR1B data can be used to assess the orbit quality. PCOs of ionosphere free carrier phase observations and SLR reflector positions in the science reference frame (SRF) are listed in Table 2. CODE’s GPS final orbits and 5 s clock products are used, and the associated OSB products are also applied to allow for SD AR [55,56].

3. Estimation of PCV Corrections

3.1. Mathematical Models

Modeling GNSS observation requires the computation of the geometric distance between the antenna phase center location of the GNSS satellite at signal emission time and the antenna phase center location of the receiving antenna at signal reception time. The phase center location usually differs from the mechanical antenna reference point (ARP). The difference vector is conventionally described by a set of phase center corrections. Such a set of corrections consists of a PCO vector, which defines the position of the mean antenna phase center, with respect to the ARP and a consistent function, which models the azimuth and zenith dependent PCVs [36].
Phase center corrections have some inherent degrees of freedom. One set of corrections consisting of a PCO vector r 0 and an azimuth and zenith dependent function ϕ ( α , z ) can be transformed into a new set, consisting of r 0 and ϕ ( α , z ) , which gives the same result:
r 0 = r 0 + Δ r ϕ ( α , z ) = ϕ ( α , z ) Δ r e + Δ ϕ
where α and z are azimuth and zenith angles, Δ ϕ is an arbitrary offset and cannot be separated from the receiver clock. The unit vector e denotes the direction from the receiver to the satellite. The offset vector Δ r can be chosen arbitrarily. Preferably, PCVs should not induce a PCO and Δ r should be zero. In that case, the mean antenna phase center is explicitly defined by the PCO. This convention is particularly important if one would only apply PCO and no PCV [14].
Taking computational burden into consideration, residual approach is employed in this research. The antenna PCV is represented as piecewise linear function with respect to the zenith and azimuth angles in the antenna fixed coordinate system. The model assumes that PCV is composed of different grids, and zenith and azimuth angles are equally divided. When the observation is at the point of P within the grid ABCD, its PCV value can be linearly interpolated:
Δ P C V , P = ( 1 γ ) ( 1 β ) Δ P C V , A + γ ( 1 β ) Δ P C V , B + γ β Δ P C V , C + ( 1 γ ) β Δ P C V , D γ = ( α α 1 ) / ( α 2 α 1 ) , β = ( z z 1 ) / ( z 2 z 1 )
where Δ P C V , P is the PCV at point P and is observation, Δ P C V , i ( i = A , B , C , D ) are the PCVs at points A, B, C, D are parameters to be estimated. γ and β are combination coefficients, α and z are azimuth and zenith angles of P, α 1 is the azimuth of point A and D, α 2 is the azimuth of B and C, z 1 is the zenith of A and B, z 2 is the zenith of C and D.

3.2. Results

The PCOs of GRACE-FO satellites provided in the VGN1B product are introduced as fixed and only the PCVs are estimated. Residual approach [36] is employed with three iterations to obtain the PCV maps. Resolutions in azimuth and zenith are 5°. Two kinds of PCV maps are created—one is generated with the phase residuals from DD AR solution and the other is derived from the SD AR solution.
Figure 2 are PCV maps of GRACE-FO satellites produced from the DD AR solution. The azimuth of 0° points into the direction of flight. A similarity exists between the PCVs of GRCC and GRCD satellites, especially in the azimuth ranges of 90°–180° and 270°–360°. The absolute maximum PCV of GRCC satellite is −14.56 mm with azimuth angle 185° and elevation angle 25°. For the GRCD satellite, the absolute maximum PCV is −14.95 mm with azimuth angle 190° and elevation angle 30°.
Figure 3 presents the PCV maps produced from the SD AR solution. The absolute maximum PCV of GRCC satellite is −12.57 mm with azimuth angle 185° and elevation angle 25°, while the GRCD satellite is −12.15 mm with azimuth angle 185° and elevation angle 30°. The PCV maps generated from DD AR and SD AR solutions have similar patterns when comparing Figure 2a with Figure 3a and comparing Figure 2b with Figure 3b. The PCVs made from the SD AR solution have less maximum values.
PCV maps in Figure 2 and Figure 3 are used to correct the carrier phase observations and then to perform the POD. All available carrier phase residuals are averaged and Figure 4 illustrates the results of GRCC satellite. Figure 4a presents the carrier phase residuals from the DD AR solution while Figure 4b presents the residuals from the SD AR solution. It should be noted that the limit of the color bar is 5 mm. With application of the PCV corrections, most of the carrier phase residuals of the DD AR solution are less than 1 mm and the mean value is 0.13 mm. The carrier phase residuals of the SD AR solution are a bit larger, especially in the range with azimuth of 0°–30° and elevation of 30°–60°, the average is 0.41 mm. Similar results are obtained for the GRCD satellite and are not shown here. Carrier phase residuals can be used to measure the consistency between the applied models and GPS observations. It can be induced that the application of PCV corrections improves the orbit accuracy. The effects of PCV correction on POD and PBD will be presented in the following sections.

4. POD with Ambiguity Resolution

The GNSS observation model and the methods employed to resolve the SD and DD ambiguities are introduced in this part. FA, DD AR, and SD AR schemes are applied to calculate the GRACE-FO orbits. The generated orbits are assessed with JPL’s orbit, SLR, and KBR data. Effects of PCV corrections on POD and PBD are also analyzed.

4.1. Mathematical Models

4.1.1. GNSS Observation Model

GNSS code and carrier phase observations between a satellite and a receiver are usually described by the following equations [35,57]:
P r , j s = ρ r s + c ( d t r d t s ) + I r , j s + b r , j b j s L r , j s = ρ r s + c ( d t r d t s ) I r , j s + λ j ( N r , j s + B r , j B j s ) + λ j ω r s
where P r , j s and L r , j s are code and carrier phase observations in meters, r and s represent the receiver and GNSS satellite, j is the signal frequency, ρ r s is the geometry distance between receiver and satellite, c is the speed of light, d t r and d t s are receiver and satellite clock errors, I r , j s is ionospheric delay [58], b r , j and b j s are receiver and satellite hardware biases of pseudorange, B r , j and B j s are hardware biases of the carrier phase observation, λ j is signal wavelength and N r , j s is the integer ambiguity. ω r s is the phase wind-up error and can be corrected by model [53], in the remaining part, this item will be omitted.
Both ionosphere free and ionosphere float models can be used for POD. The first order ionospheric delay can be eliminated by the ionosphere free model, while in the ionosphere float model, the ionospheric delay is estimated as an unknown parameter. The ionosphere free approach eliminates the ionospheric delay twice and leads to smaller redundancy. The ionosphere free method owns fewer unknown parameters and it also has the drawback of lacking flexibility for further model strengthening [59]. In this investigation, the ionosphere free model is employed.
Following the IGS convention, the IF pseudorange biases will be assimilated into the receiver and satellite clock offsets. IF observations can be formulated as:
P r , I F s = ρ r s + c ( d ¯ t r d ¯ t s ) L r , I F s = ρ r s + c ( d ¯ t r d ¯ t s ) + λ 1 N ¯ r , I F s N ¯ r , I F s = N r , I F s + ( B r , I F b r , I F / λ 1 ) ( B I F s b I F s / λ 1 )
where P r , I F s and L r , I F s are IF code and carrier phase observations, d ¯ t r and d ¯ t s are receiver and satellite clock errors, including the pseudorange hardware biases. b r , I F and b I F s are the IF pseudorange hardware biases while B r , I F and B I F s are the carrier phase biases. λ 1 is the wavelength of L1. N ¯ r , I F s is the IF ambiguity, which includes the hardware biases in both pseudorange and carrier phase observations, N r , I F s is a combination of integer L1 and L2 ambiguities.
In order to fix the ambiguities, Blewitt presents a sequential integer ambiguity resolution method to resolve the DD ambiguities of long baseline. The integer DD WL ambiguities are firstly determined with a combination of pseudorange and carrier phase observations. Then the DD narrow-lane (NL) ambiguities are fixed using the estimated IF ambiguities and already fixed integer WL ambiguities [57]. For a single receiver integer ambiguity resolution, GNSS satellites SD observations are used to cancel the receiver based phase bias. Guo et al. provide a detailed description about the model to realize single receiver AR [35]. These methods are employed in this contribution to fix the SD and DD ambiguities.

4.1.2. Single Receiver Ambiguity Resolution

The float IF ambiguity can be rewritten as the combination of integer WL and float NL ambiguities as the following [57]:
N ¯ r , I F s = f 1 f 2 f 1 2 f 2 2 N r , W L s + f 1 f 1 + f 2 N ¯ r , N L s N r , W L s = N r , 1 s N r , 2 s N ¯ r , N L s = N r , N L s + d r , N L d N L s
where N r , W L s and N r , N L s are integer WL and NL ambiguities, N ¯ r , N L s is float NL ambiguity, d r , N L and d N L s are the NL FCBs of the receiver and satellite.
The integer WL and NL ambiguities are fixed with two steps. Firstly, the WL ambiguity is resolved using the Hatch–Melbourne–Wübbena combination [60,61,62]:
N ¯ r , W L s = ( f 1 L r , 1 s f 2 L r , 2 s f 1 f 2 f 1 P r , 1 s + f 2 P r , 2 s f 1 + f 2 ) / λ W L = N r , W L s + d r , W L d W L s
where d r , W L and d W L s are receiver and satellite WL FCBs, λ W L is WL wavelength. After correction of the WL FCBs, the integer WL ambiguity can be determined. Then integer NL ambiguity can be derived with the fixed WL ambiguity, the estimated float IF ambiguity, and the associated NL FCBs using Equation (5).
Once the WL and NL ambiguities are resolved, the IF ambiguity can be reconstructed according to Equation (5). Usually satellite dependent WL and NL FCBs can be estimated with ground network stations or provided by GNSS analysis institutions [26,55,63]. However, receiver dependent FCBs are not available and cannot be removed from Equations (5) and (6). SD operation between satellites is suggested to cancel the receiver FCBs. The SD integer WL and NL ambiguities can be calculated using Equation (7):
Δ N r , W L s 1 , s 2 = Δ N ¯ r , W L s 1 , s 2 + Δ d W L s 1 , s 2 Δ N r , N L s 1 , s 2 = Δ N ¯ r , N L s 1 , s 2 + Δ d N L s 1 , s 2
where Δ is single difference operator, Δ N ¯ r , W L s 1 , s 2 and Δ N ¯ r , N L s 1 , s 2 are float SD WL and NL ambiguities of satellite s 1 and s 2 while Δ N r , W L s 1 , s 2 and Δ N r , N L s 1 , s 2 are integer ambiguities. Δ d W L s 1 , s 2 and Δ d N L s 1 , s 2 are SD hardware biases between satellites. Once SD WL and NL ambiguities are fixed, the SD IF ambiguity can be recovered as follows:
Δ N ¯ r , I F s 1 , s 2 = f 1 f 2 f 1 2 f 2 2 Δ N r , W L s 1 , s 2 + f 1 f 1 + f 2 ( Δ N r , N L s 1 , s 2 Δ d N L s 1 , s 2 )
The recovered SD IF ambiguity is taken as pseudo-observation to constrain the undifferenced IF ambiguity parameters in the POD estimation using Equation (9):
Δ N ¯ r , I F s 1 , s 2 = N ¯ r , I F s 1 N ¯ r , I F s 2 , W r s 1 , s 2
where W r s 1 , s 2 is the weight of the pseudo-observation and will be assigned with a large value to put a strong confidence.

4.1.3. Double Difference Ambiguity Resolution

With regard to the DD ambiguity, both receiver and satellite FCBs can be removed. The WL and NL DD ambiguities are theoretically integers according to Equations (5) and (6). The DD IF ambiguity can be reconstructed by:
Δ N ¯ r 1 , r 2 , I F s 1 , s 2 = f 1 f 2 f 1 2 f 2 2 Δ N r 1 , r 2 , W L s 1 , s 2 + f 1 f 1 + f 2 Δ N r 1 , r 2 , N L s 1 , s 2
where Δ is double difference operator, r 1 and r 2 denote a receiver pair. The DD IF ambiguity then is also treated as a pseudo-observation to constrain the undifferenced IF ambiguities.
Δ N ¯ r 1 , r 2 , I F s 1 , s 2 = N ¯ r 2 , I F s 2 N ¯ r 2 , I F s 1 N ¯ r 1 , I F s 2 + N ¯ r 1 , I F s 1 , W r 1 , r 2 s 1 , s 2
W r 1 , r 2 s 1 , s 2 is the weight of the pseudo-observation. The precision of the DD and SD IF ambiguities is set to be 1 × 10−4 L1 cycles (about 0.02 mm) to impose strong constraints.

4.1.4. Integer Ambiguity Validation

Several integer ambiguity resolution methods have been proposed including the integer rounding, integer bootstrapping, integer least squares, and partial ambiguity resolution methods [64]. The integer rounding approach is used in this research to resolve the WL and NL ambiguities. The ambiguity validation procedure is presented here for clearness. The reader can refer to [65] for more details.
Let b be the float WL or NL ambiguities and it is assumed that the probability density function for b is:
P ( b n ) = 1 ( 2 π ) 1 / 2 σ exp [ ( b n ) 2 2 σ 2 ]
where n is the true integer ambiguity, σ is the formal uncertainty from the adjustment. Let I be the nearest integer of b . In the ambiguity rounding procedure, it should be determined whether to fix the ambiguity to I or leave it as float. The deviation
x = b I
must fall in the interval [−0.5, 0.5]. If the true integer ambiguity n is not equal to the nearest integer I , the probability of wrong ambiguity resolution can be calculated using Equation (14):
Q 0 = m = 1 e r f c ( m x 2 σ ) e r f c ( m + x 2 σ ) α
e r f c ( x ) = 2 π x e t 2 d t
where α is the allowable rate of wrongly fixed ambiguity, e r f c is the complementary error function. In order to eliminate some extremes, a decision function is defined:
d ( x , σ ) = T / Q 0
T = 0 , x T 1 o r σ T 2 ( 1 x T 1 ) ( 3 T 2 3 σ ) , o t h e r w i s e
where T 1 is the limit for the deviation x and T 2 is the limit for σ . In this analysis, T 1 and T 2 are both 0.25 cycles for WL ambiguity resolution, and are 0.15 cycles for NL ambiguity resolution. When d ( x , σ ) is greater than 1000, the ambiguity will be fixed to the nearest integer [28,66].

4.2. Results

One year of data in 2019 collected by GRACE-FO satellites were analyzed; 38 days were excluded due to large satellite maneuvers or data gaps. The method introduced in Section 4.1.2 is used to resolve the SD WL and NL ambiguities, and the method in Section 4.1.3 is adopted to fix the DD ambiguities. Firstly, the residuals and fixing rates of WL and NL ambiguities are presented. Then the absolute orbit accuracy of each satellite and baseline accuracy between two satellites are evaluated.

4.2.1. Ambiguity Resolution Results

Ambiguity resolution performance is evaluated in terms of the residual and fixing rate. Ambiguity residual is the difference between the float ambiguity and its nearest integer. After removing the satellite and receiver FCBs, the float SD WL and NL ambiguities should be close to the integers. For DD ambiguities, the FCBs of the receiver and satellite are cancelled, which leads to an integer nature. As illustrated by Teunissen [67], the distribution of the ambiguity residuals is non-Gaussian. It is symmetric and the point of symmetry is the origin, which implies that the mean of the ambiguity residuals is zero. The distribution will become more peaked when the estimated float ambiguity is more precise. Figure 5a,b are residual distributions of SD WL and NL ambiguities while Figure 5c,d are distributions of DD ambiguities. Blue bars in the figures are percentages of the related residuals. The mean values of SD and DD WL and NL residuals are all less than 0.01 cycles and verify the symmetry of distribution. The standard deviations (STDs) of SD WL and NL residuals are 0.10 cycles and 0.07 cycles and validate the high quality of the estimated float ambiguities. For DD WL and NL ambiguities, the STDs are both 0.06 cycles. The smaller STDs of DD ambiguity residuals confirm the effectiveness of the DD operation to cancel the receiver and satellite FCBs.
The ambiguity fixing rate is the number of fixed ambiguities divided by the total number of ambiguities. Figure 6 illustrates the daily ambiguity fixing rates of SD and DD ambiguities. Figure 6a is for SD ambiguities while Figure 6b is for DD ambiguities. The upper and lower panels of Figure 6a are results of GRCC and GRCD satellites. Blue and red dots in Figure 6 are fixing rates of WL and NL ambiguities. For the GRCC satellite, one year average fixing rates of SD WL and NL ambiguities are 98.5% and 95.8%. A similar performance is observed for the GRCD satellite; the fixing rates are 98.4% and 96.5%. The excellent fixing rates confirm the good quality of applied phase bias products. For the DD case, higher fixing rates are obtained for WL and NL ambiguities, which are 99.1% and 96.7%. Gaps in day of year (DOY) 32–51, 146–149, 205–211, and 299–302 are because of data missing.

4.2.2. Single Satellite Orbit Validation

JPL’s precise orbit of GRACE-FO satellites is also included in the level 1B product. The orbit is calculated with single receiver ambiguity resolution [32,44] and can be used to assess the orbits generated in this investigation. The effect of PCV correction on POD is first evaluated. Figure 7 presents the along-track orbit differences between SD AR solution and JPL’s orbit. The upper panel is for GRCC satellite and lower part is for GRCD satellite. The red dots are daily RMSs of differences of the orbit determined with PCV corrections while blue dots are RMSs without PCV corrections. Using created PCV maps, the RMS of the GRCC satellite decreases from 7.2 to 6.0 mm, and from 6.4 to 5.9 mm for the GRCD satellite. The detailed results for three components of SD AR and DD AR orbits can be found in Table 3. The average improvements in along-track, cross-track, and radial direction are 12.2%, 1.6%, and 1.0% for SD AR solution, and are 10.3%, 1.5%, and 3.5% for DD AR solution. The improvement of along-track orbit accuracy is the most significant, which validates the necessity of PCV calibration, especially for PBD.
RMSs on DOY 100, 141, and 159 of GRCC satellite, and RMSs on DOY 204 and 290 of GRCD satellite are larger than other days. It has been found that sometimes, in these days, fewer GPS satellites are available, which leads to a decrease in orbit accuracy.
Orbit accuracy with different AR methods is also evaluated. Figure 8 presents the RMSs of along-track orbit differences of FA, DD AR, and SD AR solutions with JPL’s orbit. For GRCC satellite, the along-track RMS of the FA solution is 8.8 mm, the RMS is improved to 8.3 mm with DD AR and further achieves 6.0 mm by using the SD AR method. Results of three orbit components of the two satellites are listed in Table 4. Compared with the FA solution, the DD AR orbit has an average improvement of 6.3% while the SD AR result has an improvement of 18.3%. Similar to PCV correction, the along-track orbit accuracy has the most significant improvement, which is more than 30%.
GRACE-FO satellites are equipped with SLR retroreflectors. The independent SLR data can be used as external validation. SLR normal points from nine laser stations (Yarragadee, Greenbelt, Haleakala, Zimmerwald, Mount Stromlo, Wettzell, Graz, Herstmonceux, and Potsdam) are used for orbit evaluation. Residuals larger than 10 cm are rejected and an elevation cutoff of 10° is applied. Figure 9 shows the RMSs of SLR residuals. Red, green, and blue dots are RMSs of FA, DD AR, and SD AR orbits. The RMSs of the GRCC satellite are 11.5, 10.2, and 9.6 mm, which confirms the positive impact of AR on the orbit quality. Compared with the DD AR solution, the SD AR result has better performance, which is also validated by JPL’s orbit.

4.2.3. Baseline Validation

The GRACE-FO mission offers the possibility of validating the computed orbit, in particular, the along-track component, with the ultra-precise KBR measurements [33,68]. Figure 10 illustrates the STDs of KBR residuals (named KBR STD hereafter) calculated using FA, DD AR, and SD AR orbits. The effect of antenna PCV correction on PBD is also illustrated. Figure 10a presents the KBR STDs of the SD AR orbit with PCV corrections (blue dots), the SD AR orbit without PCV corrections (green dots), and FA orbit with PCV corrections (red dots). PCV corrections used to calculate FA orbit are generated from the SD AR solution. Figure 10b is the result of the DD AR solution.
Both SD and the DD AR lead to a significant reduction in KBR residuals. Compared with the FA result, the SD AR orbit improves the KBR STD from 5.9 to 3.1 mm, and the DD AR orbit improves the STD to 1.6 mm. With PCV corrections, the KBR STD of the SD AR solution decreases from 3.1 to 1.8 mm, and the STD of the DD AR orbit decreases from 1.6 to 0.9 mm. The baseline precision of the SD AR solution gains similar performance, as presented by Arnold et al. [33]. The DD AR orbit achieves a sub-millimeter level as reported in previous research [6,19,69]. Comparisons with JPL’s orbit and SLR data show that the SD AR solution possesses better performance, while in the KBR validation, the residual STD of the DD AR orbit is smaller. In other words, the SD AR orbit has better absolute orbit accuracy, while the DD AR orbit gets better relative performance. This is explained by two facts: (a) phase bias products for SD AR are not free of errors and these errors may degrade the KBR validation [35]; (b) both JPL’s orbit and the SD AR orbit in this investigation adopt the single receiver ambiguity resolution method and a better consistency can be inferred.

5. Discussions

Precise knowledge of the phase center location of the GNSS antenna is a prerequisite for high precision LEO orbit determination. PCO and PCV values of the LEO antenna obtained from ground calibration cannot reflect the influence of error sources, which are additionally encountered in the actual spacecraft environment. An in-flight calibration of the LEO antenna is thus mandatory. In this analysis, antenna PCV is represented as a piecewise linear function with respect to zenith and azimuth angles in the antenna fixed coordinate system. A residual approach is employed to estimate the PCV correction. Different ambiguity resolution strategies, including SD AR and DD AR, are investigated to fully exploit the precision of GPS observations for POD. The single receiver integer ambiguity resolution concept employed here makes use of the carrier phase biases and clock products provided by CODE. One year of GRACE-FO data in 2019 are used to demonstrate the effect of integer ambiguity resolution and PCV correction on POD.
PCV maps derived from SD AR and DD AR solutions have very similar patterns. With PCV corrections, carrier phase residuals are reduced from about 10 to 1–2 mm, which implies better consistency between the observations and applied models. The average improvements of the SD AR solution in the along-track, cross-track, and radial directions are 12.2%, 1.6%, and 1.0%, and for the DD AR solution, the improvements are 10.3%, 1.5%, and 3.5%. The great enhancement in the along-track orbit accuracy validates the necessity to calibrate the PCV errors in the relative POD.
The mean values of both SD and DD ambiguity residuals are less than 0.01 cycles and confirm the symmetry of distribution. The fixing rates of SD WL and NL ambiguities are more than 98% and 95%. Compared with the FA solution, DD AR orbit accuracy has an average improvement of 6.3% while the SD AR result gains an increase of 18.3%. For high-grade SLR stations, range residuals with RMS less than 10 mm are achieved for SD AR orbit, which marks a 17% improvement compared to the FA result. Independent KBR measurements are also used as external validations. KBR residuals STD of the SD AR solution is 1.8 mm while STD of the DD AR orbit is 0.9 mm and reaches the sub-millimeter level. The better baseline accuracy for the DD AR solution is explained by the errors in phase bias products.

6. Conclusions

Integer ambiguity resolution plays a crucial role in achieving the best positioning or orbit accuracy. Single receiver and double difference integer ambiguity resolution models are introduced in this contribution. One year of GRACE-FO data were analyzed to verify the orbit improvement with different ambiguity resolution strategies. In-flight antenna PCV maps were developed to further exploit the POD accuracy.
With PCV corrections, along-track orbit accuracy of the SD AR solution is improved by 12.6%, and that of the DD AR solution is improved by 10.3%, which verifies the necessity of in-flight antenna calibration. With phase biases cancelled by the double difference operation, fixing rates of DD ambiguities are more than 96.7%. Similar performance is achieved by SD AR, which confirms the consistency of theoretical models and bias/clock products.
JPL’s orbited together with independent SLR and KBR measurements are used to assess the results of different AR strategies. Both SD AR and DD AR solutions improve the orbit accuracy. SD AR solution provides the best performance of absolute orbit. Compared with JPL’s orbit, the RMS of the SD AR orbit differences is better than 6 mm and the SLR residual RMS is less than 10 mm. The DD AR solution realizes the highest baseline accuracy and the STD of KBR residuals achieves 0.9 mm. In addition, sometimes there are fewer GPS satellites available, which leads to the reduction of orbit accuracy. With the application of multi-mode GNSS receivers and multi-GNSS phase bias products [70,71], orbit accuracy based on single receiver ambiguity resolution can be further improved.

Author Contributions

Conceptualization, B.J. and Y.L.; methodology, B.J.; software, B.J. and K.J.; validation, B.J. and S.C.; writing—original draft preparation, B.J.; writing—review and editing, S.C.; supervision, Y.L. and Z.L.; funding acquisition, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 12033009.

Data Availability Statement

GRACE-FO data in this paper can be freely accessed at ftp://isdcftp.gfz-potsdam.de/grace-fo/ (accessed on 17 October 2021). The GPS final orbit, clock, and OSB products can be found at ftp://ftp.aiub.unibe.ch/CODE/ (accessed on 17 October 2021).

Acknowledgments

The authors thank GFZ for providing GRACE-FO GPS data. The authors are grateful to the CODE for providing GPS products.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bertiger, W.; Bar-Sever, Y.; Christensen, E.; Davis, E.; Guinn, J.; Haines, B.; Ibanez-Meier, R.; Jee, J.; Lichten, S.; Melbourne, W.; et al. GPS precise tracking of TOPEX/POSEIDON: Results and implications. J. Geophys. Res. 1994, 99, 24449–24464. [Google Scholar] [CrossRef]
  2. Tapley, B.D.; Ries, J.C.; Davis, G.W.; Eanes, R.J.; Schutz, B.E.; Shum, C.K.; Watkins, M.M.; Marshall, J.A.; Nerem, R.S.; Putney, B.H.; et al. Precision orbit determination for TOPEX/POSEIDON. J. Geophys. Res. 1994, 99, 24383–24404. [Google Scholar] [CrossRef]
  3. Wu, S.C.; Yunck, T.P.; Thornton, C.L. Reduced-dynamic technique for precise orbit determination of low earth satellites. J. Guid. Control Dyn. 1991, 14, 24–30. [Google Scholar] [CrossRef]
  4. Jäggi, A.; Hugentobler, U.; Beutler, G. Pseudo-Stochastic Orbit Modeling Techniques for Low-Earth Orbiters. J. Geod. 2006, 80, 47–60. [Google Scholar] [CrossRef] [Green Version]
  5. Dow, J.M.; Neilan, R.E.; Gendt, G. The International GPS Service: Celebrating the 10th anniversary and looking to the next decade. Adv. Space Res. 2005, 36, 320–326. [Google Scholar] [CrossRef]
  6. Jäggi, A.; Hugentobler, U.; Bock, H.; Beutler, G. Precise orbit determination for GRACE using undifferenced or doubly differenced GPS data. Adv. Space Res. 2007, 39, 1612–1619. [Google Scholar] [CrossRef]
  7. Hackel, S.; Montenbruck, O.; Steigenberger, P.; Balss, U.; Gisinger, C.; Eineder, M. Model improvements and validation of TerraSAR-X precise orbit determination. J. Geod. 2017, 91, 547–562. [Google Scholar] [CrossRef]
  8. Van den Ijssel, J.; Visser, P.; Patiño Rodriguez, E. Champ precise orbit determination using GPS data. Adv. Space Res. 2003, 31, 1889–1895. [Google Scholar] [CrossRef]
  9. Kang, Z.; Tapley, B.D.; Bettadpur, S.; Ries, J.C.; Nagel, P.; Pastor, R. Precise orbit determination for the GRACE mission using only GPS data. J. Geod. 2006, 80, 322–331. [Google Scholar] [CrossRef]
  10. Luthcke, S.B.; Zelensky, N.P.; Rowlands, D.D.; Lemoine, F.G.; Williams, T.A. The 1-Centimeter Orbit: Jason-1 Precision Orbit Determination Using GPS, SLR, DORIS, and Altimeter Data Special Issue: Jason-1 Calibration/Validation. Mar. Geod. 2003, 26, 399–421. [Google Scholar] [CrossRef]
  11. Bertiger, W.; Desai, S.D.; Dorsey, A.; Haines, B.J.; Harvey, N.; Kuang, D.; Sibthorpe, A.; Weiss, J.P. Sub-centimeter precision orbit determination with GPS for ocean altimetry. Mar. Geod. 2010, 33, 363–378. [Google Scholar] [CrossRef]
  12. Bock, H.; Jäggi, A.; Meyer, U.; Visser, P.; van den IJssel, J.; van Helleputte, T.; Heinze, M.; Hugentobler, U. GPS-derived orbits for the GOCE satellite. J. Geod. 2011, 85, 807–818. [Google Scholar] [CrossRef] [Green Version]
  13. Bock, H.; Jäggi, A.; Beutler, G.; Meyer, U. GOCE: Precise orbit determination for the entire mission. J. Geod. 2014, 88, 1047–1060. [Google Scholar] [CrossRef]
  14. Peter, H.; Jäggi, A.; Fernández, J.; Escobar, D.; Ayuga, F.; Arnold, D.; Wermuth, M.; Hackel, S.; Otten, M.; Simons, W.; et al. Sentinel-1A-First precise orbit determination results. Adv. Space Res. 2017, 60, 879–892. [Google Scholar] [CrossRef]
  15. Montenbruck, O.; Hackel, S.; Jäggi, A. Precise orbit determination of the Sentinel-3A altimetry satellite using ambiguity-fixed GPS carrier phase observations. J. Geod. 2018, 92, 711–726. [Google Scholar] [CrossRef] [Green Version]
  16. Guo, J.Y.; Kong, Q.; Qin, J.; Sun, Y. On precise orbit determination of HY-2 with space geodetic techniques. Acta Geophys. 2013, 61, 752–772. [Google Scholar] [CrossRef]
  17. Li, M.; Li, W.; Shi, C.; Jiang, K.; Guo, X.; Dai, X.; Meng, X.; Yang, Z.; Yang, G.; Liao, M. Precise orbit determination of the Fengyun-3C satellite using onboard GPS and BDS observations. J. Geod. 2017, 91, 1313–1327. [Google Scholar] [CrossRef] [Green Version]
  18. Li, X.; Zhang, K.; Meng, X.; Zhang, W.; Zhang, Q.; Zhang, X.; Li, X. Precise Orbit Determination for the FY-3C Satellite Using Onboard BDS and GPS Observations from 2013, 2015, and 2017. Engineering 2019, 6, 904–912. [Google Scholar] [CrossRef]
  19. Kroes, R.; Montenbruck, O.; Bertiger, W.; Visser, P. Precise GRACE baseline determination using GPS. GPS Solut. 2005, 9, 21–31. [Google Scholar] [CrossRef]
  20. Zhao, Q.; Hu, Z.; Guo, J.; Li, M.; Ge, M. Precise relative orbit determination of twin GRACE satellites. Geo Spat. Inf. Sci. 2010, 13, 221–225. [Google Scholar] [CrossRef]
  21. Zhang, B.; Teunissen, P.J.G.; Odijk, D. A novel un-differenced PPP-RTK concept. J. Navig. 2011, 64, S180–S191. [Google Scholar] [CrossRef] [Green Version]
  22. Odijk, D.; Teunissen, P.J.G.; Zhang, B. Single-frequency integer ambiguity resolution enabled GPS precise point positioning. J. Surv. Eng. 2012, 138, 193–202. [Google Scholar] [CrossRef]
  23. De Jonge, P.J. A processing strategy for the application of the GPS in networks. Ph.D. Thesis, Delft University of Technology, Delft, The Netherlands, 1998. [Google Scholar]
  24. Teunissen, P.J.G.; Odijk, D.; Zhang, B. PPP-RTK: Results of CORS Network-Based PPP with Integer Ambiguity Resolution. J. Aeronaut. Astronaut. Aviat. Ser. 2010, 42, 223–230. [Google Scholar]
  25. Laurichesse, D.; Mercier, F.; Berthias, J.P.; Broca, P.; Cerri, L. Integer ambiguity resolution on undifferenced GPS phase measurements and its application to PPP and satellite precise orbit determination. Navigation 2009, 56, 135–149. [Google Scholar] [CrossRef]
  26. Loyer, S.; Perosanz, F.; Mercier, F.; Capdeville, H.; Marty, J.C. Zero-difference GPS ambiguity resolution at CNES–CLS IGS Analysis Center. J. Geod. 2012, 86, 991–1003. [Google Scholar] [CrossRef]
  27. Collins, P.; Bisnath, S.; Lahaye, F.; Héroux, P. Undifferenced GPS Ambiguity Resolution Using the Decoupled Clock Model and Ambiguity Datum Fixing. Navigation 2010, 57, 123–135. [Google Scholar] [CrossRef] [Green Version]
  28. Ge, M.; Gendt, G.; Rothacher, M.; Shi, C.; Liu, J. Resolution of GPS carrier-phase ambiguities in precise point positioning (PPP) with daily observations. J. Geod. 2008, 82, 389–399. [Google Scholar] [CrossRef]
  29. Geng, J.; Bock, Y. Triple-frequency GPS precise point positioning with rapid ambiguity resolution. J. Geod. 2013, 87, 449–460. [Google Scholar] [CrossRef]
  30. Teunissen, P.J.G. Zero Order Design: Generalized Inverses, Adjustment, the Datum Problem and S-Transformations. In Optimization and Design of Geodetic Networks; Grafarend, E., Sanso, F., Eds.; Springer: Berlin/Heidelberg, Germany, 1985; pp. 11–55. [Google Scholar]
  31. Teunissen, P.J.G.; Khodabandeh, A. Review and principles of PPP-RTK methods. J. Geod. 2015, 89, 217–240. [Google Scholar] [CrossRef]
  32. Bertiger, W.; Desai, S.D.; Haines, B.; Harvey, N.; Moore, A.W.; Owen, S.; Weiss, J.P. Single receiver phase ambiguity resolution with GPS data. J. Geod. 2010, 84, 327–337. [Google Scholar] [CrossRef]
  33. Arnold, D.; Schaer, S.; Villiger, A.; Dach, R.; Jäggi, A. Undifferenced ambiguity resolution for GPS-based precise orbit determination of low Earth orbiters using the new CODE clock and phase bias products. In Proceedings of the International GNSS Service Workshop 2018, Wuhan, China, 29 October–2 November 2018. [Google Scholar]
  34. Li, X.; Wu, J.; Zhang, K.; Li, X.; Xiong, Y.; Zhang, Q. Real-Time Kinematic Precise Orbit Determination for LEO Satellites Using Zero-Differenced Ambiguity Resolution. Remote Sens. 2019, 11, e2815. [Google Scholar] [CrossRef] [Green Version]
  35. Guo, X.; Geng, J.; Chen, X.; Zhao, Q. Enhanced orbit determination for formation flying satellites through integrated single and double difference GPS ambiguity resolution. GPS Solut. 2020, 24, 1–12. [Google Scholar] [CrossRef]
  36. Jäggi, A.; Dach, R.; Montenbruck, O.; Hugentobler, U.; Bock, H.; Beutler, G. Phase center modeling for LEO GPS receiver antennas and its impact on precise orbit determination. J. Geod. 2009, 83, 1145–1162. [Google Scholar] [CrossRef] [Green Version]
  37. Bock, H.; Jäggi, A.; Meyer, U.; Dach, R.; Beutler, G. Impact of GPS antenna phase center variations on precise orbits of the GOCE satellite. Adv. Space Res. 2011, 47, 1885–1893. [Google Scholar] [CrossRef]
  38. Van Den IJssel, J.; Encarnação, J.; Doornbos, E.; Visser, P. Precise science orbits for the Swarm satellite constellation. Adv. Space Res. 2015, 56, 1042–1055. [Google Scholar] [CrossRef]
  39. Montenbruck, O.; Garcia-Fernandez, M.; Yoon, Y.; Schön, S.; Jäggi, A. Antenna phase center calibration for precise positioning of LEO satellites. GPS Solut. 2008, 13, 23–34. [Google Scholar] [CrossRef]
  40. Haines, B.; Bar-Sever, Y.; Bertiger, W.; Desai, S.D.; Willis, P. One-centimeter orbit determination for Jason-1: New GPS-based strategies. Mar. Geod. 2004, 27, 299–318. [Google Scholar] [CrossRef]
  41. Yi, B.; Gu, D.; Ju, B.; Shao, K.; Zhang, H. Enhanced baseline determination for formation flying LEOs by relative corrections of phase center and code residual variations. Chin. J. Aeronaut. 2019, in press. [Google Scholar]
  42. Yuan, J.; Zhou, S.; Hu, X.; Yang, L.; Cao, J.; Li, K.; Liao, M. Impact of Attitude Model, Phase Wind-Up and Phase Center Variation on Precise Orbit and Clock Offset Determination of GRACE-FO and CentiSpace-1. Remote Sens. 2021, 13, e2636. [Google Scholar] [CrossRef]
  43. Xia, Y.; Liu, X.; Guo, J.; Yang, Z.; Qi, L.; Ji, B.; Chang, X. On GPS data quality of GRACE-FO and GRACE satellites: Effects of phase center variation and satellite attitude on precise orbit determination. Acta Geod. Geophys. 2021, 56, 93–111. [Google Scholar] [CrossRef]
  44. Wen, H.; Kruizinga, G.; Paik, M.; Landerer, F.; Bertiger, W.; Sakumura, C.; Bandikova, T.; Mccullough, C. Technical Report JPL D-56935: Grace-Fo Level-1 Data Product User Handbook; Jet Propulsion Laboratory: Pasadena, CA, USA, 2019. [Google Scholar]
  45. Tapley, B.D.; Bettadpur, S.; Watkins, M.M.; Reigber, C. The gravity recovery and climate experiment: Mission overview and early results. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef] [Green Version]
  46. Dragon, K. Technical Report JPL D-71781: GRACE Follow-On Mission Plan; Jet Propulsion Laboratory: Pasadena, CA, USA, 2015. [Google Scholar]
  47. Liu, J.; Ge, M. PANDA software and its preliminary result of positioning and orbit determination. Wuhan Univ. J. Nat. Sci. 2003, 8, 603–609. [Google Scholar]
  48. Berger, C.; Biancale, R.; Ill, M.; Barlier, F. Improvement of the empirical thermospheric model DTM: DTM94-A comparative review of various temporal variations and prospects in space geodesy applications. J. Geod. 1998, 72, 161–179. [Google Scholar] [CrossRef]
  49. Förste, C.; Bruinsma, S.L.; Shako, R.; Marty, J.C.; Flechtner, F.; Abrikosov, O.; Dahle, C.; Lemoine, J.M.; Neumayer, K.H.; Biancale, R.; et al. EIGEN-6—A new combined global gravity field model including GOCE data from the collaboration of GFZ-Potsdam and GRGS Toulouse. In Proceedings of the EGU General Assembly, Vienna, Austria, 3–8 April 2011. [Google Scholar]
  50. Petit, G.; Luzum, B.; Al, E. IERS Conventions. IERS Tech. Note 2010, 36, 1–95. [Google Scholar]
  51. Lyard, F.; Lefevre, F.; Letellier, T.; Francis, O. Modeling the global ocean tides: Modern insights from FES2004. Ocean Dyn. 2006, 56, 394–415. [Google Scholar] [CrossRef]
  52. Desai, S.D. Observing the pole tide with satellite altimetry. J. Geophys. Res. 2002, 107, 1–13. [Google Scholar] [CrossRef]
  53. Wu, J.T.; Wu, S.C.; Hajj, G.A.; Bertiger, W.I.; Lichten, S.M. Effects of antenna orientation on GPS carrier phase. Manuscr. Geod. 1993, 18, 91–98. [Google Scholar]
  54. Meehan, T.; Esterhuizen, S.; Franklin, G.; Tien, J.; Young, L.; Bachman, B.; Munson, T.; Robison, D.; Stecheson, T. Development status of NASA’s TriG GNSS science instrument. In Proceedings of the International Radio Occultation Working Group (IROWG-2), Estes Park, CO, USA, 28 March–3 April 2012. [Google Scholar]
  55. Schaer, S.; Villiger, A.; Arnold, D.; Dach, R.; Prange, L.; Jäggi, A. The CODE ambiguity-fixed clock and phase bias analysis products: Generation, properties, and performance. J. Geod. 2021, 95, 1–25. [Google Scholar] [CrossRef]
  56. Villiger, A.; Schaer, S.; Dach, R.; Prange, L.; Sušnik, A.; Jäggi, A. Determination of GNSS pseudo-absolute code biases and their long-term combination. J. Geod. 2019, 93, 1487–1500. [Google Scholar] [CrossRef]
  57. Blewitt, G. Carrier phase ambiguity resolution for the Global Positioning System applied to geodetic baselines up to 2000 km. J. Geophys. Res. Solid Earth 1989, 94, 10187–10203. [Google Scholar] [CrossRef] [Green Version]
  58. Jin, B.; Chen, S.; Li, D.; Takka, E.; Li, Z.; Qu, P. Ionospheric correlation analysis and spatial threat model for SBAS in China region. Adv. Space Res. 2020, 66, 2873–2887. [Google Scholar] [CrossRef]
  59. Teunissen, P.J.G. GNSS Precise Point Positioning. In Position, Navigation, and Timing Technologies in the 21st Century: Integrated Satellite Navigation, Sensor Systems, and Civil Applications, 1st ed.; Morton, Y.J., van Diggelen, F., Spilker, J.J., Jr., Parkinson, B.W., Lo, S., Gao, G., Eds.; Wiley-IEEE Press: Hoboken, NJ, USA, 2021; Volume 1, pp. 503–528. [Google Scholar]
  60. Hatch, R. The synergism of GPS Code and carrier measurements. In Proceedings of the Third International Symposium on Satellite Doppler Positioning, Physical Sciences Laboratory of New Mexico State University, Las Cruces, NM, USA, 8–12 February 1982; pp. 1213–1231. [Google Scholar]
  61. Melbourne, W. The Case for Ranging in GPS Based Geodetic Systems. In Proceedings of the 1st International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15–19 April 1985; pp. 373–386. [Google Scholar]
  62. Wubbena, G. Software developments for geodetic positioning with GPS using TI 4100 code and carrier measurements. In Proceedings of the 1st International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15–19 April 1985; pp. 403–412. [Google Scholar]
  63. Geng, J.; Chen, X.; Pan, Y.; Zhao, Q. A modified phase clock/bias model to improve PPP ambiguity resolution at Wuhan University. J. Geod. 2019, 93, 2053–2067. [Google Scholar] [CrossRef]
  64. Teunissen, P.J. Carrier phase integer ambiguity resolution. In Springer Handbook of Global Navigation Satellite Systems, 1st ed.; Teunissen, P., Montenbruck, O., Eds.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 661–685. [Google Scholar]
  65. Dong, D.; Bock, Y. Global positioning system network analysis with phase ambiguity resolution applied to crustal deformation studies in California. J. Geophys. Res. 1989, 94, 3949–3966. [Google Scholar] [CrossRef]
  66. Ge, M.; Gendt, G.; Dick, G.; Zhang, F.P. Improving carrier-phase ambiguity resolution in global GPS network solutions. J. Geod. 2005, 79, 103–110. [Google Scholar] [CrossRef]
  67. Teunissen, P.J.G. The parameter distributions of the integer GPS model. J. Geod. 2002, 76, 41–48. [Google Scholar] [CrossRef] [Green Version]
  68. Kang, Z.; Bettadpur, S.; Nagel, P.; Save, H.; Poole, S.; Pie, N. GRACE-FO precise orbit determination and gravity recovery. J. Geod. 2020, 94, 85. [Google Scholar] [CrossRef]
  69. Gu, D.; Ju, B.; Liu, J.; Tu, J. Enhanced GPS-based GRACE baseline determination by using a new strategy for ambiguity resolution and relative phase center variation corrections. Acta Astronaut. 2017, 138, 176–184. [Google Scholar] [CrossRef]
  70. Geng, J.; Yang, S.; Guo, J. Assessing IGS GPS/Galileo/BDS-2/BDS-3 phase bias products with PRIDE PPP-AR. Satell. Navig. 2021, 2, 1–15. [Google Scholar] [CrossRef]
  71. Montenbruck, O.; Hackel, S.; Wermuth, M.; Zangerl, F. Sentinel-6A precise orbit determination using a combined GPS/Galileo receiver. J. Geod. 2021, 95, 109. [Google Scholar] [CrossRef]
Figure 1. Appearance and instrument installation position of GRACE-FO satellites.
Figure 1. Appearance and instrument installation position of GRACE-FO satellites.
Remotesensing 13 04204 g001
Figure 2. The 5° × 5° PCV corrections of (a) GRCC and (b) GRCD satellites based on ionosphere free carrier phase residuals from the DD AR solution.
Figure 2. The 5° × 5° PCV corrections of (a) GRCC and (b) GRCD satellites based on ionosphere free carrier phase residuals from the DD AR solution.
Remotesensing 13 04204 g002
Figure 3. The 5° × 5° PCV corrections of (a) GRCC and (b) GRCD satellites based on ionosphere free carrier phase residuals from SD AR solution.
Figure 3. The 5° × 5° PCV corrections of (a) GRCC and (b) GRCD satellites based on ionosphere free carrier phase residuals from SD AR solution.
Remotesensing 13 04204 g003
Figure 4. Carrier phase residuals of GRCC satellite derived from (a) DD AR solution and (b) SD AR solution with application of PCV corrections.
Figure 4. Carrier phase residuals of GRCC satellite derived from (a) DD AR solution and (b) SD AR solution with application of PCV corrections.
Remotesensing 13 04204 g004
Figure 5. Residual distributions of (a) SD WL, (b) SD NL, (c) DD WL, and (d) DD NL ambiguities of GRACE-FO satellites.
Figure 5. Residual distributions of (a) SD WL, (b) SD NL, (c) DD WL, and (d) DD NL ambiguities of GRACE-FO satellites.
Remotesensing 13 04204 g005
Figure 6. Fixing rates of (a) SD WL and NL ambiguities of GRCC and GRCD satellites and (b) DD WL and NL ambiguities.
Figure 6. Fixing rates of (a) SD WL and NL ambiguities of GRCC and GRCD satellites and (b) DD WL and NL ambiguities.
Remotesensing 13 04204 g006
Figure 7. Daily RMSs of along-track orbit differences between JPL’s orbit and SD AR solution with (red dots) and without (blue dots) PCV corrections. The upper panel is for the GRCC satellite and lower panel is for the GRCD satellite.
Figure 7. Daily RMSs of along-track orbit differences between JPL’s orbit and SD AR solution with (red dots) and without (blue dots) PCV corrections. The upper panel is for the GRCC satellite and lower panel is for the GRCD satellite.
Remotesensing 13 04204 g007
Figure 8. Daily RMSs of along-track orbit differences of FA, DD AR, and SD AR solutions with JPL’s orbit. The upper panel is for the GRCC satellite and the lower panel is for the GRCD satellite.
Figure 8. Daily RMSs of along-track orbit differences of FA, DD AR, and SD AR solutions with JPL’s orbit. The upper panel is for the GRCC satellite and the lower panel is for the GRCD satellite.
Remotesensing 13 04204 g008
Figure 9. RMSs of SLR residuals derived from FA (red dots), DD AR (green dots), and SD AR (blue dots) solutions. The upper panel is for the GRCC satellite and the lower panel is for the GRCD satellite.
Figure 9. RMSs of SLR residuals derived from FA (red dots), DD AR (green dots), and SD AR (blue dots) solutions. The upper panel is for the GRCC satellite and the lower panel is for the GRCD satellite.
Remotesensing 13 04204 g009
Figure 10. Daily STDs of KBR residuals of (a) FA orbit with PCV corrections (red dots), SD AR orbit without PCV corrections (green dots), and SD AR orbit with PCV corrections (blue dots); (b) FA orbit with PCV corrections (red dots), DD AR orbit without PCV corrections (green dots), and DD AR orbit with PCV corrections (blue dots).
Figure 10. Daily STDs of KBR residuals of (a) FA orbit with PCV corrections (red dots), SD AR orbit without PCV corrections (green dots), and SD AR orbit with PCV corrections (blue dots); (b) FA orbit with PCV corrections (red dots), DD AR orbit without PCV corrections (green dots), and DD AR orbit with PCV corrections (blue dots).
Remotesensing 13 04204 g010
Table 1. POD strategy of GRACE-FO satellites.
Table 1. POD strategy of GRACE-FO satellites.
ParameterDescription
Background force models
Static gravity field modelEIGEN_06C (130 × 130) [49]
Solid earth and pole tidesIERS Conventions 2010 [50]
Ocean tidesFES2004 (30 × 30) [51]
Ocean pole tidesDesai [52] (30 × 30)
Third body perturbationsJPL’s DE405
General relativistic effectsIERS Conventions 2010 [50]
Solar radiation pressureMacro model [44]
Atmospheric dragMacro model [44], DTM94 [48]
Input data and products
GPS ObservationsUndifferenced IF code and carrier phase
Sample interval10 s
Elevation mask
GPS orbitCODE’s final GPS orbit
GPS clock and hardware biasCODE’s 5 s clock and OSB products
GPS antenna correctionigs14.atx
GRACE-FO antenna correctionPCO applied [44]
Phase wind-upApplied [53]
Gravitational bendingIERS Conventions 2010 [50]
Relativistic correctionIERS Conventions 2010 [50]
Estimated parameters
Initial state vectorPosition and velocity per satellite per arc
Atmospheric drag coefficientOne per orbital resolution
Scale coefficient of solar radiation pressureOne per satellite per arc
Empirical along- and cross-track accelerations1-CPR accelerations per orbital resolution
Receiver clock errorsEpoch wise
Carrier phase ambiguitiesOne per satellite tracking pass
Table 2. GPS antenna PCOs and SLR reflector Coordinates of GRACE-FO satellites.
Table 2. GPS antenna PCOs and SLR reflector Coordinates of GRACE-FO satellites.
SatelliteGPS Antenna PCO/mSLR Reflector/m
XYZXYZ
GRCC0.26023−0.00128−0.476970.600000.327500.22080
GRCD0.26004−0.00107−0.476180.600000.327500.22080
Table 3. Orbit improvements with application of PCV maps for SD AR and DD AR solutions.
Table 3. Orbit improvements with application of PCV maps for SD AR and DD AR solutions.
SolutionSatelliteOrbitNo PCV/mmWith PCV/mmImprovement
SD AR
Solution
GRCCAlong-track7.26.016.7%
Cross-track6.26.20
Radial5.25.11.9%
GRCDAlong-track6.45.97.8%
Cross-track6.36.13.2%
Radial5.15.10
DD AR
Solution
GRCCAlong-track9.58.312.6%
Cross-track6.56.41.5%
Radial5.85.63.4%
GRCDAlong-track8.88.18.0%
Cross-track6.46.31.6%
Radial5.75.53.5%
Table 4. RMSs of differences between the FA, DD AR and SD AR orbits and JPL’s orbit.
Table 4. RMSs of differences between the FA, DD AR and SD AR orbits and JPL’s orbit.
SatelliteOrbitSolutions/mmImprovement
FADD ARSD ARDD ARSD AR
GRCCAlong-track8.88.36.05.7%31.8%
Cross-track6.76.46.24.5%7.5%
Radial6.05.65.16.7%15.0%
GRCDAlong-track8.88.15.98.0%33.0%
Cross-track6.76.36.16.0%9.0%
Radial5.95.55.16.8%13.6%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jin, B.; Li, Y.; Jiang, K.; Li, Z.; Chen, S. GRACE-FO Antenna Phase Center Modeling and Precise Orbit Determination with Single Receiver Ambiguity Resolution. Remote Sens. 2021, 13, 4204. https://doi.org/10.3390/rs13214204

AMA Style

Jin B, Li Y, Jiang K, Li Z, Chen S. GRACE-FO Antenna Phase Center Modeling and Precise Orbit Determination with Single Receiver Ambiguity Resolution. Remote Sensing. 2021; 13(21):4204. https://doi.org/10.3390/rs13214204

Chicago/Turabian Style

Jin, Biao, Yuqiang Li, Kecai Jiang, Zhulian Li, and Shanshan Chen. 2021. "GRACE-FO Antenna Phase Center Modeling and Precise Orbit Determination with Single Receiver Ambiguity Resolution" Remote Sensing 13, no. 21: 4204. https://doi.org/10.3390/rs13214204

APA Style

Jin, B., Li, Y., Jiang, K., Li, Z., & Chen, S. (2021). GRACE-FO Antenna Phase Center Modeling and Precise Orbit Determination with Single Receiver Ambiguity Resolution. Remote Sensing, 13(21), 4204. https://doi.org/10.3390/rs13214204

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop