Next Article in Journal
A New Approach of 3D Lightning Location Based on Pearson Correlation Combined with Empirical Mode Decomposition
Previous Article in Journal
A Weighted Mean Temperature Model with Nonlinear Elevation Correction Using China as an Example
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the BDGIM Performance in BDS Single Point Positioning

1
School of Geography and Information Engineering, China University of Geosciences, Wuhan 430078, China
2
GNSS Research Center, Wuhan University, Wuhan 430079, China
3
College of Marine Science and Technology, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(19), 3888; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193888
Submission received: 28 July 2021 / Revised: 13 September 2021 / Accepted: 18 September 2021 / Published: 28 September 2021

Abstract

:
The broadcast ionospheric model is mainly used to correct the ionospheric delay error for single-frequency users. Since the BeiDou global ionospheric delay correction model (BDGIM) is a novel broadcast ionospheric model for BDS-3, its performance was analyzed through single point positioning (SPP) in this study. Twenty-two stations simultaneously receiving B1C, B2a, B1I and B3I signals were selected from the International GNSS Service (IGS) and the International GNSS Monitoring and Assessment System (iGMAS) tracking networks for the SPP experiments. The differential code bias (DCB) parameters were used to correct the hardware delays in the signals of B1C and B2a. The results showed that the BDGIM performs the best in high-latitude areas, and can effectively improve the positioning accuracy compared with the Klobuchar model. The average 3D positioning accuracy of the four civil signals can reach 3.58 m in high-latitude areas. The positioning accuracies with the BDGIM in the northern hemisphere are better than those in the southern hemisphere, and the global average 3D positioning accuracy of the four civil signals is 4.60 m. The performance of the BDGIM also shows some seasonal differences. The BDGIM performs better than the Klobuchar model on the days of spring equinox and winter solstice, while the opposite is true on the days of summer solstice and autumn equinox. On the day of winter solstice, the average 3D accuracies with the BDGIM on the signals of B1C, B2a, B1I and B3I are 4.13 m, 5.32 m, 4.40 m and 4.49 m, respectively. Although the SPP accuracies are to some extent affected by the geomagnetic storm, the BDGIM generally performs better and are more resistant to the geomagnetic storm than the Klobuchar model.

Graphical Abstract

1. Introduction

Ionosphere is one of the main error sources affecting the signals of global navigation satellite system (GNSS), and the magnitude of ionospheric effects could reach tens of meters [1,2]. In GNSS data processing, the ionospheric error should not be negligible, and efforts have been made to avoid its adverse effects on positioning. Although the ionosphere-free combination is usually employed by dual-frequency users to eliminate the first-order ionospheric term, it always leads to significant increase in the noise level [3,4]. When only single-frequency measurements are available, the ionospheric error has to be corrected with proper ionospheric models [2]. The choice of ionospheric models is of great importance in obtaining precise positioning results. The American global positioning system (GPS) and the European Galileo navigation satellite system (Galileo) adopt the Klobuchar model and improved NeQuick-G model, respectively [1,5,6]. An approach was proposed to reduce the computational load of the NeQuick-G algorithm by increasing the update interval of the ionospheric corrections [7]. China’s BeiDou navigation satellite system (BDS) adopts both the Klobuchar model and the BeiDou global ionospheric delay correction model (BDGIM) [2,6,8]. It is worth noting that the Klobuchar models with GPS and BDS are different in the term of coordinate system. The Klobuchar model with BDS is referenced to the geographic coordinate system, while that with GPS is to the geomagnetic coordinate system [9,10].
Due to the simplicity in calculation and the performance in mid-latitude areas, the Klobuchar model has been widely used and extensively studied. By improving the nighttime term and the amplitude of the cosine term, Bi et al. [11] proposed a modified Klobuchar model without any additional coefficients and the new model can reduce the ionospheric error by 60% over the polar region. Based on the k-means clustering of ionospheric daily variations, a correction model was proposed and was able to improve the performance of Klobuchar model by 36.24% in mid-latitude areas [12]. Liu et al. [13] adjusted the number of model parameters to adapt areas of different scales, and the regional correction rates of the Klobuchar-like models exceeded 90%. Considering the applicability and the issue of broadcasting, the traditional Klobuchar model with eight parameters are still the most popular.
The information of the eight-parameter Klobuchar model used by BDS-2 is contained in D1 and D2 navigation messages and is broadcast every two hours [14]. The D1 navigation message is broadcast by the B1I signals of the medium Earth orbit (MEO) and the inclined geo synchronous orbit (IGSO) satellites, while the D2 navigation message is broadcast by the B1I signals of the geostationary Earth orbit (GEO satellites) [14]. Adopted by BDS-3, the BDGIM contains nine parameters, which are broadcast through the navigation messages by the B1C and B2a signals, respectively [15,16].
To evaluate the performance of ionospheric parameters of BDS, Zhang et al. [17] took the global ionospheric map (GIM) model released by the Center for Orbit Determination in Europe (CODE) as a reference and found that more than 70% of the ionospheric error could be corrected by the BDS Klobuchar model, which performs better in the northern hemisphere and in the mid-latitude areas [10]. Over the polar area, the eight-parameter Klobuchar model of BDS is remarkably outperformed by that of the GPS [18]. The initial performance assessment by Yuan et al. [8] shows that the BDGIM is able to mitigate the ionospheric errors by 80.9% in China and by 77.6% on the global scale, outperforming both the GPS Klobuchar model and the NeQuick Galileo model. With precise ionospheric grid data released by the CODE and slant ionospheric delay derived from dual-frequency observations, Zhou et al. [10] conducted comparison among the performances of the BDS Klobuchar model, the GPS Klobuchar Model and the BDGIM. It was found that the BDGIM in China region is comparable to the BDS Klobuchar model, and that the BDGIM exhibit an overall advantage on a global scale. Similar results were achieved by Yu et al. [2], and the global accuracy of the BDGIM is about 3.6 total electron content units (TECU). Wang et al. re-estimated the BDGIM parameters and compared the vertical total electron content (VTEC) with the VTECs from the GIM and Jason-2/3 altimetry missions to find that the overall correction capability of the BDGIM was better than 75% for over 90% of all the observed samples [19].
Since the adoption of ionospheric model could lead to an improvement in the precision of positioning with single-frequency observables, the results of single point positioning (SPP) are used to evaluate the performance of different ionospheric models [9,20,21,22]. The SPP experiments based on GPS L1 C/A code pseudoranges suggest that the positioning performance with the BDGIM is comparable to the NeQuick-C corrections and is better than GPS ionospheric correction algorithm in both quiet and disturbed conditions [19]. The SPP results of B1 and B2 code measurements with the BDS Klobuchar model are 26.3% and 28.2% better than those with the GPS Klobuchar model [23]. For the users in the China region, the accuracy of SPP with BDS Klobuchar model is better than that with GPS Klobuchar model [24,25,26]. The opposite results are also achieved due to the heavy multipath and poor orbit accuracy of the GEO satellites, and a regional satellite clock model and a regional ionospheric model are able to help improve the SPP accuracy by 50% [27]. For BDS-3 users, SPP accuracy with the BDGIM is overall better than that with the BDS Klobuchar model, especially in vertical component [10,28]. However, the performances of both the BDGIM and the BDS Klobuchar model are undermined in some extreme ionospheric environments, such as magnetic storm and polar region [6,18,29].
The existing studies with the BDGIM performance are mainly based on either the GPS code pseudoranges or the legacy BDS signals, while more comprehensive research on the BDGIM with BDS-3 satellites and the novel signals are still necessary. This work is to further investigate the performance of the BDGIM in term of SPP accuracy. From the users’ perspective, we not only compare the different ionospheric models, but also pay intensive attention on the applicability of the BDGIM to all the four civil signals of BDS. Immediately after this introduction section, the algorithm of ionospheric correction with the BDGIM is introduced, and the data processing strategies are described in detail. In the following section, the performances of SPP with the four civil signals of BDS-3 are analyzed, and the BDGIM are compared with the BDS Klobuchar model and the global ionosphere maps (GIM). In addition, the performance of SPP with the BDGIM during magnetic storm is discussed. This is followed by the conclusions.

