Next Article in Journal
The Design and Performance Evaluation of a 1550 nm All-Fiber Dual-Polarization Coherent Doppler Lidar for Atmospheric Aerosol Measurements
Next Article in Special Issue
Downlink Analysis of a Low-Earth Orbit Satellite Considering an Airborne Interference Source Moving on Various Trajectory
Previous Article in Journal
TransHSI: A Hybrid CNN-Transformer Method for Disjoint Sample-Based Hyperspectral Image Classification
Previous Article in Special Issue
LEO Satellite Clock Modeling and Its Benefits for LEO Kinematic POD
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Topside Ionosphere Modeling Parameters on Differential Code Bias (DCB) Estimation Using LEO Satellite Observations

1
Joint Laboratory of Power Remote Sensing Technology, Electric Power Research Institute, Yunnan Power Grid Co., Ltd., Kunming 650217, China
2
State Key Laboratory of Geodesy and Earth’s Dynamics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430077, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(22), 5335; https://doi.org/10.3390/rs15225335
Submission received: 19 October 2023 / Revised: 9 November 2023 / Accepted: 10 November 2023 / Published: 13 November 2023
(This article belongs to the Special Issue LEO-Augmented PNT Service)

Abstract

:
Given the potential of low-earth orbit (LEO) satellites in terms of navigation enhancement, accurately estimating the differential code bias (DCB) of GNSS satellites and LEO satellites is an important research topic. In this study, to obtain accurate DCB estimates, the effects of vertical total electron content (VTEC) modeling parameters of the topside ionosphere on DCB estimation were investigated using LEO observations for the first time. Different modeling parameters were set in the DCB estimations, encompassing modeling spacing in the dynamic temporal mode and degree and order (D&O) in spherical harmonic modeling. The DCB precisions were then evaluated, and the impacts were analyzed. Thus, a number of crucial and beneficial conclusions are drawn: (1) The maximum differences in the GPS DCB estimates after adopting different modeling spacings and different D&Os exhibit that the different modeling spacings or D&Os both affect the GPS DCB estimates and their root-mean square (RMS), and the effects of the two are at the same level. (2) The maximum differences in receiver DCBs using different modeling spacings indicate that the modeling spacing has a significant impact on the receiver DCBs, compared with GPS DCBs. Whereas, the maximum differences in receiver DCBs with different modeling D&Os are inferior to the differences in the GPS DCBs. That is, the modeling spacing has a greater impact on the LEO DCBs than those of the modeling D&O. (3) The experimental results indicate that the GPS DCB estimates using a modeling spacing of 12H have higher precisions than the others, whereas LEO receiver DCBs using a spacing of 4H or 6H obtain optimal STD. In terms of modeling D&O, adopting 8D&O in the LEO-based VTEC modeling can attain superior estimates and precisions for both GPS and LEO DCBs. The research conclusions can provide references for LEO-augmented DCB estimation.

1. Introduction

Differential code bias (DCB) is a critical error source in navigation positioning and ionosphere modeling [1,2,3,4,5,6,7,8]. Several DCB estimation methods based on ground station data are available: the first method, used by a few institutes such as the Center of Orbit Determination of Europe (CODE), uses observation data from global ground stations to estimate DCB and vertical total electron content (VTEC) simultaneously [1,4,9,10,11,12,13]; the second method proposed by the Chinese Academy of Science (CAS) applies single ground station observations to model the VTEC and estimate the DCB and VTEC parameters simultaneously [2,7,14,15]; the third method is similar to the one used by the German Aerospace Center (DLR) and estimates the DCB parameters of GNSS satellite and receiver after removing the ionosphere impact by introducing a more accurate prior ionosphere model [5,6]. The accuracy of DCB estimates depends critically on the amount of observation data and their distribution. Low-earth orbit (LEO) satellites can be used as space-based monitoring stations to estimate DCB parameters, compensating for the lack of ground stations. In addition, as LEO is above the F layer of the ionosphere and ionized electrons are derived from the topside ionosphere or plasmasphere, which is less variable in both spatial and temporal domains, the observations suffer from minor signal delays.
With the launch of tens of thousands of LEO satellites in the future, accurately estimating the DCB of GNSS and LEO satellites is an important research topic. In recent years, scholars have studied and used LEO onboard observation data to estimate DCB parameters. A few scholars estimated the DCBs of the GPS satellite and LEO receiver as unknown parameters simultaneously [16,17,18], and others [19,20,21,22,23,24,25] not only estimated the DCB parameters of GPS and LEO satellites but also estimated the corresponding VTEC parameters. Certain scholars introduced external GPS DCB products and estimated only the corresponding topside ionosphere VTEC and LEO receiver DCB parameters [26,27,28]. Presently, the aforementioned scholars have conducted research on DCB estimation employing LEO observation data by choosing different modeling parameters [20,24,25,28]. However, there is no research on the effects of topside ionosphere estimation strategies on DCB estimation, and superior modeling parameters must be determined to obtain accurate DCB estimates.
In this paper, the impacts of topside ionosphere VTEC modeling parameters on LEO-based DCB estimation are investigated and analyzed for the first time to obtain superior modeling parameters and accurate DCB estimates. The research conclusions can provide a reference for LEO-augmented DCB estimation. Section 2 presents the model and strategy for GPS and LEO DCB estimation using LEO observation data. Section 3 analyzes the results of experiments in which different modeling parameters are set in the DCB estimations, encompassing the modeling spacing in the dynamic temporal mode and the degree and order in spherical harmonic modeling, employing GRACE-FO [29] observation data. The differences in the GPS and LEO DCB estimates using different modeling parameters are shown, and the precision of the DCB estimates is evaluated. Then, the effects of the modeling parameters on DCB estimation are analyzed, and superior modeling parameters for DCB estimation are provided. Finally, the conclusions are summarized in Section 4.

2. Model and Strategy for DCB Estimation

