Improved Multi-GNSS PPP Partial Ambiguity Resolution Method Based on Two-Step Sorting Criterion

: Multi-GNSS PPP partial ambiguity resolution (PAR) can improve the ﬁxing success rate and shorten the time to ﬁrst ﬁx (TTFF). Ambiguity subset selection based on the bootstrapping success rate sorting criterion (BSSC) is widely used in PPP PAR due to its ease of computation and comprehensive evaluation of the global quality of ambiguity solutions. However, due to the inﬂuence of unmodeled errors, such as atmospheric residuals and gross errors, ambiguity parameter estimation will inevitably introduce bias. For ambiguity parameters with bias, their variance will converge incorrectly and will not accurately reﬂect the estimation accuracy. As a result, the selected ambiguity subset based on the BSSC becomes inaccurate, affecting the ﬁxing success rate and TTFF. Therefore, we proposed an improved multi-GNSS PPP PAR method based on a two-step sorting criterion (TSSC). This method aims to address the inﬂuence of inaccurate variance of ambiguity parameters, particularly those with low observation quality, on the ambiguity subset selection based on the BSSC. The ambiguity subset satisfying the preset success rate threshold is selected to reduce the inﬂuence of unconverged ambiguity on the TSSC. In the ﬁrst step of the sorting process, the observations whose elevation angle is below 30 ◦ or whose posterior residual falls into the IGG3 model reduction domain are clustered together. The posterior observation weight criterion (POWC) instead of the BSSC is adopted to sort ambiguities to overcome the false convergence of variance of ambiguity parameters. In the second step of the sorting process, the remaining ambiguities with reasonable variances are sorted based on the BSSC. Finally, the bottom ambiguity is removed one by one from the ambiguity subset sorted based on the two-step sorting criterion (TSSC) until the requirements of the ratio test for LAMBDA are met. The static data from 10 MGEX stations over a period of 30 days, along with urban kinematic data, were collected to validate the proposed method. Compared with the PAR based on the BSSC, the static experiments demonstrated a reduction of 8.7% and 16.8% in the TTFF and convergence time, respectively. Additionally, the positioning accuracy in the east, north, and up directions was improved by 20.1%, 17.1%, and 4.67%, respectively. Furthermore, the kinematic experiment revealed that the TTFF and convergence time decreased from 1.65 min and 10.5 min to 1.3 min and 1.8 min, respectively, with higher positioning accuracy.