2. The Mathematical Model of SPP with the BDGIM

Only code measurements are used for SPP, and the ionospheric delay can be significantly reduced through the linear combination of dual-frequency measurements. We focus on single-frequency scenarios so as to investigate the performance of the ionospheric models.

2.1. Observation Equation of SPP

Generally, the code measurement of BDS can be modeled as [30]
p r , j s = ρ r s + T r s + d t r d τ s + I r , j s + d r , j D j s + ε r , j s
where, s , r and j indicate the numbers or names of satellite, receiver and frequency, respectively, ρ denotes the geometric distance between satellite and receiver, T denotes the tropospheric delay, d t r and d τ s denote receiver clock and satellite clock, respectively, I is the ionospheric delay, d and D are the hardware delays at receiver and satellite ends, respectively, and ε represents the sum of the unmodeled errors such as observing noise and multipath effects. The variable of epoch time is omitted, and all the terms are in the dimension of length. When no confusion is caused, the superscripts and subscripts might be omitted in the following discussion without additional instruction.
To form the observation equation, the geometric distance is linearized through Taylor expansion at the approximate coordinate, and the corrections of coordinates will be estimated. The parameter of receiver clock is also to be estimated, and the hardware delay at receiver end will be lumped into it. Therefore, the observation equation can be rewritten as
p r , j s ρ 0 r s T r s + d τ ˜ s δ j s I r , j s = ι x v x + ι y v y + ι z v z + d t ˜ r + ε r , j s
where ρ 0 r s is the computed geometric distance between the satellite and approximate coordinate of the receiver, and d τ ˜ s is the satellite clock computed from broadcast ephemeris. Since the BDS satellite clock in the broadcast ephemeris is referenced to the B3I signal, there is d τ ˜ s = d τ s + D B 3 I s and the differential code biases (DCBs) δ j s should be considered for the other signals. ι x v x , ι y v y and ι z v z are the corrections of coordinates multiplied by their coefficients. d t ˜ r is the sum of the receiver clock and corresponding hardware delay at the receiver end. In our experiment, the tropospheric delay is computed with empirical model, and so is the ionospheric delay which will be discussed below in details. The left-hand side is actually the observed-minus-computed (OMC) values.

2.2. Correction of Hardware Delay