The geometry-free (GF) combination observation estimation method involves a simpler calculation and does not require outlier information. Therefore, we employed GF combinations of pseudo-range observations to estimate the DCB parameters. The DCB values (P1-P2) for the GPS satellites and LEO receivers were estimated daily as constant values and simultaneously with the topside ionosphere VTEC parameters.
DCB estimation using LEO observations is not affected by the troposphere. Hence, dual-frequency code observations are commonly expressed as shown in Equation (1), and the GF combination of the pseudo-range observations is formed as expressed in Equation (2) [1,6,7].
P r , 1 s = ρ r s + c ( d t r d t s ) + α 1 S T E C + c ( b s , 1 + b r , 1 ) + ε r , 1 s   ,   α 1 = 40.28 f 1 2 P r , 2 s = ρ r s + c ( d t r d t s ) + α 2 S T E C + c ( b s , 2 + b r , 2 ) + ε r , 2 s   ,   α 2 = 40.28 f 2 2
P r , G F s = P r , 1 s P r , 2 s = α S T E C + D C B s + D C B r + ε r , G F s ,   α = 40.28 1 f 1 2 1 f 2 2
Here scripts s and r denote the GPS satellite and LEO receiver; 1 and 2 denote frequency numbers; f 1   and   f 2 denote the L1 and L2 frequencies of the GPS signal, respectively; ρ r s represents the geometric distance; c is the light speed in vacuum; d t r   and   t s represent the LEO receiver and GPS satellite clock offsets, respectively; bs,i and br,i (i = 1,2) refer to the instrument delays from GPS satellite s and LEO receiver r at two frequencies, respectively; STEC represents the slant total electron content of the LEO-based ionosphere; α1, α2, and α are the coefficients of the STEC related to the L1 and L2 frequencies of the GPS signal; DCBs and DCBr refer to the DCBs of the GPS satellite s and LEO receiver r; ε r , i s and ε r , G F s represent the noise from pseudo-range observations and the GF combination observations, respectively. P r , 1 s , P r , 2 s denote the pseudo-range observations at the L1 and L2 frequencies of the GPS signals from GPS satellite s to LEO receiver r; and P r , G F s represents the GF combination of pseudo-range observations.
In this study, we screen code observation data using residuals and LEO orbits derived from phase observations, and the preprocessing for pseudo-range also draws on the onboard data preprocessing method of LEO precise orbit determination. The pseudo-range is relatively clean in this way.
The F&K mapping function [30] was originally developed for slant-path atmospheric water vapor conversion and later applied to LEO-based total electron content (TEC) conversion [31]. It is more suitable than the single-layer mapping function (SLM) for LEO-based TEC conversion under certain conditions [32]. In this study, the F&K mapping function was applied to convert the LEO-based TEC from a vertical to a slant direction. The expected height of maximum electron density is called the ionospheric effective height (IEH) [32]. In the solar-geomagnetic reference frame, the LEO-based VTEC values were modeled using spherical harmonic expansion. The temporal modeling mode was set to dynamic status. The modeling spacing means the modeling interval, and it is the calculation period of the modeling parameters. A set of modeling coefficients can be estimated within each modeling spacing. The modeling spacing of the VTEC and the degree and order (D&O) of the spherical harmonic expansion are the main research objects. The sampling interval of the observation data was set at 30 s, and the cut-off elevation angle was set at 15° to reduce the effects of multipath. The LEO-based TEC model is expressed in Equation (3) [1,30,32,33].
S T E C = m f F & K V T E C m f F & K = 1 + ( R E + h I E H ) / R L E O cos z + ( R E + h I E H ) 2 / R L E O 2 ( sin z ) 2 h I E H = ( 0.0027 F 107 + 1.79 ) h L E O 5.52 F 107 + 1350 m o d e l e d   V T E C = n = 0 n max m = 0 n P ˜ n m ( sin φ ) ( A ˜ n m cos ( m λ ) + B ˜ n m sin ( m λ ) )
where m f F & K represents the F&K mapping function; h I E H denotes IEH; R E represents the earth radius; h L E O refers to the altitude of the LEO satellite above the surface of the Earth; R L E O = R E + h L E O ; F 107 denotes the solar radio flux at 10.7 cm; z represents the zenith angle of the slant ray path. STEC and VTEC denote the slant and vertical total electron content of the LEO-based topside ionosphere, respectively, φ and λ represent the geomagnetic latitude and sun-fixed longitude in the equation of modeled VTEC, respectively; P ˜ n m denotes the normalized associated Legendre functions of degree n and order m; and A ˜ n m and B ˜ n m refer to the VTEC coefficients of the spherical harmonic function; n max represents the maximum degree of the spherical harmonic expansion.
Introduced Equations (3) into (2), unknown estimation parameters contain DCBs of GNSS satellite and LEO receiver and LEO-based VTEC model coefficients. Equation (4) is obtained from Equations (2) and (3). The GPS and LEO DCBs can be estimated and determined together with the VTEC model coefficients. Assuming that the numbers of GPS and LEO DCB estimates are μ1 and μ2, then the number of estimated VTEC model coefficients is μ3. The parameter estimation equation can be written in matrix form as:
Z n × 1 = F X ^ F n × ( μ 1 + μ 2 + μ 3 ) = A n × μ 1   B n × μ 2   C n × μ 3 X ^ ( μ 1 + μ 2 + μ 3 ) × 1 = X ^ s , d c b T 1 × μ 1   X ^ r , d c b T 1 × μ 2   X ^ v t e c , c o e f T 1 × μ 3 T
Here the vector Z with a total of n rows is related to DCB and VTEC modeling coefficients; F is the design matrix that consists of the matrix A, B, and C, related to GNSS DCB, LEO DCB, and VTEC modeling coefficients; and X ^ is the DCB estimate and VTEC modeling coefficient, which consists of the vectors of GNSS DCB ( X ^ s , d c b ), LEO DCB ( X ^ r , d c b ), and VTEC modeling coefficients ( X ^ v t e c , c o e f ). Assuming that the modeling spacing and D&O are s (in hours) and n max , respectively, the number of VTEC modeling parameters is ( n max + 1) × ( n max + 1) × (24/s + 1), the number of DCB parameters estimated daily as constant values for each GPS satellite and each LEO receiver is 32 + 1 = 33, and the number of total unknown parameters on one day is ( n max + 1) × ( n max + 1) × (24/s + 1) + 33.
Considering that the DCBs of the GPS satellite and LEO receiver are closely correlated, the DCB datum is defined by a zero-mean condition in the DCB estimation for de-correlation. Equation (5) represents the zero-mean condition [34]:
s = 1 S max D C B s = 0
Here S max denotes the total number of observed GPS satellites.
Daily DCB values were realigned by applying a shift from non-all satellite DCB values, which were calculated using a common set of satellites over the study period. After alignment [35], we evaluated the DCB estimates using their monthly stability to represent internal agreement as a standard deviation (STD) and compared them with CODE products to represent external agreement as mean differences and root-mean square (RMS). The monthly stability of the DCB values reflects the stability and reliability of the DCB estimates to a certain extent and can be expressed as follows:
S T D s = d = 1 D ( d c b d s d c b ¯ s ) 2 D 1
Here scripts s and d denote one satellite and one day, respectively; D refers to the total days of a month; d c b ¯ s represents the monthly mean of the DCB estimates of satellite s; d c b d s denotes the DCB estimate of satellite s on day d in a month; and S T D s is the monthly stability of the DCB estimates of satellite s.
The RMS values for GPS DCB estimates with respect to external reference products after conducting the alignment procedures is calculated as follows:
R M S s = d = 1 D ( d c b d s e d c b d s ) 2 D
Here d c b d s represents the DCB estimate of satellite s on day d in a month; e d c b d s denotes the DCB values of satellite s from the external reference products on day d; D refers to the total number of days in one month; R M S s is the difference RMS of the DCB estimates of satellite s relative to the reference products.

3. Experiments and Results

To obtain suitable and superior topside ionospheric VTEC modeling parameters in LEO-based DCB estimation, we investigated the effects of the different modeling parameters on DCB estimates and precision, including modeling spacings in dynamic temporal mode and degrees and orders in spherical harmonic modeling. The precision of DCB estimates was then evaluated, and the effects of the parameters were analyzed. We used onboard GPS observation data from the twin GRACE-FO satellites (GRCC and GRCD) from 1–30 January 2019 to estimate the GPS and LEO satellite DCBs. We comprehensively evaluated the internal and external consistency precisions of GPS DCB estimation and compared the monthly stability of GPS and LEO DCB estimates as internal precision, which was presented STD. The comparisons between the GPS DCB estimates and CODE DCB products represent the external precision in the form of the mean and RMS of the difference between the estimates and CODE products. The LEO DCB can provide internal precision but not external precision owing to a lack of external reference products.
Additionally, this study requires the precise orbits of the twin GRACE-FO satellites owing to the needs of Equation (3). Before DCB estimation, we conducted precise orbit determination (POD) of the GRACE-FO satellites. The main data processing of LEO POD can be found in the references [36,37,38,39,40,41]. The LEO orbit determination is not necessary. If external reference orbits are available, they can also be used to carry out experiments.