Introduction
Precise Point Positioning (PPP) is a fundamental technology in the field of global navigation satellite systems (GNSS) precise positioning, offering the advantage of achieving sub-decimeter or even centimeter-level accuracy with a single receiver [1].However, the long convergence time of float-based PPP affects real-time applications.In comparison, carrier phase integer ambiguity resolution (AR) can reduce convergence time and enhance positioning accuracy [2][3][4].However, the full ambiguity resolution (FAR) of PPP typically takes tens of min [5][6][7].With the development of the Global Positioning System (GPS), GLObal NAvigation Satellite System (GLONASS), BeiDou Navigation Satellite System (BDS), and Galileo, multi-GNSS PPP can improve the convergence time and positioning accuracy due to better spatial geometry and more observations [8][9][10][11].However, multi-GNSS PPP FAR always takes a longer time to first fix (TTFF) and a lower fixing success rate due to the frequent occurrence of low-quality observation satellites [5].
Partial ambiguity resolution (PAR) can enhance the AR success rate and reduce the TTFF [12].The selection criterion for the ambiguity subset is crucial to reducing the TTFF and improving the fixing success rate of PAR.In order to improve the fixing success rate of AR, Teunissen et al. proposed an ambiguity subset selection method based on the success rate criterion (SRC) [13].The ambiguity subset is determined when its success rate satisfies the preset success rate threshold (PSRT).Although the SRC can provide a probability of successfully fixing the ambiguity, its fixed solution reliability is difficult to guarantee due to low-quality observation satellites.Mowlam et al. proposed an advanced SRC ambiguity subset selection method combined with an elevation angle [14].This method iteratively excludes satellites with the lowest elevation angle to make the ambiguity subset meet the PSRT.Additionally, the ratio test is a commonly employed method for assessing the correctness of ambiguity estimation.To further enhance the reliability of the fixed solution, the ratio test is typically conducted following the LAMBDA search [15][16][17][18].Wang et al. proposed a method for selecting ambiguity subsets based on the bootstrapping success rate sorting criterion (BSSC).In this method, the PSRT and ratio test were used to verify the fixed solution to reduce the probability of accepting incorrect fixed solutions [19].However, the fixing success rate will be limited because the ratio test is only applied to the maximum ambiguity subset, which meets the PSRT, whereas the remaining subsets, which may satisfy the ratio test within the maximum ambiguity subset, are not validated.Li et al. improved the method by removing ambiguities with the lowest bootstrapping success rate until the requirements of both the PSRT and ratio test were met for the subsubset, which improved the fixing success rate [12].Therefore, the ambiguity subset selection based on the BSSC is widely used in PPP PAR due to its simple calculation and comprehensive evaluation of the global quality of ambiguity solutions.However, the success rate calculation of bootstrapping assumes that the ambiguities are unbiased-biases consisting of unmodeled errors with low elevation angle atmospheric residuals and gross errors in harsh environments are introduced to the ambiguity parameters.This leads to false convergence for the variance of the biased ambiguity parameters [20,21].Due to the presence of false convergence in the biased ambiguity parameters, the reliability of the ambiguity subset selection method based on the BSSC is difficult to guarantee, and the performance of multi-GNSS PPP PAR will degrade.
It is well known that unmodeled errors, such as atmospheric residuals and multipath effects, have a greater impact on satellite observations when the elevation angle of the satellite decreases.Currently, most studies divide the stochastic model of observations into two segments with a 30 • elevation angle boundary [20,22,23].It is generally believed that satellites below a 30 • elevation angle are significantly affected by atmospheric residuals and multipath effects.When the satellite elevation angle is below 30 • , the weight of observations will decrease.Otherwise, the weight values of observations are equal.The posterior observation weight based on the IGG3 model reflects the accuracy of observations more reasonably than the prior observation weight based on the elevation angle stochastic model [24].Although the robust Kalman filter based on the IGG3 model can mitigate the impact of unmodeled errors to some extent, the false convergence for the variance of the biased ambiguity parameters cannot be completely excluded [25,26].The Schmidt-Kalman filter (SKF) is the simplest filtering algorithm for stochastic systems which takes into account the uncertainty statistics of parameters but does not update them [27,28].While the SKF is well known and has been used in covariance analysis, not much attention has been given to its actual implementations in a real-time recursive estimation algorithm.
To mitigate the impact of inaccurate variance of biased ambiguity parameters caused by low-quality observations on the BSSC ambiguity subset selection, we proposed an improved multi-GNSS PPP PAR method based on the two-step sorting criterion (TSSC).In this contribution, observations under the following two circumstances will be clustered together: the first is observations with an elevation angle below 30 • , and the second is the ones whose posterior residual falls into the IGG3 model reduction domain.The posterior observation weight criterion (POWC) instead of the BSSC is adopted to sort ambiguities to overcome the false convergence of variance, and the remaining ambiguities with reasonable variances are sorted based on the BSSC.Finally, a subset of the sorted ambiguities based on TSSC is obtained, and all ambiguities are gradually excluded from the ambiguity subset until the requirements of the ratio test for LAMBDA are met.
The chapters of this paper are arranged as follows.In Section 2, we give a brief introduction to the multi-GNSS PPP AR model.In Section 3, we introduce PPP PAR based on the TSSC process.In Section 4, we validate the effectiveness of the proposed method using static data collected from 10 MGEX stations over a 30-day period, along with urban kinematic data.Finally, we present the discussion and conclusion of our study.

Multi-GNSS PPP AR Model
For a user receiver r tracking the satellite s on frequency j, the functional model of code and carrier phase observation equations can be formulated as: The symbols P and L denote code and carrier phase observations in meters, respectively, and their noises are ε P and ε L in meters, respectively; sys denotes constellation system; ρ is the geometric distance between receiver antenna phase center and satellite antenna phase center in meters; Z w is the wet part of the tropospheric delay in the zenith direction in meters; M w is the mapping function of Z w ; c is the speed of light in meters per second; dt and dt s are, respectively, the receiver clock and satellite clock errors in seconds; I is the slant ionospheric delay in meters; d r and d s are the receiver-side and the satellite-side code hardware delay in meters; λ j is the wavelength of frequency j in meters/cycles; N is the undifferenced integer ambiguity in cycles; b r is the receiver-side carrier phase hardware delay in cycles; b s is the satellite-side carrier phase hardware delay in cycles.In brief, the errors including phase wind-up, dry tropospheric delay, tidal displacements, earth rotation, and relativistic effects can be corrected using existing models [29].
Generally, carrier phase hardware delay bias of the satellite-and receiver-side will be merged into the undifferenced ambiguity.In order to recover the integer characteristic of undifferenced carrier phase ambiguity, the receiver-side hardware delay bias is usually eliminated by the between-satellite single difference (BSSD).The satellite-side hardware delays are generally eliminated by external correction information.In order to meet the needs of PPP AR users under various observation combination modes, IGS proposed observable-specific signal biases (OSBs) based on the SINEX format, which can directly correct the original observations [30,31] where the subscripts 1 and 2 represent the corresponding BDS B1I and B3I frequencies, GPS L1 and L2 frequencies, and Galileo E1 and E5a frequencies, respectively; where N s,sys r,wl denotes wide-lane integer ambiguity in cycles; Ns,sys r,nl denotes narrow-lane float ambiguity.
The float wide-lane ambiguity is obtained by Hatch-Melbourne-Wübbena (HMW) combination [32][33][34]: where N s,sys r,wl is the float wide-lane ambiguity in cycles; is the wide-line ambiguity wavelength; N s,sys r,wl is the integer wide-lane ambiguity in cycles; µ sys r,12 is the HMW combination code and carrier phase hardware biases at receiver in cycles.
In order to recover the integer characteristic of wide-lane ambiguity, the receiver-side hardware delay bias is usually eliminated by the BSSD.The nearest rounding decision function is used to fix the BSSD wide-lane ambiguity [20].When the BSSD wide-lane ambiguity is successfully fixed, the BSSD narrow-lane ambiguity with integer characteristics can be obtained by Equation (4).The LAMBDA method satisfying the ratio test is used to fix the BSSD narrow-lane ambiguity.

