Next Article in Journal
A Feasible Calibration Method for Type 1 Open Ocean Water LiDAR Data Based on Bio-Optical Models
Next Article in Special Issue
Cycle Slip Detection during High Ionospheric Activities Based on Combined Triple-Frequency GNSS Signals
Previous Article in Journal
Flash Flood Risk Analysis Based on Machine Learning Techniques in the Yunnan Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Ionospheric Delay Effects on Multi-GNSS Positioning Performance

1
Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
School of Remote Sensing and Geomatics Engineering, Nanjing University of Information Science and Technology, Nanjing 210044, China
4
Institute of Communications and Navigation, German Aerospace Center, Kalkhorstweg 53, 17235 Neustrelitz, Germany
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(2), 171; https://doi.org/10.3390/rs11020171
Submission received: 22 December 2018 / Revised: 14 January 2019 / Accepted: 16 January 2019 / Published: 17 January 2019

Abstract

:
Ionospheric delay is a significant error source in multi-GNSS positioning. We present different processing strategies to fully exploit the ionospheric delay effects on multi-frequency and multi-GNSS positioning performance, including standard point positioning (SPP) and precise point positioning (PPP) scenarios. Datasets collected from 10 stations over thirty consecutive days provided by multi-GNSS experiment (MGEX) stations were used for single-frequency SPP/PPP and dual-frequency PPP tests with quad-constellation signals. The experimental results show that for single-frequency SPP, the Global Ionosphere Maps (GIMs) correction achieves the best accuracy, and the accuracy of the Neustrelitz TEC model (NTCM) solution is better than that of the broadcast ionospheric model (BIM) in the E and U components. Eliminating ionospheric parameters by observation combination is equivalent to estimating the parameters in PPP. Compared with the single-frequency uncombined (UC) approach, the average convergence time of PPP with the external ionospheric models is reduced. The improvement in BIM-, NTCM- and GIM-constrained quad-constellation L2 single-frequency PPP was 15.2%, 24.8% and 28.6%, respectively. The improvement in convergence time of dual-frequency PPP with ionospheric models was different for different constellations and the GLONASS-only solution showed the least improvement. The improvement in the convergence time of BIM-, NTCM- and GIM-constrained quad-constellation L1/L2 dual-frequency PPP was 5.2%, 6.2% and 8.5%, respectively, compared with the UC solution. The positioning accuracy of PPP is slightly better with the ionosphere constraint and the performance of the GIM-constrained PPP is the best. The combination of multi-GNSS can effectively improve the positioning performance.

Graphical Abstract

1. Introduction

Ionospheric delay effect is one of the main errors in global navigation satellite system (GNSS) positioning, navigation and timing (PNT), and it is normally several or tens of meters although it can exceed 100 m under the severe ionospheric disturbances [1,2]. Nowadays, several ionospheric delay models can be used for real-time or post-processed applications. Among these ionospheric delay models, the Klobuchar model, which is driven by broadcast ephemeris, is widely used by GPS users [3]. The improved Klobuchar model is slightly different and has been specially adapted for Beidou satellites. Similarly, the NeQuick model has also been developed for Galileo users [4]. Additionally, the Neustrelitz TEC Model (NTCM) is a global 2D empirical ionospheric model that can be operated autonomously without external ionospheric measurement [5]. The only driving parameter of the NTCM is the solar radio flux index F10.7. Global Ionosphere Maps (GIMs) are routinely estimated by the International GNSS service (IGS) Analysis Centers (ACs) and have an accuracy of 2–8 TECU (total electron content unit).
With the rapid development of current navigation constellations and construction of more stations, GNSS represents a new era and has wide applications all over the world. In this new multi-frequency and multi-GNSS environment, observations from different signals and constellations are processed together, however, a consideration of the error sources is necessary and mandatory to ensure that the positioning model is consistent. Standard point positioning (SPP) and precise point positioning (PPP) are two types of positioning modes with a stand-alone GNSS receiver that have been used to fulfill different requirements. SPP has been widely used in many fields such as land surveying and vehicle navigation. In the case of the single-frequency SPP approach, the GNSS observations are corrected with existing ionospheric models to ensure a positioning accuracy of meters. The PPP technique has also proved to be an efficient technique due to its high accuracy, global coverage and cost-efficiency. The ionosphere-free (IF) combination is an effective and popular way to eliminate the first order of ionosphere delay in dual-frequency PPP [6]. Another approach is to use the undifferenced and uncombined (UC) observations in dual-frequency PPP processing to extract ionospheric delays and avoid noise amplification [7]. The ionospheric delay is estimated by employing a prior ionosphere model and has the potential to improve positioning performance for wider applications [8,9]. For instance, the experimental tests of [10] showed that the convergence time of dual-frequency PPP was reduced by 30% with the IGS GIM product constraint. Along with dual-frequency PPP, single-frequency PPP has also attracted increasing attention due to its low cast and high accuracy. Similarly, IF combination on a single frequency, known as the Group and Phase Ionospheric Correction (GRAPHIC) can be utilized to eliminate the ionospheric delay in single-frequency PPP [11]. The rank-deficient mathematical problem exists in the GRAPHIC approach, which utilizes the average of the carrier phase and pseudo-range observations. Ionospheric delays can also be estimated simultaneously along with other parameters when proper constraints are given by the ionospheric models. Several studies have been done to improve the convergence time of single-frequency PPP with the GIM [8,12,13]. The quality of ionospheric delay derived from the models will affect the positioning performances. Thus, ionospheric delay is a key issue to resolve in multi-GNSS combined positioning.
In this study, we first present the ionospheric processing strategies for current multi-GNSS single point positioning and summarize the positioning models. Multi-GNSS combined positioning models consist of single-frequency SPP/PPP and dual-frequency PPP. Thereafter, three types of ionospheric correction models, including the broadcast ionospheric model (BIM), NTCM and GIM are introduced for multi-GNSS positioning scenarios. Then, the data and processing strategy is given. A unified time, varying weight, ionospheric scheme for PPP is presented to adapt three ionospheric models. Finally, a comprehensive analysis of the impact of ionospheric delay correction on multi-GNSS performances is performed.

2. Materials and Methods

2.1. Multi-GNSS Positioning Models