3.1. Different Modeling Spacings

To investigate the effects of modeling spacings on DCB estimate and precision and to obtain superior parameter settings, we set different modeling spacings of 2H, 4H, 6H, and 12H in the LEO-based VTEC dynamic temporal modeling mode. The modeling degree and order were set to 8. Then, the effects of modeling spacing on DCB estimation were analyzed by evaluating the precision of the GPS and LEO DCB estimates.

3.1.1. GPS Satellite DCB Estimates

The mean DCB can exhibit the averaged ranges of DCBs in the form of specific and intuitive numbers; thereby, Table 1 presents the mean values of GPS satellite DCB estimates after adopting different modeling spacings of 2H, 4H, 6H, and 12H using the GRCC (a) and GRCD (b) observations. The GPS DCB values for G04 were not estimated owing to a lack of observations. The mean GPS DCB estimates range between −10 and 10 ns and are stable. There is no significant difference in mean GPS DCBs with different modeling spacings. Considering the limitation of mean statistics in Table 1 in terms of reflecting all DCB values, Figure 1 displays the maximum differences in the GPS satellite DCB estimates after adopting different modeling spacings using GRCC and GRCD observations during the experiment period. The differences in the GPS DCB estimates using different modeling spacings are in the range of 0.05 ns. Specifically, in Figure 1, the difference in G15 DCB estimates adopting different modeling spacings is the maximum value, 0.04 ns, when using GRCC data, whereas for GRCD data, the difference in G12 DCB is the largest, approximately 0.05 ns. Therefore, corresponding to Figure 1, the maximum, minimum, and differences in G15 and G12 DCB estimated based on different modeling spacings using GRCC and GRCD observations, respectively, are presented in Table 2. In Table 2a using GRCC data, the difference in the G15 DCB estimates based on different modeling spacings on DOY 5 is the largest, and the maximum, minimum, and difference are 2.6815, 2.6415, and 0.0400 ns, respectively. In Table 2b using GRCD data, the difference in G12 DCB estimated on DOY 4 is the largest, and the maximum, minimum, and difference are 4.3340, 4.2874, and 0.0466 ns, respectively. Therefore, there is a difference in the GPS DCB estimated using different modeling spacings, within 0.05 ns. The different modeling spacings affect the GPS DCB estimates.

3.1.2. Precision Evaluation of GPS Satellite DCB Estimates

Figure 2 presents the monthly stabilities (STD) of the GPS satellite DCB estimates adopting different modeling spacings of 2H, 4H, 6H, and 12H using GRCC (top) and GRCD (bottom) observation data. The STDs of the GPS DCB estimates using GRCC and GRCD data are in the range of 0.11 ns with high stability. Meanwhile, the maximum differences in the STD of the GPS DCB estimates adopting different modeling spacings are shown in Figure 3. The maximum differences in the STD of the GPS DCB using different modeling spacings are within 0.006 ns. The STD values of the GPS DCB estimates adopting different modeling spacings exhibit no significant differences.
Table 3 presents the mean STDs of the GPS DCB estimates adopting the different modeling spacings. Statistically, the mean STD values of GPS DCBs with modeling spacings of 2H, 4H, 6H, and 12H using GRCC and GRCD data are 0.0624, 0.0621, 0.0620, and 0.0619 ns and 0.0648, 0.0643, 0.0641, and 0.0639 ns, respectively. The GPS DCB estimates with modeling spacings of 12H have the optimal mean STD results. The GPS DCB estimates using twin LEO satellite data attain similar and high stability.
Figure 4 depicts the mean differences for GPS DCB estimates with different modeling spacings relative to CODE products. The mean differences in Figure 4 vary between −0.6 and 0.8 ns, which indicates that the GPS DCBs estimated by LEO satellite and ground station data have good consistency. Figure 5 showcases the RMS of the differences between the GPS DCB estimates adopting different modeling spacings (2H, 4H, 6H, and 12H) and the CODE DCB products. To visually demonstrate the differences in RMS of GPS DCBs with different modeling spacings, Figure 6 displays the maximum differences in RMS of GPS DCB estimates using different modeling spacings. The maximum differences in RMS of GPS DCBs using different modeling spacings are in the range of 0.03 ns. The different modeling spacings affect the RMS of the GPS DCB estimates. The corresponding mean RMS values are listed in Table 4. The RMS results of the GPS DCBs with different modeling spacings relative to the CODE products are in the range of 0.8 ns. The mean RMSs of the GPS DCBs are in the range of 0.3 ns, which indicates that the GPS DCBs estimated by LEO satellite and ground station data have good consistency. Statistically, the mean RMS results in Table 4 indicate that the GPS DCBs with a modeling spacing of 12H had slightly smaller RMS values than the others. Additionally, the GPS DCB estimates using twin satellite data have similar precisions, and the RMS results using GRCC observations are slightly poorer than those obtained using GRCD data.

3.1.3. LEO Receiver DCB Estimates and Stability Evaluation

Figure 7 presents the time series of the receiver DCB estimates with modeling spacings of 2H, 4H, 6H, and 12H for the twin GRACE-FO satellites. The GRCC receiver DCB estimates are located at approximately −1.0 ns, whereas the GRCD DCB estimates fluctuate at approximately 4.6 ns, and their receiver DCB estimates are stable. In Figure 7, the GRCC and GRCD receiver DCB estimates are sorted in descending order as follows: receiver DCBs with modeling spacings of 12H, 6H, 4H, and 2H, where the receiver DCBs decrease as the modeling spacing decreases.
In order to visually demonstrate the differences in receiver DCBs with different modeling spacings, Figure 8 exhibits the maximum differences in the LEO receiver DCB estimates using different modeling spacings. Figure 9 shows the maximum differences in the STD of the receiver DCB using different modeling spacings. The maximum difference in receiver DCBs using different modeling spacings is 0.22 ns, whereas the maximum differences in STD of receiver DCBs are within 0.005 ns. This indicates that different modeling spacings have a significant impact on receiver DCB estimates.
The mean values and STD results for the LEO receiver DCB estimates with different modeling spacings are listed in Table 5. Statistically, the maximum receiver DCB values for both the GRCC and GRCD are the DCBs with modeling spacings of 12H, whereas the minimum receiver DCBs are those with modeling spacings of 2H. The differences in the mean receiver DCBs with different modeling spacings for the GRCC and GRCD are in the ranges of 0.1569 and 0.1684 ns, respectively. The GRCC and GRCD receiver DCBs achieve optimal STD results when modeling spacings of 6H and 4H are applied. The receiver DCBs of the twin GRACE-FO satellites are different; however, their DCB estimates have similar STD results.
In summary, the GPS DCB estimates using a modeling spacing of 12H have higher precision than the others, whereas LEO receiver DCBs applying the modeling spacings of 4H or 6H obtain optimal STD.

3.2. Different Modeling Degrees and Orders

This section investigates the impact of the D&O of spherical harmonic modeling on DCB estimates and precision, thereby obtaining suitable modeling D&O parameters. The 6, 8, and 10D&Os of spherical harmonic modeling were introduced into topside ionosphere modeling and DCB estimation. Subsequently, the GPS and LEO DCB estimates were analyzed and evaluated. The modeling spacing was fixed at 4H to improve the temporal resolution.

3.2.1. GPS Satellite DCB Estimates

