Next Article in Journal
Multitemporal Analysis of Slow-Moving Landslides and Channel Dynamics through Integrated Remote Sensing and In Situ Techniques
Next Article in Special Issue
Tropospheric Delay Parameter Estimation Strategy in BDS Precise Point Positioning
Previous Article in Journal
The SOLAR-HRS New High-Resolution Solar Spectra for Disk-Integrated, Disk-Center, and Intermediate Cases
Previous Article in Special Issue
Region-Specific and Weather-Dependent Characteristics of the Relation between GNSS-Weighted Mean Temperature and Surface Temperature over China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Epoch Decimeter-Level Precise Point Positioning with a Galileo Five-Frequency Ionosphere-Reduced Combination

1
School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China
2
Key Laboratory of Micro-Inertial Instrument and Advanced Navigation Technology, Southeast University, Nanjing 210096, China
3
Laboratory of Science and Technology on Marine Navigation and Control, China State Shipbuilding Corporation, Tianjin 300131, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(14), 3562; https://0-doi-org.brum.beds.ac.uk/10.3390/rs15143562
Submission received: 20 June 2023 / Revised: 12 July 2023 / Accepted: 14 July 2023 / Published: 15 July 2023
(This article belongs to the Special Issue New Progress in GNSS Data Processing Technology and Modeling)

Abstract

:
Currently, there are two main methods for single-epoch decimeter-level precise point positioning (PPP); one is a model based on ambiguity-fixed ionosphere-free (AFIF) observations, and the other is based on uncombined (UC) PPP. The implementation of these two models requires both extra-wide-lane (EWL) and wide-lane (WL) ambiguity fixing. Different from the existing methods, this paper proposes a multi-frequency ionosphere-reduced (IR) PPP model suitable for single-epoch decimeter-level positioning. Based on Galileo five-frequency data, the optimal selection strategy of IR combinations is first studied with ionosphere, noise level and wavelength factors considered. Then, based on the selected IR combination, two IR PPP models, namely IR(EST) and IR(IGN), are established according to whether ionosphere parameters are estimated or ignored. Finally, the proposed models are verified with real tracked data from globally distributed stations, and further compared with the existing AFIF/UC models in terms of positioning performance and time consumption. The relationship between the ionosphere equivalent ranging error and satellite elevation in the IR models is analyzed. The lower the elevation is, the larger the residual ranging error is, and its impact on positioning is weakened by downweighting its observations and adjusting the cut-off elevation during the partial ambiguity fixing (PAF) process. The results show that the performance of the two IR models is basically the same, and both can achieve horizontal and vertical accuracies better than 20 cm and 40 cm, respectively. Compared with the existing AFIF/UC models, the proposed IR models can achieve similar decimeter-level accuracy with only one step of EWL ambiguity fixing, and at the same time, the IR models have varying degrees of improvement in time consumption: 38% shorter than the AFIF model and 97% shorter than the UC model.

1. Introduction

The global navigation satellite system (GNSS) precise point positioning (PPP) technology, known for its ability to provide positioning, navigation and timing (PNT) services with high accuracy on a global scale, has been widely used in both military and civilian fields [1,2]. However, traditional float PPP without ambiguity fixing often requires a rather long convergence time to achieve centimeter-level accuracy [3,4], and this also limits the development of PPP technology, which cannot be widely used like real-time kinematic (RTK). By performing ambiguity resolution, the convergence time can be slightly shortened, and the positioning accuracy can also be improved to a certain extent [5,6]. Especially with the development of multi-GNSS, the initialization time is significantly shortened, but even so, the time to first fix (TTFF) is still about 10 min. Improving its real-time performance has always been the pursuit of scholars in the field. Currently, there are two commonly used strategies to improve PPP real-time performance. One is to adopt the similar idea of network real-time kinematic (NRTK) to extract and model the regional undifferenced atmosphere and then achieve fast (or even single-epoch) centimeter-level positioning by means of atmosphere enhancement. This technology is also called PPP-RTK by some scholars [7,8]. Another strategy is to achieve high real-time performance at the cost of degrading positioning accuracy to a certain extent (e.g., from centimeter level to decimeter level) without external atmosphere constraints, and the typical representative is the single-epoch wide-lane ambiguity resolution (WAR) positioning with decimeter-level accuracy [9,10], which takes advantage of the fact that the extra-wide-lane (EWL) and wide-lane (WL) ambiguity could be fixed reliably in single-epoch mode. In general, these two modes have their own advantages and applicable scenarios, and the content of this paper mainly focuses on the latter.
The modernization of GNSS can greatly improve the efficiency and reliability of ambiguity fixing, and especially in the case of good satellite geometry, an instantaneous accuracy of about 20 cm can be achieved [11]. With the fixed EWL and WL ambiguity, the so-called ambiguity-fixed ionosphere-free (AFIF) observation is reconstructed, since its noise and multipath effect outperform the raw pseudorange; this observation can be used to perform instantaneously decimeter-level positioning or to assist narrow-lane (NL) ambiguity fixing [12]. Based on Galileo/BDS multi-frequency data, Li et al. [13,14] studied the performance of various AFIF combinations composed of triple-/quad-/five-frequency, and the results show that decimeter-level accuracy can be obtained within 0.5 min with the instantaneous fixing of EWL/WL ambiguity, and the positioning accuracy of the five-frequency combination is better than that of the triple-/quad-frequency combinations. With Galileo multi-frequency data, Guo and Xin [15] found that compared with other combinations, the E1/E5a/E6 combination has a lower noise level, and the resulting AFIF observation is theoretically expected to improve the single-epoch accuracy to about 10 cm. In addition to the above-mentioned AFIF-based model, some scholars have also achieved similar positioning accuracy based on uncombined (UC) PPP (UCPPP) with EWL/WL/NL ambiguities sequentially fixed [16]. The UCPPP model can make full use of pseudoranges and carrier observations on all frequencies and is an ideal model for processing multi-frequency data. A major advantage of using the UCPPP model is that EWL/WL ambiguities can be estimated by a geometry-based (GB) model, so the ambiguity fixing is theoretically more reliable in single-epoch mode compared with a conventional geometry-free (GF) model, such as the Melbourne–Wübbena (MW) combination [17,18]. The results of Gu et al.’s research [19] based on BDS-2 triple-frequency data show that the EWL and WL ambiguity can be fixed reliably within 2 min, and a horizontal and vertical accuracy of about 0.5 m can be obtained. Gao et al. [20] also analyzed the performance of single-epoch ambiguity-fixed PPP with regional atmosphere enhancement, and the WL ambiguity-fixed solution can achieve an accuracy of several centimeters. Qu et al. [21] also obtained similar positioning results and found that compared with triple-frequency WAR, the performance improvement of quad-/five-frequency cases is not obvious. Generally speaking, UCPPP-based single-epoch positioning can achieve an accuracy of about 20 cm in the horizontal direction on a global scale [9,10].
Compared with a dual-frequency case, multi-frequency can provide more combined observations with excellent characteristics, such as those demonstrated by a ionosphere-reduced (IR) combination [22]. Usually in the ionosphere-free (IF) combination case, in order to eliminate the influence of ionosphere, the combined ambiguity lose its integer nature; in contrast, the IR combination based on integer coefficients can still guarantee the integer solvability of ambiguity under the condition of almost eliminating ionosphere, which provides new opportunities for GNSS data processing. Both Feng [22] and Guo et al. [23] systematically studied the selection criteria of IR combinations. Based on BDS-3 multi-frequency signals, Li et al. [24,25] studied an IR combination suitable for medium-/long-baseline RTK ambiguity resolution, whose ambiguity maintains integer characteristics, and the noise level is comparable to the NL combination. At present, IR-based research mainly focus on differential positioning, such as RTK, and there is little targeted research on PPP.
Based on the above background, this paper mainly studies the theory of a multi-frequency IR combination in single-epoch PPP positioning and further compares it with the existing AFIF and UC models. The article is structured as follows. Section 2 mainly introduces and compares different positioning models, namely the AFIF model, UC model and IR model. Section 3 presents the experiment’s results and analysis. A summary of our results is given in Section 4 and Section 5.

2. Materials and Methods