Posterior Observation Weight Criterion
The IGG3 model can suppress abnormal observations by constructing equivalent weight factors.The equivalent weight function is as follows [35]: where p i is the weight reduction factor; v i is the standardized posterior residual of the code or carrier phase; k 0 and k 1 are constants.In this paper, k 0 takes 1.5 and k 1 takes 3.0.In the IGG3 model, once the satellite code or carrier phase standardized posterior residuals fall into the rejection domain, they will be excluded.If multiple satellites fall into the rejection region at the same time, only the observation with the largest variance is excluded [36].In order to ensure computational efficiency, the maximum number of iterations is set to 6.The posterior observation weight matrix based on the IGG3 model can be written as follows: where P L k is the prior observation weight matrix using elevation angle stochastic model in dimension n × n; P L k is the posterior observation weight matrix based on the IGG3 model in dimension n × n.
Due to the influence of unmodeled errors, such as atmospheric residuals and gross errors, the variances of ambiguity parameters exist in false convergence, which cannot accurately reflect its estimation accuracy.The POWC is sorted from small to large according to the posterior weight of the carrier phase observation.To avoid the adverse effects of ambiguity variance inaccuracy on the ambiguity subset selection when using the BSSC, the ambiguity subset based on POWC can be determined as follows: (1) The ambiguity subset âm that meets the PSRT to reduce the influence of unconverged ambiguities is selected.
(2) The ambiguities whose observation elevation angle is below 30 • or whose posterior residual falls into the IGG3 model reduction domain are clustered together, denoted as â1 .The POWC is adopted to sort ambiguities to overcome the false convergence of variance of ambiguity parameters.

Bootstrapping Success Rate Sorting Criterion
Due to the lower bound of the success rate of integer least squares estimation, the success rate of bootstrapping has been proven to be a solution to the success rate of integer least squares with a very high degree of approximation.The success rate of bootstrapping can be calculated as follows [13]: where Ni|I is the abbreviation of Ni|i−1,••• ,1 , which is the conditional estimation of the ith ambiguity with the previous i−1 ambiguity fixed to the integer; σ is the standard deviation of Ni|I ; Φ(x) is a normal distribution function.
The BSSC method sorts the ambiguity from small to large according to the bootstrapping success rate of the ambiguity; then, partial ambiguity screening is realized.The remaining subset of ambiguities with reasonable variances except â1 in âm are denoted as â2 .The â2 are sorted based on the BSSC.
Combined with the POWC and BSSC, we proposed an improved multi-GNSS PPP PAR method based on TSSC, abbreviated as TSSC-PAR.Figure 1 shows the processing of the TSSC-PAR.The basic process is as follows: 1.
Stage 1.The ambiguity subset satisfying the PSRT is selected to reduce the influence of unconverged ambiguity on the TSSC.2.
Stage 2. In the first step of sorting based on TSSC, the POWC method is used to sort the ambiguities in the ambiguities subset â1 , and the sorted ambiguities subset is denoted as â1 .This process helps avoid selecting inaccurate variance of the ambiguity based on the BSSC method.3.
Stage 3. In the second step of sorting based on TSSC, the BSSC method is used to sort the ambiguities in the ambiguities subset â2 , and the sorted ambiguities subset is denoted as â2 .â1 and â2 are combined to form âm = â2 â1 T .

4.
Stage 4. Finally, the bottom ambiguity from âm is one by one excluded until either the ratio test is met or the number of participating fixed ambiguities falls below 5.
Otherwise, the float solution is output.
the ambiguities in the ambiguities subset ˆ2