The mean values of GPS DCB estimates adopting different modeling parameters of 6, 8, and 10D&Os using the GRCC (a) and GRCD (b) observations are listed in Table 6. The GPS DCB values for G04 were not estimated owing to a lack of observations. The mean GPS DCB estimates range between −10 and 10 ns and are stable. There is no significant difference in the GPS DCB estimates obtained using different modeling D&Os. Considering the limitation of mean statistics in Table 6, Figure 10 exhibits the maximum differences in GPS DCB estimates adopting different modeling D&Os using GRCC and GRCD data. The differences in the GPS DCB estimates using different modeling D&Os are in the range of 0.05 ns. It indicates that the effects of the modeling spacing and D&O on GPS DCB estimates are at the same level. Specifically, the differences in G05 DCB estimates adopting different modeling D&Os using GRCC and GRCD data are both the largest, within 0.04 and 0.05 ns, respectively. The maximum, minimum, and differences in G05 DCB estimated with different modeling D&Os using GRCC and GRCD observations are presented in Table 7. In Table 7a, using GRCC data, the difference in G05 DCB estimates between different modeling D&Os on DOY 13 is the largest, and the maximum, minimum, and difference are 3.0099, 2.9776, and 0.0323 ns, respectively. In Table 7b, using GRCD data, the difference in G12 DCB estimated on DOY 4 is the largest, and the maximum, minimum, and difference are 3.0515, 3.0048, and 0.0467 ns, respectively. Therefore, there is a certain difference in the GPS DCB estimates using different modeling D&Os, all within 0.05 ns. The different modeling D&Os have certain effects on the GPS DCB estimates.

3.2.2. Precision Evaluation of GPS DCB Estimates

Figure 11 showcases the STD of GPS DCB estimates adopting different modeling D&Os of 6, 8, and 10 using GRCC and GRCD data. The STDs of GPS DCB estimates based on GRCC and GRCD data are within 0.11 ns. Additionally, Figure 12 displays the maximum differences in the STD of the GPS DCB estimates by applying different modeling D&Os. The maximum differences are within 0.005 ns. The STDs of the GPS DCB estimates adopting different modeling D&Os show no marked differences.
Table 8 presents the mean STDs of GPS satellite DCB estimates by applying different modeling D&Os. Statistically, the mean STD values of GPS DCBs with modeling D&Os of 6, 8, and 10 using GRCC and GRCD data are 0.0620, 0.0621, and 0.0623 ns and 0.0643, 0.0643, and 0.0644 ns, respectively. The GPS DCB estimates, with modeling D&Os of 6 and 8, attain the optimal STD results.
Figure 13 exhibits the mean differences for GPS DCB estimates with different D&Os relative to CODE products. The mean differences in Figure 13 vary between −0.6 and 0.7 ns, which indicates that the GPS DCBs estimated by LEO satellite and ground station data have good consistency. Figure 14 presents the RMS of the differences between the GPS DCB estimates adopting different modeling D&Os of 6, 8, and 10 and the CODE DCB products. To visually exhibit the differences in the RMS, Figure 15 displays the maximum differences in the RMS of the GPS DCB estimates using different modeling D&Os. The maximum differences reach 0.02 ns. The different modeling D&Os have certain effects on the RMS of the GPS DCB estimates. The corresponding mean RMS values are listed in Table 9. These RMS results are in the range of 0.8 ns. The mean RMS values are in the range of 0.3 ns, which indicates that the GPS DCBs estimated by LEO satellite and ground station data have good consistency. Statistically, Table 9 shows that the GPS DCBs with modeling D&Os of 8 and 10 have slightly higher RMS than the others. The RMS results using the GRCD observations are slightly higher than those obtained using the GRCC data. Considering the STD and RMS results of the GPS DCBs, the GPS DCB estimates with modeling D&Os of 8 attain superior precision.

3.2.3. LEO Receiver DCB Estimates and Stability Evaluation

Figure 16 exhibits the time series of the receiver DCB estimates with different modeling D&Os of 6, 8, and 10 for twin GRACE-FO satellites. The GRCC receiver DCB estimates are located at approximately −1.0 ns, whereas the GRCD DCB estimates fluctuate at approximately 4.55 ns. The GRCC and GRCD receiver DCBs with the 6 D&O in spherical harmonic modeling are greater than the others. In order to visually present the differences in receiver DCBs with different modeling D&Os, Figure 17 shows the maximum differences in LEO DCB estimates using different modeling D&Os. Meanwhile, Figure 18 displays the maximum differences in the STD of the receiver DCB using different modeling D&Os. The maximum differences are within 0.02 ns, inferior to the differences in the GPS DCBs, while the maximum differences in STDs are within 0.002 ns.
The mean values and STD results for the receiver DCB estimates with different modeling D&Os are listed in Table 10. Statistically, the mean receiver DCBs for GRCC and GRCD are both sorted in descending order as follows: receiver DCBs with 6, 8, and 10D&Os modeling, where the receiver DCBs decrease as the modeling D&Os increase. The GRCC receiver DCBs with 10D&O modeling gain the minimum STD result, whereas the GRCD DCBs applying 8 and 10D&Os have a better STD than the others. The receiver DCBs of the GRCC and GRCD differ; however, their DCB estimates have similar STD results.
In summary, adopting 8D&O in the LEO-based VTEC modeling exhibits superior estimates and precisions for both GPS and LEO DCBs. Additionally, compared to the results in Section 3.1, this demonstrates that the impact of modeling spacing on DCB estimates is greater than that of modeling D&O.

4. Conclusions

The effects of topside ionosphere VTEC modeling parameters on LEO-based DCB estimates and precision were investigated using GRACE-FO observation data to obtain superior modeling parameters. Different modeling parameters were set into the DCB estimations, encompassing the modeling spacing in the dynamic temporal mode and D&O in LEO-based VTEC modeling. The differences in GPS and LEO DCB estimates with the different modeling parameters were showcased, the precision of the DCB estimates was evaluated, and the effects of these parameters on DCB estimation were analyzed. The beneficial conclusions are drawn as follows:
(1) The mean GPS DCB estimates range between −10 and 10 ns. Although there is no significant difference in the mean GPS DCB estimates obtained using different modeling spacings, the mean DCB may have a limitation in reflecting all real DCBs. The maximum differences in the GPS DCB estimates after adopting different modeling spacings reach 0.05 ns. The different modeling spacings affect the GPS DCB estimates. The maximum differences in the RMS with different modeling spacings are in the range of 0.03 ns. The different modeling spacings have certain effects on the RMS of the GPS DCBs. The GPS DCB estimates with a modeling spacing of 12H have slightly higher precisions than the others.
The maximum difference in receiver DCBs adopting different modeling spacings is 0.22 ns, which indicates that modeling spacing has a significant impact on the receiver DCBs compared with GPS DCBs. The GRCC and GRCD receiver DCBs gain more optimal STD than the others when applying modeling spacings of 6H and 4H, respectively.
In summary, the GPS DCB estimates using a modeling spacing of 12H have higher precisions than the others, whereas LEO DCBs applying the modeling spacings of 4H or 6H obtain optimal STD results.
(2) The differences in the GPS DCB estimates using different modeling D&Os of 6, 8, and 10 are also in the range of 0.05 ns. The different modeling D&Os affect the GPS DCB estimates. The maximum difference in the RMS with different modeling spacings is 0.02 ns. The GPS DCB estimates, with a modeling D&O of 8, attain superior STD and RMS values.
In terms of receiver DCBs, the maximum differences in receiver DCBs using different modeling D&Os are within 0.02 ns, inferior to the differences in the GPS DCBs. The GRCC and GRCD DCBs achieve the optimal STD when adopting modeling D&Os of 10 and 8, respectively.
In summary, adopting 8D&O in the LEO-based VTEC modeling can obtain superior estimates and precisions for both GPS and LEO DCBs.
(3) The modeling spacing has a greater effect on the LEO receiver DCBs than on the GPS DCBs. The effects of the different modeling spacings and D&Os on the GPS DCB estimates and their RMS are at the same level. The modeling spacing has a greater impact on the LEO receiver DCBs than those of the modeling D&O. The research conclusions in this study can provide references for estimating the GNSS and LEO satellite DCBs using LEO onboard observation data.