The zero-differenced observations of pseudorange P r , j s and carrier phase L r , j s between the receiver r and satellite s on frequency j can be expressed as:
P r , j s = ρ r s + c · t r c · t s + M r , w s · T r , w + γ j · I r , 1 s + d r , j d j s + e r , j s
L r , j s = ρ r s + c · t r c · t s + M r , w s · T r , w γ j · I r , 1 s + λ j · ( N r , j s + b r , j b j s ) + ε r , j s
where ρ r s denotes the geometric distance between the satellite and receiver; c is the speed of light; t r and t s denote the receiver and satellite clock offsets; T r , w is the residual zenith wet-tropospheric delay (ZWD) with mapping function M r , w s ; I r , 1 s is the slant ionosphere delay at the first frequency f 1 , with an amplification factor γ j = f 1 2 / f j 2 depending on the frequency f j ; d r , j and d j s denote the receiver and satellite code biases; b r , j and b j s denote the receiver and satellite phase biases; λ j and N r , j s are the wavelength and integer ambiguity of frequency f j , respectively; and e r , j s and ε r , j s denote observation noises. It should be noted that some error items that could be accurately corrected by the existing models are not listed, such as the phase center offset/variation (PCO/PCV), earth tides loading, phase wind-up, etc. [26,27].
Since the zero-differenced PPP data processing cannot eliminate the receiver and satellite hardware bias through differential means used in relative positioning, it is not possible to estimate each parameter accurately, but only the combination of multiple parameters, i.e., re-parameterization. This strategy is widely adopted for each analysis center or client end. In order to simplify the description of subsequent sections, some re-parameterization strategies commonly used in current PPP data processing are explained first.
At present, dual-frequency IF observations are usually adopted for satellite precise orbit determination (POD) and precise clock estimation (PCE); some frequency-related hardware biases are correspondingly absorbed by satellite clock offsets.
c · t ˜ s = c · t s + α · d 1 s + β · d 2 s
where t ˜ s denotes the satellite clock offset provided by each analysis center and α and β denote the coefficients of a dual-frequency IF combination, with specific forms as:
{ α = f 1 2 / ( f 1 2 f 2 2 ) β = f 2 2 / ( f 1 2 f 2 2 )
Similarly, the receiver code bias in the form of the IF combination is also lumped with the re-parameterized receiver clock offset t ˜ r .
c · t ˜ r = c · t r + α · d r , 1 + β · d r , 2
For an uncombined PPP model with estimated ionosphere parameters, the frequency dependent code biases are merged with ionosphere parameters.
I ˜ r , 1 s = I r , 1 s β · ( d r , 2 d r , 1 d 2 s + d 1 s )
The re-parameterization process described in the above Equations (3), (5) and (6) has been widely used in current dual-/multi-frequency PPP data processing, and the subsequent function models are also based on these re-parameterizations. Based on Galileo five-frequency data, as shown in Table 1, the three positioning models will be described below, namely the AFIF model, UC model and IR model. For the convenience of description, the following sections use f i ( i = 1 , 2 , 3 , 4 , 5 ) to represent the E1/E5a/E6/E5b/E5 frequencies of Galileo, respectively.

2.1. AFIF Model

The AFIF observation is composed of an IF combination of specific EWL and WL observations. Taking E6/E5a and E1/E5a as the selected EWL and WL combinations, the corresponding observation equations are,
L r , E 23 s = f 3 f 3 f 2 · L r , 3 s f 2 f 3 f 2 · L r , 2 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w f 1 2 f 2 · f 3 · I ˜ r , 1 s + λ E 23 · ( N r , E 23 s + B r , E 23 B E 23 s ) + ε r , E 23 s
L r , W 12 s = f 1 f 1 f 2 · L r , 1 s f 2 f 1 f 2 · L r , 2 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w f 1 f 2 · I ˜ r , 1 s + λ W 12 · ( N r , W 12 s + B r , W 12 B W 12 s ) + ε r , W 12 s
where L r , E 2 3 s denotes the EWL observation composed of frequency f 2 and f 3 ; L r , W 12 s denotes the WL observation composed of frequency f 1 and f 2 ; λ E 2 3 and λ W 12 are EWL and WL wavelengths, with λ E 2 3 = 2.93   m and λ W 12 = 0.75   m ; N r , E 23 s and N r , W 12 s are EWL and WL integer ambiguities; B r , E 2 3 and B E 23 s are receiver and satellite EWL fractional cycle biases (FCBs); B r , W 12 and B W 12 s are receiver and satellite WL FCBs; and ε r , E 23 s and ε r , W 12 s are observation noises of EWL/WL observations.
Due to the long wavelengths of EWL combinations, typically several meters, the EWL ambiguity can be reliably fixed in single-epoch mode. Specifically, the float EWL ambiguity N ¯ r , E 23 s is usually generated by using the well-known MW combination with both geometry and ionosphere terms eliminated. After receiver and satellite FCBs corrections are applied, it can be fixed by integer rounding with a typical threshold of 0.2~0.3 cycles. It should be noted that in order to be consistent with the PCE estimation strategy and avoid the introduction of additional code bias by multi-frequency pseudoranges, when constructing MW combinations, only E1/E5a dual-frequency pseudoranges consistent with the PCE reference are used for various EWL combinations.
N ¯ r , E 23 s = L r , E 23 s [ η · P r , 1 s + ( 1 η ) · P r , 2 s ] λ E 23 = N r , E 23 s + B r , E 23 B E 23 s
with η = f 1 2 · ( f 3 f 2 ) f 3 · ( f 1 2 f 2 2 ) .
Once the EWL ambiguity is fixed, the unambiguous EWL observation L ¯ r , E 23 s can be obtained.
L ¯ r , E 23 s = L r , E 23 s ( N r , E 23 s + B r , E 23 B E 23 s )
where N r , E 23 s denotes fixed integer EWL ambiguity.
Then, this unambiguous EWL observation is combined with the WL observation in Equation (8) to form the AFIF observation.
L r , AFIF 123 s = f 1 f 1 f 3 · L r , W 12 s f 3 f 1 f 3 · L ¯ r , E 23 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + f 1 f 1 f 3 · λ W 12 · ( N r , W 12 s + B r , W 12 B W 12 s ) + ε r , AFIF 123 s
where L r , AFIF 123 s denotes the AFIF observation and ε r , AFIF 123 s denotes corresponding observation noise.
Assume that the accuracy of carrier observations at different frequencies is equal, according to the law of error propagation, the prior accuracy of AFIF observations L r , AFIF 123 s can be obtained as:
σ AFIF 123 2 = κ AFIF 123 2 · σ L 2 = ( ( f 1 f 1 f 3 ) 2 · f 1 2 + f 2 2 ( f 1 f 2 ) 2 + ( f 3 f 1 f 3 ) 2 · f 3 2 + f 2 2 ( f 3 f 2 ) 2 ) · σ L 2
where κ AFIF 123 denotes the noise amplification factor and σ L and σ AFIF 123 denote the accuracy of raw carrier phase and AFIF observations, respectively. The σ L is usually set to 3 mm.
Based on Equation (11), if the WL ambiguity is successfully fixed, the coordinate parameters can be constrained to about decimeter-level accuracy [12,14]. For WL ambiguity resolution, there are generally two approaches. One is to use the MW combination similarly to EWL ambiguity fixing, that is, the GF method; the other one is to calculate the float WL ambiguity and its variance–covariance based on the GB model and then use the Least-squares AMBiguity Decorrelation Adjustment (LAMBDA) algorithm to fix it [28].
If the GF method is adopted for WL ambiguity resolution, the AFIF observation becomes unambiguous, which means that it can be treated as a virtual pseudorange.
L ¯ r , AFIF 123 s = L r , AFIF 123 s f 1 f 1 f 3 · λ W 12 · ( N r , W 12 s + B r , W 12 B W 12 s )
where L ¯ r , AFIF 123 s denotes the unambiguous AFIF observation, and N r , W 12 s denotes fixed integer WL ambiguity.
In case of sufficient visible satellites, theoretically, absolute positioning can be achieved by only using these unambiguous AFIF observations without the assistance of pseudoranges. However, in order to ensure the stability of the observation model, it can still be combined with pseudorange observations for parameter estimation. On this basis, further combining various AFIF observations composed of other frequencies, such as E1/E5a/E5b and E1/E5a/E5, the following observation model can be obtained:
{ P r , IF s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + e r , IF s L ¯ r , AFIF 123 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + ε r , AFIF 123 s L ¯ r , AFIF 124 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + ε r , AFIF 124 s L ¯ r , AFIF 125 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + ε r , AFIF 125 s
where P r , IF s = α · P r , 1 s + β · P r , 2 s denotes dual-frequency IF pseudoranges. The estimated parameter vector in Equation (14) includes:
X ^ = [ x , y , z t ˜ r T r , w ]
where ( x , y , z ) denotes receiver coordinates.
If the GB model is adopted, due to the existence of ambiguity parameters, pseudoranges must be introduced to ensure that the model is full rank:
{ P r , IF s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + e r , IF s L r , AFIF 123 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + f 1 f 1 f 3 · λ W 12 · N ¯ r , W 12 s + ε r , AFIF 123 s L r , AFIF 124 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + f 1 f 1 f 4 · λ W 12 · N ¯ r , W 12 s + ε r , AFIF 124 s L r , AFIF 125 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + f 1 f 1 f 5 · λ W 12 · N ¯ r , W 12 s + ε r , AFIF 125 s
where N ¯ r , W 12 s denotes float WL ambiguity and can be expressed as:
N ¯ r , W 12 s = N r , W 12 s + B r , W 12 B W 12 s
In case of n visible satellites, the estimated vector corresponding to Equation (16) includes:
X ^ = [ x , y , z t ˜ r T r , w N ¯ r , W 12 1 , N ¯ r , W 12 2 , , N ¯ r , W 12 n n ]
It is not difficult to find from Equation (16) that although three different AFIF combinations are used, only one common WL ambiguity parameter needs to be estimated. For different AFIF combinations, the corresponding WL ambiguity wavelengths are amplified to varying degrees, with magnifications of f 1 f 1 f 3 , f 1 f 1 f 4 and f 1 f 1 f 5 , respectively [12]. As shown in Table 2, the equivalent wavelength can reach about 3~4 m, which is obviously larger than the original 0.751 m. The larger equivalent wavelength is also more conducive to efficient WL ambiguity resolution.
For convenience, AFIF(GF) and AFIF(GB) are used in subsequent sections to represent the models in Equation (14) and Equation (16), respectively.

2.2. Uncombined Model

The uncombined PPP model has been demonstrated to be the most ideal model for multi-frequency data processing. To avoid additional code bias correction for multi-frequency pseudoranges, only E1/E5a dual-frequency pseudoranges are used for parameter estimation, with the observation model as:
{ P r , 1 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + I ˜ r , 1 s + e r , 1 s P r , 2 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + γ 2 · I ˜ r , 1 s + e r , 2 s L r , 1 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w I ˜ r , 1 s + λ 1 · N ¯ r , 1 s + ε r , 1 s L r , 2 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w γ 2 · I ˜ r , 1 s + λ 2 · N ¯ r , 2 s + ε r , 2 s L r , j s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w γ j · I ˜ r , 1 s + λ j · N ¯ r , j s + ε r , j s
where N ¯ r , j s denotes float ambiguity at frequency j , which can be expressed as:
N ¯ r , j s = N r , j s + B r , j B j s
where B r , j and B j s denote receiver and satellite FCBs at frequency j .
In the case of n visible satellites with j frequencies, the estimated vector contains the receiver position, receiver clock offset, ZWD, slant ionosphere delay and multi-frequency ambiguities for each satellite.
X ^ = [ x , y , z t ˜ r T r , w I ˜ r , 1 1 , I ˜ r , 1 2 , , I ˜ r , 1 n n N ¯ r , 1 , N ¯ r , 2 , , N ¯ r , n n × j ]
with N ¯ r , s = [ N ¯ r , 1 s , N ¯ r , 2 s , , N ¯ r , j s j ] .
After the float solution X ^ and its variance–covariance Q X ^ X ^ are obtained, the raw ambiguity of each frequency can be mapped to EWL/WL ambiguity with a properly selected transformation matrix. In this paper, the same EWL/WL combinations as shown in Table 2 are adopted, that is, E5–E5a, E5b–E5a, E6–E5a and E1–E5a, with corresponding wavelengths of 19.54 m, 9.77 m, 2.93 m and 0.75 m. According to the error propagation law, the transformed vector X ^ LC and its variance–covariance Q X ^ LC X ^ LC can be derived as:
X ^ LC = [ I D W D E 1 D E 2 D E 3 ] · X ^ ,       Q X ^ LC X ^ LC = [ I D W D E 1 D E 2 D E 3 ] · Q X ^ X ^ · [ I D W D E 1 D E 2 D E 3 ] T
where I denotes the identity matrix of dimension ( m + 5 ) ; D W denotes the E1–E5a WL ambiguity mapping matrix; and D E 1 , D E 2 and D E 3 denote the E6–E5a, E5b–E5a and E5–E5a EWL ambiguity mapping matrix, respectively.
To simplify illustration, the first satellite is chosen as the reference satellite to eliminate the effect of receiver FCBs, and the transformed vector X ^ LC contains:
X ^ LC = [ x , y , z t ˜ r T r , w I ˜ r , 1 1 , I ˜ r , 1 2 , , I ˜ r , 1 n n N ¯ r , W 1 , N ¯ r , E 1 1 , N ¯ r , E 2 1 , N ¯ r , E 3 1 ( n 1 ) × 4 ]
where N ¯ r , W 1 , N ¯ r , E 1 1 , N ¯ r , E 2 1 , N ¯ r , E 3 1 denote the between-satellite single-differenced (BSSD) E1–E5a, E6–E5a, E5b–E5a and E5–E5a ambiguity of each satellite, and their specific expressions are as follows:
{ N ¯ r , W 1 = [ N ¯ r , W 12 , N ¯ r , W 13 , N ¯ r , W 1 n ] N ¯ r , E 1 1 = [ N ¯ r , E 1 12 , N ¯ r , E 1 13 , , N ¯ r , E 1 1 n ] N ¯ r , E 2 1 = [ N ¯ r , E 2 12 , N ¯ r , E 2 13 , , N ¯ r , E 2 1 n ] N ¯ r , E 3 1 = [ N ¯ r , E 3 12 , N ¯ r , E 3 13 , , N ¯ r , E 3 1 n ]
The transformation matrix D in Equation (22) can be generically expressed as:
D = [ C C C C C C ]
where the symbol represents W, E1, E2 and E3; C denotes a specific combination coefficient; and the expression is as follows:
{ C W = [ 1 1 0 0 0 ] C E 1 = [ 0 1 1 0 0 ] C E 2 = [ 0 1 0 1 0 ] C E 3 = [ 0 1 0 0 1 ]
Based on the transformed vector X ^ LC and its variance–covariance Q X ^ LC X ^ LC , after satellite FCBs are corrected, a stepwise ambiguity resolution process is carried out with EWL/WL ambiguities fixed sequentially. It should be noted that the three selected EWL ambiguities are all fixed in one step, instead of being divided into three steps [16]. Regardless of whether the EWL or WL ambiguity is fixed, the precision of the position parameter can be improved theoretically. When both EWL and WL ambiguities are successfully fixed, decimeter-level positioning accuracy can generally be achieved. Subsequently, the term UC is used to represent the model in Equation (19).

2.3. Ionosphere-Reduced Model

2.3.1. Multi-Frequency Observation Combination

For the E1/E5a/E6/E5b/E5 five-frequency case, the linear combination of observation based on integer coefficients can be expressed as [22,29]:
L r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s = k = 1 5 ( i k · f k · L r , k s ) k = 1 5 ( i k · f k )
with the specific observation equation as
L r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · I ˜ r , 1 s + λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s + μ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · ε r , 1 s
where γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) , λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) , N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s and μ ( i 1 , i 2 , i 3 , i 4 , i 5 ) denote the ionosphere scalar factor, wavelength, float ambiguity and noise amplification factor of the combined signal, respectively. Their specific forms are as follows:
{ γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) = f 1 2 · ( i 1 / f 1 + i 2 / f 2 + + i 5 / f 5 ) f ( i 1 , i 2 , i 3 , i 4 , i 5 ) λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) = c f ( i 1 , i 2 , i 3 , i 4 , i 5 ) N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s = N r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s + B r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) B ( i 1 , i 2 , i 3 , i 4 , i 5 ) s μ ( i 1 , i 2 , i 3 , i 4 , i 5 ) = ( i 1 · f 1 ) 2 + ( i 2 · f 2 ) 2 + + ( i 5 · f 5 ) 2 f ( i 1 , i 2 , i 3 , i 4 , i 5 )
where B r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) and B ( i 1 , i 2 , i 3 , i 4 , i 5 ) s denote receiver and satellites FCBs of the combined ambiguity, and f ( i 1 , i 2 , i 3 , i 4 , i 5 ) denotes the frequency of the combined signal and can be expressed as:
f ( i 1 , i 2 , i 3 , i 4 , i 5 ) = i 1 · f 1 + i 2 · f 2 + + i 5 · f 5
One thing to note is that γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) only reflects the influence of the ionosphere on the ranging accuracy. To further analyze its impacts on ambiguity resolution, it should be converted to β ( i 1 , i 2 , i 3 , i 4 , i 5 ) in the unit of c y c l e · m 1 , as shown below:
β ( i 1 , i 2 , i 3 , i 4 , i 5 ) = f 1 2 · ( i 1 / f 1 + i 2 / f 2 + + i 5 / f 5 ) c