As is mentioned above, the hardware delay of the B3I code measurement has already been included in the broadcast satellite clock, and therefore no additional correction is needed, meaning δ B 3 I s equals to zero. For the other signals, the code hardware delays at satellite end have to be corrected either by the timing group delay (TGD) parameters in broadcast ephemeris or post-time DCB products.
According to the specification of the format RINEX 3.04, the TGD parameter for B1I signal is recorded in the broadcast ephemeris of BDS, while the corrections for B1C and B2a signals are not. DCB products can be used as the substitute of the TGD parameters, since their equivalence has been verified through both theoretic induction and experimental and analysis [30,31,32]. Taking the pilot components of the B1C and B2a signals as examples, the DCB parameters δ B 1 C s and δ B 2 a s in Equation (2) are calculated as
{ δ B 1 C s = Δ C 1 P , C 6 I s δ B 2 a s = Δ C 1 P , C 6 I s Δ C 1 P , C 5 P s
where Δ C 1 P , C 6 I s is the DCB between the C1P and C6I signals for satellite s , and Δ C 1 P , C 5 P s is the DCB between the C1P and C5P signals. These two DCB parameters are those which can be looked up in the DCB products.

2.3. Algorithm of the BDGIM

Similar to the tropospheric delay, the ionospheric delay for single-frequency SPP users is also computed with empirical model, and it is dependent on the total electron content (TEC) and the signal frequency. The spherical harmonic functions are able to effectively describe the subtle variations in TEC, and are widely used in the modeling and prediction of global TEC in solar geomagnetic frame [33].
Based on the modified spherical harmonics method, the algorithm of the BDGIM is detailed in the interface control documents (ICDs) for open service signals B1C and B2a, released by China Satellite Navigation Office (CSNO) in December 2017 [15,16]. The BDGIM coefficients are divided into two groups [2,10]. The nine low-order spherical harmonic function coefficients are broadcast through the navigation messages, while the other 17 coefficients, along with the calculation method, are provided in the ICDs. The nine broadcast coefficients, playing a dominant role in the prediction of TEC, are calculated with the global ionospheric background prediction model and the ground tracking data of BDS master control stations and A/B monitoring stations in China [34]. The 17 non-broadcast coefficients are supposed to describe the small-scale variations of TEC [2].
With the BDGIM, the line-of-sight ionospheric delay (in meters) of the signal with frequency f can be calculated as
I = M F 40.28 × 10 16 f [ A 0 + i = 1 9 α i A i ]
where α i ( i = 1 9 ) are the nine broadcast coefficients, and f is the carrier frequency of the signal. The M F is the ionospheric mapping function for the conversion from vertical to slant TEC. The values of A i are calculated as [15,16]
A i = { N n i , m i P | n i | , | m i | sin φ cos ( m i λ ) m i 0 N n i , m i P | n i | , | m i | sin φ sin ( m i λ ) m i < 0
while A 0 , the predictive vertical TEC (VTEC), are calculated as [15,16]
{ A 0 = j = 1 17 β j B j β j = a 0 , j + k = 1 12 ( a k , j cos ( ω k t p ) + b k , j sin ( ω k t p ) ) ω k = 2 π / T k B j = { N n j , m j P | n j | , | m j | sin φ cos ( m j λ ) m j 0 N n j , m j P | n j | , | m j | sin φ sin ( m j λ ) m j < 0
In Equations (5) and (6), φ and λ are the geomagnetic latitude and longitude of ionospheric pierce point (IPP) in the solar-fixed reference frame. a k , j and b k , j are the non-broadcast coefficients of BDGIM, and T k are the corresponding periods for prediction. t p is the odd hour of the modified Julian day. P | n i | , | m i | are the un-normalized Legendre functions with degree | n i | and order | m i | , and N n i , m i are the normalization functions. According to the ICDs, the values of n i and m i are listed in Table 1 [15,16].
The mapping function in Equation (4) can be expressed as
M F = 1 / 1 ( R e R e + H ion cos E ) 2
where, E is the elevation angle of the satellite, R e is the mean radius of the Earth, and H ion is the altitude of the ionospheric single-layer shell. Shown in Figure 1 is the flow chart for computing the ionospheric delay with the BDGIM.

3. Processing Strategies

At present, the available signals of BDS include B1I, B3I, B1C and B2a, and it is necessary to investigate the adaptability of BDGIM on different signals. Since the performance of BDGIM can be reflected in single-frequency positioning results, experiments of SPP with different civil signals are conducted.

3.1. Data Selection

Data was collected during the days of year (DOYs) 23–46, 2020 (i.e., from 23 January to 15 February 2020) at 22 global tracking stations of the International GNSS Service (IGS) multi-GNSS experiment (MGEX) and the international GNSS Monitoring and Assessment System (iGMAS). During the 24 days, the Kp index stayed around 2.0, the 3-h forecast value never exceeded 4.0. The distribution of the stations is shown in Figure 2, and the detailed information of the 22 tracking stations is listed in Table 2. All of the four civil signals are available at each station. The broadcast ephemeris with RINEX format Version 3.04 are obtained from the website of the IGS, and the coefficients of both the BDS Klobuchar model and the BDGIM are provided by the Test and Assessment Research Center (TARC) of CSNO.
The numbers of visible BDS satellites and the position dilution of precision (PDOP) at each station are shown in Figure 3. The average numbers of BDS-3 and BDS-2/3 satellites per epoch are displayed as blue squares and green circles, respectively, while the top and bottom of the vertical lines represent the corresponding maximum and minimum values. The average PDOP with BDS-2/3 satellites at each station is indicated by a grey triangle. For all the stations, an average of no less than six BDS-3 satellites can be tracked, although the numbers of visible BDS-3 satellites might sometimes be less than four at the stations REDU, PADO, CEBR, NNOR and MGUE. When the BDS-2 satellites are also counted, the numbers of visible satellites increase, especially for the stations in Asia-Pacifica area. The average PDOP values are better than 2.2 for all the stations except CNYR, RDJN and MGUE, and the worst average PDOP value is 4.0 at the station MGUE.

3.2. SPP Strategies

When code measurements are used for SPP, the satellite equipment group delay ought to be considered. The equipment group delay of B3I is regarded as the reference and is included in the clock parameter. Therefore, no additional correction is needed for B3I code measurement. The differential equipment group delay between B1I and B3I signals is given as the TGD parameters broadcast in the navigation message. The signals of B1C and B2a contain pilot components and data components. For the pilot components of B1C and B2a signals, the equipment group delays also need to be compensated by the TGD parameters. If data component is used, the group delay differential between the data component and the corresponding pilot component has to be additionally compensated by the inter-signal correction (ISC) parameters besides the TGD parameter. Both the TGD and the ISC parameters are broadcast in the navigation message.
Since the TGD and the ISC parameters for B1C and B2a signals are not recorded in the current ephemeris files, the differential code bias (DCB) parameters were used as the substitute to correct the equipment delays in B1C and B2a code measurements. In the study, we adopted the multi-GNSS DCB products released by China Academy of Sciences (CAS).
The accuracy of the GIMs provided by the IGS analysis centers could reach 2–5 TECU, equivalent to a zenith delay of 0.32–0.80 m in GPS code measurements of L1 frequency [35,36]. Therefore, the GIM released by CODE were employed as the control. The GIM and broadcast ionospheric models were used in SPP, respectively, and the positioning results were compared and analyzed. The data processing strategies are listed in Table 3, and the 95th percentile of the absolute values of positioning errors was adopted to measure the daily positioning accuracy of each station. According to the user range accuracy (URA) information in the broadcast ephemeris and the specifications in the ICDs, the a priori precisions of BDS code measurements were set as 2.40 m and they were weighted with the elevation angles [14,15,16]. For the selected stations, the coordinates extracted from the officially issued software independent exchange (SINEX) files were taken as the reference true position for comparison.

4. Results

The global performances of SPP with the four BDS civil signals were investigated, and the SPP results with the BDGIM, the BDS Klobuchar model and the GIM were compared. To further study the adaptability of the BDGIM in different seasons, the positioning performances on four typical days were analyzed. In addition, the performance of the BDGIM during the geomagnetic storm was also analyzed.

4.1. Global Positioning Accuracies

Shown in Figure 4 and Figure 5 are the accuracies in each direction (i.e., east, north and up) and the three-dimensional (3D) accuracies of SPP with different models and different signals. In each panel, arranged in the far left is the northernmost station, and the other stations are sorted rightwards by latitude from north to south. In Figure 4, the results of SPP with the BDGIM, the Klobuchar model and the GIM are displayed in the panels of left, middle and right columns, and each row of panels represent the results with each of the four BDS signals. The SPP accuracies in east, north and up directions are displayed as red, green and blue bars, respectively. In Figure 5, the 3D accuracies of SPP with the BDGIM, the Klobuchar model and the GIM are displayed as red circles, green stars and blue squares, respectively. Generally, the horizontal accuracies are significantly better than the vertical for all the stations except CANB, and the east components are usually the most accurate among all the three directions. The SPP accuracies of most of the stations in the northern hemisphere are better than their counterparts in the southern hemisphere, which might be explained by the relatively abundant monitoring stations in the northern hemisphere. With more stations in a certain area, a larger amount of data can be used to monitor the ionosphere, and the model parameters are theoretically able to describe the ionospheric variations much better. It is also worth noting that around the polar regions both the BDGIM and the GIM perform better than the Klobuchar model.
To investigate the performance of the BDGIM in the regions of different latitudes, we categorized all the 22 stations into groups of low-latitude (0°–30°), mid-latitude (30°–60°) and high-latitude (60°–90°), and for each group the average accuracies were calculated. Listed in Table 4 and Table 5 are the 3D accuracies of SPP with the novel and the legacy signals, as well as the improvement percentages of the accuracies with the BDGIM relative to those with the Klobuchar model. On the global scale, the average accuracies of SPP with the BDGIM for the signals of B1C, B2a, B1I and B3I reach 4.10 m, 5.35 m, 4.39 m and 4.57 m, respectively, improving by 24.9%, 37.0%, 15.4% and 29.0% compared with the corresponding results using the Klobuchar model. For all of the four signals, the best average accuracies of SPP with the BDGIM are achieved at the high-latitude stations, while the opposite is true with the Klobuchar model.
From the perspective of average SPP accuracies, the Klobuchar model performs better at the low-latitude stations, while the BDGIM performs significantly better at the high-latitude stations. At the mid-latitude stations, the BDGIM performs better than the Klobuchar model for all the signals except B1I. Among the four signals, the BDGIM almost always performs the best for the B1C signal, either regionally or globally. At the high-latitude stations, the highest average positioning accuracy, reaching 2.98 m, is obtained with the BDGIM on the B1C signal, and the average positioning accuracy of the four signals reaches 3.58 m. Compared with the Klobuchar model, the average SPP accuracies at the high-latitude stations with the BDGIM on the signals of B1C, B2a, B1I and B3I are improved by 78.2%, 83.2%, 70.5% and 79.4%, respectively. When the Klobuchar model is adopted, the average SPP accuracy at the high-latitude stations seems the worst, especially at the two arctic stations, CNYR and KIRU, whose accuracies exceed 10 m, as are shown in Figure 4 and Figure 5.
At the mid-latitude stations, the BDGIM performs 5.5%, 19.8% and 7.3% better than the Klobuchar model for the signals of B1C, B2a and B3I, respectively, while is 5.7% outperformed by the Klobuchar model only for the signal of B1I. For all of the four civil signals, the Klobuchar model outperforms the BDGIM at the low-latitude stations by around 30%. For each signal, the best accuracy of SPP with the BDGIM is obtained at the high-latitude stations, and the SPP accuracy seems to increase with the latitude.
To investigate the difference between the two hemispheres, each of the abovementioned station groups is divided into the northern and the southern subgroups. Shown in Figure 6 are the average accuracies with the BDGIM, the Klobuchar model and the GIM in the six different regions of the world. The results show that the global average 3D positioning accuracy of BDGIM reaches 4.87 m, and that of GIM reaches 3.43 m. In the northern hemisphere, the positioning accuracies with the BDGIM are similar to those with the GIM, and the differences between them are 0.20 m, 0.63 m and 0.78 m in high-, mid- and low-latitude areas, respectively. The SPP accuracies with the BDGIM in the southern hemisphere are not as good as their counterparts in the northern hemisphere. The differences between the accuracies with the BDGIM and the GIM are 1.11 m, 2.16 m and 3.79 m in high-, mid- and low-latitude areas, respectively.
The sharp deterioration of the positioning accuracy with the Klobuchar model in high-latitude areas can be explained from the following two aspects. In the BeiDou Klobuchar model, the value of the parameter α is relatively large, and the latitude of pierce point is not constrained. The combined effect of the two factors leads to significant increase in the calculated amplitude in high-latitude areas. To avoid these adverse impacts, the BDGIM is based on the method of modified spherical harmonics and calculates the amplitude with normalized Legendre function.
Figure 7 presents the typical time series of positioning errors at the station CLGY on 25 January 2020 (DOY 25), and the three models are used to correct the ionospheric delays in all of the four signals. For the signals of B2a and B3I, the BDGIM performs significantly better than the Klobuchar model, while for the signals of B1C and B1I the performances of the two models are at the same level. Although the positioning accuracies with the GIM are the highest among the three models, the GIM cannot be obtained in real time. Since the GIM is a post-time high-precision product with a certain time delay, it is not as convenient as the broadcast ionospheric models.

4.2. Positioning Accuracies in Different Seasons

To investigate the performance of the BDGIM in different seasons, four typical days (i.e., 20 March 2020, 20 June 2020, 22 September 2020 and 22 December 2019) were selected to represent spring, summer, autumn and winter. The average SPP accuracies of the 22 stations with the four signals on the four typical days are shown in Figure 8. Note that the values larger than 10 m are not displayed. It can be seen that the BDGIM generally performs the best on the day of winter solstice and the worst on the day of autumn equinox. On the day of winter solstice, the average 3D accuracies of the signals of B1C, B2a and B3I are 4.13 m, 5.32 m and 4.49 m, respectively, and they are better than those of all the other three days. For the signal of B1I, the 3D accuracies on the three days other than autumn equinox are almost the same, i.e., around 4.40 m, and that on the day of autumn reaches 6.31 m.
Listed in Table 6 and Table 7 are the average 3D accuracies of the 22 stations with the four signals and the three ionospheric models in different seasons. Since the GIM is post-time ionospheric products with high precision, it performs the best in all seasons for the four signals. The BDGIM performs better than the Klobuchar model on the days of spring equinox and winter solstice, while the opposite is true on the days of summer solstice and autumn equinox. On the day of spring equinox, the average 3D accuracies with the BDGIM on the signals of B1C, B2a, B1I and B3I signals are 8.23 m, 5.78 m, 4.43 m and 4.74 m, respectively, i.e., 44.1%, 71.8%, 59.7% and 70.1% better than those with the Klobuchar model. On the day of winter solstice, the average 3D accuracies with the BDGIM on the signals of B1C, B2a, B1I and B3I are 4.13 m, 5.32 m, 4.40 m and 4.49 m, respectively, i.e., 16.1%, 27.9%, 6.2% and 20.5% better than those with the Klobuchar model.

4.3. Positioning Accuracies during the Geomagnetic Storm

A geomagnetic storm can not only cause strong disturbances in the geomagnetic field, but also influence the ionosphere. According to the alerts issued by the Space Environment Prediction Center of CAS, the 3-h forecast Kp index reached 5.0 twice on 20 April, 2020 (DOY 111), indicating an event of a geomagnetic storm. To investigate the performance of the BDGIM around the geomagnetic storm, we conducted SPP experiments with the BDS data collected on the DOYs 110–112, 2020.
Shown in Figure 9 are the average positioning accuracies of the four signals before, during and after the occurrence of the geomagnetic storm. Generally, the geomagnetic storm has perceptible but slight influence on SPP with the BDGIM, and the two novel signals and the zenith direction are affected more obviously. On the very day of the geomagnetic storm, the positioning accuracies of B1C and B2a signals decrease by 5.6% and 6.7% with respect to the previous day. In addition, the positioning accuracies of the two days before and after the occurrence of the geomagnetic storm are almost the same. For the signals of B1I and B3I, the 3D accuracies stay around 5 m all through the three days, and the horizontal accuracies are around 2 m. It seems that the positioning accuracies with the two legacy signals are resistant to the influence of the geomagnetic storm.
Shown in Figure 10 are the time series of positioning errors at the station GUA1 during the three days around the geomagnetic storm. On the day of the geomagnetic storm, the vertical accuracies of SPP with B1C, B2a, B1I and B3I signals decreases by 0.34 m, 0.64 m, 0.17 m and 0.06 m, respectively, while the horizontal accuracies keep stable for all of the four signals.
The averages of accuracies with the four civil signals at each station on the day of the geomagnetic storm (DOY 111) are listed in Table 8, and the positioning accuracies with the BDGIM and the Klobuchar model are compared. The mean accuracies of the 22 stations are also calculated. It can be seen that the mean 3D accuracy with the BDGIM is 6.22 m, much better than that with the Klobuchar model. Considering the poor performance of the Klobuchar model at the three high-latitude stations (CNYR, KIRU and ZHON), the averages without the three stations are also calculated and compared. Even after the three high-latitude stations are excluded, the mean SPP accuracies with the BDGIM and the Klobuchar model are 6.44 m and 7.66 m, respectively. The results show that the BDGIM performs better than the Klobuchar model even during the geomagnetic storm.

5. Conclusions and Discussion

In the study, the performance of the BDGIM in SPP was assessed, and the DCB parameters were used to correct the hardware delays in the signals of B1C and B2a. To quantitatively analyze the SPP accuracies, 24 days of data collected at 22 global stations were processed, and the positioning accuracies were investigated and compared.
On the global scale, the average accuracy of SPP with the BDGIM can reach 4.60 m, which is close to that with the GIM. Compared with Klobuchar model, the positioning accuracies of the four signals are improved by 15.4%–37.0% when the BDGIM is adopted. In high-latitude areas, the BDGIM performs much higher than the Klobuchar model. With the BDGIM, the positioning accuracy of the B1C signal can reach 2.98 m in high-latitude areas, and the average positioning accuracy of the four civil signals is 3.58 m.
In mid-latitude areas, the BDGIM outperforms the Klobuchar model by 5.5%, 19.8% and 7.3% for the signals of B1C, B2a and B3I, respectively. In low-latitude areas, the Klobuchar model performs better than the BDGIM by around 30%. Seeming to increase with the latitude, the accuracies of SPP with BDGIM are similar to those with the GIM in the northern hemisphere, and their counterparts in the southern hemisphere are not as good.
The positioning performance of the BDGIM show seasonal differences and the general SPP accuracy is the best in winter and the worst in autumn. The BDGIM performs better than the Klobuchar model on the days of spring equinox and winter solstice, while the opposite is true on the days of summer solstice and autumn equinox. On the day of winter solstice, the average 3D accuracies with the BDGIM on the signals of B1C, B2a, B1I and B3I are 4.13 m, 5.32 m, 4.40 m and 4.49 m, respectively, i.e., 16.1%, 27.9%, 6.2% and 20.5% better than those with the Klobuchar model.
Although the geomagnetic storm can to some extent affect the accuracies of SPP with the BDGIM, the influences are more perceptible only for the two novel signals and the zenith direction. The accuracies of SPP with the two legacy signals are less affected. On the very day of the geomagnetic storm, the average 3D accuracy of SPP with the BDGIM is 6.22 m, and is much better than that with the Klobuchar model.
In general, the BDGIM performs better in the northern hemisphere and in high-latitude areas than in the southern hemisphere and in low-latitude areas. Seasonal differences are also seen in the performance of the BDGIM, and it outperforms the Klobuchar model on the days of spring equinox and winter solstice. Despite some slight influences, the BDGIM are more resistant to the geomagnetic storm than the Klobuchar model. The advantages of the BDGIM are probably attributed to the spherical harmonics and more favorable distribution of the monitoring stations used in the generation of model parameters. It is also worth noting that the above results and analysis are derived from a relatively limited data set during the time of solar minimum. To evaluate the long-term performance of the BDGIM, further investigation will be conducted in the future.

Author Contributions

G.W. and Z.H. provided the initial idea; Z.Y. and Y.B. designed and performed the experiments; G.C. and W.L. reviewed the data processing methods and helped with the analysis; G.W. and Z.Y. drafted and refined the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the funds of Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University (Grant No. 18-01-06), the National Science Foundation of China (Grant No. 41804033 and 42004073), and the Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan) (Grant No. CUGL180831).

Data Availability Statement

The observations of the MGEX stations and corresponding broadcast ephemeris were obtained from ftp://igs.gnsswhu.cn/pub/gps/data/daily/2020 (accessed on 28 October 2020). The observations of the iGMAS stations were obtained from Wuhan University. The BDGIM and Klobuchar model parameters were obtained from http://www.csno-tarc.cn/datacenter/ephemeris (accessed on 28 October 2020). The GIM products were obtained from ftp://igs.gnsswhu.cn/pub/gps/products/ionex/2020 (accessed on 28 October 2020), while the DCB products were obtained from ftp://igs.ign.fr/pub/igs/products/mgex/dcb/2020 (accessed on 28 October 2020).

Acknowledgments

The authors would like to thank the IGS Multi-GNSS Experiment (MGEX), iGMAS, CAS and CSNO for providing GNSS observations, relevant products and models, all of which enable this study. The authors gratefully acknowledge X. Liu at Fugro Intersite B.V. and H. Lyu at UPC-IonSAT research group, both of whom provided valuable advice and generous encouragement.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Klobuchar, J.A. Ionospheric Time-Delay Algorithm for Single-Frequency GPS Users. IEEE Trans. Aerosp. Electron. Syst. 1987, AES-23, 325–331. [Google Scholar] [CrossRef]
  2. Yu, S.; Ma, F.J.; Wang, W.; Wu, R.P.; Zeng, Q. Preliminary Evaluation on the Precision of the BDS-3 Global Ionospheric Model. In China Satellite Navigation Conference; Sun, J., Yang, C., Yang, Y., Eds.; Lecture Notes in Electrical Engineering; Springer Singapore Pte Ltd.: Singapore, 2019; Volume 563, pp. 131–139. [Google Scholar]
  3. Gioia, C.; Borio, D.; Realini, E.; Gatti, A.; Tagliaferro, G. From Single to Precise Point Positioning: The Impact on Time Retrieval. In Proceedings of the 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS), Portland, OR, USA, 20–23 April 2020; pp. 1104–1110. [Google Scholar]
  4. Paziewski, J.; Sieradzki, R.; Baryla, R. Multi-GNSS high-rate RTK, PPP and novel direct phase observation processing method: Application to precise dynamic displacement detection. Meas. Sci. Technol. 2018, 29, 035002. [Google Scholar] [CrossRef]
  5. Wang, N.; Yuan, Y.; Li, Z.; Huo, X. Improvement of Klobuchar model for GNSS single-frequency ionospheric delay corrections. Adv. Space Res. 2016, 57, 1555–1569. [Google Scholar] [CrossRef]
  6. Zhu, Y.; Tan, S.; Zhang, Q.; Ren, X.; Jia, X. Accuracy evaluation of the latest BDGIM for BDS-3 satellites. Adv. Space Res. 2019, 64, 1217–1224. [Google Scholar] [CrossRef]
  7. Gioia, C.; Borio, D. NeQuick-G and Android Devices: A Compromise between Computational Burden and Accuracy. Sensors 2020, 20, 5908. [Google Scholar] [CrossRef]
  8. Yuan, Y.; Wang, N.; Li, Z.; Huo, X. The BeiDou global broadcast ionospheric delay correction model (BDGIM) and its preliminary performance evaluation results. J. Inst. Navig. 2019, 66, 55–69. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, X.; Hu, X.; Wang, G.; Zhong, H.; Tang, C. Evaluation of COMPASS ionospheric model in GNSS positioning. Adv. Space Res. 2013, 51, 959–968. [Google Scholar] [CrossRef]
  10. Zhou, R.Y.; Hu, Z.G.; Su, M.D.; Li, J.Z.; Li, P.B.; Zhao, Q.L. Preliminary Performance Evaluation of BeiDou Global Ionospheric Delay Correction Model; Geomatics and Information Science of Wuhan University: Wuhan, China, 2019; Volume 44, pp. 1457–1464. (In Chinese) [Google Scholar]
  11. Bi, T.; An, J.; Yang, J.; Liu, S. A modified Klobuchar model for single-frequency GNSS users over the polar region. Adv. Space Res. 2017, 59, 833–842. [Google Scholar] [CrossRef]
  12. Pongracic, B.; Wu, F.; Fathollahi, L.; Brčić, D. Midlatitude Klobuchar correction model based on the k-means clustering of ionospheric daily variations. GPS Solut. 2019, 23, 80. [Google Scholar] [CrossRef]
  13. Liu, C.; Liu, C.J.; Feng, X.; Xu, L.F.; Du, Y. Klobuchar-like Ionospheric Model for Different Scales Areas. Acta Geod. Cartogr. Sin. 2016, 45, 54–63. [Google Scholar] [CrossRef]
  14. China Satellite Navigation Office. BeiDou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal B1I (Version 3.0). Available online: http://www.beidou.gov.cn/xt/gfxz/201902/P020190227593621142475.pdf (accessed on 8 June 2021).
  15. China Satellite Navigation Office. BeiDou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal B1C (Version 1.0). Available online: http://www.beidou.gov.cn/xt/gfxz/201712/P020171226741342013031.pdf (accessed on 8 June 2021).
  16. China Satellite Navigation Office. BeiDou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal B2a (Version 1.0). Available online: http://www.beidou.gov.cn/xt/gfxz/201712/P020171226742357364174.pdf (accessed on 8 June 2021).
  17. Zhang, Q.; Zhao, Q.L.; Zhang, H.P.; Hu, Z.G.; Wu, Y. Evaluation on the Precision of Klobuchar Model for BeiDou Navigation Satellite System; Geomatics and Information Science of Wuhan University: Wuhan, China, 2014; Volume 39, pp. 142–146. (In Chinese) [Google Scholar]
  18. Wang, Z.M.; Bi, T.; Sun, W.; An, J.C.; Liu, S.L. Analysis of GNSS broadcast ionospheric model correction in polar regions. Chin. J. Polar Res. 2016, 28, 235–242. (In Chinese) [Google Scholar]
  19. Wang, N.; Li, Z.; Yuan, Y.; Huo, X. BeiDou Global Ionospheric delay correction Model (BDGIM): Performance analysis during different levels of solar conditions. GPS Solut. 2021, 25, 1–13. [Google Scholar] [CrossRef]
  20. Hoque, M.M.; Jakowski, N.; Berdermann, J.; Inst, N. An Ionosphere Broadcast Model for Next Generation GNSS. In Proceedings of the 28th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2015), Tampa, FL, USA, 14–18 September 2015; pp. 3755–3765. [Google Scholar]
  21. Ye, Y.L.; Liu, C.; Liu, C.J. Evaluation of ionospheric model based on standard point positioning technique. Sci. Surv. Mapp. 2017, 42, 126–131. (In Chinese) [Google Scholar]
  22. Rovira-Garcia, A.; Ibáñez-Segura, D.; Orús-Perez, R.; Juan, J.M.; Sanz, J.; González-Casado, G. Assessing the quality of ionospheric models through GNSS positioning error: Methodology and results. GPS Solut. 2020, 24, 1–12. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, A.H.; Chen, J.P.; Zhang, Y.Z.; Meng, L.D.; Wang, J.X. Performance of Selected Ionospheric Models in Multi-Global Navigation Satellite System Single-Frequency Positioning over China. Remote Sens. 2019, 11, 2070. [Google Scholar] [CrossRef] [Green Version]
  24. Li, W.W.; Li, M.; Hu, Z.G.; Zhao, Q.L. Comparative Analysis of COMPASS and GPS Ionospheric Model on Positioning and Navigation Precision. In Proceedings of the 4th CSNC, Wuhan, China, 15–17 May 2013; pp. 95–99. (In Chinese). [Google Scholar]
  25. Wang, L.J.; Cai, C.L.; Li, Z.B.; Wei, Z.C.; Li, T.S. COMPASS User Positioning Accuracy Evaluation Based on Single Point Pseudo-Range Positioning Principle. In Proceedings of the 4th CSNC, Wuhan, China, 15–17 May 2013; pp. 168–173. (In Chinese). [Google Scholar]
  26. Zhang, Q.; Zhao, Q.L.; Zhang, H.P.; Hu, Z.G.; Liu, Z.Y. Research on BeiDou Navigation Satellite System Ionospheric Model Accuracy. In Proceedings of the 4th CSNC, Wuhan, China, 15–17 May 2013; pp. 271–277. (In Chinese). [Google Scholar]
  27. Zheng, F.; Gu, S.; Gong, X.; Lou, Y.; Fan, L.; Shi, C. Real-time single-frequency pseudorange positioning in China based on regional satellite clock and ionospheric models. GPS Solut. 2020, 24, 6. [Google Scholar] [CrossRef]
  28. Guo, S.R.; Cai, H.L.; Meng, Y.N.; Geng, C.J.; Jia, X.L.; Mao, Y.; Geng, T.; Rao, Y.N.; Zhang, H.J.; Xie, X. BDS-3 RNSS technical characteristics and service performance. Acta Geod. Cartogr. Sin. 2019, 48, 810–821. (In Chinese) [Google Scholar]
  29. Wang, W.; Liu, C.; Wang, Y.P.; Jiang, Y.; Xu, F. Analysis on the influence of geomagnetic storm of BDS pseudorange positioning. GNSS World China 2019, 44, 82–88. (In Chinese) [Google Scholar]
  30. Guo, F.; Zhang, X.; Wang, J. Timing group delay and differential code bias corrections for BeiDou positioning. J. Geod. 2015, 89, 427–445. [Google Scholar] [CrossRef]
  31. Ge, Y.; Zhou, F.; Sun, B.; Wang, S.; Shi, B. The Impact of Satellite Time Group Delay and Inter-Frequency Differential Code Bias Corrections on Multi-GNSS Combined Positioning. Sensors 2017, 17, 602. [Google Scholar] [CrossRef]
  32. Zhang, Y.; Chen, J.; Gong, X.; Chen, Q. The update of BDS-2 TGD and its impact on positioning. Adv. Space Res. 2020, 65, 2645–2661. [Google Scholar] [CrossRef]
  33. Li, Z.; Yuan, Y.; Wang, N.; Hernandez-Pajares, M.; Huo, X. SHPTS: Towards a new method for generating precise global ionospheric TEC map based on spherical harmonic and generalized trigonometric series functions. J. Geod. 2015, 89, 331–345. [Google Scholar] [CrossRef]
  34. Yang, Y.; Tang, J.; Montenbruck, O. Chinese Navigation Satellite Systems. In Springer Handbook of Global Navigation Satellite Systems; Teunissen, P.J.G., Montenbruck, O., Eds.; Springer: Cham, Switzerland, 2017; pp. 273–304. [Google Scholar]
  35. Hernández-Pajares, M.; Juan, J.M.; Sanz, J.; Orus, R.; Garcia-Rigo, A.; Feltens, J.; Komjathy, A.; Schaer, S.C.; Krankowski, A. The IGS VTEC maps: A reliable source of ionospheric information since 1998. J. Geod. 2009, 83, 263–275. [Google Scholar] [CrossRef]
  36. Li, Z.S.; Wang, N.B.; Li, M.; Zhou, K.; Yuan, Y.B.; Yuan, H. Evaluation and analysis of the global ionospheric TEC map in the frame of international GNSS services. Chin. J. Geophys. 2017, 60, 3718–3729. [Google Scholar] [CrossRef]
Figure 1. Computation of the ionospheric delay with the BDGIM.
Figure 1. Computation of the ionospheric delay with the BDGIM.
Remotesensing 13 03888 g001
Figure 2. Distribution of the 22 MGEX (blue triangle) and iGMAS (red square) stations from which the data was used in this study.
Figure 2. Distribution of the 22 MGEX (blue triangle) and iGMAS (red square) stations from which the data was used in this study.
Remotesensing 13 03888 g002
Figure 3. The average numbers of visible BDS-3 (blue square) and BDS-2/3 (green circle) satellites with the vertical lines indicating the corresponding extreme values, as well as the average PDOPs (grey triangle).
Figure 3. The average numbers of visible BDS-3 (blue square) and BDS-2/3 (green circle) satellites with the vertical lines indicating the corresponding extreme values, as well as the average PDOPs (grey triangle).
Remotesensing 13 03888 g003
Figure 4. The east (red), north (green) and up (blue) accuracies of SPP with the BDGIM (left), the Klobuchar model (middle) and the GIM (right) using the signals of B1C, B2a, B1I and B3I at the selected 22 stations, which are sorted by latitude from north to south.
Figure 4. The east (red), north (green) and up (blue) accuracies of SPP with the BDGIM (left), the Klobuchar model (middle) and the GIM (right) using the signals of B1C, B2a, B1I and B3I at the selected 22 stations, which are sorted by latitude from north to south.
Remotesensing 13 03888 g004
Figure 5. The 3D positioning accuracies of SPP with the BDGIM (red circle), the Klobuchar model (green star) and the GIM (blue square) at the 22 stations. The four panels represent the signals of B1C, B2a, B1I and B3I.
Figure 5. The 3D positioning accuracies of SPP with the BDGIM (red circle), the Klobuchar model (green star) and the GIM (blue square) at the 22 stations. The four panels represent the signals of B1C, B2a, B1I and B3I.
Remotesensing 13 03888 g005
Figure 6. The average 3D positioning accuracies with the BDGIM, the Klobuchar model and the GIM in different regions of the world. The letters H, M and L represent high, middle and low latitudes, respectively, while N and S represent the northern and southern hemispheres, respectively.
Figure 6. The average 3D positioning accuracies with the BDGIM, the Klobuchar model and the GIM in different regions of the world. The letters H, M and L represent high, middle and low latitudes, respectively, while N and S represent the northern and southern hemispheres, respectively.
Remotesensing 13 03888 g006
Figure 7. The time series of SPP errors at the station CLGY on January 25, 2020. The four rows represent the signals of B1C, B2a, B1I and B3I from top to bottom, and in each row the left, the middle and the right panels represent the BDGIM, the Klobuchar model and the GIM. In each panel, the time series and statistics corresponding to the directions of east (E), north (N) and up (U) are displayed in red, green and blue, respectively.
Figure 7. The time series of SPP errors at the station CLGY on January 25, 2020. The four rows represent the signals of B1C, B2a, B1I and B3I from top to bottom, and in each row the left, the middle and the right panels represent the BDGIM, the Klobuchar model and the GIM. In each panel, the time series and statistics corresponding to the directions of east (E), north (N) and up (U) are displayed in red, green and blue, respectively.
Remotesensing 13 03888 g007aRemotesensing 13 03888 g007b
Figure 8. The average SPP accuracies with the BDGIM on the days of spring equinox (red), summer solstice (green), autumn equinox (blue) and winter solstice (orange). The four panels represent the signals of B1C (upper left), B2a (upper right), B1I (lower left) and B3I (lower right).
Figure 8. The average SPP accuracies with the BDGIM on the days of spring equinox (red), summer solstice (green), autumn equinox (blue) and winter solstice (orange). The four panels represent the signals of B1C (upper left), B2a (upper right), B1I (lower left) and B3I (lower right).
Remotesensing 13 03888 g008
Figure 9. The average positioning accuracies of SPP with the novel (top) and the legacy (bottom) signals before (red), during (green) and after (blue) the day on which geomagnetic storm occurred.
Figure 9. The average positioning accuracies of SPP with the novel (top) and the legacy (bottom) signals before (red), during (green) and after (blue) the day on which geomagnetic storm occurred.
Remotesensing 13 03888 g009
Figure 10. The time series of SPP errors at the station GUA1 before (left column), during (middle column) and after (right column) the day on which geomagnetic storm occurred. The four rows of panels represent the signals of B1C, B2a, B1I and B3I from top to bottom. In each panel, the time series and statistics corresponding to the directions of east (E), north (N) and up (U) are displayed in red, green and blue, respectively.
Figure 10. The time series of SPP errors at the station GUA1 before (left column), during (middle column) and after (right column) the day on which geomagnetic storm occurred. The four rows of panels represent the signals of B1C, B2a, B1I and B3I from top to bottom. In each panel, the time series and statistics corresponding to the directions of east (E), north (N) and up (U) are displayed in red, green and blue, respectively.
Remotesensing 13 03888 g010aRemotesensing 13 03888 g010b
Table 1. The values of n i and m i .
Table 1. The values of n i and m i .
i 123456789
n i / m i 0/01/01/11/−12/02/12/−12/22/−2
Table 2. Detailed information of the 22 GNSS tracking stations.
Table 2. Detailed information of the 22 GNSS tracking stations.
StationLatitude (°)Longitude (°)LocationReceiver TypeSource
CNYR78.9211.93ArcticUNICORE UB4B0IiGMAS
KIRU67.8620.97Kiruna, SwedenSEPT POLARX5IGS
BRCH52.2610.53Braunschweig, GermanyCETC-54-GMR-4016iGMAS
ICUK51.51−0.13London, EnglandCETC-54-GMR-4016iGMAS
CLGY51.05−114.07Calgary, CanadaCETC-54-GMR-4016iGMAS
REDU50.005.14Redu, BelgiumSEPT POLARX5IGS
PADO45.4111.90Padova, ItalySEPT POLARX5IGS
CHU143.79125.44Changchun, ChinaBD070iGMAS
GUA143.4787.18Urumqi, ChinaBD070iGMAS
CEBR40.45−4.37Cebreros, SpainSEPT POLARX5TRIGS
BJF139.61115.89Beijing, ChinaCETC-54-GMR-4016iGMAS
XIA134.37109.22Xian, ChinaGNSS_GGRiGMAS
SHA131.10121.20Shanghai, ChinaUNICORE UB4B0IiGMAS
WUH130.52114.49Wuhan, ChinaCETC-54-GMR-4016iGMAS
LHA129.6691.10Lhasa, ChinaCETC-54-GMR-4016iGMAS
KNDY7.2980.63Kandy, Sri LankaCETC-54-GMR-4016iGMAS
FAA1−17.56−149.61Papeete, SOCISEPT POLARX5IGS
RDJN−22.91−43.20Rio, BrazilUNICORE UB4B0IiGMAS
NNOR−31.05116.19New Norcia, AustraliaSEPT POLARX5TRIGS
CANB−35.30149.20Canberra, AustraliaCETC-54-GMR-4011iGMAS
MGUE−35.78−69.40Malargue, ArgentinaSEPT POLARX5TRIGS
ZHON−69.3776.37AntarcticCETC-54-GMR-4011iGMAS
Table 3. The data processing strategies of single-frequency SPP.
Table 3. The data processing strategies of single-frequency SPP.
ItemsSettings
Sampling interval30s
Elevation cutoff
Orbits and clocksBroadcast ephemeris
Tropospheric delaySaastamoinen model
Ionospheric delayKlobuchar/BDGIM/GIM
Hardware delayB1I: TGD; B1C/B2a: DCB
WeightElevation-based
EstimatorWeighted least squares
Table 4. Regional and global averages of the 3D accuracies of SPP with the BDGIM and the Klobuchar model for the BDS novel signals B1C and B2a.
Table 4. Regional and global averages of the 3D accuracies of SPP with the BDGIM and the Klobuchar model for the BDS novel signals B1C and B2a.
RegionsB1CB2a
BDGIM
(m)
Klobuchar
(m)
Improvement Percentage
(%)
BDGIM
(m)
Klobuchar
(m)
Improvement Percentage
(%)
High-latitude2.9813.6878.24.0023.7783.2
Mid-latitude3.954.185.54.896.1019.8
Low-latitude5.504.12−33.58.055.99−34.4
Global4.105.4624.95.358.4937.0
Table 5. Regional and global averages of the 3D accuracies of SPP with the BDGIM and the Klobuchar model for the BDS legacy signals B1I and B3I.
Table 5. Regional and global averages of the 3D accuracies of SPP with the BDGIM and the Klobuchar model for the BDS legacy signals B1I and B3I.
RegionB1IB3I
BDGIM
(m)
Klobuchar
(m)
Improvement Percentage
(%)
BDGIM
(m)
Klobuchar
(m)
Improvement Percentage
(%)
High-latitude3.6612.4070.53.6617.7779.4
Mid-latitude4.254.02−5.74.224.557.3
Low-latitude5.484.17−31.46.575.00−31.4
Global4.395.1915.44.576.4429.0
Table 6. The average 3D positioning accuracies of different models for novel BDS signals in different seasons (unit: m).
Table 6. The average 3D positioning accuracies of different models for novel BDS signals in different seasons (unit: m).
SeasonB1CB2a
BDGIMKlobucharGIMBDGIMKlobucharGIM
Spring8.2314.717.665.7820.513.88
Summer8.848.238.165.564.993.47
Autumn10.288.458.368.775.703.71
Winter4.134.923.295.327.383.64
Table 7. The average 3D positioning accuracies of different models for old BDS signals in different seasons (unit: m).
Table 7. The average 3D positioning accuracies of different models for old BDS signals in different seasons (unit: m).
SeasonB1IB3I
BDGIMKlobucharGIMBDGIMKlobucharGIM
Spring4.4311.003.794.7415.833.56
Summer4.383.563.335.134.233.50
Autumn6.313.623.608.094.803.85
Winter4.404.693.504.495.652.99
Table 8. The average accuracies of each station during the geomagnetic storm (unit: m).
Table 8. The average accuracies of each station during the geomagnetic storm (unit: m).
StationBDGIMKlobuchar
ENU3DENU3D
CNYR1.711.605.265.7619.0724.4159.2066.81
KIRU1.181.984.194.7811.8225.3436.6846.12
BRCH1.342.745.616.392.088.038.7712.07
ICUK1.293.035.916.771.868.288.4812.00
CLGY1.331.815.075.551.856.8410.8512.96
REDU1.412.715.065.912.878.646.1610.99
PADO1.142.314.295.001.785.934.597.71
CHU10.961.434.544.861.222.524.765.52
GUA10.721.092.332.670.751.863.303.86
CEBR1.122.355.756.311.243.004.145.26
BJF11.682.074.445.181.552.394.935.69
XIA11.021.183.103.470.951.313.463.82
SHA12.153.026.307.312.052.996.127.11
WUH11.973.536.767.881.883.156.787.71
LHA12.834.579.4510.872.794.189.1810.47
KNDY1.722.716.196.971.572.125.045.69
FAA13.673.5711.0212.153.523.706.758.46
RDJN2.572.846.677.692.562.898.108.97
NNOR1.051.714.354.791.112.065.976.41
CANB9.143.485.3811.169.104.474.9711.29
MGUE1.191.503.784.241.241.514.464.87
ZHON1.742.003.144.117.2713.4117.9123.53
Mean1.952.425.396.223.646.3210.4812.77
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, G.; Yin, Z.; Hu, Z.; Chen, G.; Li, W.; Bo, Y. Analysis of the BDGIM Performance in BDS Single Point Positioning. Remote Sens. 2021, 13, 3888. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193888

AMA Style

Wang G, Yin Z, Hu Z, Chen G, Li W, Bo Y. Analysis of the BDGIM Performance in BDS Single Point Positioning. Remote Sensing. 2021; 13(19):3888. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193888

Chicago/Turabian Style

Wang, Guangxing, Zhihao Yin, Zhigang Hu, Gang Chen, Wei Li, and Yadong Bo. 2021. "Analysis of the BDGIM Performance in BDS Single Point Positioning" Remote Sensing 13, no. 19: 3888. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193888

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