Improving the Triple-Carrier Ambiguity Resolution with a New Ionosphere-Free and Variance-Restricted Method

The ionospheric bias and the combined observation noise are two crucial factors affecting the reliability of the triple-carrier ambiguity resolution (TCAR). In order to obtain a better reliability of TCAR, a new ionosphere-free and variance-restricted TCAR method is proposed through exploring the ambiguity link between each step of TCAR. The method constructs an ionosphere-free combination and simultaneously restricts the combined observation noise with respect to the wavelength to a sufficiently low level for each step of TCAR. The performance of the proposed method is tested by the datasets from the BeiDou navigation satellite system (BDS), with the baseline varying from 7.7 km to 68.8 km. Comparing with the state-of-the-art TCAR methods, the experimental results indicate that the proposed method can obtain a better performance of ambiguity resolution, even though the double-differenced ionospheric delay increases up to 72.4 cm at the baseline of 68.8 km.


Introduction
Integer ambiguity resolution is the key technique for high-precision relative positioning using carrier phase observation such as real-time kinematic (RTK) [1].With the modernization of global navigation satellite system (GNSS), significant efforts have been made to achieve reliable ambiguity resolution by using the triple-frequency signals.The triple-carrier ambiguity resolution (TCAR) method proposed by the authors of [2,3] is one of the widely-used approaches for resolving the triple-frequency ambiguity in a short baseline case so that the ionospheric error can be neglected.When the ionospheric bias cannot be neglected, especially in medium and long baseline cases, much research has also been done to improve its performance [4][5][6][7].A common conclusion is that both the ionospheric bias and combined observation noise are two crucial factors for the reliability of TCAR [4,8].Therefore, it is very necessary to eliminate the ionospheric bias and restrict the combined observation noise level simultaneously to improve the reliability of TCAR.
Considering these two crucial factors for TCAR, the ionosphere-estimated method can be used to compensate the effect of ionospheric bias without amplifying the observation noise [6,9].However, the performance of the ionosphere-estimated method depends on the precision of the ionospheric model, which includes a functional and stochastic model.Moreover, the increased ionosphere-related parameters will reduce the strength of the estimated-model [10].The rate of total electrical content (TEC) corrections can be utilized to improve the strength of the ionosphere-estimated model, and it has proven to be efficient for different latitudes and ionospheric conditions [11,12].
In addition to the ionosphere-estimated method, most studies of TCAR focus on choosing the optimal combination coefficients among triple-frequency observations.These studies aim at suppressing the effect of ionospheric bias and combined observation noise [4,5,8,13,14].The coefficients can be determined based on different categories to suppress these two kinds of observation errors.The widely-used category includes the ionosphere-free combination [13,[15][16][17] and the ionosphere-reduced combination [4].The ionosphere-free combination ensures the absence of ionospheric bias, but the combined observation noise is largely amplified, which will worsen ambiguity resolution performance [18,19].Therefore, a moving average method has to be used to restrict the effect of combined observation noise [7,13,17,20,21].Furthermore, the ionosphere-reduced combination can also be applied because it is in essence a trade-off between the ionospheric effect reduction and observation variance restriction.It aims at suppressing the effect of ionospheric bias and combined observation noise simultaneously [4].Nevertheless, the residual of ionospheric delay is still a threat for the reliability of TCAR.
The ionosphere is difficult to be predicted or modeled precisely because it is subjected to undetermined variation in temporal and spatial dimensions [22].Therefore, the reliability of the ionosphere-estimated method has its limitations.In contrast, the observation combination method has enormous potential in mitigating the effect of ionospheric bias so that the reliability of TCAR can be improved [5,8,13].Typically, the ionosphere-free combination is usually used for eliminating the effect of ionospheric bias on TCAR.However, the current ionosphere-free combined methods still suffer from the contradiction between eliminating ionospheric bias and suppressing combined observation noise.The contradiction is reflected by the fact that, if the ionospheric bias is eliminated, the combined observation noise will be amplified.If the combined observation noise is restricted, the ionospheric bias is difficult to be eliminated for the traditional TCAR methods.Although the traditional ionosphere-free combination method takes the advantage of ambiguity-corrected observations from previous steps to restrict the combined observation noise with respect to the wavelength, the corresponding combination still does not successfully restrict the noise with respect to the wavelength to a sufficiently low level [13,17].To restrict the combined observation noise level under the ionosphere-free combination, we propose an ionosphere-free and variance-restricted (IFVR) TCAR method in this contribution."Variance-restricted" means the combined observation noise with respect to the wavelength can be restricted to a sufficiently low level so as to achieve reliable ambiguity resolution.The method takes into account the available information of the ambiguity link between different steps, which can be used to build a relationship between observations with respect to different frequencies.This relationship can provide an opportunity to form an optimal combination method so that the ionospheric bias can be eliminated and the combined observation noise can be restricted simultaneously.Taking advantages of ambiguity links between each step of TCAR for the proposed method, a better ambiguity resolution performance can therefore be anticipated.
The remaining parts of the paper are organized as follows: the detail derivation of the novel TCAR method is introduced in Section 2. The performance of the proposed TCAR method is extensively examined in Section 3. Finally, the conclusion remarks are summarized in Section 4.

Methodology
In order to demonstrate the proposed TCAR method, we will briefly review the general linear combinations between the double-differenced code and carrier phase observations, shown by [4].
where combination coefficients (i, j, k) ∈ Z are the arbitrary integer and f t (t = 1, 2, 3) denote BDS three frequency signals, i.e., f 1 = 1561.098MHz, f 2 = 1207.140MHz and f 3 = 1268.520MHz, respectively; p and φ represent code and phase observations in meters, respectively; ρ describes the non-dispersive terms including geometric distance and the tropospheric delay in meters; I 1 is the first-order ionospheric delay on frequency f 1 in meters; λ (i,j,k) and N (i,j,k) are wavelength in meters and ambiguity in cycles for different coefficients combination; ε p and ε φ represent noise of double-differenced code and phase observations, respectively.The ionospheric scale factor β (i,j,k) , the noise scale factor γ (i,j,k) , wavelength λ (i,j,k) , and ambiguity N (i,j,k) based on the combination, are defined as (Feng 2008) where c is the speed of light.
Similar to the traditional TCAR methods, the proposed TCAR method comprises three steps: extra-wide-lane (EWL), wide-lane (WL) and narrow-lane (NL) ambiguity resolution.EWL ambiguity can be reliably resolved by previous studies [4,6,7], but rapid and reliable WL and NL ambiguity resolution is still challenging because of both ionospheric bias and the combined observation noise [6].Therefore, the proposed TCAR method focuses on finding a better WL and NL combination method to eliminate ionospheric bias and restrict combined observation noise simultaneously.To simplify the quantitative analysis below, we assume all un-differenced observations are independent with each other, the observations among different frequencies have the same precision, i.e., un-differenced code and phase precision take 30 cm and 3 mm, respectively.Thus, the standard deviations (STD) of corresponding double-differenced observations are σ p = 60 cm and σ φ = 6 mm, respectively.The proposed TCAR method can be implemented by the following three-step procedure.

EWL Ambiguity Resolution
The widely-used EWL combination can be expressed as [4,6] where φ EWL represents EWL combination: the corresponding wavelength is λ EWL ≈ 4.88 m.Based on the assumed observation precision, the ratio of the noise to the wavelength is σ EWL /λ EWL ≈ 0.09 cycles.The sufficiently low σ EWL /λ EWL indicates the success rate of EWL ambiguity resolution by using the constructed EWL combination is close to 100%.
In the new WL combination, it can be found that the ratio of the noise to the wavelength σ WL1 /λ WL1 is only 0.15 cycles, although the noise has been amplified by about 114.4 times.Therefore, the new WL combination method can sufficiently restrict the combined observation noise under ionosphere-free combination, which is able to obtain a better ambiguity resolution performance.
Because of the presence of geometry terms, the rank deficiency issue exists for (5).Therefore, additional ionosphere-free combination should be constructed.Furthermore, the combination to be constructed and φ WL1 must be mutually independent.One option is to construct an ionosphere-free combination by using code observations [19].In view of the lower precision of code observations, the STD of combined observation noise using the code observation will be enlarged to meter level, e.g., the noise of ionosphere-free combination would reach up to 1.74 m by using f 1 and f 2 code observations [4,19].In order to avoid the enlarged observation noise under the ionosphere-free combination, another new WL ionosphere-free combination contained WL ambiguity N (1,0,−1) is constructed taking full advantage of EWL ambiguity-corrected observation φ (0,−1,1) = φ (0,−1,1) − λ (0,−1,1) N (0,−1,1) , written as where φ WL2 denotes another WL ionosphere-free combination; (l,m,n) ∈ Z represents combined coefficients of code observation.In order to construct ionosphere-free combination and retain geometry terms, the coefficients set b 1 = . The corresponding wavelength and noise are The new φ WL2 eliminates the rank deficiency issue of ( 5), which is independent with φ WL1 because of the introduction of code observations.Meanwhile, in order to obtain the best reliability of WL ambiguity resolution, the minimization criterion of σ WL2 /λ WL2 should be satisfied.In this paper, the integer coefficients of the code observation are selected as (l,m,n) = (0,0,1), the corresponding coefficients are b 1 ≈ −4.20 and b 2 ≈ 4.20, respectively.The integer coefficient selection is detailed in Appendix A. As a result, the wavelength λ WL2 is 4.30 m and the ratio of the noise to the wavelength is σ WL2 /λ WL2 ≈ 0.23 cycles based on the assumed observation precision.
These two WL ionosphere-free combinations from ( 5) and ( 6) take full advantage of ambiguity links and ambiguity-corrected observations to restrict the combined observation noise with respect to wavelength sufficiently.Moreover, because the two WL observations contain the same WL ambiguity N (1,0,−1) depending on the ambiguity link between the EWL and WL ambiguity, the strength of the WL ambiguity model can be improved by these advantages.Using these two WL ionosphere-free combinations, we construct a new linear model for WL ambiguity resolution.In view of presence of non-dispersive terms, the residual tropospheric bias also needs to be taken into account as estimated parameter to compensate the tropospheric bias, which is detailed in Section 3. Therefore, the WL linear model writes where Φ represents the vector of modelled combined observations; H denotes the geometry matrix; m is tropospheric mapping function; b includes baseline components; τ is zenith tropospheric wet delay (ZWD); N (1,0,−1) is vector of ambiguities; I is identity matrix.After applying least-square estimation, the WL ambiguities N (1,0,−1) can be determined using the least-squares ambiguity decorrelation adjustment (LAMBDA) algorithm [1].

NL Ambiguity Resolution
The reliability of NL ambiguity resolution determines the performance of precise positioning.Compared with WL and EWL ambiguity resolution, the NL ambiguity resolution is more sensitive to the ionospheric bias and the combined observation noise [17,24].Therefore, it is very necessary to eliminate the ionospheric bias and restrict the combined observation noise simultaneously for improving the reliability of NL ambiguity resolution.
Assuming that the EWL and WL ambiguity are correctly resolved above, the N (0,−1,1) and N (1,0,−1) can be used to form the ambiguity link between EWL, WL and NL ambiguities, shown by Due to the ambiguity link, we can construct two new ionosphere-free combinations contained the same ambiguity N (1,0,0) for NL ambiguity resolution, which can be expressed as where φ NL1 and φ NL2 denote two independent ionosphere-free combinations by using different phase observations in NL ambiguity resolution.In order to construct ionosphere-free combination and retain the geometry terms, the coefficients set ≈ −1.94, respectively.The wavelength and the noise are σ φ ,respectively.The corresponding ratio of the noise to the wavelength are σ NL1 /λ NL1 ≈ 0.16 cycles and σ NL2 /λ NL2 ≈ 0.20 cycles, respectively.
Similar to WL ambiguity resolution, depending on the advantage of the ambiguity link (8), we construct two new ionosphere-free combinations from (9) with sufficiently combined observation noise restriction.Using these two ionosphere-free combinations, the linear model for NL ambiguity resolution can be expressed as Similar to the WL ambiguity resolution in (7), it can be seen that ambiguities N (1,0,0) can be resolved by these two combinations.After applying a least-square estimation, the NL ambiguities N (1,0,0) can be determined using the LAMBDA algorithm.Then, the baseline parameters are adjusted by the integer constraints of the corrected NL ambiguities to obtain high-precision positioning solutions [1].In order to clearly understand the proposed method, the coefficients for the EWL, WL, NL observations of the proposed IFVR method are listed in Table 1.

Observations Combination
According to the description of the proposed TCAR method above, the advantage of the ambiguity link between each step of TCAR is demonstrated for WL and NL ambiguity resolution.The proposed IFVR method precisely eliminates the ionospheric bias and sufficiently restricts the variance of combined observation noise by the ambiguity link.From ( 5) and ( 6) for WL ambiguity resolution and (9) for NL ambiguity resolution, it can be seen that each new combination restricts the ratio of the noise to the wavelength within 0.25 cycles based on the assumed observation precision.Furthermore, the strength of WL and NL ambiguity resolution is also enhanced by the ambiguity link, which is beneficial for improving the reliability of ambiguity resolution.It can be found that the combined observation in (7) and (10) are independent from each other, and, due to the corresponding design matrix in (7) and (10) are full row rank, it can also be deduced that (7) and (10) are independent from each other.Therefore, the model strength of proposed TCAR is sufficiently strong.
In order to get an insight into the essence of the proposed method, we use the equality lemma to compare with the geometry-based un-combined TCAR (GBUC) [25], as shown by Appendix B. Although the GBUC is theoretically comparable with the proposed method, the effectiveness of GBUC is limited by the model robustness of ionospheric bias.

Experiments Setup
In order to test the proposed TCAR method, a series of actual triple-frequency BDS data were collected.The baseline distances of the selected stations vary from 7.7 km to 68.8 km.Table 2 gives a summary of the data including the baseline distance, the time of data collection and the sampling rate.The distribution of these stations is shown in Figure 1.The number of common satellites for four baselines is also given in Figure 2. To evaluate the ambiguity resolution performance later, we correctly resolve the NL ambiguities among three frequencies beforehand using the whole span of data.The stations of CUT0, PERT and JFNG belong to the international GNSS service (IGS) network, whose locations are being precisely known.In contrast, the precisely locations of WUDH, WHHP, and WHHN cannot be publicly obtained.We use the commercial post-processing software (NovAtel GrafNav) to obtain the coordinate of these stations, i.e., the precise point positioning (PPP) mode using the IGS precise products.Moreover, because the WL and NL model from the proposed TCAR method contain the non-dispersive terms, the tropospheric bias needs to be taken into account.The hydrostatic part of tropospheric delay contained in the un-differenced observations is compensated by the Saastamoinen model [26].The wet part of troposphere delay is estimated by zenith wet parameter along with Niell mapping function [27].The parameter is assumed to be a random-walk process [6,28].The corresponding state transition between two consecutive epochs is τ k = τ k−1 + ω, in which ω denotes the process noise with the corresponding variance σ ω 2 = q τ ∆t, where q τ and ∆t is the spectrum density coefficients and the sampling internal, respectively.We take q τ = 3 cm 2 /h for tropospheric estimated model [6].In addition, the cutoff elevation angle is set to 15 • .The elevation-dependent weighting is applied, and the observation STD at elevation θ is σ(θ) = σ o (1 + 1/sin(θ))with σ o being the STD in zenith, which takes 30 cm and 3 mm for triple-frequency un-differenced code and phase observation, respectively.Moreover, in order to further improve estimated precision and efficiency, we use the sequential method of the extended Kalman filter as introduced in [6].tropospheric estimated model [6].In addition, the cutoff elevation angle is set to 15°.The elevationdependent weighting is applied, and the observation STD at elevation θ is σ(θ) = σo(1 + 1/sin(θ))with σo being the STD in zenith, which takes 30 cm and 3 mm for triple-frequency un-differenced code and phase observation, respectively.Moreover, in order to further improve estimated precision and efficiency, we use the sequential method of the extended Kalman filter as introduced in [6].

Error Analysis of the Combined Observation
In order to validate the performance of the proposed TCAR method with respect to different ionospheric cases, the residual double-differenced ionospheric delay is extracted using ambiguityfixed phase observations from the proposed IFVR method.tropospheric estimated model [6].In addition, the cutoff elevation angle is set to 15°.The elevationdependent weighting is applied, and the observation STD at elevation θ is σ(θ) = σo(1 + 1/sin(θ))with σo being the STD in zenith, which takes 30 cm and 3 mm for triple-frequency un-differenced code and phase observation, respectively.Moreover, in order to further improve estimated precision and efficiency, we use the sequential method of the extended Kalman filter as introduced in [6].

Error Analysis of the Combined Observation
In order to validate the performance of the proposed TCAR method with respect to different ionospheric cases, the residual double-differenced ionospheric delay is extracted using ambiguityfixed phase observations from the proposed IFVR method.

Error Analysis of the Combined Observation
In order to validate the performance of the proposed TCAR method with respect to different ionospheric cases, the residual double-differenced ionospheric delay is extracted using ambiguity-fixed phase observations from the proposed IFVR method.
According to the relationship of ionospheric delay between frequencies, I 1 , the first-order double-differenced ionospheric delay on frequency f 1 , can be calculated using ambiguity-corrected observation φ (1,0,0) and φ (0,1,0) .It can be expressed as Figure 3 shows the calculated results of the ionospheric delay on frequency f 1 , in which the different color codes represent the ionospheric delay from each satellite-pairs.It can be found from Figure 3 that, the peaks of the ionospheric delay for the baseline distance of 7.7 km, 22.4 km, 42.5 km, and 68.8 km can reach up to 11.6 cm, 18.1 cm, 58.0 cm, and 72.4 cm, respectively.It should be pointed out that the ionospheric delay for 68.8 km can be up to 3.8 cycles with respect to the wavelength of f 1 .It is because the ionospheric delay does not only depend on the baseline distance but also the ionospheric condition.The peak of ionospheric delay fluctuations for the 68.8 km is at the local time from 18:00 to 22:00, which correspond to the local sunset and midnight.This relatively large ionospheric delay is probably related to the ionospheric disturbance, i.e., the so-called traveling ionospheric disturbances (TID) [29].Because the new ionosphere-free combinations are used in the proposed TCAR method, the ambiguity resolution will not be affected by the ionospheric bias.
Remote Sens. 2017, 9, 1108 8 of 15 According to the relationship of ionospheric delay between frequencies, I1, the first-order double-differenced ionospheric delay on frequency f1, can be calculated using ambiguity-corrected observation φ  (1,0,0) and φ Figure 3 shows the calculated results of the ionospheric delay on frequency f1, in which the different color codes represent the ionospheric delay from each satellite-pairs.It can be found from Figure 3 that, the peaks of the ionospheric delay for the baseline distance of 7.7 km, 22.4 km, 42.5 km, and 68.8 km can reach up to 11.6 cm, 18.1 cm, 58.0 cm, and 72.4 cm, respectively.It should be pointed out that the ionospheric delay for 68.8 km can be up to 3.8 cycles with respect to the wavelength of f1.It is because the ionospheric delay does not only depend on the baseline distance but also the ionospheric condition.The peak of ionospheric delay fluctuations for the 68.8 km is at the local time from 18:00 to 22:00, which correspond to the local sunset and midnight.This relatively large ionospheric delay is probably related to the ionospheric disturbance, i.e., the so-called traveling ionospheric disturbances (TID) [29].Because the new ionosphere-free combinations are used in the proposed TCAR method, the ambiguity resolution will not be affected by the ionospheric bias.In addition to the effect of ionospheric bias, the combined observation noise is another factor impacting the reliability of ambiguity resolution.Because the residuals of observation reflect the noise level of observation, we therefore use Figure 4 to show different combined observation residuals with respect to each visible satellites-pair for the baseline of 68.8 km.It can be seen that the majority of observation residuals can be restricted to within 0.25 cycles and the mean is close to 0 cycles.This demonstrates that the proposed TCAR method is not only capable of eliminating the effect of atmospheric bias, but also can restrict the noise level, which has a positive effect on the fast and reliable ambiguity resolution.The statistical results of residuals for all baselines are listed in Table 3.It can be seen that the maximum of STD of observation residuals are less than 0.15 cycles.Moreover, it can be found that the residual level of φWL2 and φEWL are larger than that of the others.This is because these two combined observations are constructed by using code observation with low precision.In addition to the effect of ionospheric bias, the combined observation noise is another factor impacting the reliability of ambiguity resolution.Because the residuals of observation reflect the noise level of observation, we therefore use Figure 4 to show different combined observation residuals with respect to each visible satellites-pair for the baseline of 68.8 km.It can be seen that the majority of observation residuals can be restricted to within 0.25 cycles and the mean is close to 0 cycles.This demonstrates that the proposed TCAR method is not only capable of eliminating the effect of atmospheric bias, but also can restrict the noise level, which has a positive effect on the fast and reliable ambiguity resolution.The statistical results of residuals for all baselines are listed in Table 3.It can be seen that the maximum of STD of observation residuals are less than 0.15 cycles.Moreover, it can be found that the residual level of φ WL2 and φ EWL are larger than that of the others.This is because these two combined observations are constructed by using code observation with low precision.However, the STD of the residuals of these two combinations can still be restricted within 0.15 cycles, which means the reliable ambiguity resolution can be achieved.
Remote Sens. 2017, 9, 1108 9 of 15 However, the STD of the residuals of these two combinations can still be restricted within 0.15 cycles, which means the reliable ambiguity resolution can be achieved.
Table 3. Statistics of residuals for each combined observation from all of data (STD: cycles).
The experimental results are evaluated by the ambiguity resolution performance and positioning accuracy.The ambiguity resolution performance of EWL, WL and NL steps are shown by the metrics of the empirical fixed rate (P fix ) and the conditional success rate (P cf ), in which P fix are defined as the ratio between the epochs of the ambiguity being fixed and total epochs, and P cf is defined as the ratio between the epochs of the ambiguity being fixed correctly and the epochs of the ambiguity being fixed, in which the ambiguity can been accepted as fixed only if the ratio-test is passed.Herein, P fix reflects the continuity of ambiguity resolution, while P cf reflects the correctness of ambiguity resolution.As discussed in [21,31], the empirical fixed rate and conditional success rate reflects the reliability of ambiguity resolution, which are listed in Table 4.The positioning errors for the three methods in the ambiguity-fixed mode are shown in Figure 5 and the corresponding STD values of positioning errors are listed in Table 5.The experimental results are evaluated by the ambiguity resolution performance and positioning accuracy.The ambiguity resolution performance of EWL, WL and NL steps are shown by the metrics of the empirical fixed rate (Pfix) and the conditional success rate (Pcf), in which Pfix are defined as the ratio between the epochs of the ambiguity being fixed and total epochs, and Pcf is defined as the ratio between the epochs of the ambiguity being fixed correctly and the epochs of the ambiguity being fixed, in which the ambiguity can been accepted as fixed only if the ratio-test is passed.Herein, Pfix reflects the continuity of ambiguity resolution, while Pcf reflects the correctness of ambiguity resolution.As discussed in [21,31], the empirical fixed rate and conditional success rate reflects the reliability of ambiguity resolution, which are listed in Table 4.The positioning errors for the three methods in the ambiguity-fixed mode are shown in Figure 5 and the corresponding STD values of positioning errors are listed in Table 5.As seen from Table 4, the empirical and correctly fixed rates of EWL and WL ambiguity resolution for three methods achieve 100% regardless of the baseline distances, which reveals the effectiveness of long wavelength by observation combination.By comparing the NL ambiguity resolution result between the GBIR and the IF-based TCAR methods, i.e., the GFIF and the proposed IFVR, the necessity of the ionosphere-free combination can be demonstrated because the continuity and the reliability of GBIR is deteriorated when a larger ionospheric bias results from the increased baselines.Particularly, when the baseline increases to 68.8 km, the empirical fixed rate of GBIR is 22.5%, and the correctly fixed rate is 30.5%,respectively.Compared with the GFIF and the IFVR, it can be found that the empirical fixed rate performances are comparable to each other, which can be also be reflected by the comparable positioning accuracy performance as listed in Table 5.Note that the position error of first column is shifted from its reference value in the Figure 5.This is probably the bias from PPP-based positioning reference.However, the ambiguity resolution correctness performance of IFVR can achieve an improvement by 5.7% at most in the 68.8 km baseline distance case, which means the proposed IFVR has stronger ambiguity resolution strength than the GFIF.

Concluding Remarks
The purpose of this study is to improve the reliability of the TCAR method by eliminating ionospheric bias and restricting the combined observation noise simultaneously.We proposed a new ionosphere-free and combined observation variance-restricted TCAR method by utilizing ambiguity links between each step of TCAR.Since the advantage of the ambiguity link is explored sufficiently, the contradiction between the ionosphere elimination and combined observation noise suppression can be removed, so that the new ionosphere-free combination method with combined observation noise restriction can be constructed especially for WL and NL ambiguity resolution.Moreover, the strength of ambiguity resolution can also be improved by the ambiguity link, which is beneficial for improving the reliability of ambiguity resolution.The proposed TCAR method is tested using real data from various baselines.As a result, the STD of residuals of combination from each step of TCAR can be restricted within 0.15 cycles, even when the baseline is increased to 68.8 km.By comparing the NL ambiguity resolution performance of GBIR with the GFIF and the IFVR, it is found that the elimination of the effect of ionospheric bias is essential to improve the ambiguity resolution, particularly in the NL step.Meanwhile, the continuity performance is significantly improved by effectively restricting the observational noise level with respect to the wavelength, as indicated by the comparison with the GFIF and the IFVR.
According to the equality lemma, the CA 2 can be derived as Therefore, the condition of CA 2 = 0 and rank (C) = n − t is fulfilled between the GBUC and the IFVR method, which means the method strength of these two combinations are equivalent.However, the effectiveness of the GBUC depends on accuracy of ionospheric model, particularly, ionospheric stochastic model.In contrast, the IFVR method can be used without regard to inaccuracy of ionospheric model, which means the proposed IFVR have better ambiguity resolution performance in the presence of ionosphere.

Figure 1 .
Figure 1.Distribution of the selected stations in our experiment.

Figure 1 .
Figure 1.Distribution of the selected stations in our experiment.

Figure 1 .
Figure 1.Distribution of the selected stations in our experiment.

Figure 4 .
Figure 4. Residuals of each combined observations from the distance of 68.8 km.The panels from top to bottom represent residuals of φEWL, φWL2, φWL1, φNL2 and φNL1, respectively.The gray dash line indicates the mean value.

Figure 4 .
Figure 4. Residuals of each combined observations from the distance of 68.8 km.The panels from top to bottom represent residuals of φ EWL , φ WL2 , φ WL1 , φ NL2 and φ NL1 , respectively.The gray dash line indicates the mean value.

Figure 5 .
Figure 5. Positioning error comparison for different methods in the ambiguity-fixed mode.The row panels from top to bottom denote GBIR, GFIF, and IFVR, respectively.The column panels from left to right denote the distances of 7.7 km, 22.4 km, 42.5 km, and 68.8 km.The colors of purple, green and red represent the positioning errors in the east, north, and up components, respectively.

Figure 5 .
Figure 5. Positioning error comparison for different methods in the ambiguity-fixed mode.The row panels from top to bottom denote GBIR, GFIF, and IFVR, respectively.The column panels from left to right denote the distances of 7.7 km, 22.4 km, 42.5 km, and 68.8 km.The colors of purple, green and red represent the positioning errors in the east, north, and up components, respectively.

Table 3 .
Statistics of residuals for each combined observation from all of data (STD: cycles).

Table 4 .
Ambiguity validation performances for different methods (×100%).value for the ratio test is set to 3. The epoch-by-epoch ambiguity resolution processing is used in order to investigate the performance of ambiguity resolution ultimately. threshold