2.3.2. Ionosphere-Reduced Combination Selection

For Equation (28), the number of linear combinations is theoretically infinite. In order to find a combination with a sufficiently small ionosphere effect (e.g., IR combination), and at the same time, guarantee the performance of ambiguity resolution and positioning (e.g., decimeter-level), the following criteria are identified:
  • The impact of the ionosphere is small. Taking the 100 m ionosphere delay as an example, its impact on ambiguity resolution and ranging accuracy is less than 0.1 cycle and 0.15 m. The first criterion is β ( i 1 , i 2 , i 3 , i 4 , i 5 ) < 0.001 and γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) < 0.0015 .
  • The amplified noise level supports decimeter-level positioning. Considering that the accuracy ratio of carrier phase and pseudorange observations is usually 100:1, and the noise amplification factor of AFIF observations ranges from tens to hundreds, as shown in Table 2, this paper takes μ ( i 1 , i 2 , i 3 , i 4 , i 5 ) < 100 as the second criterion.
  • The combination should have a long wavelength to ensure efficient ambiguity resolution and to be able to resist the influence of residual geometry errors. This paper takes the equivalent wavelength in Table 2 as a reference and takes λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) > 3   m as the third criterion.
For coefficients i k ( k = 1 , 2 , 3 , 4 , 5 ) within range [–10, 10], the IR combination list in Table 3 is finally selected after traversing each combination. Although the noise amplification factor of the selected IR combination is close to 100, it is derived from raw carrier phase observations, which means that in addition to noise level, the multipath effects are theoretically weaker than the pseudorange.