Experiments and Analysis
To demonstrate the effectiveness of the multi-GNSS PPP PAR method based on TSSC, the static observation data from 10 MGEX stations during DOY 116-DOY 145 in 2022, and the urban kinematic data on DOY 258, 2021 were collected.The data were compared with the PAR method based on the BSSC, abbreviated as BSSC-PAR.The BSSC-PAR accepts fixed solution evaluation criteria consistent with TSSC-PAR.
Table 1 outlines the processing strategy for the static and kinematic data.The positioning accuracy of the TSSC-PAR and BSSC-PAR shown in this paper presents fixed solutions.BDS's Geostationary Earth Orbit (GEO) satellites are excluded from this process.This is due to the poor accuracy of their observations and orbit products.The performance of the multi-GNSS PPP PAR is evaluated, including the TTFF, positioning accuracy, and convergence time.The "true" coordinates of the stations are obtained from the IGS weekly solution file and are used as the reference for assessing positioning accuracy.The TTFF in the ambiguity-fixed solution is defined as the time to get the first correctly fixed solution, which means the ambiguity validation succeeds and positioning errors of the fixed solution are smaller than that of the float solution in the same epoch.The positioning errors are defined as the differences between positioning solutions and reference solutions.The positioning accuracy is the RMS of converged positioning errors.Convergence time is defined as the time required for a solution to achieve a horizontal positioning error of less than 10 cm and remains for at least 10 epochs.

Experiments and Analysis
To demonstrate the effectiveness of the multi-GNSS PPP PAR method based on TSSC, the static observation data from 10 MGEX stations during DOY 116-DOY 145 in 2022, and the urban kinematic data on DOY 258, 2021 were collected.The data were compared with the PAR method based on the BSSC, abbreviated as BSSC-PAR.The BSSC-PAR accepts fixed solution evaluation criteria consistent with TSSC-PAR.
Table 1 outlines the processing strategy for the static and kinematic data.The positioning accuracy of the TSSC-PAR and BSSC-PAR shown in this paper presents fixed solutions.BDS's Geostationary Earth Orbit (GEO) satellites are excluded from this process.This is due to the poor accuracy of their observations and orbit products.The performance of the multi-GNSS PPP PAR is evaluated, including the TTFF, positioning accuracy, and convergence time.The "true" coordinates of the stations are obtained from the IGS weekly solution file and are used as the reference for assessing positioning accuracy.The TTFF in the ambiguity-fixed solution is defined as the time to get the first correctly fixed solution, which means the ambiguity validation succeeds and positioning errors of the fixed solution are smaller than that of the float solution in the same epoch.The positioning errors are defined as the differences between positioning solutions and reference solutions.The positioning accuracy is the RMS of converged positioning errors.Convergence time is defined as the time required for a solution to achieve a horizontal positioning error of less than 10 cm and remains for at least 10 epochs.

The Static Experiment
The distribution of the stations is presented in Figure 2. Due to differences in the satellite visibility from site to site, the selection of stations across the globe allows us to get a clear insight into the expected performance that is globally applicable and not locationrestricted.The sampling frequency of static experimental data collected by the MGEX stations is 30 s.The static experimental data were processed using the kinematic model.To simulate PPP PAR hourly solutions, restarts were performed every hour.

The Static Experiment
The distribution of the stations is presented in Figure 2. Due to differences in th satellite visibility from site to site, the selection of stations across the globe allows us to ge a clear insight into the expected performance that is globally applicable and not location restricted.The sampling frequency of static experimental data collected by the MGEX sta tions is 30 s.The static experimental data were processed using the kinematic model.T simulate PPP PAR hourly solutions, restarts were performed every hour.

