Differential code bias (DCB) is one of the main errors involved in ionospheric total electron content (TEC) retrieval using a global navigation satellite system (GNSS). It is typically assumed to be constant over time. However, this assumption is not always valid because receiver DCBs have long been known to exhibit apparent intraday variations. In this paper, a combined method is introduced to estimate the epoch-wise receiver DCB, which is divided into two parts: the receiver DCB at the initial epoch and its change with respect to the initial value. In the study, this method was proved feasible by subsequent experiments and was applied to analyze the possible reason for the intraday receiver DCB characteristics of 200 International GNSS Service (IGS) stations in 2014 (high solar activity) and 2017 (low solar activity). The results show that the proportion of intraday receiver DCB stability less than 1 ns increased from 72.5% in 2014 to 87% in 2017, mainly owing to the replacement of the receiver hardware in stations. Meanwhile, the intraday receiver DCB estimates in summer generally exhibited more instability than those in other seasons. Although more than 90% of the stations maintained an intraday receiver DCB stability within 2 ns, substantial variations with a peak-to-peak range of 5.78 ns were detected for certain stations, yielding an impact of almost 13.84 TECU on the TEC estimates. Moreover, the intraday variability of the receiver DCBs is related to the receiver environment temperature. Their correlation coefficient (greater than 0.5 in our analyzed case) increases with the temperature. By contrast, the receiver firmware version does not exert a great impact on the intraday variation characteristics of the receiver DCB in this case.
Keywords: differential code bias variation; modified carrier-to-code leveling; ionospheric total electron content; solar activity
Due to high temporal continuity and spatial coverage, global navigation satellite systems (GNSSs) have become a key technology, leading to significant advances in ionosphere sounding [[
The existing studies that have examined the time variability of receiver DCBs using real GNSS data can be classified into two groups. The first group considers receiver DCB estimates as byproducts of ionospheric modeling. Their methods consist of either estimating the DCBs and the parameters of the ionospheric model in a common adjustment process or determining the DCB estimates using pre-computed global ionosphere maps (GIMs) [[
To eliminate the need for either ionospheric modeling or an external reference such as a short-baseline setup or GNSS simulator, we introduced an effective method called modified carrier-to-code leveling (MCCL) for detecting receiver DCB variations (relative values with respect to the first epoch) in an epoch-wise manner using a standalone receiver [[
This paper proceeds as follows. First, we review the derivation procedure for MCCL and describe how to address the rank deficiencies in detail. Then, we introduce the approach to retrieving the absolute receiver DCBs based on MCCL and validate this approach by means of simulated DCB variations and an experiment using co-located receivers. Following this validation, we use the proposed approach to assess the intraday stability of the receiver DCBs for 200 IGS stations in two years (2014 and 2017) with different solar activities. Hence, the characteristics of the intraday variations in these receiver DCB estimates are investigated in detail. Furthermore, the significant intraday variations in the receiver DCBs and their adverse effect on STEC estimates are presented. The impacts of the receiver firmware version and the surrounding temperature on the intraday variability of the receiver DCBs are also analyzed. In the last section, we summarize the results and conclusions.
This section briefly reviews the basic principles and technological aspects of modified carrier-to-code leveling (MCCL), which are used to estimate receiver DCB variations with respect to its value at the initial epoch, and a thin-layer ionospheric model used for isolating the initial value of the receiver DCB from ionosphere observations. We primarily focus on the development of a functional model and the process of epoch-wise receiver DCB estimates.
The starting point for CCL is the geometry-free (GF) combination of carrier phase
(
where
(
with
Unfortunately, Equation (
(
The second type of rank deficiency appears among
(
Considering Equations (
(
It should be noted that
As described above, to obtain the absolute receiver DCB variation, the reference datum
(
where
(
where
By inserting Equation (
(
where
To test the performance of the receiver DCB estimation scheme presented in Section 2, the following performance metrics are employed in this section: (
The experimental data were collected by two sets of collocated GPS receivers with a sampling rate of 30 s and a cut-off elevation angle of 5 degrees. The first set of data were also used by the authors of [[
Table 1 presents an overview of the basic characteristics of the receiver pairs, including the receiver and antenna types, the approximate receiver locations and the periods of time covered by the observations. DLF4 and DLF5, which are two identical receivers, are separated by a zero baseline because they are connected to the same antenna via an antenna-splitter device located 10 m from the antenna of the DLFT. By contrast, ALGO, ALG2 and ALG3 are independently connected to three different antennas, creating three short baselines of distances between 70 and 150 m. During data processing, we used a cut-off elevation angle of 20 degrees to discard particularly noisy GPS data and adopted an elevation-dependent weighting of the observations. The elevation angles were calculated on the basis of the satellite positions computed using broadcast ephemeris and the known receiver positions. The zenith-referenced standard deviations were empirically set to 30 cm for the code and 0.3 cm for the phase. First, simulations were carried out to illustrate the efficiency of the proposed approach. The variations in the receiver DCBs simulated using the given mathematical functions were incorporated into the GPS C1/P1 data, whereas the remaining P2, L1 and L2 data were left as the original data. Then, the changed GPS observables, including the receiver DCBs known a priori, could be obtained. Second, based on [[
To verify the effectiveness of the proposed method, we modified the code observation by adding the DCB variation simulated from three types of mathematical representations. This was also used as the reference value for method valuation. In Figure 1, the red lines indicate the simulated values of the receiver DCB obtained from the three given mathematical representations, namely, a constant function (Figure 1a), a cosine function (Figure 1b) and a sine function (Figure 1c), which were also incorporated into the raw code observables collected by ALG3 on day 07 of 2017. The scattered blue dots represent the corresponding estimates retrieved from the modified code observables using the proposed method. It can be seen that an intraday variation of as small as 2 ns in the receiver DCB can successfully be detected, and good agreement between the estimated and simulated values is evident. This agreement demonstrates the feasibility of the proposed approach for retrieving absolute receiver DCB variations using a standalone receiver. The standard deviation of these estimates is at the level of approximately 0.5 ns, which is similar to the noise level of the measurements.
In Figure 2, the BR-DCBs are plotted for three baselines, namely, DLFT-DLF4, DLFT-DLF5 and DLF4-DLF5, on day 172 in 2010. The red lines depict the BR-DCB estimates obtained using the SD-GF method [[
Figure 3 shows the variations in the single-receiver DCBs estimated using the MCCL plus datum method for the same three receivers (DLFT, DLF4 and DLF5), where the epoch-wise variations in the DCBs were obtained via MCCL and the datum was isolated from the ionospheric observables based on the ionospheric model. Among the considered receivers, none of their DCBs remained stable throughout the analyzed time period. In particular, the receiver DCB of DLFT showed a substantial variation, with a peak-to-peak range of almost 8 ns. This finding also confirms the previous BR-DCB estimates, which indicate that a significant variation in the form of an inverted "W-shape" can mainly be attributed to the instability of the DLFT receiver. At the same time, both DLF4 and DLF5 present less pronounced variations in their receiver DCBs, following a similar trend with a peak-to-peak value of almost 3 ns, since they are of the same receiver type and are connected to the same antenna.
Figure 4 and Figure 5 shows results that are analogous to those in Figure 2 and Figure 3 but based on the second set of GPS data listed in Table 1. It can be seen from Figure 5 that the DCB estimates for station ALGO clearly vary during the day, with a peak-to-peak range of approximately 4 ns, while those from stations ALG2 and ALG3 are rather constant. These characteristics of the receiver DCB variations, as retrieved with the MCCL plus datum method, show strong agreement with those observed in the time series of BR-DCBs shown in Figure 4. In conclusion, compared with the BR-DCBs, the analysis of single-receiver DCB estimates allows a receiver with significant intraday DCB variations to be directly distinguished. On the other hand, these results further prove the validity of the proposed method for retrieving absolute receiver DCB estimates with an epoch-wise temporal resolution.
Many previous studies have analyzed the annual stability of, or day-to-day variations in, receiver DCBs. However, the intraday variations in single-receiver DCBs have not yet been systematically studied at an epoch-wise temporal resolution. Therefore, our work mainly focused on the intraday stability of GPS receiver DCB estimates based on the MCCL method for 200 IGS stations, including 6 main types of receivers. The distribution of the stations is shown in Figure 6. The observation period covered 2 years, 2014 (high solar activity) and 2017 (low solar activity). The data sampling rate was 30 s. The intraday stability of epoch-wise receiver DCB estimates can be expressed in terms of their daily standard deviation (DSTD), which is computed using the following equation:
(
where
Figure 7 demonstrates the maximum intraday stability index DSTD of receiver DCB estimates derived using the MCCL method for 200 IGS stations in 2014 (red dots) and 2017 (blue dots), which are aligned according to their geomagnetic latitudes ranging from 90°S to 90°N. In the upper scatter plot, the red dots are more scattered than the blue ones, implying that, in general, the receiver DCBs in 2017 are more stable than those in 2014. As the bottom histogram shows, the DSTD is smaller than 1 ns (or 2ns) in 72% (or 93%) of cases for all the receivers in 2014, while these variations in 2017 are less than 1 ns (or 2ns) in 87% (or 96%) of cases. However, it can be seen from the scatter plot that the maximum DSTD values in the low-latitude region (from 30°S to 30°N, with an average DSTD of 0.74 ns) are no more scattered than those in the middle- and high-latitude regions (from 30°S (N) to 80°S (N), with an average DSTD of 0.76 ns) in either 2014 or 2017. These results suggest that the intraday stability of the receiver DCBs is not associated with the geomagnetic latitudes of the receiver locations.
Figure 8 shows the DSTD distribution of the receiver DCBs in different seasons. In the northern hemisphere, the three months of June, July and August are regarded as summer and the three months of December, January and February are winter. The opposite trend is observed in the southern hemisphere. The other months with relatively moderate temperatures are classified as spring and autumn. It can be seen that the DSTD of the receiver DCB estimates within 1 ns (or 2 ns) in winter accounts for 89.5% (or 99.5%) of cases for all of the receivers in 2014, while for those in summer, it accounts for 72.5% (or 99.5%) of all the receivers. In spring and autumn, the ratio of the DSTD within 1 ns (or 2 ns) reaches 84.0% (95.5%). These results indicate that the intraday stability of receiver DCBs may be related to the seasons, especially meteorological parameters such as the temperature of the receiver environment. Moreover, compared with those in 2014, the number of receivers with a DSTD less than 1 ns (or 2 ns) in different seasons of 2017 is correspondingly increased. Even though the solar activity in 2017 is lower than that in 2014, the increase may not be attributed to this factor. As Figure 7 shows, there is no obvious discrepancy in the intraday stabilities of the receiver DCBs between the stations located at different geomagnetic latitudes.
To further investigate the difference in the intraday stabilities of the receiver DCBs between 2014 and 2017, we conducted a comparison of the intraday variation in the receiver DCBs in the two years for some specific stations where receiver DCBs estimates exhibited more stability in 2017 than in 2014. The information regarding the receiver types and antenna types is presented in Table 2. For these stations, the intraday variations in the receiver DCBs over one day in two different years are shown in Figure 9. It can be seen that the stations equipped with the same types of receivers and antennas have the same intraday variations in their receiver DCBs, as Figure 9a–c shows, while for the stations equipped with different receivers and antennas, the intraday variation in the receiver DCBs is more stable, with different patterns (Figure 9d–i) in 2017 compared to 2014. These results demonstrate that the improvement in the intraday stability of receiver DCBs is mainly attributed to the replacement of the receivers and antennas of IGS stations.
To represent this information in greater detail, the time series of the epoch-wise receiver DCBs on representative days are shown in Figure 10, where the yellow lines represent polynomial fits to the estimated receiver DCBs and the red lines represent the corresponding receiver DCB products provided by CODE. Since the MCCL-based estimates have an epoch-wise temporal resolution, the receiver DCBs change gradually rather than abruptly, as is the case for the CODE-based estimates, between these two consecutive days; however, similar variation trends between these two DCB products for DOY 016–017 are evident. We can take this similarity as an indication that the obvious jumps in the CODE-based DCBs between the different days are related, in part, to the intraday variations in the receiver DCBs. A closer look at Figure 10 reveals three types of intraday variation characteristics. The first type is a quasi-periodic variation, as shown in Figure 10a–d. These DCB series show recovery to the initial state over a one-day period, most likely due to a thermally induced variation, such as a variation in the ambient temperature. Therefore, these DCB series may exhibit daily repeatability. The second type is an almost quasi-linear variation, as shown in Figure 10e,f, but with a magnitude smaller than 1 ns. Due to the lack of knowledge about the receiver design, the exact reason for this variation is not clear and requires further study. The third type of variation in the DCB series, as shown in Figure 10g,h, appears to be similar to random noise. We can conclude that these DCBs remain stable at the intraday level relative to their noise level.
It has been shown that variations in receiver DCBs may be associated with environmental temperature variations and the receiver firmware version [[
As seen in Figure 11, the receiver DCB variation is, in principle, consistent with the temperature variation. Meanwhile, all the absolute values of the PCC on different days are greater than 0.5, indicating a strong correlation between the receiver DCB and the temperature. Moreover, the slope of the corresponding linear regression equation is presented to visualize the extent of the impact of temperature variation in the receiver DCB variation, as shown in Figure 12.
Table 3 presents the statistical results, from which it can be seen that the absolute value of the PCC and the slope of the linear regression equation become larger as the daily temperature increases, and vice versa. Thus, it is confirmed that temperature is a major contributor to the intraday variation in the receiver DCBs. In addition, during the period 2017–2018, the receiver firmware version at station ALIC was changed three times; however, there are no obvious corresponding differences in the intraday variations in the receiver DCB estimates. The PCC gaps between consecutive days corresponding to firmware version changes (from 30 January 2017 to 31 January 2017 and from 21 May 2018 to 22 May 2018) are smaller than 0.07. These observations indicate that the receiver firmware version has no great impact on receiver DCB variation in this case 3.5. Effects of intraday variations in receiver DCBs on STEC observables.
Although most receiver DCBs remain relatively stable at the daily timescale, large intraday variations in receiver DCBs have been detected at certain stations. To provide a better understanding of the impact of abnormal intraday variations in receiver DCBs on the extraction of ionospheric observables, let us consider the example of DOY 016, for which the receiver DCB at station ADIS varies significantly, whereas the receiver DCB at station ZECK remains stable. For this example, in Figure 13, we compare the intraday variations in the receiver DCBs with the STEC biases caused by these variations. The STEC biases caused by the intraday variations in the receiver DCBs can be interpreted as differences between the CCL-derived and MCCL-derived ionospheric observables, because the MCCL method compensates for possible between-epoch variations in the receiver DCBs, which, consequently, are not reflected in the ionospheric STEC observables derived using this method. When the receiver DCBs exhibit a significant variation, the STEC biases become correspondingly apparent. As the panels on the left of Figure 13 show, when the receiver DCBs at ADIS vary over time with a peak-to-peak range of almost −5.78 ns, the STEC biases reach −13.84 TECU (peak-to-peak value), while the right panels show that the receiver DCBs at ZECK remain stable throughout the time period and the corresponding STEC biases are smaller than 1 TECU. Thus, under the assumption that the receiver DCBs are constant over a one-day period, intraday variations in the receiver DCBs will introduce errors to the process of ionospheric TEC extraction. Therefore, considering the receiver DCB as a time-varying parameter is helpful for enhancing the reliability of TEC estimation.
This study presented the epoch-wise estimation of single-receiver DCBs consisting of DCB variations estimated based on the MCCL, together with the absolute DCB datum isolated from MCCL-derived ionospheric measurements. The feasibility and effectiveness of this scheme for estimating single-receiver DCBs with an epoch-wise temporal resolution were validated through numerical simulation and zero-/short-baseline experiments. Compared with typical daily receiver DCB estimates, the estimates obtained using the proposed method have two advantages: (
Based on epoch-wise receiver DCB estimates, the intraday stability of a GPS receiver DCBs for 200 IGS stations, most of which contribute to the GIM product, was assessed for the two years of 2014 (high solar activity) and 2017 (low solar activity). The maximum DSTD values of the receiver DCBs in 2014 and 2017 were less than 1 ns in 72.5 and 87% of all cases, respectively. The difference in the intraday stability of the receiver DCBs between the two years is mainly attributed to the replacement of receivers and antennas at IGS stations. In addition, the intraday stability of the receiver DCBs was generally lower in summer than in the other seasons in the two years. Although more than 90% of the investigated receiver DCBs showed an intraday stability of less than 2 ns, significant variations over a one-day period were observed with a maximum range of approximately 5.78 ns (peak-to-peak value), corresponding to an impact of almost 13.85 TECU on the extraction of ionospheric STEC values in this case.
Furthermore, the epoch-wise receiver DCB estimates exhibit three representative patterns of intraday variation: quasi-periodic variation, quasi-linear variation and random noise fluctuation. In contrast to the CODE-based DCB product, the epoch-wise receiver DCB estimates do not show day-boundary jumps; instead, they vary gradually over time. However, they also present trends of day-to-day variation similar to those of the CODE-based DCB product, indicating that the abrupt change in the CODE-based DCB between two consecutive days is partially due to the intraday variation in the receiver DCB.
We also found that the significant variability in the estimates of the receiver DCBs over a one-day period showed good consistency with the intraday temperature variation. The maximum correlation coefficient between the receiver DCB and the temperature was as high as 0.9, and the minimum value was still greater than 0.5. By contrast, when the receiver firmware version was updated, the intraday variation characteristics of the receiver DCBs showed no significant change. It can be concluded that temperature is one of the main factors affecting the intraday variability of receiver DCBs.
Graph: Figure 1 Simulated values (red line) versus estimates (blue line) of receiver DCB variations at site ALG3 on DOY 07, 2017: (a) DCB(i)=0 , (b) DCB(i)=Cos (2πn⋅i)−1 and (c) DCB(i)=2⋅sin2πn⋅i for epochs i=1,2...n , where n=2880 and the interval is 30 s.
Graph: Figure 2 Epoch-by-epoch estimates of BR-DCB using SD-GF method (red line) and MCCL plus datum method (blue line) for three pairs of co-located receivers: DLFT-DLF4, DLFT-DLF5, DLF4-DLF5 on DOY 172, 2010.
Graph: Figure 3 Epoch-by-epoch estimates of single-receiver DCBs retrieved using the MCCL plus datum method for three co-located receivers: DLFT, DLF4, DLF5 on DOY 172, 2010.
Graph: Figure 4 Epoch-by-epoch estimates of BR-DCB using SD-GF method (red line) and MCCL plus datum method (blue line) for three pairs of co-located receivers: ALG2-ALG3, ALGO-ALG3, ALGO-ALG2 on DOY 07, 2017.
Graph: Figure 5 Epoch-by-epoch estimates of single-receiver DCBs retrieved with the MCCL plus datum method for three co-located receivers: ALGO, ALG2, ALG3 on DOY 07, 2017.
Graph: Figure 6 Geographic locations of receivers (of different types) that provide the second set of data analyzed in this work.
Graph: Figure 7 The upper scatter plot displays the maximum DSTD of the estimated receiver DCBs in 2014 and 2017, where the stations are aligned in accordance with their geomagnetic latitudes along the horizontal axis. The bottom histograms show the distribution of the DSTD values, where the vertical axis represents the number of receivers.
Graph: Figure 8 The distribution of the intraday stability DSTD of the estimated receiver DCBs determined using MCCL in different seasons of 2014 (red) and 2017 (blue). The vertical axis stands for the number of selected receivers and the horizontal axis stands for the DSTD of receiver DCB estimates.
Graph: Figure 9 Intraday variations in the receiver DCBs estimates obtained using MCCL for the analyzed stations over one day (DOY 06) in 2014 (red lines) and 2017 (blue lines). The lines are shifted to the same datum for ease of comparison.
Graph: Figure 10 Intraday variation characteristics of receiver DCBs spanning consecutive days. The blue scattered dots represent the epoch-wise estimates of the receiver DCBs, the yellow lines represent fits to these values, and the red solid lines represent the CODE DCB products.
Graph: Figure 11 Receiver DCB variations (blue lines) at station ALIC, as estimated with the MCCL method for various receiver firmware versions and various temperature conditions. The intraday temperature values (red lines) were extracted from IGS meteorological data. PCC denotes the Pearson correlation coefficient between the blue line and the red line.
Graph: Figure 12 Receiver DCB variations (vertical axis) plotted as a function of the intraday temperature values (horizontal axis). The red line represents the results obtained by fitting a linear regression equation to the same dataset used in Figure 11.
Graph: Figure 13 Biases of the ionospheric STEC (upper panels, where different colors correspond to different satellite arcs) stemming from the intraday variations in the receiver DCBs (bottom panels) for stations ADIS and ZECK on DOY 16, 2017.
Table 1 A general overview of two groups of co-located receivers.
Receiver Name Receiver Type Antenna Type Location Observation Period DLFT JPS LEGACY JPSREGANT_DD_E 51.98°E, 4.39°N 2010, day 172 DLF4 SEPT POLARX2 LEIAT504 DLF5 SEPT POLARX2 LEIAT504 ALGO AOA BENCHMARK ACT AOAD/M_T 78.07°W, 45.95°N 2017, day 07 ALG2 TPS NET-G3A NOV750.R4 ALG3 TPS NETG3 TPSCR.G3
Table 2 Information on receiver and antenna types used by the analyzed stations.
Station 2014 2017 Receiver Type Antenna Type Receiver Type Antenna Type CHIL TPS NET-G3A TPSCR.G3 TPS NET-G3A TPSCR.G3 COYQ TRIMBLE NETRS ASH700936D_M TRIMBLE NETRS ASH700936D_M ZECK JAVAD TRE_G3TH JAVRINGANT_DM JAVAD TRE_G3TH JAVRINGANT_DM POVE TRIMBLE NETR5 TRM29659.00 TRIMBLE NETRS ASH701945B_M COCO TRIMBLE NETR8 AOAD/M_T SEPT POLARXS AOAD/M_T SFER TRIMBLE NETRS TRM29659.00 LEICA GR25 LEIAR20 PRDS TPS NET-G3A AOAD/M_T JAVAD TRE_G3TH AOAD/M_T PALV JPS EGGDT ASH700936D_M JAVAD TRE_G3TH ASH700936D_M YAR2 ASHTECH UZ-12 AOAD/M_T SEPT POLARX4TR AOAD/M_T
Table 3 Factors influencing the intraday variability of the receiver DCBs at station ALIC, with a LEICA GR25 receiver and a LEIAR25.R3 antenna.
Time Firmware Version Temperature Span (°C) PCC Slope 30 January 2017 4.02.386 22.9~44.2 −0.914 −0.11 31 January 2017 4.10.598 21.6~42.6 −0.927 −0.16 17 June 2017 4.11.606 2.5~24.7 −0.791 −0.092 29 June 2017 4.11.606 1.7~17.5 −0.693 −0.078 21 May 2018 4.11.606 2.3~22.5 −0.52 −0.036 22 May 2018 4.20.232 7.3~22.92 −0.65 −0.051
Conceptualization, X.Z. and L.X.; methodology, X.Z.; software, X.Z.; validation, X.Z., H.L. and Q.L.; writing—original draft preparation, X.Z.; writing—review and editing, X.Z. and L.X.; visualization, X.Z.; supervision, L.X.; project administration, L.X. All authors have read and agreed to the published version of the manuscript.
Not applicable.
The authors declare no conflict of interest.
By Xiao Zhang; Linyuan Xia; Hong Lin and Qianxia Li
Reported by Author; Author; Author; Author