2.3.3. Ionosphere-Reduced Observation Model

On the basis of the above selected IR combination, pseudorange observations are also needed to build a full-rank model. Like the aforementioned AFIF model and uncombined model, in order to avoid the additional code bias caused by the inconsistency between multi-frequency pseudoranges and precise clock references, only the E1/E5 dual-frequency pseudoranges are used.
For the selected IR observation Equation (28), if the influence of the residual ionosphere is small, it can be ignored directly, thereby reducing the dimension of the estimated parameters and improving the model strength to a certain extent. Depending on whether ionosphere parameters are considered, two models can be established.
If the ionosphere parameters are estimated, raw dual-frequency observations are integrated with Equation (28), and the corresponding observation model is as follows:
{ P r , 1 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + I ˜ r , 1 s + e r , 1 s P r , 2 s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + γ 2 · I ˜ r , 1 s + e r , 2 s L r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w     γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · I ˜ r , 1 s + λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s + μ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · ε r , 1 s
In case of n visible satellites, the estimated vector in Equation (32) includes:
X ^ = [ x , y , z t ˜ r T r , w I ˜ r , 1 1 , I ˜ r , 1 2 , , I ˜ r , 1 n n N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) 1 , N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) 2 , , N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) n n ]
If the ionosphere parameters in Equation (28) are directly ignored, a dual-frequency pseudorange in the form of an IF combination is used for parameter estimation. The observation model with the ionosphere ignored can be expressed as:
{ P r , IF s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + e r , IF s L r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s = ρ r s + c · t ˜ r c · t ˜ s + M r , w s · T r , w + λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) s + μ ( i 1 , i 2 , i 3 , i 4 , i 5 ) · ε r , 1 s
The corresponding estimated vector in Equation (34) contains:
X ^ = [ x , y , z t ˜ r T r , w N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) 1 , N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) 2 , , N ¯ r , ( i 1 , i 2 , i 3 , i 4 , i 5 ) n n ]
To simplify description, the terms IR(EST) and IR(IGN) are used to represent the models described in Equation (32) and Equation (34), respectively.

2.4. Comparison of Different Models

Table 4 briefly compares the above models in the case of n visible satellites per epoch. Since it is difficult to accurately estimate troposphere delay using only a single epoch of data, and also considering that the accuracy of the prior model can basically meet the decimeter-level positioning requirements, especially since zenith hydrostatic delay (ZHD) can be corrected with high accuracy, therefore, referring to the existing related research [10], the troposphere delay is directly corrected by a prior model instead of being estimated as an unknown, so the ZWD parameter does not appear in the table. The comparisons between various models adopted the same observation configurations, namely E1/E5a dual-frequency pseudoranges and E1/E5a/E6/E5b/E5 multi-frequency carrier phase observations.
In terms of ionosphere processing strategies, except for the UC and IR(EST) models that need to estimate n-dimensional ionosphere parameters, the other three models (i.e., AFIF(GF), AFIF(GB) and IR(IGN)) do not need to estimate the ionosphere. Among them, the two AFIF models use an IF combination to eliminate its influence, while the IR(IGN) model adopts the strategy of directly ignoring it.
In terms of the number of observation equations, the UC model requires the most equations for parameter estimation, up to 7n; the two AFIF models have the same number of equations, both at 4n. In contrast, the IR model requires fewer equations, and the number of IR(EST) and IR(IGN) models are 3n and 2n, respectively.
In terms of the ambiguity number involved in each model, the AFIF(GF) model resolves EWL/WL ambiguities in advance, and there is no need to estimate them in the observation model, while the other four models all need to estimate ambiguities of different dimensions. In particular, the ambiguity dimension of the UC model can reach up to 5n dimensions. In contrast, the AFIF(GB), IR(EST) and IR(IGN) models only need to estimate a small number of ambiguity parameters in n dimensions.
In terms of the number of fixed ambiguities, the UC and two AFIF models all need to fix 3n-dimensional EWL ambiguities and n-dimensional WL ambiguities, while the two IR models only need to fix n-dimensional EWL ambiguities without fixing WL ambiguities.
In terms of ambiguity fixing strategy, the UC and two IR models all adopt the strategy of GB+LAMBDA to fix all ambiguities [28]. For the two AFIF models, the strategies of EWL ambiguity fixing are the same, with both using GF+ROUND, while for WL ambiguity fixing, AFIF(GF) and AFIF(GB) adopt the different fixing strategies of GF+ROUND and GB+LAMBDA, respectively.
Compared with the UC model, the IR model has a smaller dimension of both observation equations and estimated parameters, so the calculation efficiency is higher; compared with the AFIF model, the IR model also has fewer observation equations and uses a more rigorous GB+LAMBDA strategy to resolve all ambiguities, so the performance of ambiguity resolution is better in theory. In summary, the IR model takes into account both performance and efficiency factors.

3. Experiments and Results

About 230 stations from the International GNSS Service (IGS) Multi-GNSS Experiment (MGEX) on 2023 DOY061 with sample rate of 30 s were selected for experiments, and their distribution is shown in Figure 1. All stations could receive Galileo E1/E5a/E6/E5b/E5ab five-frequency data, of which the stations marked by blue solid circles were used for FCBs estimation, and the stations marked with red stars were used for PPP verification. The selected 15 user stations are distributed in different latitudes and equipped with various types of receivers, as shown in Table 5. In data processing, the elevation cut-off angle was set to 10°, and the elevation dependent stochastic model was adopted with the priori precision of 0.003 m and 0.3 m for raw carrier and pseudorange, respectively. The precise orbit and clock products from the Center for Orbit Determination in Europe (CODE) were used. The troposphere delay was corrected by the Saastamoinen model with the pressure, temperature and water vapor information provided by the Global Pressure and Temperature 2 wet (GPT2w) [30,31], and the Vienna mapping function (VMF) was adopted [32]. For ambiguity fixing, a threshold of 0.25 cycles was set for integer rounding, and a ratio threshold of 2.0 was set for the LAMBDA algorithm. In addition, a partial ambiguity fixing (PAF) strategy was used simultaneously to improve the fixing success rate, and the ambiguity subset was determined according to successively increased elevation, until the number of satellites was less than five or the number of ambiguities was less than four [33]. Since there are not as many Galileo satellites as GPS, in order to ensure the redundancy of the model, only epochs with no less than five satellites were calculated. In order to avoid the large interpolation errors caused by orbit day boundary bias at the beginning and end of each day, only observation data from 1:00 to 23:00 were processed [34]. Figure 2 shows the sequence of Kp indexes, and it can be seen that the ionosphere was relatively quiet, and there was no obvious ionosphere scintillation phenomenon during the experiment period. The experimental data were processed by self-designed software, and the operating platform was a desktop computer equipped with an Intel i7-11700 CPU.