The AMC4 Station Experiment
Figure 3 shows the positioning accuracy of the BSSC-PAR and TSSC-PAR for the AMC4 station in the east, north, and up directions.It can be clearly observed that the positioning accuracy of TSSC-PAR is significantly improved in the east, north, and up directions, compared with the BSSC-APR.As shown in Table 2, the TTFF and convergence time of TSSC-PAR are 4.3 min and 7.7 min, respectively, which is an improvement of 28.3% and 18.9%, respectively, compared with the BSSC-PAR.Furthermore, the positioning error RMS of TSSC-PAR are 13.1 mm 11.5 mm, and 42.5 mm in the east, north, and up directions, respectively, with an improvement of about 23.4%, 5.7%, and 8.2% compared with the BSSC-PAR.Figure 4 shows the number of ambiguity-fixed for the AMC4 station on the 116th day of 2022.The inaccurate variance of the ambiguity parameter due to the lower quality of the observations affects the accuracy of the ambiguity selection by the BSSC, while the TSSC-PAR uses the POWC to select a subset of ambiguities with inaccurate variances.As shown in Figure 4, during the convergence phase of the PPP, the number of ambiguity-fixed for the TSSC-PAR is higher than that for the BSSC-PAR.Combined with Figure 3 and Table 2, TSSC-PAR shows a more significant improvement in positioning performance during the convergence phase of PPP compared with the BSSC-PAR.It means that the ambiguity subset selected by TSSC-PAR is more accurate than the one selected by BSSC-PAR.Furthermore, compared with BSSC-PAR, the average number of ambiguity-fixed for TSSC-PAR increased from 14.4 to 15.1.
with the BSSC-PAR.Figure 4 shows the number of ambiguity-fixed for the AMC4 station on the 116th day of 2022.The inaccurate variance of the ambiguity parameter due to the lower quality of the observations affects the accuracy of the ambiguity selection by the BSSC, while the TSSC-PAR uses the POWC to select a subset of ambiguities with inaccurate variances.As shown in Figure 4, during the convergence phase of the PPP, the number of ambiguity-fixed for the TSSC-PAR is higher than that for the BSSC-PAR.Combined with Figure 3 and Table 2, TSSC-PAR shows a more significant improvement in positioning performance during the convergence phase of PPP compared with the BSSC-PAR.It means that the ambiguity subset selected by TSSC-PAR is more accurate than the one selected by BSSC-PAR.Furthermore, compared with BSSC-PAR, the average number of ambiguity-fixed for TSSC-PAR increased from 14.4 to 15.1.To further investigate the performance of the BSSC-PAR and TSSC-PAR at the AMC4 station, the data from DOY 116-DOY 145 in 2022 were analyzed.Figure 5 shows the convergence series of the BSSC-PAR and TSSC-PAR for 68% quantile in the east, north, and To further investigate the performance of the BSSC-PAR and TSSC-PAR at the AMC4 station, the data from DOY 116-DOY 145 in 2022 were analyzed.Figure 5 shows the convergence series of the BSSC-PAR and TSSC-PAR for 68% quantile in the east, north, and up directions.From Figure 5, we can see that the convergence series for the TSSC-PAR is better than that for the BSSC-PAR.Additionally, compared with the BSSC-PAR, the east convergence trend of TSSC-PAR is significantly improved.As is shown in Figure 6, the cumulative frequency of the TTFF and convergence time for the TSSC-PAR are better than the BSSC-PAR.As shown in Table 3, the RMS of positioning error for TSSC-PAR are 15.1 mm, 13.9 mm, and 47.5 mm in the east, north, and up directions, respectively.There is an improvement of 25.2%, 22.3%, and 10.8%, respectively, compared with the BSSC-PAR.Furthermore, the average TTFF and convergence time of the TSSC-PAR are 4.8 min and 7.3 min, respectively.Compared with the BSSC-PAR, the TSSC-PAR has an improvement of 11.1% and 21.5%, respectively.Figure 7 shows the distribution of the ambiguities-fixed number for BSSC-PAR and TSSC-PAR.From Figure 7, we can see that the distribution of the ambiguity-fixed number for TSSC-PAR is better than that for BSSC-PAR.Combined with Figures 5 and 6 and Table 3, with the increase of the ambiguity-fixed number, the positioning performance of the TSSC-PAR is better than the BSSC-PAR.Based on the multiday data in the AMC4 station, it can be concluded that the TTFF, convergence time, and positioning accuracy of the TSSC-PAR are superior to that of the BSSC-PAR.In order to assess the global applicability of the TSSC-PAR performance, static data from 10 MGEX stations were analyzed.Figure 8 shows the RMS of positioning error for BSSC-PAR and TSSC-PAR in the east, north, and up directions.The positioning accuracy for the TSSC-PAR is better than that for the BSSC-PAR at each station as shown in Figure 8.In Figure 9, the TTFF and convergence time for the TSSC-PAR are better than that for the BSSC-PAR.Table 4 gives the statistic indicators of the BSSC-PAR and TSSC-PAR.Compared with the BSSC-PAR, the positioning accuracy of the TSSC-PAR has improved by 20.1%, 17.1%, and 4.67%, and it has eventually reached 20.7 mm, 14.8 mm, and 58.8 mm in the east, north, and up directions, respectively.Furthermore, compared with the BSSC-PAR, the average TTFF and convergence time for the TSSC-PAR are 3.15 min and 5.2 min, which is an improvement of 8.7% and 16.8%, respectively.As shown in Figure 10, the distribution of the ambiguity-fixed number for the TSSC-PAR is better than that for the BSSC-PAR.The average number of the ambiguity-fixed solutions for TSSC-PAR and BSSC-PAR are 17.0 and 16.1, respectively.Therefore, compared with the BSSC-PAR, the ambiguity subset selected by the TSSC-PAR is more accurate.To sum up, the positioning performance of TSSC-PAR is still better than that of BSSC-PAR under the verification of  In order to assess the global applicability of the TSSC-PAR performance, static data from 10 MGEX stations were analyzed.Figure 8 shows the RMS of positioning error for BSSC-PAR and TSSC-PAR in the east, north, and up directions.The positioning accuracy for the TSSC-PAR is better than that for the BSSC-PAR at each station as shown in Figure 8.In Figure 9, the TTFF and convergence time for the TSSC-PAR are better than that for the BSSC-PAR.Table 4 gives the statistic indicators of the BSSC-PAR and TSSC-PAR.Compared with the BSSC-PAR, the positioning accuracy of the TSSC-PAR has improved by 20.1%, 17.1%, and 4.67%, and it has eventually reached 20.7 mm, 14.8 mm, and 58.8 mm in the east, north, and up directions, respectively.Furthermore, compared with the BSSC-PAR, the average TTFF and convergence time for the TSSC-PAR are 3.15 min and 5.2 min, which is an improvement of 8.7% and 16.8%, respectively.As shown in Figure 10, the distribution of the ambiguity-fixed number for the TSSC-PAR is better than that for the BSSC-PAR.The average number of the ambiguity-fixed solutions for TSSC-PAR and BSSC-PAR are 17.0 and 16.1, respectively.Therefore, compared with the BSSC-PAR, the ambiguity subset selected by the TSSC-PAR is more accurate.To sum up, the positioning performance of TSSC-PAR is still better than that of BSSC-PAR under the verification of large amounts of static data.