Author Contributions

In this work, M.L. proposed the idea, designed the experiment results, and provided the software; Y.W. collected the experimental data; Y.W. and M.L. analyzed the experiment results and completed the paper. Y.Y., G.W. and H.G. supervised its analysis and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the China Postdoctoral Science Foundation (2021MD703896), supported by the Science and technology project of China Southern Power Grid Yunnan Power Grid Co., Ltd. (YNKJXM20220033), the China Natural Science Funds (No. 42074045, 42374043), and the Independent Project of the State Key Laboratory of China.

Data Availability Statement

CODE precise products are available at ftp://ftp.aiub.unibe.ch/ (accessed on 27 February 2023); onboard GPS observation data of the GRACE-FO satellites are available at ftp://rz-vm152.gfz-potsdam.de/ (accessed on 27 February 2023).

Acknowledgments

The authors are grateful to CODE for providing the precise products and onboard GPS observation data of GRACE-FO satellites, respectively. Meanwhile, we would like to thank the anonymous reviewers for their valuable comments.

Conflicts of Interest

The authors declare that this study received funding from China Southern Power Grid Yunnan Power Grid Co., Ltd. (YNKJXM20220033). The funder had the following involvement with the study: Yifan Wang and Mingming Liu analyzed the experiment results and completed the paper, Yunbin Yuan, Guofang Wang, and Hao Geng supervised its analysis and edited the manuscript.

References

  1. Schaer, S. Mapping and Predicting the Earth’s Ionosphere Using the Global Positioning System. Ph.D. Dissertation, University of Berne, Bern, Switzerland, 1999. [Google Scholar]
  2. Yuan, Y. Study on Theories and Methods of Correcting Ionospheric Delay and Monitoring Ionosphere Based on GPS. Ph.D. Dissertation, Institute of Geodesy and Geophysics, Chinese Academy of Sciences, Wuhan, China, 2002. [Google Scholar]
  3. Dach, R.; Schaer, S.; Hugentobler, U.; Schildknecht, T.; Gäde, A. Combined multi-system GNSS analysis for time and frequency transfer. In Proceedings of the 20th European Frequency and Time Forum, Braunschweig, Germany, 27–30 March 2006; pp. 530–537. [Google Scholar]
  4. Hernández-Pajares, M.; Juan, J.M.; Sanz, J.; Aragón-Àngel, À.; García-Rigo, A.; Salazar, D.; Escudero, M. The ionosphere: Effects, GPS modeling and the benefits for space geodetic techniques. J. Geod. 2011, 85, 887–907. [Google Scholar] [CrossRef]
  5. Montenbruck, O.; Hauschild, A. Code biases in multi-GNSS point positioning. In Proceedings of the 2013 International Technical Meeting of the Institute of Navigation, San Diego, CA, USA, 27–29 January 2013. [Google Scholar]
  6. Montenbruck, O.; Hauschild, A.; Steigenberger, P. Diferential code bias estimation using Multi-GNSS observations and global ionosphere maps. Navigation 2014, 61, 191–201. [Google Scholar] [CrossRef]
  7. Wang, N.; Yuan, Y.; Li, Z.; Montenbruck, O.; Tan, B. Determination of differential code biases with multi-GNSS observations. J. Geod. 2016, 90, 209–228. [Google Scholar] [CrossRef]
  8. Liu, T.; Zhang, B. Estimation of code observation-specific biases (OSBs) for the modernized multi-frequency and multi-GNSS signals: An undifferenced and uncombined approach. J. Geod. 2021, 95, 97. [Google Scholar] [CrossRef]
  9. Mannucci, A.; Wilson, B.; Yuan, D.; Ho, C.; Lindqwister, U.; Runge, T. A global mapping technique for GPS-derived ionospheric total electron content measurements. Radio Sci. 1998, 33, 565–582. [Google Scholar] [CrossRef]
  10. Schaer, S.; Gurtner, W.; Feltens, J. IONEX: The ionosphere map exchange format version 1. In Proceedings of the IGS AC Workshop, Darmstadt, Germany, 9–11 February 1998. [Google Scholar]
  11. Schaer, S. Differential code biases (DCB) in GNSS analysis. In Proceedings of the IGS Workshop, Miami Beach, FL, USA, 2–6 June 2008. [Google Scholar]
  12. Schaer, S. Activities of IGS Bias and Calibration Working Group. In IGS Technical Report 2011; Meindl, M., Dach, R., Jean, Y., Eds.; Astronomical Institute, University of Berne: Berne, Switzerland, 2012; pp. 139–154. [Google Scholar]
  13. Krankowski, A.; Shagimuratov, I.I.; Ephishov, I.I.; Krypiak-Gregorczyk, A.; Yakimova, G. The occurrence of themid-latitude ionospheric trough in GPS-TEC measurements. Adv. Space Res. 2009, 43, 1721–1731. [Google Scholar] [CrossRef]
  14. Yuan, Y.; Ou, J. A generalized trigonometric series function model for determining ionospheric delay. Prog. Nat. Sci. 2004, 14, 1010–1014. [Google Scholar] [CrossRef]
  15. Li, Z.; Yuan, Y.; Li, H.; Ou, J.; Huo, X. Two-step method for the determination of the differential code biases of COMPASS satellites. J. Geod. 2012, 86, 1059–1076. [Google Scholar] [CrossRef]
  16. Lin, J.; Yue, X.; Zhao, S. Estimation and analysis of GPS satellite DCB based on LEO observations. GPS Solut. 2014, 20, 251–258. [Google Scholar] [CrossRef]
  17. Zhong, J.; Lei, J.; Dou, X.; Yue, X. Is the long-term variation of the estimated GPS differential code biases associated with ionospheric variability? GPS Solut. 2015, 20, 313–319. [Google Scholar] [CrossRef]
  18. Jin, R.; Jin, S.; Feng, G. M_DCB: Matlab code for estimating GNSS satellite and receiver differential code biases. GPS Solut. 2012, 16, 541–548. [Google Scholar] [CrossRef]
  19. Wautelet, G.; Loyer, S.; Mercier, F.; Perosanz, F. Computation of GPS P1–P2 Differential Code Biases with JASON-2. GPS Solut. 2017, 21, 1619–1631. [Google Scholar] [CrossRef]
  20. Chen, P.; Yao, Y.; Li, Q.; Yao, W. Modeling the plasmasphere based on LEO satellites onboard GPS measurements. J. Geophys. Res. Space Phys. 2017, 122, 1221–1233. [Google Scholar] [CrossRef]
  21. Li, W.; Li, M.; Shi, C.; Fang, R.; Zhao, Q.; Meng, X.; Bai, W. GPS and BeiDou Differential Code Bias Estimation Using Fengyun-3C Satellite Onboard GNSS Observations. Remote Sens. 2017, 9, 1239. [Google Scholar] [CrossRef]
  22. Li, X.; Ma, T.; Xie, W.; Zhang, K.; Huang, J.; Ren, X. FY-3D and FY-3C onboard observations for differential code biases estimation. GPS Solut. 2019, 23, 57. [Google Scholar] [CrossRef]
  23. Zhou, P.; Nie, Z.; Xiang, Y.; Wang, J.; Du, L.; Gao, Y. Differential code bias estimation based on uncombined PPP with LEO onboard GPS observations. Adv. Space Res. 2019, 65, 541–551. [Google Scholar] [CrossRef]
  24. Ren, X.; Chen, J.; Zhang, X.; Schmidt, M.; Li, X.; Zhang, J. Mapping topside ionospheric vertical electron content from multiple LEO satellites at different orbital altitudes. J. Geod. 2020, 94, 86. [Google Scholar] [CrossRef]
  25. Liu, M.; Yuan, Y.; Huo, X.; Li, M.; Chai, Y. Simultaneous estimation of GPS P1–P2 differential code biases using low earth orbit satellites data from two different orbit heights. J. Geod. 2020, 94, 121. [Google Scholar] [CrossRef]
  26. Zakharenkova, I.; Cherniak, I. How can GOCE and TerraSAR-X contribute to the topside ionosphere and plasmasphere research? Space Weather 2015, 13, 271–285. [Google Scholar] [CrossRef]
  27. Zhong, J.; Lei, J.; Yue, X.; Dou, X. Determination of differential code bias of GNSS receiver onboard low earth orbit satellite. IEEE Trans. Geosci. Remote Sens. 2016, 54, 4896–4905. [Google Scholar] [CrossRef]
  28. Zhang, X.; Tang, L. Estimation of COSMIC LEO satellite GPS receiver differential code bias. Chin. J. Geophys. 2014, 57, 377–383. [Google Scholar]
  29. Wen, H.; Kruizingga, G.; Paik, M.; Landerer, F.; Bertiger, W.; Sakumura, C.; Bandikova, T.; Mccullough, C. Gravity Recovery and Climate Experiment Follow-On (GRACE-FO) Level-1 Data Product User Handbook; NASA: Silicon Valley, CA, USA, 2019. [Google Scholar]
  30. Foelsche, U.; Kirchengast, G. A simple “geometric” mapping function for the hydrostatic delay at radio frequencies and assessment of its performance. Geophys. Res. Lett. 2002, 29, 111-1–111-4. [Google Scholar] [CrossRef]
  31. Yue, X.; Schreiner, W.S.; Hunt, D.C.; Rocken, C.; Kuo, Y.H. Quantitative evaluation of the low Earth orbit satellite based slant total electron content determination. Space Weather 2011, 9, S09001. [Google Scholar] [CrossRef]
  32. Zhong, J.; Lei, J.; Dou, X.; Yue, X. Assessment of vertical TEC mapping functions for space-based GNSS observations. GPS Solut. 2015, 20, 353–362. [Google Scholar] [CrossRef]
  33. Klobuchar, J.A. A First-Order World-Wide Ionospheric Time-Delay Algorithm; Air Force Cambridge Research Labs.: Bedford, MA, USA, 1975; No. 324. [Google Scholar]
  34. IGS. IGS Technical Report 2011; Astronomical Institute, University of Berne: Berne, Switzerland, 2012. [Google Scholar]
  35. Sanz, J.; Miguel-Juan, J.; Rovira-Garcia, A.; Gonza’lez-Casado, G. GPS differential code biases determination: Methodology and analysis. GPS Solut. 2017, 21, 1549–1561. [Google Scholar] [CrossRef]
  36. 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]
  37. Kang, Z.; Tapley, B.; Bettadpur, S.; Ries, J.; Nagel, P.; Pastor, R. Precise orbit determination for the GRACE mission using only GPS data. J. Geod. 2006, 80, 322–331. [Google Scholar] [CrossRef]
  38. Montenbruck, O.; Andres, Y.; Bock, H.; van Helleputte, T.; van den Ijssel, J.; Loiselet, M.; Marquardt, C.; Silvestrin, P.; Visser, P.; Yoon, Y. Tracking and orbit determination performance of the GRAS instrument on MetOp-A. GPS Solut. 2008, 12, 289–299. [Google Scholar] [CrossRef]
  39. Hwang, C.; Tseng, T.-P.; Lin, T.; Švehla, D.; Schreiner, B. Precise orbit determination for the FORMOSAT-3/COSMIC satellite mission using GPS. J. Geod. 2008, 83, 477–489. [Google Scholar] [CrossRef]
  40. Cerri, L.; Berthias, J.P.; Bertiger, W.I.; Haines, B.J.; Lemoine, F.G.; Mercier, F.; Ries, J.C.; Willis, P.; Zelensky, N.P.; Ziebart, M. Precision Orbit Determination Standards for the Jason Series of Altimeter Missions. Mar. Geod. 2010, 33, 379–418. [Google Scholar] [CrossRef]
  41. Liu, M.; Yuan, Y.; Ou, J.; Chai, Y. Research on Attitude Models and Antenna Phase Center Correction for Jason-3 Satellite Orbit Determination. Sensors 2019, 19, 2408. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Maximum differences in GPS DCB estimates adopting different modeling spacings using twin LEO observations.