3.1. Performance Analysis of Ionosphere-Reduced Model

For the ionosphere-reduced model, although the IR(IGN) model with ionosphere estimated can be used, if the equivalent ranging errors are small and will not affect the ambiguity fixing, it can be ignored, and the IR(IGN) model is used instead. In this way, the dimensions of both estimated parameters and observation equations can be reduced, thereby improving model redundancy and data processing efficiency. Taking the BRUX station as an example, Figure 3 shows the slant ionosphere delay and the corresponding ranging errors of each satellite during the whole period, and the elevation information is also given in the figure. It can be seen from the figure that the variation range of the slant ionosphere was about −5~25 m, and its magnitude was relatively large at around 12:00. The corresponding equivalent ranging error of each satellite did not exceed 3 cm. Considering that the wavelength of the selected ionosphere-reduced combination is about 3.9 m, even if these residual errors are ignored, reliable EWL ambiguity resolution can still be achieved theoretically. Figure 4 further gives the positioning results of the two IR models. It can be seen that ignoring the ionosphere does not have a significant impact on the results, and there are only a few millimeters in difference in the RMS values of the two models. The positioning accuracy of the IR(IGN) and IR(EST) models in the north (N), east (E) and up (U) directions are (8.3, 5.7, 18.9) cm and (8.3, 5.7, 18.5) cm, respectively.
It can be preliminarily found from Figure 3 that the ionosphere delay and its corresponding ranging errors are related to satellite elevation. Therefore, Figure 5 shows the distribution relationship between the ionosphere ranging error and satellite elevation of all 15 user stations. It can be clearly found that low-elevation satellites often have larger residual ranging errors, and the ranging errors gradually decrease as the satellite elevation increases. At intervals of five degrees, the 95th, 90th and 80th percentile values within the range of 10–90 degrees are counted, as given in Figure 5 and Table 6. Even for low-elevation satellites (below 20 degrees), the 95th percentile value does not exceed 4 cm, which theoretically still supports EWL ambiguity resolution with a wavelength of several meters and positioning with decimeter-level accuracy. In practical applications, the impact of low-elevation satellites can be avoided or weakened by increasing the cut-off elevation (for example, in the case of multi-GNSS with sufficient satellites) or downweighting their observations.
Figure 6 shows the horizontal and vertical positioning accuracy of each station. The accuracies of the two IR models are almost the same. The horizontal and vertical accuracies of IR(IGN) and IR(EST) are (16.3, 35.8) cm and (16.2, 35.4) cm, respectively, and both can achieve decimeter-level positioning horizontally. Although the large residual ranging errors of low-elevation satellites will theoretically affect the positioning performance, the two models do not show significant differences, mainly for the following two reasons. First, the elevation-dependent stochastic model is adopted, and the weight of low-elevation satellites is small; secondly, in the iterative process of PAF, the low-elevation satellites are excluded, and the selected satellites are all high-elevation with small ionosphere effects.

3.2. Comparison between Ionosphere-Reduced and AFIF/Uncombined Models

In GNSS positioning, in addition to accuracy, real-time performance is also an important evaluation index. For this reason, this section will compare the IR model and the AFIF/UC model from the two aspects of positioning performance and time consumption.

3.2.1. Positioning Performance Comparison

Since the positioning performance of the two IR models is basically the same, the IR(IGN) model is selected in this section as a representative. Figure 7 shows the positioning results of the four models for the KITG station. For the two AFIF models, the main difference between the AFIF(GF) and AFIF(GB) models lies in the ambiguity fixing strategy used. The AFIF(GF) model adopts an integer rounding strategy based on the MW combination, while the AFIF(GB) model adopts an integer least square (ILS) method, namely the LAMBDA algorithm. Although the MW combination could eliminate the geometry and ionosphere terms, due to the enlarged noise of pseudorange observations, multiple epoch smoothing is often required to ensure reliable fixing. However, in single-epoch mode, direct rounding is likely to cause incorrect fixing, especially for combinations with short wavelengths. Figure 8 shows the fractional part distribution of different EWL/WL ambiguities in the AFIF(GF) model, where the FCBs are corrected in advance. It can be clearly seen that compared with the WL combination, the distribution of the three EWL combinations (i.e., E6–E5a/E5b–E5a/E5–E5a) is more concentrated, and the percentages within 0.25 cycles are close to 100%, reaching 99.94%, 100% and 100%, respectively. In contrast, the distribution of WL ambiguity is scattered, and the proportion within 0.25 cycles is only 90.80%. This means that in single-epoch mode, direct rounding is likely to cause incorrect ambiguity fixing, which will lead to larger positioning errors. As shown in Figure 7, the results of the AFIF(GF) model jump at multiple epochs, and the accuracy in the N, E and U directions is only (16.5, 14.0, 47.3) cm, respectively. Except for the AFIF(GF) model, the positioning results of the other three models (i.e., AFIF(GB), UC and IR(IGN)) are basically the same, and the accuracy in the N, E and U directions are (10.6, 10.1, 24.3) cm, (11.1, 9.2, 22.0) cm and (9.4, 8.2, 19.8) cm, respectively.
While the AFIF(GF) model uses an integer rounding strategy, the AFIF(GB), UC and IR(IGN) models all use the LAMBDA method with a ratio test to fix ambiguities. Figure 9 shows the ratio values of the three models. The ratio distributions of the different models are similar, and all three models are able to successfully fix ambiguities with large ratio values.
Figure 10 shows the epoch fixing rate of different models for each station. Due to geographical differences, the fixing rates vary among different stations. For the same station, the fixing rates of different models are very close. Among them, the AFIF(GF) model had the highest fixing rate of 98.9%, but it may include multiple wrongly fixed epochs. The fixing rates of the other three models (i.e., AFIF(GB), UC and IR(IGN)) were 98.6%, 98.7% and 98.7%, respectively.
Figure 11 shows the horizontal and vertical accuracy of different models for all stations. Since the results of the AFIF(GF) model have larger errors caused by incorrectly fixed ambiguities, its accuracy is the worst, with an average horizontal and vertical accuracy of (34.7, 71.5) cm. Except for the AFIF(GF) model, the accuracy of the other three models at different stations is basically the same, and the horizontal and vertical accuracy can basically be maintained within 20 cm and 50 cm, respectively. The horizontal and vertical average accuracies of the AFIF(GB), UC and IR(IGN) models are (17.8, 35.2) cm, (17.2, 33.8) cm and (16.3, 35.8) cm, respectively.

3.2.2. Time Consumption Comparison

Figure 12 shows the time consumption comparison per epoch of different models for each station. For the UC model, due to the high dimensions of both estimated parameters and observation equations, the processing efficiency is the lowest, and the results are not completely shown in Figure 12. The detailed data can be found in Table 7. Except for the UC model, the time consumption of the other four models is basically at the same level, maintaining around 0.5 ms. At the same time, it can be seen intuitively from the figure that the IR(IGN) model of each station has the shortest time consumption.
It can be seen from Table 7 that the UC model takes an average of 14.64ms to process one epoch data, which is significantly longer than other models. In contrast, the average time consumption of the AFIF(GF), AFIF(GB), IR(EST) and IR(IGN) models is only 0.50 ms, 0.69 ms, 0.56 ms and 0.42 ms, respectively. For the two AFIF models, previous analysis has shown that the positioning performance of the AFIF(GF) model is poor. Therefore, on the premise of ensuring positioning performance, the results of the AFIF(GB) model are selected for subsequent comparison. It can be found that the time consumption of both IR models is shorter than that of the AFIF(GB) model. For the two IR models, since the IR(EST) model estimates the ionosphere parameters and has more observation equations than the IR(IGN) model, its time consumption is also slightly longer.
Overall, the IR(IGN) model takes the least time consumption, being 38% shorter than the AFIF(GB) model and 97% shorter than the UC model. Compared with the AFIF/UC models, the IR model shows varying degrees of advantages in terms of process efficiency under the premise of ensuring positioning performance.