The Kinematic Experiment
Table 5 summarizes the basic information from kinematic experiments.Figure 11 shows the reference station position and vehicle trajectory, where the longest baseline between the reference station and the vehicle is 2 km.RTK can obtain centimeter-level or even millimeter-level positioning accuracy.Therefore, we use fixed solution RTK of BDS/GPS combination by RTKLIB software as the reference true value, which is used to evaluate the performance of the TSSC-PAR and BSSC-PAR.

The Kinematic Experiment
Table 5 summarizes the basic information from kinematic experiments.Figure 11 shows the reference station position and vehicle trajectory, where the longest baseline between the reference station and the vehicle is 2 km.RTK can obtain centimeter-level or even millimeter-level positioning accuracy.Therefore, we use fixed solution RTK of BDS/GPS combination by RTKLIB software as the reference true value, which is used to evaluate the performance of the TSSC-PAR and BSSC-PAR.Figure 12 shows the value of the PDOP and HDOP of the kinematic experiment.From Figure 12, we can see that the value of the PDOP and HDOP varies rapidly.Figure 13 shows that the positioning accuracy of the TSSC-PAR is significantly better than that of the BSSC-PAR, in the east, north, and up directions.Figure 14 shows the number of ambiguity-fixed solutions for the BSSC-PAR and TSSC-PAR.Compared with the BSSC-PAR, the average number of ambiguity-fixed solutions for the TSSC-PAR increased from 12.9 to 13.6.Combined with Figures 13 and 14, it can be seen that the number of ambiguity-fixed solutions by the TSSC-PAR is better than that by the BSSC-PAR in the PPP initialization phase.The positioning accuracy of the TSSC-PAR is better than that of the BSSC-PAR in the PPP initialization phase.Although there is little difference in the TTFF between the TSSC-PAR and BSSC-PAR, the convergence time of TSSC-PAR is significantly shorter than that of the BSSC-PAR.Table 6 gives the statistic indicators of the kinematic experiment for the BSSC-PAR and TSSC-PAR.The results are evaluated, including positioning accuracy, TTFF, and convergence time.Compared with BSSC-PAR, the positioning accuracy of TSSC-PAR improved by 38.1%, 1.4%, and 2.2% in the east, north, and up directions, respectively.Additionally, the TTFF and convergence time decreased from 1.65 min and 10.5 min to 1.3 min and 1.8 min, respectively.In conclusion, the experimental results have demonstrated that the TSSC-PAR has a better performance than the BSSC-PAR in kinematic scenes.Figure 12 shows the value of the PDOP and HDOP of the kinematic experiment.From Figure 12, we can see that the value of the PDOP and HDOP varies rapidly.Figure 13 show that the positioning accuracy of the TSSC-PAR is significantly better than that of the BSSC PAR, in the east, north, and up directions.Figure 14 shows the number of ambiguity-fixe solutions for the BSSC-PAR and TSSC-PAR.Compared with the BSSC-PAR, the averag number of ambiguity-fixed solutions for the TSSC-PAR increased from 12.9 to 13.6.Com bined with Figures 13 and 14, it can be seen that the number of ambiguity-fixed solution by the TSSC-PAR is better than that by the BSSC-PAR in the PPP initialization phase.Th positioning accuracy of the TSSC-PAR is better than that of the BSSC-PAR in the PPP initia ization phase.Although there is little difference in the TTFF between the TSSC-PAR an BSSC-PAR, the convergence time of TSSC-PAR is significantly shorter than that of the BSSC PAR.Table 6 gives the statistic indicators of the kinematic experiment for the BSSC-PAR an TSSC-PAR.The results are evaluated, including positioning accuracy, TTFF, and conver gence time.Compared with BSSC-PAR, the positioning accuracy of TSSC-PAR improved b 38.1%, 1.4%, and 2.2% in the east, north, and up directions, respectively.Additionally, th TTFF and convergence time decreased from 1.65 min and 10.5 min to 1.3 min and 1.8 min respectively.In conclusion, the experimental results have demonstrated that the TSSC-PAR has a better performance than the BSSC-PAR in kinematic scenes.Figure 12 shows the value of the PDOP and HDOP of the kinematic experiment.From Figure 12, we can see that the value of the PDOP and HDOP varies rapidly.Figure 13 show that the positioning accuracy of the TSSC-PAR is significantly better than that of the BSSC PAR, in the east, north, and up directions.Figure 14 shows the number of ambiguity-fixe solutions for the BSSC-PAR and TSSC-PAR.Compared with the BSSC-PAR, the averag number of ambiguity-fixed solutions for the TSSC-PAR increased from 12.9 to 13.6.Com bined with Figures 13 and 14, it can be seen that the number of ambiguity-fixed solution by the TSSC-PAR is better than that by the BSSC-PAR in the PPP initialization phase.Th positioning accuracy of the TSSC-PAR is better than that of the BSSC-PAR in the PPP initia ization phase.Although there is little difference in the TTFF between the TSSC-PAR an BSSC-PAR, the convergence time of TSSC-PAR is significantly shorter than that of the BSSC PAR.Table 6 gives the statistic indicators of the kinematic experiment for the BSSC-PAR an TSSC-PAR.The results are evaluated, including positioning accuracy, TTFF, and conver gence time.Compared with BSSC-PAR, the positioning accuracy of TSSC-PAR improved b 38.1%, 1.4%, and 2.2% in the east, north, and up directions, respectively.Additionally, th TTFF and convergence time decreased from 1.65 min and 10.5 min to 1.3 min and 1.8 min respectively.In conclusion, the experimental results have demonstrated that the TSSC-PAR has a better performance than the BSSC-PAR in kinematic scenes.