Figure 1. Maximum differences in GPS DCB estimates adopting different modeling spacings using twin LEO observations.
Remotesensing 15 05335 g001
Figure 2. Monthly stability (STD) of GPS DCB estimates adopting different modeling spacings ((top): GRCC and (below): GRCD data).
Figure 2. Monthly stability (STD) of GPS DCB estimates adopting different modeling spacings ((top): GRCC and (below): GRCD data).
Remotesensing 15 05335 g002
Figure 3. Maximum differences in STD of GPS DCB estimates adopting different modeling spacings.
Figure 3. Maximum differences in STD of GPS DCB estimates adopting different modeling spacings.
Remotesensing 15 05335 g003
Figure 4. Mean differences for GPS DCBs using different spacings relative to CODE products.
Figure 4. Mean differences for GPS DCBs using different spacings relative to CODE products.
Remotesensing 15 05335 g004aRemotesensing 15 05335 g004b
Figure 5. Root-mean square (RMS) for differences between the GPS DCB estimates using different modeling spacings and CODE DCB products ((top): GRCC and (below): GRCD data).
Figure 5. Root-mean square (RMS) for differences between the GPS DCB estimates using different modeling spacings and CODE DCB products ((top): GRCC and (below): GRCD data).
Remotesensing 15 05335 g005
Figure 6. Maximum differences in RMS of the GPS DCB estimates using different modeling spacings.
Figure 6. Maximum differences in RMS of the GPS DCB estimates using different modeling spacings.
Remotesensing 15 05335 g006
Figure 7. GRACE-FO receiver DCB estimates with different modeling spacings ((top): GRCC and (below): GRCD).
Figure 7. GRACE-FO receiver DCB estimates with different modeling spacings ((top): GRCC and (below): GRCD).
Remotesensing 15 05335 g007
Figure 8. Maximum differences in LEO DCB estimates using different modeling spacings.
Figure 8. Maximum differences in LEO DCB estimates using different modeling spacings.
Remotesensing 15 05335 g008
Figure 9. Maximum differences in STD of LEO DCB estimates using different spacings.
Figure 9. Maximum differences in STD of LEO DCB estimates using different spacings.
Remotesensing 15 05335 g009
Figure 10. Maximum differences in GPS DCB estimates adopting different D&Os using twin LEO observations.
Figure 10. Maximum differences in GPS DCB estimates adopting different D&Os using twin LEO observations.
Remotesensing 15 05335 g010
Figure 11. Monthly stability of GPS satellite DCB estimates using different modeling D&Os ((top): GRCC obs. and (below): GRCD obs.).
Figure 11. Monthly stability of GPS satellite DCB estimates using different modeling D&Os ((top): GRCC obs. and (below): GRCD obs.).
Remotesensing 15 05335 g011
Figure 12. Maximum differences in STD of GPS DCB estimates using different D&Os.
Figure 12. Maximum differences in STD of GPS DCB estimates using different D&Os.
Remotesensing 15 05335 g012
Figure 13. Mean differences for GPS DCBs using different D&Os relative to CODE products.
Figure 13. Mean differences for GPS DCBs using different D&Os relative to CODE products.
Remotesensing 15 05335 g013
Figure 14. RMS of differences between the GPS DCB estimates using different modeling D&Os and CODE DCB products ((top): GRCC data and (below): GRCD data).
Figure 14. RMS of differences between the GPS DCB estimates using different modeling D&Os and CODE DCB products ((top): GRCC data and (below): GRCD data).
Remotesensing 15 05335 g014
Figure 15. Maximum differences in RMS of the GPS DCB estimates using different D&Os.
Figure 15. Maximum differences in RMS of the GPS DCB estimates using different D&Os.
Remotesensing 15 05335 g015
Figure 16. GRACE-FO receiver DCB estimates adopting different modeling D&Os ((top): GRCC and (below): GRCD).
Figure 16. GRACE-FO receiver DCB estimates adopting different modeling D&Os ((top): GRCC and (below): GRCD).
Remotesensing 15 05335 g016
Figure 17. Maximum differences in LEO DCB estimates using different modeling D&Os.
Figure 17. Maximum differences in LEO DCB estimates using different modeling D&Os.
Remotesensing 15 05335 g017
Figure 18. Maximum differences in STD of LEO DCB estimates using different D&Os.
Figure 18. Maximum differences in STD of LEO DCB estimates using different D&Os.
Remotesensing 15 05335 g018
Table 1. Mean GPS differential code bias (DCB) estimates based on different modeling spacings using GRCC (a) and GRCD (b) observations [ns].
Table 1. Mean GPS differential code bias (DCB) estimates based on different modeling spacings using GRCC (a) and GRCD (b) observations [ns].
(a) GRCC
GPS2H4H6H12H
G01−7.4145−7.4174−7.4183−7.4191
G029.12439.12989.13309.1375
G03−5.0039−5.0161−5.0208−5.0262
G052.98262.97822.97752.9772
G06−6.6772−6.6762−6.6750−6.6726
G073.35693.35363.35323.3534
G08−7.5298−7.5226−7.5205−7.5182
G09−4.9578−4.9469−4.9431−4.9388
G10−5.0343−5.0449−5.0487−5.0529
G113.72133.73323.73743.7420
G124.16004.14804.14364.1385
G133.41753.42773.43073.4336
G141.91241.92051.92231.9231
G152.63212.64672.65232.6588
G162.79662.78442.78012.7753
G172.93212.94062.94292.9452
G18−0.2903−0.2915−0.2913−0.2908
G195.63295.64055.64255.6440
G201.47071.46401.46211.4597
G212.27432.27912.28212.2857
G227.27497.26487.26127.2572
G239.13299.14349.14699.1511
G24−5.5552−5.5542−5.5536−5.5526
G25−7.1077−7.1219−7.1273−7.1336
G26−8.5698−8.5834−8.5883−8.5941
G27−5.2432−5.2366−5.2348−5.2330
G283.20783.19533.19103.1861
G292.20612.21352.21542.2167
G30−6.0024−6.0065−6.0074−6.0079
G314.96354.95884.95734.9560
G32−4.3827−4.3735−4.3709−4.3688
(b) GRCD
GPS2H4H6H12H
G01−7.2909−7.2886−7.2910−7.2904
G029.30209.30949.31259.3178
G03−4.8126−4.8174−4.8270−4.8316
G053.02883.02633.02883.0305
G06−6.6031−6.6008−6.5987−6.5960
G073.41993.41943.41543.4150
G08−7.2286−7.2240−7.2212−7.2195
G09−4.7793−4.7650−4.7637−4.7593
G10−4.9828−4.9938−4.9933−4.9963
G113.78893.80413.80553.8099
G124.24064.22844.21844.2111
G133.47163.48043.48543.4884
G142.01122.01142.02172.0232
G152.64362.65802.66372.6704
G162.94512.93382.92752.9224
G172.98432.99412.99562.9973
G18−0.2142−0.2115−0.2117−0.2098
G195.73025.73935.74115.7430
G201.59281.58631.58731.5864
G212.34552.34982.35542.3602
G227.50317.49497.48977.4856
G239.28149.29489.29589.3000
G24−5.4638−5.4579−5.4583−5.4559
G25−7.1005−7.1194−7.1222−7.1292
G26−8.5873−8.5992−8.6059−8.6114
G27−5.1256−5.1196−5.1166−5.1145
G283.26583.25263.24593.2399
G292.25012.25472.26102.2630
G30−5.9075−5.9087−5.9144−5.9158
G315.03445.02645.02645.0244
G32−4.4046−4.4053−4.3915−4.3886
Table 2. Maximum, minimum, and differences in G15 and G12 DCB estimates based on different spacings using GRCC (a) and GRCD (b) observations [ns].
Table 2. Maximum, minimum, and differences in G15 and G12 DCB estimates based on different spacings using GRCC (a) and GRCD (b) observations [ns].
(a) GRCC-G15
DOYMaxMinDiff.
12.65762.63210.0255
22.65442.63000.0244
32.71292.67960.0333
42.62712.59600.0311
52.68152.64150.0400
62.63312.59470.0384
72.74362.72340.0202
82.61722.58130.0359
92.65172.62260.0291
102.64602.61300.0330
112.59142.56750.0239
122.73382.71490.0189
132.65732.63110.0262
142.65952.63030.0292
152.62392.59470.0292
162.77302.74690.0261
172.55072.52570.0250
182.73132.71190.0194
192.61082.58380.0270
202.54932.52690.0224
212.60502.57890.0261
222.58172.56080.0209
232.56382.54290.0209
242.67312.65080.0223
252.64412.61130.0328
262.70222.67260.0296
272.73572.71160.0241
282.74912.72880.0203
292.67882.65580.0230
302.72422.70110.0231
(b) GRCD-G12
DOYMaxMinDiff.
14.19294.17040.0225
24.31654.28220.0343
34.18624.14280.0434
44.33404.28740.0466
54.22774.19500.0327
64.20234.17310.0292
74.09694.07640.0205
84.18864.14540.0432
94.23704.19400.0430
104.26234.22940.0329
114.27044.23740.0330
124.25134.22830.0230
134.30404.27590.0281
144.14814.11610.0320
154.29784.26380.0340
164.31964.29970.0199
174.24994.23160.0183
184.26844.24360.0248
194.25864.22210.0365
204.22464.18790.0367
214.31874.29140.0273
224.23454.21270.0218
234.25934.23330.0260
244.28214.25260.0295
254.21354.18110.0324
264.20524.17900.0262
274.21334.19310.0202
284.20854.18960.0189
294.17184.15060.0212
304.27384.24840.0254
Table 3. Mean STD of GPS DCBs adopting different modeling spacings [ns].
Table 3. Mean STD of GPS DCBs adopting different modeling spacings [ns].
LEOVar.STDLEOVar.STD
GRCC2H0.0624GRCD2H0.0648
4H0.06214H0.0643
6H0.06206H0.0641
12H0.061912H0.0639
Table 4. Mean RMS statistics for GPS DCBs using different modeling spacings [ns].
Table 4. Mean RMS statistics for GPS DCBs using different modeling spacings [ns].
LEOVar.RMSLEOVar.RMS
GRCC2H0.2801GRCD2H0.2763
4H0.27694H0.2732
6H0.27596H0.2722
12H0.274812H0.2711
Table 5. Mean values and STD statistics for receiver DCBs with different modeling spacings [ns].
Table 5. Mean values and STD statistics for receiver DCBs with different modeling spacings [ns].
LEOVar.MeanSTD
GRCC2H−1.06840.0465
4H−1.02020.0459
6H−0.98330.0456
12H−0.91150.0458
GRCD2H4.48470.0439
4H4.53650.0434
6H4.57610.0450
12H4.65310.0468
Table 6. Mean GPS DCB estimates based on different modeling D&Os using GRCC (a) and GRCD (b) data [ns].
Table 6. Mean GPS DCB estimates based on different modeling D&Os using GRCC (a) and GRCD (b) data [ns].
(a) GRCC
GPS6D&O8D&O10D&O
G01−7.4142−7.4174−7.4178
G029.12619.12989.1342
G03−5.0144−5.0161−5.0181
G052.96862.97822.9862
G06−6.6785−6.6762−6.6719
G073.35603.35363.3555
G08−7.5223−7.5226−7.5236
G09−4.9485−4.9469−4.9473
G10−5.0425−5.0449−5.0479
G113.73383.73323.7316
G124.14864.14804.1491
G133.42543.42773.4261
G141.92231.92051.9152
G152.64272.64672.6493
G162.78122.78442.7874
G172.94112.94062.9395
G18−0.2882−0.2915−0.2940
G195.64205.64055.6387
G201.46151.46401.4668
G212.27612.27912.2807
G227.26707.26487.2637
G239.14259.14349.1420
G24−5.5488−5.5542−5.5550
G25−7.1203−7.1219−7.1223
G26−8.5842−8.5834−8.5830
G27−5.2361−5.2366−5.2396
G283.19813.19533.1957
G292.21582.21352.2079
G30−6.0024−6.0065−6.0056
G314.96134.95884.9594
G32−4.3729−4.3735−4.3770
(b) GRCD
GPS6D&O8D&O10D&O
G01−7.2891−7.2886−7.2899
G029.30469.30949.3138
G03−4.8226−4.8174−4.8238
G053.01693.02633.0372
G06−6.6027−6.6008−6.5959
G073.41773.41943.4179
G08−7.2217−7.2240−7.2245
G09−4.7683−4.7650−4.7695
G10−4.9894−4.9938−4.9927
G113.80183.80413.7994
G124.22734.22844.2246
G133.48023.48043.4802
G142.02132.01142.0155
G152.65492.65802.6601
G162.92952.93382.9345
G172.99452.99412.9910
G18−0.2099−0.2115−0.2140
G195.74065.73935.7359
G201.58541.58631.5916
G212.34902.34982.3540
G227.49437.49497.4922
G239.29229.29489.2901
G24−5.4554−5.4579−5.4609
G25−7.1148−7.1194−7.1166
G26−8.6021−8.5992−8.6001
G27−5.1179−5.1196−5.1212
G283.25433.25263.2514
G292.26052.25472.2539
G30−5.9095−5.9087−5.9124
G315.03115.02645.0287
G32−4.3943−4.4053−4.3967
Table 7. Maximum, minimum, and differences in G05 DCB estimates with different D&Os using GRCC (a) and GRCD (b) data [ns].
Table 7. Maximum, minimum, and differences in G05 DCB estimates with different D&Os using GRCC (a) and GRCD (b) data [ns].
(a) GRCC-G05
DayMaxMinDiff.
12.97532.96940.0059
22.93702.90830.0287
33.05353.02790.0256
43.07843.05840.0200
53.04443.03160.0128
63.02643.02020.0062
72.97242.95340.0190
82.97802.95070.0273
92.99062.96380.0268
102.94402.94320.0008
112.94632.94550.0008
123.04053.02690.0136
133.00992.97760.0323
143.03112.99850.0326
152.94262.92010.0225
163.01093.00310.0078
172.93052.92310.0074
182.96712.94540.0217
192.96062.93280.0278
202.98212.95930.0228
212.89372.88200.0117
223.01973.01480.0049
232.97962.96830.0113
243.01982.99330.0265
252.97092.94030.0306
263.00222.98460.0176
272.93352.92780.0057
282.97002.96400.0060
292.91862.89800.0206
303.05513.02520.0299
(b) GRCD-G05
DayMaxMinDiff.
12.98492.97190.0130
22.98262.95220.0304
33.07703.04000.0370
43.05513.02590.0292
53.05233.03510.0172
63.01773.00860.0091
73.06113.03900.0221
83.02893.00090.0280
93.03172.99640.0353
102.96612.94880.0173
113.00903.00800.0010
123.12873.11310.0156
133.09023.05690.0333
143.05153.00480.0467
152.98742.96430.0231
163.04143.03520.0062
172.97142.96720.0042
183.03203.01240.0196
193.05033.01950.0308
203.02533.00440.0209
212.95722.94510.0121
223.09513.09010.0050
233.06273.05070.0120
243.10393.07700.0269
253.03413.00540.0287
263.06613.04660.0195
272.98572.97890.0068
283.03393.02630.0076
293.00922.98640.0228
303.12413.09640.0277
Table 8. Mean STD of GPS DCBs using different modeling D&Os [ns].
Table 8. Mean STD of GPS DCBs using different modeling D&Os [ns].
LEOVar.STDLEOVar.STD
GRCC6D&O0.0620GRCD6D&O0.0643
8D&O0.06218D&O0.0643
10D&O0.062310D&O0.0644
Table 9. Mean RMS statistics for GPS DCBs using different modeling D&Os [ns].
Table 9. Mean RMS statistics for GPS DCBs using different modeling D&Os [ns].
LEOVar.RMSLEOVar.RMS
GRCC6D&O0.2780GRCD6D&O0.2742
8D&O0.27698D&O0.2732
10D&O0.276810D&O0.2732
Table 10. Mean values and STD statistics for receiver DCBs with different modeling D&Os [ns].
Table 10. Mean values and STD statistics for receiver DCBs with different modeling D&Os [ns].
LEOVar.MeanSTD
GRCC6D&O−1.01320.0460
8D&O−1.02020.0459
10D&O−1.02150.0457
GRCD6D&O4.54390.0449
8D&O4.53650.0434
10D&O4.53530.0441
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Y.; Liu, M.; Yuan, Y.; Wang, G.; Geng, H. Effects of Topside Ionosphere Modeling Parameters on Differential Code Bias (DCB) Estimation Using LEO Satellite Observations. Remote Sens. 2023, 15, 5335. https://doi.org/10.3390/rs15225335

AMA Style

Wang Y, Liu M, Yuan Y, Wang G, Geng H. Effects of Topside Ionosphere Modeling Parameters on Differential Code Bias (DCB) Estimation Using LEO Satellite Observations. Remote Sensing. 2023; 15(22):5335. https://doi.org/10.3390/rs15225335

Chicago/Turabian Style

Wang, Yifan, Mingming Liu, Yunbin Yuan, Guofang Wang, and Hao Geng. 2023. "Effects of Topside Ionosphere Modeling Parameters on Differential Code Bias (DCB) Estimation Using LEO Satellite Observations" Remote Sensing 15, no. 22: 5335. https://doi.org/10.3390/rs15225335

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