4. Discussion

With the development of multi-frequency GNSS, considering the accuracy, reliability, real-time calculation and other indicators comprehensively, instantaneous decimeter-level positioning with PPP has become a promising positioning mode. Based on existing research, two positioning models, namely AFIF and UC, were derived using different ionosphere processing strategies.
In the AFIF model, the AFIF observations, constructed by the IF combination of EWL and WL observations, are used for positioning. Since the coefficients of the IF combination are not integers, the combined ambiguity does not have integer characteristics, and the EWL and WL ambiguities need to be fixed separately during the implementation process. In the UC model, raw frequency ambiguity is mapped to EWL and WL ambiguity, and then decimeter-level positioning is achieved through the cascading ambiguity resolution. This model can theoretically process observation of any number of frequencies, but as it needs to deal with high-dimensional matrix operations, the computational load is relatively large.
In contrast, the IR model based on integer coefficients greatly eliminates the influence of ionosphere delay under the condition of ensuring the integer characteristics of the ambiguity. By selecting an appropriate wavelength and controlling noise level, in addition to similar decimeter-level accuracy comparable to the AFIF/UC models, the calculation efficiency is also improved to varying degrees due to the avoidance of high-dimensional matrix operations. The proposed IR model takes into account both positioning performance and calculation efficiency. Especially for the IR(IGN) model, at the expense of losing the unbiasedness of the observation model, the redundancy is improved, and the calculation efficiency is also greatly improved. After all, for the practical application of PPP, in addition to positioning accuracy, real-time performance is also an important indicator that users are concerned about.
Of course, the IR model also has certain limitations. First of all, the model is suitable for multi-frequency signal situations. In practical applications, when multi-frequency observations cannot be fully received, the model will not be applicable. Nevertheless, since the model is based on single-epoch mode, that is, there is no parameter transfer between previous and subsequent epochs, it can be flexibly switched to the AFIF or UC models that use fewer frequency signals between different epochs, when a certain frequency is missing. Secondly, the proposed IR models are based on the assumption that the residual ionosphere effect is small, which means that under extreme ionosphere conditions, the applicability of the IR models, especially the IR(IGN) model that ignores ionosphere effect, still needs further study. In a disturbed ionosphere case, as the ionosphere delay increases sharply, the corresponding positioning performance will be greatly reduced [35]. In such case, it is recommended to use the IR(EST) model in Equation (34) instead of blindly using the IR(IGN) model. Since the IR(EST) model properly estimates the ionosphere parameters, the model is theoretically still applicable.
This paper mainly evaluates the positioning accuracy and time-consuming performance of Galileo-only five-frequency IR PPP. On this basis, the next research work will be carried out from the following two aspects: first, to further study the performance and applicability of the proposed IR model in active ionosphere conditions; secondly, to further explore the use of IR combinations from other multi-frequency GNSSs, such as BDS-3, and to evaluate the performance of multi-GNSS IR PPP.

5. Conclusions

This paper investigates the theory and feasibility of a multi-frequency ionosphere-reduced (IR) combination in the application of single-epoch PPP with decimeter-level accuracy. Firstly, the existing AFIF and UC models are introduced. Then, based on multi-frequency observation combination theory, the selection method of an ionosphere-reduced (IR) combination is studied by taking Galileo five-frequency signals as an example. Based on the selected IR combination, two IR positioning models, namely IR(EST) and IR(IGN), are proposed according to whether ionosphere parameters are estimated or ignored, and the proposed models are further compared with the existing AFIF/UC models in various aspects.
Real tracked observations from globally distributed stations were used for model validation. The magnitude of the ionosphere equivalent ranging errors in the IR models and its relationship with satellite elevation were analyzed. The higher the satellite elevation, the smaller the residual ranging errors. The 95th percentile value of low-elevation satellites generally does not exceed 4 cm, which theoretically still supports EWL ambiguity resolution with a wavelength of several meters. In practical applications, the impact of low-elevation satellites can be avoided or weakened by increasing the cut-off elevation or downweighting their observations. The positioning results of the two IR models are almost the same, and the horizontal and vertical accuracies of IR(IGN) and IR(EST) are (16.3, 35.8) cm and (16.2, 35.4) cm, respectively.
The proposed IR models are further compared with the existing AFIF/UC models in terms of positioning performance and time consumption. In terms of positioning performance, the AFIF(GF) model suffers from incorrect ambiguity fixing and has the worst accuracy. In contrast, the IR and AFIF(GB)/UC models all use the LAMBDA method to determine the optimal candidate and can successfully fix ambiguities with large ratio values. The corresponding positioning performance of different models is basically the same, and all can achieve horizontal and vertical accuracy better than 20 cm and 40 cm, respectively. In terms of time consumption, the UC model suffers from high dimensions of both estimated parameters and observation equations, with the longest average processing time of one epoch data, reaching 14.64 ms; meanwhile, the other models take a relatively shorter time, maintaining around 0.5 ms, and the IR(IGN) model takes the shortest time of 0.42 ms.
Overall, the proposed IR models can achieve single-epoch decimeter-level positioning, and this performance is comparable to that of the existing AFIF/UC models. On the basis of guaranteeing positioning performance, the IR model shows varying degrees of advantages in terms of time consumption, compared with the AFIF/UC models. In particular, the IR(IGN) model, with no ionosphere parameters estimated, has the shortest time consumption, being 38% faster than the AFIF(GB) model and 97% faster than the UC model.

Author Contributions

Conceptualization, Q.Z., S.P., J.L. and W.G.; Investigation, Q.Z. and S.P.; Methodology, Q.Z. and W.G.; Validation, S.P., Y.L. and P.Z.; Writing—original draft, Q.Z., S.P., J.L. and W.G.; Writing—review and editing, Y.L. and P.Z.; Supervision, S.P. and W.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by the National Key Research and Development Program of China (Grant No. 2021YFC3000502), the National Natural Science Foundation of China (Grant No. 42204027) and the Foundation of Laboratory of Science and Technology on Marine Navigation and Control, China State Shipbuilding Corporation (Grant No. 2021010104).

Data Availability Statement

The multi-GNSS precise product provided by CODE is available at https://cddis.nasa.gov/archive/gnss/products/ (accessed on 20 June 2023). The MGEX multi-frequency observation is available at https://cddis.nasa.gov/archive/gnss/data/daily/ (accessed on 20 June 2023).

Acknowledgments