The simplest positioning method, SPP is where pseudo-range observations made by a receiver are processed to acquire its position with the orbits of GNSS satellites and their clock offsets, as computed from the broadcast navigation message. Point positioning based on the pseudo-range and phase observations as well as precise orbit and clock products is known as PPP. The quad-constellation GNSS raw pseudo-range and carrier phase observation on the j frequency can be expressed as [14]:
p r , j s = ρ r s + d t r d t s + T r s + I r , j s + d r , j d j s + ε p
φ r , j s = ρ r s + d t r d t s + T r s I r , j s + λ j w r s + λ j N r , j s + λ j ( b r , j b j s ) + ε φ
where the superscript s denotes a GNSS satellite; the subscript r denotes the receiver; p r , j s denotes the observed pseudo-range; φ r , j s is the corresponding carrier phase; ρ r s denotes the computed geometrical range; d t r is the receiver clock offset; dts is the satellite clock offset; T r s is the slant tropospheric delay; λ j is the wavelength of carrier phase; and I r , j s is the slant ionospheric delay on jth frequency. It is possible to convert the first order ionospheric delays at different frequencies by I r , j s = ( λ j 2 / λ 1 2 ) I r , 1 s ; d r , j and d j s are the uncalibrated code delays (UCDs) of the receiver and the satellite on jth frequency; w r s is the phase wind-up delay; N r , j s is the integer ambiguity; b r , j and b j s are the uncalibrated phase delays (UPDs) for the receiver and satellites, respectively. ε p and ε φ are the pseudo-range and carrier phase observation noises including multipath, respectively.
For the following statement, we define the notations as:
{ α 1 , 2 s = ( f 1 s ) 2 / ( ( f 1 s ) 2 ( f 2 s ) 2 ) , β 1 , 2 s = ( f 2 s ) 2 / ( ( f 1 s ) 2 ( f 2 s ) 2 ) , DCB 1 , 2 s = d 1 s d 2 s , DCB r , 1 , 2 = d r , 1 d r , 2 , d I F s = α 1 , 2 s d 1 s + β 1 , 2 s d 2 s , d r , I F = α 1 , 2 s d r , 1 + β 1 , 2 s d r , 2 , b I F s = α 1 , 2 s λ 1 b 1 s + β 1 , 2 s λ 2 b 2 s , b r , I F = α 1 , 2 s λ 1 b r , 1 + β 1 , 2 s λ 2 b r , 2 ,
where f 1 s and f 2 s are the signal frequency; α 1 , 2 s and β 1 , 2 s are the frequency-dependent factors; DCB 1 , 2 s and DCB r , 1 , 2 are the satellite and receiver differential code bias (DCB).

2.1.1. Single-Frequency SPP

It is possible to find a SPP solution based on at least 4 pseudo-range observations and by utilizing dual- or multi-frequency IF combinations, which eliminate the ionospheric delay. For a single-frequency SPP, ionospheric corrections are usually computed by an external model. In general, with m observed satellites and considering the satellite orbits and clocks, satellite timing group delay and atmosphere delays, the single-frequency SPP model can be given for a single epoch as [15]:
[ p r , 1 s , 1 ρ ¯ r s p r , 1 s , m ρ ¯ r s ] = [ u r s 1 ] [ x d t ¯ r ] + [ ε p ] , Q L
d t ¯ r = d t r + d r , 1
where ρ ¯ r s denotes the geometrical range including the measurement errors; u r s is the unit vector of the component from the receiver to the satellite; 1 is a vector with corresponding rows and one column; x is the vector of the receiver position increments relative to the a priori position; Q L denotes the stochastic model of observed minus computed values of observations; The estimable receiver clock d t ¯ r is a combination of true clock offset plus the receiver hardware delay.
The redundancy of the SPP model is the number of observables minus the number of estimable parameters: m s ( 3 + 1 ) = m s 4 (satellites), m s denotes the number of observables. Hence, the SPP model by least square (LS) is solvable for m s 4 .

2.1.2. Single-Frequency PPP

One crucial problem in single-frequency PPP, is that the ionospheric delay cannot be mitigated by the combination of observations on different frequencies. In the GRAPHIC approach, the arithmetic mean of the pseudo-range and carrier phase is considered as the observation where the ionospheric delays are canceled out since the pseudo-range and carrier phase delays have the same values but opposite signs. Thus, the observations with m observed satellites can be expressed as
[ ( p r , 1 s , 1 + λ 1 φ r , 1 s , 1 ) / 2 ρ ¯ r s ( p r , 1 s , m + λ 1 φ r , 1 s , m ) / 2 ρ ¯ r s ] = [ u r s 1 M W R 1 ] [ x d t ¯ r Z w N ¯ r , 1 s ] + [ ε p , φ ] , Q L
d t ¯ r = d t r + d r , 1 / 2 + λ 1 b r , 1 / 2
N ¯ r , 1 s = ( λ 1 ( N r , 1 s b 1 s ) + d I F s ) / λ 1 / 2
where M W is the wet mapping function of tropospheric delay; R 1 denotes the corresponding wavelength of carrier phase to the ambiguity parameters on the first frequency; Z w is the tropospheric zenith wet delay (ZWD); and N ¯ r , 1 s denotes the estimated float ambiguity values on the first frequency.
For the single-frequency UC PPP, the pseudo-range and carrier phase observations with m observed satellites for a single epoch can be written as
[ p r , 1 s , 1 ρ ¯ r s φ r , 1 s , 1 ρ ¯ r s p r , 1 s , m ρ ¯ r s φ r , 1 s , m ρ ¯ r s ] = [ u r s 1 M W K R 1 ] [ x d t ¯ r Z w I r , 1 s N ¯ r , 1 s ] + [ ε p ε φ ] , Q L
d t ¯ r = d t r + d r , 1
N ¯ r , 1 s = N r , 1 s + b r , 1 b 1 s d r , 1 / λ 1 + d I F s / λ 1
where the matrix K denotes the corresponding coefficient of slant ionospheric delay, the element for pseudo-range observation is 1 while the element for carrier phase is −1.
Compared to single-frequency UC PPP, the ionosphere-constrained single-frequency PPP adds virtual ionospheric observations from the external model with their corresponding constraints, which can be expressed as
[ p r , 1 s , 1 ρ ¯ r s φ r , 1 s , 1 ρ ¯ r s p r , 1 s , m ρ ¯ r s φ r , 1 s , m ρ ¯ r s I ˜ r , 1 s , 1 I ˜ r , 1 s , m ] = [ u r s 1 M W K R 1 O O O I O ] [ x d t ¯ r Z w I r , 1 s N ¯ r , 1 s ] + [ ε p ε φ ε i o n ] , Q L , Q I
where I ˜ r , 1 s is the virtual ionospheric observation derived from external ionospheric models; O denotes zero matrix; I is the identity matrix; and Q I denotes the stochastic model of virtual ionospheric observations.

2.1.3. Dual-Frequency PPP

The IF combination is often utilized for dual-frequency PPP to remove the first order ionospheric delay, which can be expressed as
[ α 1 , 2 s p r , 1 s , 1 + β 1 , 2 s p r , 2 s , 1 ρ ¯ r s α 1 , 2 s φ r , 1 s , 1 + β 1 , 2 s φ r , 2 s , 1 ρ ¯ r s α 1 , 2 s p r , 1 s , m + β 1 , 2 s p r , 2 s , m ρ ¯ r s α 1 , 2 s φ r , 1 s , m + β 1 , 2 s φ r , 2 s , m ρ ¯ r s ] = [ u r s 1 M W R 1 ] [ x d t ¯ r Z w N ¯ r , I F s ] + [ ε p ε φ ] , Q L
d t ¯ r = d t r + d r , I F
N ¯ r , I F s = ( α 1 , 2 s λ 1 N r , 1 s + β 1 , 2 s λ 2 N r , 2 s + b r , I F b I F s d r , I F + d I F s ) / λ 1
where N ¯ r , I F s denotes the estimated float ambiguity values.
For the dual-frequency UC PPP, the receiver UCD will be absorbed by the receiver clock offset and slant ionospheric delay and the model can be expressed as
[ p r , 1 s , 1 ρ ¯ r s φ r , 1 s , 1 ρ ¯ r s p r , 2 s , m ρ ¯ r s φ r , 2 s , m ρ ¯ r s ] = [ u r s 1 M W K R 1 R 2 ] [ x d t ¯ r Z w I ¯ r , 1 s N ¯ r , 1 s N ¯ r , 2 s ] + [ ε p ε φ ] , Q L
{ d t ¯ r = d t r + d r , I F I ¯ r , 1 s = β 1 , 2 s ( DCB 1 , 2 s DCB r , 1 , 2 ) λ 1 N ¯ r , 1 s = λ 1 ( N r , 1 s + b r , 1 b 1 s ) + d I F s d r , I F + β 1 , 2 s ( DCB 1 , 2 s DCB r , 1 , 2 ) λ 2 N ¯ r , 2 s = λ 2 ( N r , 2 s + b r , 2 b 2 s ) + d I F s d r , I F + ( f 1 s / f 2 s ) 2 β 1 , 2 s ( DCB 1 , 2 s DCB r , 1 , 2 )
For the ionosphere-constrained dual-frequency PPP, virtual ionospheric observations with their constraints are added and the model can be expressed as
[ p r , 1 s , 1 ρ ¯ r s φ r , 1 s , 1 ρ ¯ r s p r , 2 s , m ρ ¯ r s φ r , 2 s , m ρ ¯ r s I ˜ r , 1 s , 1 I ˜ r , 1 s , m ] = [ u r s 1 J M W K R 1 R 2 O O O O I O O ] [ x d t ¯ r DCB r , 1 , 2 Z w I r , 1 s N ¯ r , 1 s N ¯ r , 2 s ] + [ ε p ε φ ε i o n ] , Q L , Q I
{ λ 1 N ¯ r , 1 s = λ 1 ( N r , 1 s + b r , 1 b 1 s ) + d I F s d r , I F λ 2 N ¯ r , 2 s = λ 2 ( N r , 2 s + b r , 2 b 2 s ) + d I F s d r , I F
In matrix J , the element for p r , 1 s ρ ¯ r s is set to β 1 , 2 s , and the element for p r , 2 s ρ ¯ r s is set to α 1 , 2 s , to acquire the corresponding receiver DCB.

2.2. Ionospheric Correction Models

Ionospheric delay models are generally categorized as an empirical model or a mathematical functions model. The former is based on data observed in long-term records and represents the characteristic variation patterns while the latter is fitted by mathematical functions and utilize the actual measured ionospheric delay of a certain area for a period of time.

2.2.1. BIM

The Klobuchar model corrects the ionospheric delay by broadcast ephemeris, which has the advantages of convenient calculation and simple structure [3]. The model considers the periodic and amplitude variations of the ionosphere on a daily scale when setting the parameters, and reflects the characteristic variations of the ionosphere to ensure the reliability of large-scale ionospheric forecasts. The drawbacks of the model are that it is only suitable for the mid-latitude area and it has limited accuracy on the ionospheric delay correction. The input parameters of the Klobuchar model are the eight model coefficients, geodetic latitude and longitude of the GNSS antenna, GPS time, as well as the azimuth and the elevation of the observed satellite.
The improved Klobuchar model is slightly different from the GPS Klobuchar model as it is especially used for Beidou to adapt to the Chinese region. The broadcast parameters of the Beidou Klobuchar are derived based on the monitoring stations in China. The method of Beidou Klobuchar can be found in the Beidou navigation satellite system signal in space interface control document [16].
The NeQuick model is a time dependent ionospheric electron density and three-dimensional model and is applicable for real-time single-frequency correction for the Galileo satellite navigation system. It provides the electron concentration in the ionosphere and total electronic content (TEC) along the ground to satellite ray-path by means of numerical integration [17]. The input parameters of the NeQuick model include the geodetic position of the GNSS antenna, GPS observing time and effective ionospheric level factor AZ determined by three coefficients broadcasted by Galileo ephemeris. The description of the Galileo NeQuick model is available at https://www.gsc-europa.eu/system/files/galileo_documents/Galileo_Ionospheric_Model.pdf.
In the BIM of this study, the GPS and GLONASS satellite ionospheric delay is corrected by the Klobuchar model, the Beidou satellite ionospheric delay is corrected by the improved Klobuchar model, and the NeQuick model is used for Galileo ionospheric delay correction.

2.2.2. NTCM

The NTCM is a user-friendly TEC model that covers all levels and a global scale of solar activity [18]. This empirical approach describes the TEC/ionospheric delay dependencies on the geographic location, local time, and solar irradiance and activity. Considering the irradiation angle of the sun, the TEC variation with the local time is divided into diurnal, semi-diurnal and ter-diurnal harmonic components in the NTCM. The solar activity level is provided by the solar radio flux index F10.7. The 12 model coefficients, together with the covariance matrix in NTCM are fixed and determined by LS adjustment with observation data sets. The standard deviations for the coefficients are utilized to compute 95% confidence intervals. The input parameters of the NTCM are the F10.7 index, day of year, local time in hours, and the geographic and geomagnetic latitude in degrees [5].

2.2.3. GIM

The GIM assumes that the ionosphere is composed of a thin spherical layer at a height of 450 km above the earth’s surface [19]. The GIM provides vertical TEC (VTEC) at certain grid points over the globe. The GIM reconstructed from the predicted and actual measured datasets can be utilized for real-time and post-processed applications, respectively. The vertical electron content at the ionospheric pierce point (IPP) location is generally obtained by interpolation of VTEC values from four neighboring grid points. Hence, space and time interpolations are necessary to obtain the VTEC with high accuracy. The VTEC at an IPP needs to be firstly interpolated between two consecutive maps and then mapped to the slant ionospheric delay. The method for the interpolation and mapping of GIM is explained in detail in [20].

3. Data and Processing Strategy

To evaluate the ionospheric delay effects on multi-GNSS combined positioning performances, datasets from 10 stations of the MGEX network for 1–30 September 2018 (day of year (DOY), from 244 to 273) were selected and utilized for numerical analysis. Figure 1 shows the geographical distribution of the selected MGEX stations. Figure 2 provides the F10.7 values and geomagnetic Kp index for 1–30 September 2018. All of the stations receive quad-system observations from GPS, Beidou, GLONASS and Galileo constellations. The quad-constellation satellite orbits and clock offsets are corrected by the broadcast ephemris provided by MGEX, which is available at ftp://cddis.gsfc.nasa.gov/pub/gps/data/campaign/mgex/daily/rinex3/2018/brdm/, or the precise orbit and clock offset products provided by Deutsches GeoForschungsZentrum (GFZ). The sampling rate of orbit and clock products are 5 min and 30 s, respectively. The detailed data processing strategies for multi-GNSS positioning are shown in Table 1. When integrated multi-GNSS combined positioning is carried out, the additional system time difference is needed for each newly added GNSS system. To maintain the consistency of the pseudo-range and carrier phase observations of the quad constellation, we use b1 to denote observations on the first frequency, b2 to denote observations on the second frequency, and b3 to denote observations on the third frequency. The cut-off angle is set to 7° in this study. The stochastic models of the observations represented by a sinusoidal function and based on the elevation angle of the satellite are utilized to reduce the effect of these elevation angle related errors [21]. The GPS and GLONASS phase observation precision is set to 0.003 m and the GPS code observation precision is set to 0.3 m. To reduce the effect of the inter-frequency bias (IFB) in GLONASS, the GLONASS code observation precision is set to 0.9 m. The Beidou and Galileo code observation precision is set to 0.9 m and 0.4 m, respectively, and the phase observation precision is set to 0.004 m. The LS method is used for SPP data processing and the Kalman filter is utilized in the PPP process of parameter estimation. The tropospheric ZWD is estimated as a random walk process and the receiver clock is estimated as white noise. The time system difference parameters are estimated as constant. The ambiguities are estimated as constant for each epoch.
The determination of the weight of virtual ionospheric parameters will affect the ionosphere-constrained PPP performances. In this study, we define the variance of virtual ionospheric parameters as:
σ 2 = ( 1 + 0.2 Δ t ) σ i o n 2 m f ( z ) 2
m f ( z ) = [ 1 sin 2 z / ( 1 + H i o n / R ) 2 ] 1 / 2
where Δ t denotes the epoch number in the filter processing. σ i o n 2 is the mean VTEC spatial variance. Hion is the height of the ionosphere single layer. R is the average radius of the Earth. Equation (20) indicates that the virtual ionospheric observations are weighed at a relatively larger weight at the beginning, and then the weight gradually decreases to get a better positioning accuracy after convergence [12]. The a priori variance of ionospheric delay is determined by the precision of the ionospheric models. Since the Klobuchar model in BIM can remove the ionospheric delay by more than 50%, the variance for BIM is uniformly set to the square of the corresponding BIM ionospheric correction value [3]. The variance of NTCM is set to 1.22 m2 since the model has a root mean square (RMS) deviation of 7.5 TECU [5]. With regard to the GIM, a temporal and spatial correlation function is utilized, which can be expressed as in [22].
σ i o n 2 = { σ i o n , 0 2 ,     t < 8   o r   t > 20   o r   B > π / 3 σ i o n , 0 2 + σ i o n , 1 2 cos ( B ) cos ( t 14 12 π ) , o t h e r
where t is the local time at IPP in hours, B denotes the satellite elevation angle. Since the GIM products have an accuracy of 2–8 TECU, the variance of the zenith ionospheric delay σ i o n , 0 2 and the zenith ionospheric delay variation σ i o n , 1 2 are set to 0.09 m2 and 0.09 m2, respectively.
To obtain the single-frequency SPP solution, three different schemes, ionosphere delay correction with BIM, NTCM and GIM were applied. The rapid GIM products for real-time ionospheric correction were provided by the Center for Orbit Determination in Europe (CODE). For single-frequency PPP, five different schemes including GRAPHIC approach, UC approach, BIM-, NTCM- and GIM-constrained approaches were utilized. For dual-frequency PPP, five schemes including the IF approach, UC approach, BIM-, NTCM- and GIM-constrained approaches, were also applied. The ionosphere processing strategy is summarized in Table 2.

4. Results and Discussion

4.1. Performance of Single-Frequency SPP

Figure 3 depicts the positioning errors of the CUT0 station on 1 September 2018 for the north (N), east (E) and up (U) component of the b1 and b2 based nine different constellation combinations SPP. It can be seen that the SPP results with b1 and b2 code observations behave similarly. The horizonal error of single-frequency SPP with broadcast ephemeris can reach meter-level, while the vertical error is relatively larger. It is obvious that the accuracy of single-frequency SPP can be improved with multi-GNSS combination. The single-frequency SPP with GIM based ionospheric correction shows the best accuracy, especially in the U direction.
The RMS of different single-frequency SPP solutions were calculated, and the mean values of all the tests are summarized in Table 3. As shown in Table 3, the positioning accuracy of multi-GNSS single-frequency SPP with GIM correction performs the best because the GIM is based on the interpolation of the TEC data observed on the selected day and has the highest model accuracy. The solutions from NTCM correction performed better than the BIM correction in the E and U components. Only Beidou, Galileo, GPS/Beidou, and GPS/Galileo b3 SPP solutions were tested, considering that part of the GPS satellites and GLONASS have no triple-frequency signal. For all the single-system single-frequency SPP solutions, Beidou SPP performed the worst, mainly because the satellite orbit and clock offset of the Beidou Geostationary Earth Orbit (GEO) have a relatively lower accuracy. The GPS/Beidou single-frequency SPP achieved slightly better positioning accuracy than the GPS SPP because the Beidou pseudo range observations were assigned less weight than the GPS observations. For the GPS/Galileo b3-based SPP solution, the performance was mainly influenced by the Galileo satellites because few GPS satellites have the b3 signal. For GPS/Beidou/GLONASS/Galileo b1-based SPP, the positioning accuracy was improved by 17.6%, 6.3% and 12.6%, respectively, in the N, E and U components, after the GIM ionospheric correction and compared with the BIM correction. For GPS/Beidou/GLONASS/Galileo b2-based SPP, the positioning accuracy was improved by 28.0%, 12.3% and 31.0%, respectively, in the N, E and U components. For GPS/Galileo b3-based SPP, the positioning accuracy was improved by 30.9%, 21.3% and 16.2%, respectively, in the N, E and U components. For all frequency cases, the multi-GNSS single-frequency SPP performances are improved for more satellites and constellations and the subsequently improved satellite geometry.

4.2. Performance of Single-Frequency PPP

Table 4 shows the parameter estimation method in single-frequency quad-constellation PPP models, where m, n, p and q denote the satellite number of the GPS, Beidou, GLONASS and Galileo, respectively. The method of eliminating ionospheric parameters is actually equivalent to estimating the parameters [28]. The variance-covariance matrix and the solution vector are identical though they may apply different algorithms [29]. Hence, the GRAPHIC and UC approach are equivalent provided the variance-covariance matrix is transformed according to the law of covariance propagation. The ionosphere-constrained approaches will produce different results for the priori ionospheric information provided by different models.
Figure 4 shows the comparison of b2 single-frequency PPP with different schemes for GPS, GLONASS and GPS/Beidou/GLONASS/Galileo solutions at station MR01 on 2 September 2018. The value errors of respective solutions are shifted by 1 m in order to avoid overlapping. For the single-frequency PPP approach, the positioning filter is considered to have converged when the horizontal components of the error reach 0.3 m and keep within 0.3 m. Because the vertical positioning errors are larger than the horizontal components, the criterion is enlarged to 0.5 m in the vertical direction. The convergence time for N, E and U components are also given and marked in Figure 4. It is obvious that the GRAPHIC and UC approaches have the same convergence time for GPS-only, GLONASS-only and quad-constellation PPP. Besides, the convergence time of ionosphere-constrained PPP solutions is less than the UC solutions at MRO1 for three schemes. The convergence time of the GIM-constrained solutions is the least for its highest model accuracy. The convergence time of quad-constellation single-frequency is reduced for the improvement of spatial geometry with more GNSS satellite compared with the GPS-only solution.
Figure 5 shows the average convergence time of b1, b2 and b3 single-frequency PPP with different schemes on all days over all the test stations. Table 5 summarizes the RMS errors of the single-frequency PPP in multi-constellation combinations, which is calculated from the convergence epoch to the end epoch of a day. The GLONASS-only and Beidou-only solutions clearly perform worse than other solutions, which is caused by the strong correlation between IFBs and slant ionospheric delays in GLONASS and the current distribution and orbit determination accuracy of the Beidou constellation [14]. The results also indicate that the integration of the multi-GNSS can improve the single-frequency PPP performance in terms of convergence time and positioning accuracy. As it is shown in Figure 5 and Table 5, the GRAPHIC and UC approaches have much the same performance for the equivalence of two solutions. By introducing external ionosphere information, the average convergence time of PPP solutions is reduced. For instance, the improvement in convergence time of BIM-, NTCM- and GIM-constrained quad-constellation b2 single-frequency PPP is 15.2%, 24.8% and 28.6%, respectively, compared with the UC solution. The positioning accuracy of the quad-constellation single-frequency PPP can reach 1–3 cm horizontally and 8–9 cm vertically. The improvement in the positioning accuracy of the ionosphere-constrained PPP, for giving higher weighting after convergence, was not significant. The positioning accuracy of the GIM-constrained PPP solutions is the best and NTCM-constrained PPP solutions perform better than the BIM-constrained solutions in the E and U components.

4.3. Performance of Dual-Frequency PPP

Table 6 shows the parameter estimation method in dual-frequency quad-constellation PPP models. The first order of the ionospheric effect has been eliminated in the IF model by the observation combination. However, the combined ambiguities are not integers any more, and the combined observations have higher standard deviations. Similarly, the IF and UC approaches are equivalent as the variance-covariance matrix follows the law of covariance propagation and the ionosphere constrained approaches also leads to different results.
Figure 6 provides an intuitive comparison of b1/b2 dual-frequency PPP with different schemes for GPS, GLONASS and GPS/Beidou/GLONASS/Galileo solutions at station SEYG on 3 September 2018. The value errors of the respective solutions were shifted by 0.4 m in order to avoid overlapping. The convergence criterion is defined when the component of the error is less than 0.1 m and stays within 0.1 m in the subsequent epochs. During the period 00:00 to 02:00 UTC, approximately 8–10 and 6–8 satellites can be seen at SEYG in each epoch for GPS-only and GLONASS-only PPP and the corresponding variation position dilution of precision (PDOP) is 1.5–2.6 and 1.7–3.0, respectively. The observed satellite number of quad-constellation PPP are between 27 and 32 and the PDOP values are below 1.0 and very stable. It is clear that the IF and UC approaches have almost the same performance at SEYG. The convergence time of ionosphere constraint GPS-only PPP is less than the UC approach thanks to the external ionospheric information at SEYG. The convergence time of ionosphere constrained GLONASS-only PPP is little longer compared with the UC solution, which is mainly caused by the existing IFBs in GLONASS though the pseudo-range observations are assigned a relatively lower weighting. Similarly, the IFBs of GLONASS can also explain why the convergence time of quad-constellation PPP in E and U components is not effectively reduced with the external ionospheric information at SEYG.
The receiver DCB is estimated as an unknown parameter in the ionosphere-constrained approaches. Figure 7 shows the time series of the estimated GPS, Beidou, GLONASS and Galileo receiver DCB at SEYG. It should be noted that the b1/b3 DCB of GLONASS is not given because the GLONASS has no triple-frequency signal. The multi-GNSS MGEX DCB provided by the Chinese Academy of Sciences (CAS) (available at: ftp://igs.ign.fr/pub/igs/products/mgex/dcb/) is also shown in Figure 7 for reference and comparison, together with the standard deviation (STD) and RMS of corresponding DCB (see top of each panel). It can be seen that the estimated DCB of GIM-constrained approach is much swifter than the other two schemes. Owing to growing weighting for virtual ionospheric observation in the filter process, the existing errors of the ionospheric models have little influence on the positioning accuracy after convergence and are absorbed in the estimable receiver DCB [30]. Hence, as it is shown in Figure 7, this also reflects that the GIM model has the highest accuracy in three models at SEYG.
Figure 8 shows the average convergence time of b1/b2 and b1/b3 dual-frequency PPP with different schemes. Table 7 summarizes the RMS errors of dual-frequency PPP in multi-constellation combinations, the average values of all PPP tests after convergence are calculated. It can be seen that the b1/b2 based PPP performs better than the b1/b3 solutions. The performance of b1/b2 GPS/Beidou and GPS/Galileo PPP differs greatly from the b1/b3 solutions in that only GPS BLOCK-IIF can transmit signals on three frequencies. The positioning accuracy reaches approximately 0.5–2.5 cm in the horizontal and 2–5 cm in the vertical for dual-frequency PPP. Similar to single-frequency PPP, the combination of multi-GNSS speeds up the convergence process and improves the positioning accuracy. For instance, the convergence of quad-constellation b1/b2 GIM-constrained PPP is obviously improved by 40.6% from 13.88 to 8.21 min, compared with the GPS-only solution. The performance of the IF and UC PPP approaches are generally consistent for the equivalence of two models. Introducing different ionospheric models as constraints will reduce the convergence time of PPP schemes in different levels. We can see that the improving average convergence performance for GLONASS-only PPP with ionosphere constraints is the lowest compared with other solutions, where the performances of part stations like the results of Figure 6 may even deteriorate, which is primarily attributed to the existing GLONASS pseudo-range IFBs. The results for dual-frequency ionosphere-constrained GLONASS PPP achieved by [13] is slightly different, which is mainly attributed to different GLONASS code observation precision. In three ionosphere-constrained approaches, the GIM-constrained PPP has the best convergence performance and NTCM-constrained PPP performs better than the BIM-constrained solution. For quad-constellation b1/b2 PPP schemes, the improvement in the convergence time of BIM-, NTCM- and GIM-constrained is 5.2%, 6.2% and 8.5%, respectively, compared with the UC solution. A small but not obvious improvement can be seen in the positioning accuracy with the ionosphere constraint, mainly due to the relatively higher weighting after convergence. The GIM-constrained PPP has the best positioning accuracy in the different schemes.

5. Conclusions

To fully evaluate the ionospheric delay effects on multi-GNSS positioning performances, multi-GNSS positioning models for single (b1, b2 and b3) and dual-frequency (b1/b2 and b1/b3) GNSS signals were assessed by different ionosphere processing schemes. The positioning models were extended to SPP or PPP processing with observations from single-, dual-, triple- and quad-constellations. Static datasets collected at ten stations over thirty consecutive days were utilized to evaluate the impact on positioning performance with ionospheric delay models.
Comparative analysis showed that the positioning accuracy of single-frequency SPP with the external ionospheric model correction can obtain meter-level accuracy, and the vertical error is relatively larger than the horizontal components. For single-frequency SPP, the GIM correction solution achieves the best accuracy and the accuracy of SPP with NTCM correction is better than the solution with the BIM correction in the E and U components. The multi-GNSS single-frequency PPP results indicate that the GRAPHIC and UC approaches have the same performance levels. The average convergence time of single-frequency PPP with external ionosphere correction is reduced. Compared with the quad-constellation b2 single-frequency UC approach, the improvement in the convergence time of BIM-, NTCM- and GIM-constrained solutions is 15.2%, 24.8% and 28.6%, respectively. The improvement in positioning accuracy of the ionosphere-constrained PPP is not obvious for the given higher weighting after convergence. The positioning accuracy of the GIM-constrained PPP is the best and NTCM-constrained PPP solutions perform better than the BIM-constrained solutions in the E and U components. For multi-GNSS dual-frequency PPP, the performance of IF and UC solutions are the same. The improvement in convergence time of dual-frequency PPP varies for different constellations. The existing GLONASS pseudo-range IFBs cause the improvement in the GLONASS-only solutions to be the lowest of all schemes. Compared with the quad-constellation b1/b2 UC approach, the improvement in convergence time of the BIM-, NTCM- and GIM-constrained solution is 5.2%, 6.2% and 8.5%, respectively. The positioning accuracy of the GIM-constrained PPP is the best in all positioning schemes. Also, the combination of multi-GNSS can improve the positioning performance for all three positioning schemes.

Author Contributions

Conceptualization of the Manuscript Idea: K.S. and S.J.; Methodology and Software: K.S. and M.M.H.; Writing-Original Draft Preparation: K.S.; Writing-Review and Editing: S.J. and M.M.H.; Supervision and Funding Acquisition: S.J.

Funding

This research was funded by the National Natural Science Foundation of China-German Science Foundation (NSFC-DFG) Project, grant number 41761134092.

Acknowledgments

The authors gratefully acknowledged the GFZ for providing products and IGS for the data from MGEX. Many thanks go to Dr. Guoqiang Jiao for his valuable suggestions. We also thank the anonymous reviewers for their detailed and constructive comments that significantly improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jin, S.G.; Jin, R.; Kutoglu, H. Positive and negative ionospheric responses to the March 2015 geomagnetic storm from BDS observations. J. Geod. 2017, 91, 613–626. [Google Scholar] [CrossRef]
  2. Chen, L.; Yi, W.; Song, W.; Shi, C.; Lou, Y.; Cao, C. Evaluation of three ionospheric delay computation methods for ground-based GNSS receivers. GPS Solut. 2018, 22, 125. [Google Scholar] [CrossRef]
  3. 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]
  4. Radicella, S.M. The NeQuick model genesis, uses and evolution. Ann. Geophys. 2009, 52, 417–422. [Google Scholar]
  5. Jakowski, N.; Hoque, M.M.; Mayer, C. A new global TEC model for estimating transionospheric radio wave propagation errors. J. Geod. 2011, 85, 965–974. [Google Scholar] [CrossRef]
  6. Odijk, D. Ionosphere-free phase combinations for modernized GPS. J. Surv. Eng. 2003, 129, 165–173. [Google Scholar] [CrossRef]
  7. Pengfei, C.; Wei, L.; Jinzhong, B.; Hanjiang, W.; Yanhui, C.; Hua, W. Performance of precise point positioning (PPP) based on uncombined dual-frequency GPS observables. Survey Rev. 2011, 43, 343–350. [Google Scholar] [CrossRef]
  8. Lou, Y.; Zheng, F.; Gu, S.; Wang, C.; Guo, H.; Feng, Y. Multi-GNSS precise point positioning with raw single-frequency and dual-frequency measurement models. GPS Solut. 2016, 20, 849–862. [Google Scholar] [CrossRef]
  9. Gu, S.; Shi, C.; Lou, Y.; Liu, J. Ionospheric effects in uncalibrated phase delay estimation and ambiguity-fixed PPP based on raw observable model. J. Geod. 2015, 89, 447–457. [Google Scholar] [CrossRef]
  10. Zhang, H.; Gao, Z.; Ge, M.; Niu, X.; Huang, L.; Tu, R.; Li, X. On the convergence of ionospheric constrained precise point positioning (IC-PPP) based on undifferential uncombined raw GNSS observations. Sensors 2013, 13, 15708–15725. [Google Scholar] [CrossRef]
  11. Cai, C.; Liu, Z.; Luo, X. Single-frequency ionosphere-free precise point positioning using combined GPS and GLONASS observations. J. Navig. 2013, 66, 417–434. [Google Scholar] [CrossRef]
  12. Cai, C.; Gong, Y.; Gao, Y.; Kuang, C. An Approach to Speed up Single-Frequency PPP Convergence with Quad-Constellation GNSS and GIM. Sensors 2017, 17, 1302. [Google Scholar] [CrossRef] [PubMed]
  13. Zhou, F.; Dong, D.; Li, W.; Jiang, X.; Wickert, J.; Schuh, H. GAMP: An open-source software of multi-GNSS precise point positioning using undifferenced and uncombined observations. GPS Solut. 2018, 22, 33. [Google Scholar] [CrossRef]
  14. Schönemann, E.S.; Becker, M.; Springer, T. A new approach for GNSS analysis in a multi-GNSS and multi-signal environment. J. Geod. Sci. 2011, 1, 204–214. [Google Scholar] [CrossRef]
  15. 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]
  16. CSNO. BeiDou Navigation Satellite System Signal in Space Interface Control Document—Open Service Signal B3I. Version 1.0; China Satellite Navigation Office: Beijing, China, 2018. [Google Scholar]
  17. Nava, B.; Coisson, P.; Radicella, S.M. A new version of the NeQuick ionosphere electron density model. J. Atmos. Sol. Terrest. Phys. 2008, 70, 1856–1862. [Google Scholar] [CrossRef]
  18. Hoque, M.M.; Jakowski, N.; Berdermann, J. Ionospheric Correction using NTCM driven by GPS Klobuchar coefficients for GNSS applications. GPS Solut. 2017, 21, 1563–1572. [Google Scholar] [CrossRef]
  19. Otsuka, Y.; Ogawa, T.; Saito, A.; Tsugawa, T.; Fukao, S.; Miyazaki, S. A new technique for mapping of total electron content using GPS network in Japan. Earth Planets Space 2002, 54, 63–70. [Google Scholar] [CrossRef] [Green Version]
  20. Schaer, S.; Gartner, 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]
  21. Witchayangkoon, B. Elements of GPS Precise Point Positioning. Ph.D. Dissertation, Department of Spatial Information Science and Engineering, University of Maine, Orono, ME, USA, 2000. [Google Scholar]
  22. Gao, Z.; Ge, M.; Shen, W.; Zhang, H.; Niu, X. Ionospheric and receiver DCB-constrained multi-GNSS single-frequency PPP integrated with MEMS inertial measurements. J. Geod. 2017, 91, 1351–1366. [Google Scholar] [CrossRef]
  23. Kalman, R. A new approach to linear filtering and prediction problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef]
  24. Su, K.; Jin, S. Improvement of Multi-GNSS Precise Point Positioning Performances with Real Meteorological Data. J. Navig. 2018, 71, 1363–1380. [Google Scholar] [CrossRef]
  25. Kouba, J. A Guide to Using International GNSS Service (IGS) Products; IGS Central Bureau, Jet Propulsion Laboratory: Pasadena, CA, USA, 2009; p. 34. [Google Scholar]
  26. Gérard, P.; Luzum, B. IERS Conventions; No. IERS-TN-36; Bureau International des Poids et Mesures: Sevres, France, 2010. [Google Scholar]
  27. Wu, J.T.; Wu, S.C.; Hajj, G.A.; Bertiger, W.I.; Lichten, S.M. Effects of antenna orientation on GPS carrier phase. Astrodynamics 1992, 18, 1647–1660. [Google Scholar]
  28. Schaffrin, B.; Grafarend, E. Generating classes of equivalent linear models by nuisance parameter. Manuscr. Geod. 1986, 11, 262–271. [Google Scholar]
  29. Xu, G.; Xu, Y. GPS: Theory, Algorithms and Applications, 3nd ed.; Springer: Berlin, Germany, 2016. [Google Scholar]
  30. Jin, S.G.; Jin, R.; Li, D. Assessment of BeiDou differential code bias variations from multi-GNSS network observations. Ann. Geophys. 2016, 34, 259–269. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geographical distribution of the selective MGEX stations.
Figure 1. Geographical distribution of the selective MGEX stations.
Remotesensing 11 00171 g001
Figure 2. F10.7 values and geomagnetic Kp index on 1–30 September 2018.
Figure 2. F10.7 values and geomagnetic Kp index on 1–30 September 2018.
Remotesensing 11 00171 g002
Figure 3. Positioning error scatters of b1 or b2 SPP with different schemes in nine different constellation combinations at station CUT0 (DOY 244/2018, Kp ≈ 0.92) (a) Horizontal error of b1 SPP (b) Vertical error of b1 SPP (c) Horizontal error of b2 SPP (d) Vertical error of b2 SPP. Horizontal error: the horizontal and vertical axes represent the error of the N and E component, respectively (unit: m). Vertical error: the horizontal and vertical axes represent the universal time (unit: hour) and the error of the U component (unit: m), respectively.
Figure 3. Positioning error scatters of b1 or b2 SPP with different schemes in nine different constellation combinations at station CUT0 (DOY 244/2018, Kp ≈ 0.92) (a) Horizontal error of b1 SPP (b) Vertical error of b1 SPP (c) Horizontal error of b2 SPP (d) Vertical error of b2 SPP. Horizontal error: the horizontal and vertical axes represent the error of the N and E component, respectively (unit: m). Vertical error: the horizontal and vertical axes represent the universal time (unit: hour) and the error of the U component (unit: m), respectively.
Remotesensing 11 00171 g003aRemotesensing 11 00171 g003b
Figure 4. Comparison of positioning error of b2 single-frequency PPP with different schemes for GPS, GLONASS and GPS/Beidou/GLONASS/Galileo solutions at station MRO1 (DOY 245/2018, Kp ≈ 1.17). The corresponding satellite numbers and PDOP values are also shown.
Figure 4. Comparison of positioning error of b2 single-frequency PPP with different schemes for GPS, GLONASS and GPS/Beidou/GLONASS/Galileo solutions at station MRO1 (DOY 245/2018, Kp ≈ 1.17). The corresponding satellite numbers and PDOP values are also shown.
Remotesensing 11 00171 g004
Figure 5. Average convergence time for GPS-only, Beidou-only, GLONASS-only, Galileo-only, GPS/Beidou, GPS/GLONASS, GPS/Galileo, GPS/Beidou/GLONASS and GPS/Beidou/GLONASS/Galileo single-frequency PPP collected at ten stations over thirty days.
Figure 5. Average convergence time for GPS-only, Beidou-only, GLONASS-only, Galileo-only, GPS/Beidou, GPS/GLONASS, GPS/Galileo, GPS/Beidou/GLONASS and GPS/Beidou/GLONASS/Galileo single-frequency PPP collected at ten stations over thirty days.
Remotesensing 11 00171 g005
Figure 6. Comparison of positioning error of b1/b2 dual-frequency PPP with different schemes for GPS, GLONASS and GPS/Beidou/GLONASS/Galileo solutions at station SEYG (DOY 246/2018, Kp ≈ 1.92). The corresponding satellite numbers and PDOP values are also shown.
Figure 6. Comparison of positioning error of b1/b2 dual-frequency PPP with different schemes for GPS, GLONASS and GPS/Beidou/GLONASS/Galileo solutions at station SEYG (DOY 246/2018, Kp ≈ 1.92). The corresponding satellite numbers and PDOP values are also shown.
Remotesensing 11 00171 g006
Figure 7. Time series of the estimated receiver DCB at station SEYG (DOY 246/2018).
Figure 7. Time series of the estimated receiver DCB at station SEYG (DOY 246/2018).
Remotesensing 11 00171 g007
Figure 8. Average convergence time for GPS-only, Beidou-only, GLONASS-only, Galileo-only, GPS/Beidou, GPS/GLONASS, GPS/Galileo, GPS/Beidou/GLONASS and GPS/Beidou/GLONASS/Galileo dual-frequency PPP collected at ten stations over thirty days.
Figure 8. Average convergence time for GPS-only, Beidou-only, GLONASS-only, Galileo-only, GPS/Beidou, GPS/GLONASS, GPS/Galileo, GPS/Beidou/GLONASS and GPS/Beidou/GLONASS/Galileo dual-frequency PPP collected at ten stations over thirty days.
Remotesensing 11 00171 g008
Table 1. Data processing strategies for multi-GNSS combined positioning.
Table 1. Data processing strategies for multi-GNSS combined positioning.
ItemsStrategies
Data span1–30 September 2018
Signal SelectionGPS: L1/L2/L5; BDS: B1/B2/B3; GLONASS: G1/G2; Galileo: E1/E5a/E5b
EstimatorSPP: LS PPP: Kalman filter [23]
Observation sampling rate30 s
Elevation cutoff
Satellite orbit and clockSPP: broadcast ephemeris provided by MGEX
PPP: Fixed by MGEX orbit and clock products
Satellite TGDCorrect using broadcast ephemeris for SPP
Satellite DCBCorrect using MGEX DCB products for PPP
Tropospheric delaySPP: Modified Hopfield PPP: Modified Hopfield for dry part and estimated for wet part as random-walk process (10−9 m2/s) [24]
Ionospheric delayEliminated or estimated as white noise Strategies in Table 2
Relativistic effectGeneral relativistic models [25]
Sagnac effectCorrected model [26]
Phase windup effectCorrected model [27]
Satellite and receiver antennaCorrected with the values from MGEX and IGS [26]
Tide displacementCorrected, including solid Earth, pole and ocean tide [26]
Station reference coordinatesIGS SINEX solutions
Station coordinateEstimated as constants
Receiver clockEstimated as white noise process (1×105 m2)
Receiver IFBEstimated as constants or absorbed by receiver clock biases
Phase ambiguitiesEstimated as constants for float values
Table 2. Summaries of the ionosphere processing strategy in multi-GNSS positioning.
Table 2. Summaries of the ionosphere processing strategy in multi-GNSS positioning.
ConstellationModeSchemesComment
G, C, R, ESingle-frequency SPP: b1, b2, b3BIM-corr; NTCM-corr; GIM-corr;In the BIM, different ionospheric delay models are utilized for their corresponding systems.
Klobuchar model (G/R);
improved Klobuchar model(C);
NeQuick model(E)
G/C
G/RSingle-frequency PPP: b1, b2, b3GRAPHIC; UC; BIM-cons; NTCM-cons; GIM-cons;
G/E
G/C/RDual-frequency PPP: b1-b2;b1-b3IF; UC; BIM-cons; NTCM-cons; GIM-cons;
G/C/R/E
The abbreviation G, C, R and E represent GPS, Beidou, GLONASS and Galileo, respectively. The abbreviation “corr” and “cons” represent “correction” and “constrained”, respectively.
Table 3. RMS of single-frequency SPP with BIM, NTCM and GIM based ionospheric correction (unit: m).
Table 3. RMS of single-frequency SPP with BIM, NTCM and GIM based ionospheric correction (unit: m).
TypeSchemesb1b2b3
NEUNEUNEU
GBIM2.0841.6804.0922.4911.8065.338---
GNTCM2.1521.6343.6732.6341.7504.135---
GGIM1.7041.5703.5891.7681.5743.649---
CBIM4.4556.5168.1954.6856.1668.2264.0966.2037.364
CNTCM4.3486.1127.6214.6835.8197.4234.2195.8107.055
CGIM4.0526.0937.7174.0374.8017.4753.6815.7216.990
RBIM2.4102.1685.4222.6072.2746.382---
RNTCM2.3562.0615.0402.6182.1565.593---
RGIM2.1692.0224.9802.2042.0805.401---
EBIM3.1172.8025.1493.4662.9356.3323.5063.0336.112
ENTCM3.1902.7694.7543.8043.0265.5983.7673.0415.443
EGIM2.7732.7044.6042.9362.9105.1663.0253.0112.113
GCBIM2.0731.6674.0772.4801.7935.3242.7183.2916.770
GCNTCM2.1411.6253.6602.6241.7394.1222.8233.1035.857
GCGIM1.6901.5593.5781.7581.5623.6362.1362.9975.413
GRBIM1.7491.3963.4192.0671.4954.428---
GRNTCM1.8031.3623.1082.1831.4573.498---
GRGIM1.4531.3123.0061.5071.3183.092---
GEBIM1.6971.3613.3572.0521.4734.3383.0812.4335.706
GENTCM1.7581.3262.9922.1741.4293.4052.9022.1875.579
GEGIM1.3681.2702.9121.4331.2772.9872.1291.9154.781
GCRBIM1.7481.3943.4142.0661.4934.425---
GCRNTCM1.8021.3613.1062.1831.4553.497---
GCRGIM1.4511.3103.0051.5061.3173.091---
GCREBIM1.7281.3763.3812.0621.4834.429---
GCRENTCM1.7861.3443.0612.1811.4453.472---
GCREGIM1.4231.2902.9561.4841.3013.055---
“-” means no corresponding combination of results, the same below.
Table 4. Parameter estimation method in single-frequency quad-constellation PPP schemes.
Table 4. Parameter estimation method in single-frequency quad-constellation PPP schemes.
GRAPHICUCIonosphere Constraint
Observed quantitym+n+p+q2m+2n+2p+2q3m+3n+3p+3q
Parameter numberm+n+p+q+82m+2n+2p+2q+82m+2n+2p+2q+8
Redundancy−8−8m+n+p+q−8
Type of the parametersPosition, receiver clock, ZWD, ambiguityPosition, receiver clock, ZWD, slant ionospheric delay, ambiguityPosition, receiver clock, ZWD, slant ionospheric delay, ambiguity
Table 5. RMS of single-frequency PPP with IF approach, UC approach, BIM-, NTCM- and GIM-constrained approaches (unit: cm).
Table 5. RMS of single-frequency PPP with IF approach, UC approach, BIM-, NTCM- and GIM-constrained approaches (unit: cm).
TypeSchemesb1b2b3
NEUNEUNEU
GGRAPHIC1.6643.4199.2671.8183.4686.899---
GUC1.6643.4199.2661.8183.4686.900---
GBIM-cons1.6333.9909.2721.7753.9246.906---
GNTCM-cons1.6673.1669.2901.8263.2556.867---
GGIM-cons1.6422.8959.2861.8243.1966.825---
CGRAPHIC7.84410.15326.7827.2239.82625.3747.0379.87724.515
CUC7.85910.16426.7557.2179.82625.3517.0409.88724.494
CBIM-cons6.3188.63525.6455.6727.27124.6066.5408.03623.508
CNTCM-cons6.9898.04524.2886.0656.70024.3026.7997.80122.575
CGIM-cons6.9047.74322.0115.9626.46922.3146.6447.65821.717
RGRAPHIC4.4048.00320.4404.3918.58217.139---
RUC4.4038.00320.5034.3848.55917.142---
RBIM-cons4.3478.09519.2674.2688.44216.998---
RNTCM-cons4.3568.09018.7134.2798.44916.873---
RGIM-cons4.2487.94118.3133.9988.37316.668---
EGRAPHIC2.7173.7517.8632.4863.9018.0482.9394.8058.760
EUC2.7173.7537.8622.4883.9018.0532.9394.8058.759
EBIM-cons2.6583.7777.8322.3693.9327.8722.8324.5928.681
ENTCM-cons2.7583.6687.7862.5303.7427.0042.9174.4968.190
EGIM-cons2.5213.5327.7502.3153.2366.9062.6664.4247.906
GCGRAPHIC1.6523.3879.1681.8083.4306.7853.6998.94213.284
GCUC1.6523.3879.1651.8083.4296.7873.6998.94113.282
GCBIM-cons1.6283.9089.2421.7463.9076.8393.3047.86012.912
GCNTCM-cons1.6593.1629.2451.7963.1666.7593.5847.43012.257
GCGIM-cons1.6323.0669.2081.7933.1416.7403.1066.68112.244
GRGRAPHIC1.6623.4189.2601.8103.4616.890---
GRUC1.6623.4179.2591.8103.4616.890---
GRBIM-cons1.6183.9829.2461.7653.9196.853---
GRNTCM-cons1.6623.1599.2681.8193.1726.795---
GRGIM-cons1.6312.8769.2161.8213.1696.779---
GEGRAPHIC1.5682.9058.6571.6063.1236.8792.9034.4898.089
GEUC1.5682.9058.6561.6063.1226.8802.9034.4898.089
GEBIM-cons1.5423.2428.6551.5683.6756.8902.8064.3908.074
GENTCM-cons1.5853.0078.6341.6202.8996.8322.8644.1647.837
GEGIM-cons1.5432.4808.6201.6182.8616.8192.8293.2357.715
GCRGRAPHIC1.6513.3859.1611.7113.3426.775---
GCRUC1.6513.3859.1591.7103.3416.777---
GCRBIM-cons1.5503.0989.0131.6643.8406.762---
GCRNTCM-cons1.5912.2538.9321.7343.1456.688---
GCRGIM-cons1.5502.0478.8441.6993.1296.597---
GCREGRAPHIC1.5572.8858.5791.5993.0956.670---
GCREUC1.5572.8858.5761.5993.0946.673---
GCREBIM-cons1.5383.0358.4881.5623.5076.669---
GCRENTCM-cons1.5782.2408.4081.6102.8026.575---
GCREGIM-cons1.5411.9628.3161.5992.7536.517---
Table 6. Parameter estimation method in dual-frequency quad-constellation PPP schemes.
Table 6. Parameter estimation method in dual-frequency quad-constellation PPP schemes.
IFUCIonosphere Constraint
Observed quantity2m+2n+2p+2q4m+4n+4p+4q5m+5n+5p+5q
Parameter numberm+n+p+q+83m+3n+3p+3q+83m+3n+3p+3q+12
Redundancym+n+p+q−8m+n+p+q−82m+2n+2p+2q−12
Type of the parametersPosition, receiver clock, ZWD, ambiguityPosition, receiver clock, ZWD, slant ionospheric delay, ambiguityPosition, receiver clock, ZWD, DCB, slant ionospheric delay, ambiguity
Table 7. RMS of dual-frequency PPP with IF approach, UC approach, and BIM-, NTCM- and GIM-constrained approaches (unit: cm).
Table 7. RMS of dual-frequency PPP with IF approach, UC approach, and BIM-, NTCM- and GIM-constrained approaches (unit: cm).
TypeSchemesb1/b2b1/b3
NEUNEU
GIF0.9051.6133.147---
GUC0.9071.6093.135---
GBIM-cons0.8781.5953.162---
GNTCM-cons0.9051.5873.159---
GGIM-cons0.8981.5973.148---
CIF1.4352.4905.6791.5602.9545.886
CUC1.4402.5265.7261.5722.9565.968
CBIM-cons1.3602.3745.5351.3902.4645.828
CNTCM-cons1.3872.3095.4881.3982.4445.692
CGIM-cons1.3772.2695.4661.3402.4385.507
RIF0.9171.8733.303---
RUC0.9171.8733.302---
RBIM-cons0.9441.8753.472---
RNTCM-cons0.9441.8893.434---
RGIM-cons0.9241.8783.428---
EIF1.0281.8093.5661.1092.0993.364
EUC1.0251.8133.5601.0802.0893.357
EBIM-cons1.0401.8053.4861.0661.9183.306
ENTCM-cons1.0301.7713.4041.0971.8473.318
EGIM-cons1.0391.7423.1111.0381.7993.308
GCIF0.8861.7103.0661.4272.7895.405
GCUC0.8921.7083.0561.4252.7745.393
GCBIM-cons0.8701.5812.9461.3302.3675.426
GCNTCM-cons0.8891.5632.9571.3552.3535.238
GCGIM-cons0.8271.5292.9221.2932.3075.007
GRIF0.8111.5252.903---
GRUC0.8071.5232.907---
GRBIM-cons0.7901.4962.874---
GRNTCM-cons0.8171.4872.757---
GRGIM-cons0.7711.3722.718---
GEIF0.8711.7763.1071.0532.0303.362
GEUC0.8641.7783.0921.0541.9973.355
GEBIM-cons0.8531.6072.9201.0711.9143.445
GENTCM-cons0.8541.6042.8841.0771.8263.275
GEGIM-cons0.8341.5872.8680.9941.7713.247
GCRIF0.7941.5132.886---
GCRUC0.7871.5112.884---
GCRBIM-cons0.7821.4282.818---
GCRNTCM-cons0.7951.4082.742---
GCRGIM-cons0.7551.3052.701---
GCREIF0.7771.4262.702---
GCREUC0.7601.3932.695---
GCREBIM-cons0.7741.3452.476---
GCRENTCM-cons0.7781.3242.366---
GCREGIM-cons0.6761.2862.347---

Share and Cite

MDPI and ACS Style

Su, K.; Jin, S.; Hoque, M.M. Evaluation of Ionospheric Delay Effects on Multi-GNSS Positioning Performance. Remote Sens. 2019, 11, 171. https://doi.org/10.3390/rs11020171

AMA Style

Su K, Jin S, Hoque MM. Evaluation of Ionospheric Delay Effects on Multi-GNSS Positioning Performance. Remote Sensing. 2019; 11(2):171. https://doi.org/10.3390/rs11020171

Chicago/Turabian Style

Su, Ke, Shuanggen Jin, and M. M. Hoque. 2019. "Evaluation of Ionospheric Delay Effects on Multi-GNSS Positioning Performance" Remote Sensing 11, no. 2: 171. https://doi.org/10.3390/rs11020171

APA Style

Su, K., Jin, S., & Hoque, M. M. (2019). Evaluation of Ionospheric Delay Effects on Multi-GNSS Positioning Performance. Remote Sensing, 11(2), 171. https://doi.org/10.3390/rs11020171

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