Discussion
The PAR can enhance the AR success rate and reduce the TTFF.The core of PAR is the selection of the ambiguity subset.In this study, we aim to address the influence of inaccurate variance of ambiguity parameters, particularly those with low observation quality, on the ambiguity subset selection based on the BSSC.In this contribution, the

Discussion
The PAR can enhance the AR success rate and reduce the TTFF.The core of PAR i the selection of the ambiguity subset.In this study, we aim to address the influence o inaccurate variance of ambiguity parameters, particularly those with low observatio quality, on the ambiguity subset selection based on the BSSC.In this contribution, th

Discussion
The PAR can enhance the AR success rate and reduce the TTFF.The core of PAR is the selection of the ambiguity subset.In this study, we aim to address the influence of inaccurate variance of ambiguity parameters, particularly those with low observation quality, on the ambiguity subset selection based on the BSSC.In this contribution, the observations whose elevation angle is below 30 • or whose posterior residual falls into the IGG3 model reduction domain are clustered together.The POWC instead of the BSSC is adopted to sort ambiguities to overcome the false convergence of variance.In static and kinematic experiments and analysis, we can see that the performance of the TSSC-PAR is better than that of the BSSC-PAR.
The present study evaluates kinematic performance using limited real data; future research will involve additional kinematic experiments conducted in various urban environments.Furthermore, the method presented in this paper conducts post-processing analysis with fast precise ephemeris.Future research will focus on multi-GNSS real-time PPP PAR.Due to the potential influence on the accuracy of real-time satellite clock and orbit service, experimental verification is necessary to confirm the effectiveness of low-precision ambiguity screening in real-time PPP PAR.We hope that future research will leverage the multi-GNSS PPP PAR method based on the two-step sorting criterion to further shorten convergence time in real-time data and enhance the practicality of the proposed method.

Conclusions
The selection criterion for the ambiguity subset is crucial for improving the positioning performance of PAR.To avoid the adverse effects of ambiguity variance inaccuracy on the PAR when using the BSSC, we proposed an improved multi-GNSS PPP PAR method based on the TSSC.The observation data from 10 stations in the MGEX during the period from DOY 116 to DOY 145 in 2022, along with urban kinematic data from DOY 258 in 2021, were used to verify the performance of the proposed method.
The static experiment results demonstrate that the probability distribution of the ambiguity-fixed number for the TSSC-PAR is better than that for the BSSC-PAR.Compared with the BSSC-PAR, the positioning accuracy improved by 20.1%, 17.1%, and 4.67% in the east, north, and up directions, respectively.The average TTFF and convergence time improved by 8.7% and 16.8%, respectively.Compared with the BSSC-PAR, the urban kinematic experiment results demonstrate that the positioning accuracy for the TSSC-PAR enhanced by 38.1%, 1.4%, and 2.2%, in the east, north, and up directions, respectively.The TTFF and convergence time decreased from 1.65 min and 10.5 min to 1.3 min and 1.8 min, respectively.In summary, the positioning performance of the TSSC-PAR is better than that of the BSSC-PAR for static and kinematic data.