The authors sincerely thank IGS and CODE for providing multi-GNSS data and products. The authors would like to thank four anonymous reviewers for their constructive comments and useful suggestions to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zumberge, J.F.; Heflin, M.B.; Jefferson, D.C.; Watkins, M.M.; Webb, F.H. Precise point positioning for the efficient and robust analysis of GPS data from large networks. J. Geophys. Res. Solid. Earth. 1997, 102, 5005–5017. [Google Scholar] [CrossRef] [Green Version]
  2. Kouba, J.; Héroux, P. GPS precise point positioning using IGS orbit products. GPS Solut. 2001, 5, 12–28. [Google Scholar] [CrossRef]
  3. Cai, C.; Gao, Y. Modeling and assessment of combined GPS/GLONASS precise point positioning. GPS Solut. 2013, 17, 223–236. [Google Scholar] [CrossRef]
  4. Li, P.; Zhang, X. Integrating GPS and GLONASS to accelerate convergence and initialization times of precise point positioning. GPS Solut. 2014, 18, 461–471. [Google Scholar] [CrossRef]
  5. Ge, M.; Gendt, G.; Rothacher, M.; Shi, C.; Liu, J. Resolution of GPS carrier-phase ambiguities in precise point positioning (PPP) with daily observations. J. Geod. 2008, 82, 389–399. [Google Scholar] [CrossRef]
  6. Zhang, X.; Li, P.; Guo, F. Ambiguity resolution in precise point positioning with hourly data for global single receiver. Adv. Space Res. 2013, 51, 153–161. [Google Scholar] [CrossRef]
  7. Wabbena, G.; Schmitz, M.; Bagge, A. PPP-RTK: Precise point positioning using state-space representation in RTK networks. In Proceedings of the 18th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2005), Long Beach, CA, USA, 16 September 2005; pp. 2584–2594. [Google Scholar]
  8. Zhang, B.; Chen, Y.; Yuan, Y. PPP-RTK based on undifferenced and uncombined observations: Theoretical and practical aspects. J. Geod. 2019, 93, 1011–1024. [Google Scholar] [CrossRef]
  9. Geng, J.; Guo, J.; Chang, H.; Li, X. Toward global instantaneous decimeter-level positioning using tightly coupled multi-constellation and multi-frequency GNSS. J. Geod. 2019, 93, 977–991. [Google Scholar] [CrossRef] [Green Version]
  10. Geng, J.; Guo, J. Beyond three frequencies: An extendable model for single-epoch decimeter-level point positioning by exploiting Galileo and BeiDou-3 signals. J. Geod. 2020, 94, 14. [Google Scholar] [CrossRef]
  11. Laurichesse, D.; Banville, S. Innovation: Instantaneous Centimeter-Level Multi-frequency Precise Point Positioning. GPS World 2018, 29, 43–47. [Google Scholar]
  12. Geng, J.; Bock, Y. Triple-frequency GPS precise point positioning with rapid ambiguity resolution. J. Geod. 2013, 87, 449–460. [Google Scholar] [CrossRef]
  13. Li, X.; Li, X.; Liu, G.; Feng, G.; Yuan, Y.; Zhang, K.; Ren, X. Triple-frequency PPP ambiguity resolution with multi-constellation GNSS: BDS and Galileo. J. Geod. 2019, 93, 1105–1122. [Google Scholar] [CrossRef]
  14. Li, X.; Liu, G.; Li, X.; Zhou, F.; Feng, G.; Yuan, Y.; Zhang, K. Galileo PPP rapid ambiguity resolution with five-frequency observations. GPS Solut. 2020, 24, 1–13. [Google Scholar] [CrossRef]
  15. Guo, J.; Xin, S. Toward single-epoch 10-centimeter precise point positioning using Galileo E1/E5a and E6 signals. In Proceedings of the 32nd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2019), Miami, FL, USA, 20 September 2019; pp. 2870–2887. [Google Scholar]
  16. Zhao, Q.; Pan, S.; Gao, W.; Wu, B. Multi-GNSS fast precise point positioning with multi-frequency uncombined model and cascading ambiguity resolution. Math. Probl. Eng. 2022, 2022, 3478332. [Google Scholar] [CrossRef]
  17. Melbourne, W.G. The case for ranging in GPS-based geodetic systems. In Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15 April 1985; pp. 373–386. [Google Scholar]
  18. Wubbena, G. Software developments for geodetic positioning with GPS using TI 4100 code and carrier measurements. In Proceedings of the 1st International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15 April 1985; pp. 403–412. [Google Scholar]
  19. Gu, S.; Lou, Y.; Shi, C.; Liu, J. BeiDou phase bias estimation and its application in precise point positioning with triple-frequency observable. J. Geod. 2015, 89, 979–992. [Google Scholar] [CrossRef]
  20. Gao, W.; Zhao, Q.; Meng, X.; Pan, S. Performance of single-epoch EWL/WL/NL ambiguity-fixed precise point positioning with regional atmosphere modelling. Remote Sens. 2021, 13, 3758. [Google Scholar] [CrossRef]
  21. Qu, L.; Wang, L.; Acharya, T.D.; Du, Y.; Wang, H.; Jiang, W. Global single-epoch narrow-lane ambiguity resolution with multi-constellation and multi-frequency precise point positioning. GPS Solut. 2023, 27, 29. [Google Scholar] [CrossRef]
  22. Feng, Y. GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals. J. Geod. 2008, 82, 847–862. [Google Scholar] [CrossRef]
  23. Guo, Z.; Yu, X.; Hu, C.; Jiang, C.; Zhu, M. Research on linear combination models of BDS multi-frequency observations and their characteristics. Sustainability 2022, 14, 8644. [Google Scholar] [CrossRef]
  24. Li, J.; Yang, Y.; He, H.; Guo, H. An analytical study on the carrier-phase linear combinations for triple-frequency GNSS. J. Geod. 2017, 91, 151–166. [Google Scholar] [CrossRef]
  25. Li, J.; Yang, Y.; He, H.; Guo, H. Benefits of BDS-3 B1C/B1I/B2a triple-frequency signals on precise positioning and ambiguity resolution. GPS Solut. 2020, 24, 100. [Google Scholar] [CrossRef]
  26. Schmid, R.; Steigenberger, P.; Gendt, G.; Ge, M.; Rothacher, M. Generation of a consistent absolute phase-center correction model for GPS receiver and satellite antennas. J. Geod. 2007, 81, 781–798. [Google Scholar] [CrossRef] [Green Version]
  27. Wu, J.T.; Wu, S.C.; Hajj, G.A.; Bertiger, W.I.; Lichten, S.M. Effects of antenna orientation on GPS carrier phase. Manuscr. Geod. 1993, 18, 91–98. [Google Scholar]
  28. Teunnissen, P.J.G. The least-square ambiguity decorrelation adjustment: A method for fast GPS integer ambiguity estimation. J. Geod. 1995, 70, 65–82. [Google Scholar] [CrossRef]
  29. Liu, L.; Pan, S.; Gao, W.; Ma, C.; Tao, J.; Zhao, Q. Assessment of quad-frequency long-baseline positioning with BeiDou-3 and Galileo observations. Remote Sens. 2021, 13, 1551. [Google Scholar] [CrossRef]
  30. Saastamoinen, J. Atmospheric correction for the troposphere and stratosphere in radio ranging satellites. Use Artif. Satell. Geod. 1972, 15, 247–251. [Google Scholar]
  31. Böhm, J.; Möller, G.; Schindelegger, M.; Pain, G.; Weber, R. Development of an improved empirical model for slant delays in the troposphere (GPT2w). GPS Solut. 2015, 19, 433–441. [Google Scholar] [CrossRef] [Green Version]
  32. Landskron, D.; Böhm, J. VMF3/GPT3: Refined discrete and empirical troposphere mapping functions. J. Geod. 2018, 92, 349–360. [Google Scholar] [CrossRef]
  33. Gao, W.; Gao, C.; Pan, S. A method of GPS/BDS/GLONASS combined RTK positioning for middle-long baseline with partial ambiguity resolution. Surv. Rev. 2017, 49, 212–220. [Google Scholar] [CrossRef]
  34. Griffiths, J.; Ray, J.R. On the precision and accuracy of IGS orbits. J. Geod. 2009, 83, 277–287. [Google Scholar] [CrossRef]
  35. Kuverova, V.V.; Adamson, S.O.; Berlin, A.A.; Bychkov, V.L.; Dmitriev, A.V.; Dyakov, Y.A.; Eppelbaum, L.V.; Golubkov, G.V.; Lushnikov, A.A.; Manzhelii, M.I.; et al. Chemical physics of D and E layers of the ionosphere. Adv. Space Res. 2019, 64, 1876–1886. [Google Scholar] [CrossRef]
Figure 1. IGS MGEX station distribution: blue circles denote stations used for FCBs estimation; red stars denote user stations.
Figure 1. IGS MGEX station distribution: blue circles denote stations used for FCBs estimation; red stars denote user stations.
Remotesensing 15 03562 g001
Figure 2. The Kp index series provided by GFZ on 2023 DOY 061 (https://kp.gfz-potsdam.de/en/figures/kp-daily-plots, accessed on 20 June 2023).
Figure 2. The Kp index series provided by GFZ on 2023 DOY 061 (https://kp.gfz-potsdam.de/en/figures/kp-daily-plots, accessed on 20 June 2023).
Remotesensing 15 03562 g002
Figure 3. Satellite elevation (up), slant ionosphere (middle) and the corresponding ranging errors (down) at the BRUX station.
Figure 3. Satellite elevation (up), slant ionosphere (middle) and the corresponding ranging errors (down) at the BRUX station.
Remotesensing 15 03562 g003
Figure 4. The positioning results of the IR(IGN) and IR(EST) models at the BRUX station.
Figure 4. The positioning results of the IR(IGN) and IR(EST) models at the BRUX station.
Remotesensing 15 03562 g004
Figure 5. Distribution of ionosphere equivalent ranging errors with satellite elevation for all stations.
Figure 5. Distribution of ionosphere equivalent ranging errors with satellite elevation for all stations.
Remotesensing 15 03562 g005
Figure 6. Horizontal and vertical accuracy of IR(IGN) and IR(EST) models for all 15 stations.
Figure 6. Horizontal and vertical accuracy of IR(IGN) and IR(EST) models for all 15 stations.
Remotesensing 15 03562 g006
Figure 7. Positioning results of the AFIF(GF), AFIF(GB), UC and IR(IGN) models at the KITG station.
Figure 7. Positioning results of the AFIF(GF), AFIF(GB), UC and IR(IGN) models at the KITG station.
Remotesensing 15 03562 g007
Figure 8. Fractional part distribution of EWL/WL ambiguities in the AFIF(GF) model at the KITG station.
Figure 8. Fractional part distribution of EWL/WL ambiguities in the AFIF(GF) model at the KITG station.
Remotesensing 15 03562 g008
Figure 9. Ratio values of the AFIF(GB), UC and IR(IGN) models at the KITG station.
Figure 9. Ratio values of the AFIF(GB), UC and IR(IGN) models at the KITG station.
Remotesensing 15 03562 g009
Figure 10. Epoch fixing rates of the AFIF(GF), AFIF(GB), UC and IR(IGN) models.
Figure 10. Epoch fixing rates of the AFIF(GF), AFIF(GB), UC and IR(IGN) models.
Remotesensing 15 03562 g010
Figure 11. Horizontal and vertical accuracy of the AFIF(GF), AFIF(GB), UC and IR(IGN) models.
Figure 11. Horizontal and vertical accuracy of the AFIF(GF), AFIF(GB), UC and IR(IGN) models.
Remotesensing 15 03562 g011
Figure 12. Time consumption comparison per epoch of different models for each station.
Figure 12. Time consumption comparison per epoch of different models for each station.
Remotesensing 15 03562 g012
Table 1. Specific frequency of Galileo E1/E5a/E6/E5b/E5.
Table 1. Specific frequency of Galileo E1/E5a/E6/E5b/E5.
Frequency NotationFrequency IDFrequency/MHz
f 1 E11575.42
f 2 E5a1176.45
f 3 E61278.75
f 4 E5b1207.14
f 5 E51191.795
Table 2. Noise amplification factors (AF) and equivalent wavelengths for different AFIF combinations.
Table 2. Noise amplification factors (AF) and equivalent wavelengths for different AFIF combinations.
AFIF
Observation
EWL/WL
Combination
Noise AFPrior AccuracyWavelength AFEquivalent
Wavelength
L r , AFIF 123 s E6-E5a/E1-E5a77.70.233 m5.33.990 m
L r , AFIF 124 s E5b-E5a/E1-E5a181.30.544 m4.33.214 m
L r , AFIF 125 s E5-E5a/E1-E5a339.61.019 m4.13.086 m
Table 3. Galileo five-frequency ionosphere-reduced (IR) combination.
Table 3. Galileo five-frequency ionosphere-reduced (IR) combination.
E1E5aE6E5bE5ab λ ( i 1 , i 2 , i 3 , i 4 , i 5 ) γ ( i 1 , i 2 , i 3 , i 4 , i 5 ) β ( i 1 , i 2 , i 3 , i 4 , i 5 ) μ ( i 1 , i 2 , i 3 , i 4 , i 5 )
i 1 i 2 i 3 i 4 i 5
14−31−33.907 m−0.0012−0.000395.407
Table 4. Comparison of AFIF, Uncombined (UC) and Ionosphere-Reduced (IR) models.
Table 4. Comparison of AFIF, Uncombined (UC) and Ionosphere-Reduced (IR) models.
ItemAFIFUncombinedIonosphere-Reduced
AFIF(GF)AFIF(GB)UCIR(EST)IR(IGN)
Observation configurationAll models use E1/E5a dual-frequency pseudoranges and E1/E5a/E6/E5b/E5 five-frequency carrier phase observations.
Observation model equation(14)(16)(19)(32)(34)
Estimated vector equation(15)(18)(21)(33)(35)
Number of different types of estimated parameters
Coordinates
33333
Receiver clock
11111
Ionosphere
00nn0
Ambiguity
0n5nnn
Ionosphere delayEliminateEliminateEstimateEstimateIgnore
Number of observation equation 1(3 + 1)n(3 + 1)n(5 + 2)n(1 + 2)n(1 + 1)n
Number of fixed ambiguities 2(3 + 1)n(3 + 1)n(3 + 1)n(1 + 0)n(1 + 0)n
EWL ambiguity fixingGF+ROUNDGF+ROUNDGB+LAMBDAGB+LAMBDAGB+LAMBDA
WL ambiguity fixingGF+ROUNDGB+LAMBDAGB+LAMBDA--
1 The first and second numbers in brackets indicate the number of carrier phase and pseudorange observation equations for each satellite, respectively. 2 The first and second numbers in brackets indicate the number of fixed EWL and WL ambiguities for each satellite, respectively.
Table 5. Location, receiver and visible satellites of user stations (10 deg cutoff elevation).
Table 5. Location, receiver and visible satellites of user stations (10 deg cutoff elevation).
SiteReceiverLatitude/degLongitude/degAverage Number of Satellites per Epoch
AMC4SEPT POLARX5TR38.8−104.56.4
AREGSEPT POLARX5−16.5−71.56.9
BRUXSEPT POLARX5TR50.84.46.8
GAMGSEPT POLARX5TR35.6127.97.3
HARBSEPT POLARX5TR−25.927.77.2
KAT1SEPT POLARX5−14.4132.28.0
KIR8TRIMBLE ALLOY67.921.17.2
KIRISEPT POLARX51.4172.98.3
KITGSEPT POLARX539.166.97.2
KOUGSEPT POLARX5TR5.1−52.67.3
MCHLTRIMBLE ALLOY−26.4148.17.5
MET3JAVAD TRE_3 DELTA60.224.47.5
RGDGTRIMBLE ALLOY−53.8−67.86.5
SEYGSEPT POLARX5−4.755.58.6
YEL2SEPT POLARX5TR62.5−114.57.3
Table 6. The 95th, 90th and 80th percentile ranging errors in different elevation intervals.
Table 6. The 95th, 90th and 80th percentile ranging errors in different elevation intervals.
Elevation Interval
[deg]
Equivalent Ranging Error [cm]
95th Percentile Value90th Percentile Value80th Percentile Value
10–154.03.42.6
15–203.93.32.5
20–253.42.92.2
25–303.12.82.1
30–352.92.51.9
35–402.62.01.6
40–452.31.91.5
45–502.21.81.4
50–552.11.71.3
55–602.01.61.2
60–652.01.61.2
65–701.81.51.1
70–751.51.41.1
75–801.71.41.1
80–851.51.20.9
85–901.51.41.0
Table 7. Time consumption per epoch of different models for each station.
Table 7. Time consumption per epoch of different models for each station.
SiteTime Consumption per Epoch [ms]
AFIF(GF)AFIF(GB)UCIR(EST)IR(IGN)
AMC40.430.6313.760.510.38
AREG0.460.6414.250.530.42
BRUX0.520.6714.210.520.40
GAMG0.490.7014.720.580.45
HARB0.510.6914.580.570.44
KAT10.520.7115.380.600.43
KIR80.460.6714.570.530.43
KIRI0.510.7115.280.580.43
KITG0.520.6814.380.570.42
KOUG0.510.6814.740.550.42
MCHL0.490.7014.830.560.42
MET30.480.7114.730.560.42
RGDG0.460.6513.620.520.41
SEYG0.580.7315.720.610.45
YEL20.510.7314.780.560.43
Average0.500.6914.640.560.42
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, Q.; Pan, S.; Liu, J.; Lu, Y.; Zhang, P.; Gao, W. Single-Epoch Decimeter-Level Precise Point Positioning with a Galileo Five-Frequency Ionosphere-Reduced Combination. Remote Sens. 2023, 15, 3562. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15143562

AMA Style

Zhao Q, Pan S, Liu J, Lu Y, Zhang P, Gao W. Single-Epoch Decimeter-Level Precise Point Positioning with a Galileo Five-Frequency Ionosphere-Reduced Combination. Remote Sensing. 2023; 15(14):3562. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15143562

Chicago/Turabian Style

Zhao, Qing, Shuguo Pan, Ji Liu, Yin Lu, Peng Zhang, and Wang Gao. 2023. "Single-Epoch Decimeter-Level Precise Point Positioning with a Galileo Five-Frequency Ionosphere-Reduced Combination" Remote Sensing 15, no. 14: 3562. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15143562

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