Figure 1 .
Figure 1.Flowchart of multi-GNSS PPP PAR method based on TSSC.

Figure 1 .
Figure 1.Flowchart of multi-GNSS PPP PAR method based on TSSC.

Figure 3 .
Figure 3.The positioning accuracy of the BSSC-PAR and TSSC-PAR for AMC4 station on DOY 116, 2022.

Figure 3 .
Figure 3.The positioning accuracy of the BSSC-PAR and TSSC-PAR for AMC4 station on DOY 116, 2022.

Figure 4 .
Figure 4.The number of ambiguity-fixed for the BSSC-PAR and TSSC-PAR for AMC4 station on DOY 116, 2022.

Figure 4 .
Figure 4.The number of ambiguity-fixed for the BSSC-PAR and TSSC-PAR for AMC4 station on DOY 116, 2022.

Figure 5 .Figure 5 .
Figure 5. Convergence series of the BSSC-PAR and TSSC-PAR for 68% quantile of the east, north, and up in 1 h pass.

Figure 5 .
Figure 5. Convergence series of the BSSC-PAR and TSSC-PAR for 68% quantile of the east, and up in 1 h pass.

Figure 6 .
Figure 6.Cumulative frequency of the TTFF and convergence time for BSSC-PAR and TSSC for AMC4 station DOY 116 to DOY 145, 2022.

18 Figure 8 .
Figure 8.The positioning accuracy of the east, north, and up for the BSSC-PAR and TSSC-PAR.

Figure 9 .
Figure 9.The average TTFF and average convergence time for the BSSC-PAR and TSSC-PAR.

Figure 8 . 18 Figure 8 .
Figure 8.The positioning accuracy of the east, north, and up for the BSSC-PAR and TSSC-PAR.

Figure 9 .
Figure 9.The average TTFF and average convergence time for the BSSC-PAR and TSSC-PAR.

Figure 9 .
Figure 9.The average TTFF and average convergence time for the BSSC-PAR and TSSC-PAR.

Figure 10 .
Figure 10.Distribution of ambiguity-fixed number for the BSSC-PAR and TSSC-PAR for selection stations from DOY 116 to DOY 145, 2022.

Figure 10 .
Figure 10.Distribution of ambiguity-fixed number for the BSSC-PAR and TSSC-PAR for selection stations from DOY 116 to DOY 145, 2022.

Figure 11 .
Figure 11.Reference station position and vehicle tracking in the kinematic experiment.

Figure 11 .
Figure 11.Reference station position and vehicle tracking in the kinematic experiment.

Figure 11 .
Figure 11.Reference station position and vehicle tracking in the kinematic experiment.

Figure 13 .
Figure 13.The positioning accuracy of the BSSC-PAR and TSSC-PAR in the kinematic experiment.

Figure 14 .
Figure 14.The number of ambiguity-fixed for the BSSC-PAR and TSSC-PAR in the kinematic experiment.

Figure 14 .
Figure 14.The number of ambiguity-fixed for the BSSC-PAR and TSSC-PAR in the kinematic experiment.
Stage 4.Finally, the bottom ambiguity from ˆm a is one by one excluded until either the ratio test is met or the number of participating fixed ambiguities falls below 5. Otherwise, the float solution is output.

Table 2 .
Statistics of the TTFF, convergence time, and positioning accuracy for BSSC-PAR and TSSC-PAR for AMC4 station on DOY 116, 2022.

Table 2 .
Statistics of the TTFF, convergence time, and positioning accuracy for BSSC-PAR and TSSC-PAR for AMC4 station on DOY 116, 2022.

Table 3 .
The statistics of the TTFF, convergence time, and positioning accuracy for BSSC-PAR and TSSC-PAR for AMC4 station DOY 116 to DOY 145, 2022.

Table 3 .
The statistics of the TTFF, convergence time, and positioning accuracy for BSSC-PA TSSC-PAR for AMC4 station DOY 116 to DOY 145, 2022.

Table 4 .
Statistics of TTFF, convergence time, and positioning accuracy of the BSSC-PAR and TSSC-PAR.

Table 4 .
Statistics of TTFF, convergence time, and positioning accuracy of the BSSC-PAR and TSSC-PAR.

Table 4 .
Statistics of TTFF, convergence time, and positioning accuracy of the BSSC-PAR and TSSC-PAR.

Table 5 .
Basic information of the experiment in urban areas.

Table 5 .
Basic information of the experiment in urban areas.

Table 6 .
Statistics of TTFF, convergence time, and positioning accuracy of the BSSC-PAR and TSSC-PAR in the kinematic experiment.

Table 6 .
Statistics of TTFF, convergence time, and positioning accuracy of the BSSC-PAR and TSSC-PAR in the kinematic experiment.