Evaluation and analysis on positioning performance of BDS/QZSS satellite navigation systems in Asian-Pacific region

By using the observation data and products of precise obit and clock offset from Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS) and GNSS Research Centre, Curtin University in this paper, the positioning performance of BDS/QZSS satellite navigation system has been analyzed and evaluated in aspects of the quantity of visible satellites, DOP value, multipath effect, signal-to-noise ratio, static PPP and kinematic PPP. The analysis results show that compared to BDS single system when the cutoff angle are 30°and 40°, the DOP value of BDS/QZSS combined system has decreased above 20%, and the quantity of visible satellites increased about 16–30% respectively, because of the improved spatial geometric configuration. The magnitude of satellite multipath effect of BDS system shows the trend of MEO?>?IGSO?>?GEO, which is consistent with that of QZSS satellite system, as the constellation structure of the two systems is similar. The variation tendencies of signal-to-noise ratio with respect to elevation angle of the two systems are almost the same at all frequencies, showing that at the same elevation angle the signal-to-noise ratio of MEO satellites is higher than that of IGSO satellites, as the higher obit is the lower transmitting power is obtained. For having a specially designed obit, the variation of signal-to-noise ratio of BDS system is more stable. However, the magnitude of signal-to-noise ratio of QZSS system appears the trend of frequency 3?>?frequency 2?>?frequency 1. The static PPP performance of the BDS/QZSS combination system has been improved more significantly than the BDS single system in E, N and U directions. When the cutoff angle are at 7°, 15° and 30°, the PPP accuracy is increased about 25–34% in U direction, 10–13% and 23–34% in E and N directions respectively. When the elevation angle is large (40°), compared to BDS single system at lower elevation angles (7° and 15°) the PPP accuracy of the BDS/QZSS combination system is improved above 30% in U direction. In kinematic PPP performance, compared to BDS single system, the accuracy, availability and reliability of the BDS/QZSS combination system has been improved too, especially at large elevation angles (30° and 40°), the kinematic PPP accuracy in E and U directions has been improved about 10–50%, and above 50% in U direction. It can be concluded that the combination with QZSS system can improve the positioning accuracy, reliability and stability of BDS system. In the future, with the improvement of the satellite construction of Japan’s QZSS system and the global networking of China’s BDS satellites, the QZSS satellites will contribute greatly to improve the positioning accuracy, reliability, availability and stability of GNSS systems in areas such as cities, mountains, densely-packed buildings and severely covered areas in Asian-Pacific region.     

Performance assessment of real-time multi-GNSS integrated PPP with uncombined and ionospheric-free combined observables

For precise position services, the real-time precise point positioning (PPP) is a promising technology. The real-time PPP performance is expected to be improved by multi-system combination. The performance of real-time multi-system PPP needs to be periodically investigated, with the increasing number of available satellites and the continuously improved quality of real-time precise products of satellite clocks and orbits. In this study, a comprehensive performance assessment is conducted for the four-system integrated real-time PPP (FSIRT-PPP) with GPS, BDS, Galileo and GLONASS in both static and kinematic modes. The datasets from 118 stations spanning approximately a month are used for analysis, and the real-time stream CLK93 is employed. The superior performance of FSIRT-PPP is validated by comparing with the results of GPS/BDS, GPS/Galileo, GPS/GLONASS, GPS-only, BDS-only, Galileo-only and GLONASS-only cases. The FSIRT-PPP using ionospheric-free (IF) combined observables can achieve a convergence time of 10.9, 4.8 and 11.8 min and a positioning accuracy of 0.4, 0.5 and 0.7 cm in the static mode in the east, north and up directions, respectively, while the derived statistic is 15.4, 7.0 and 16.4 min, and 1.6, 1.2 and 3.4 cm in the kinematic mode in the three directions, respectively. Moreover, we also compare the position solutions of real-time PPP adopting IF combined and uncombined (UC) observables, and prove the mathematical equivalence between the two PPP models in the converged stage, provided that there are no external ionospheric corrections or constraints given to the estimated ionospheric delays in the UC model. The difference between the fully converged positioning accuracy of IF-based and UC-based real-time PPP is marginal, but the UC-based real-time PPP has longer convergence time due to the influence of the significant unmodeled time-varying errors in the real-time precise products as well as the different parameterization between them. For completeness, the real-time kinematic PPP results in harsh environments and the post-processed PPP results are also presented.     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

Comparison of convergence time and positioning accuracy among BDS,GPS and BDS/GPS precise point positioning with ambiguity resolution

Precise point positioning (PPP) usually takes about 30?min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60?days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0?cm, 1.5?cm, and 5.9?cm respectively, while those of the Gfix are 0.8?cm, 0.5?cm, and 1.7?cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7?cm, 0.5?cm, and 1.9?cm for the E, N, U components, respectively.     

Assessing the performance of BDS-3 for multi-GNSS static and kinematic PPP-AR

As of 2021, a total of four different GNSS constellations – namely, GPS, GLONASS, Galileo, and BDS-3 – can be used with Full Operational Capability (FOC). In this work, the contribution of BDS-3 FOC to GPS + GLONASS + Galileo (GRE) PPP-AR is investigated, considering the three different cut-off angles (7°, 30°, and 45°) and different lengths of static observation sessions (24-, 12-, 6-, 3-, 1-, 0.5-, 0.25-hour). The data of 31 IGS-MGEX stations is processed with GRE PPP-AR and GREC3 (GPS + GLONASS (using float mode) + Galileo + BDS-3) PPP-AR modes. The results showed that BDS-3 degraded the horizontal (except for 24-h sessions) and vertical accuracy of static GRE PPP-AR solutions regardless of the elevation cutoff angle and observation time. The kinematic results showed that BDS-3 significantly contributed to the accuracy of GRE kinematic PPP-AR for 30° and 45° cutoff angles. The convergence time analysis showed that BDS-3 only contributes to GRE kinematic PPP-AR for the vertical component.     

Performance improvement of real-time PPP ambiguity resolution using a regional integer clock

Obtaining reliable GNSS uncalibrated phase delay (UPD) or integer clock products is the key to achieving ambiguity-fixed solutions for real-time (RT) precise point positioning (PPP) users. However, due to the influence of RT orbit errors, the quality of UPD/integer clock products estimated with a globally distributed GNSS network is difficult to ensure, thereby affecting the ambiguity resolution (AR) performance of RT-PPP. In this contribution, by fully utilising the consistency of orbital errors in regional GNSS network coverage areas, a method is proposed for estimating regional integer clock products to compensate for RT orbit errors. Based on Centre National d’Études Spatiales (CNES) RT precise products, the regional GPS/BDS integer clock was estimated with a CORS network in the west of China. Results showed that the difference between the estimated regional and CNES global integer clock/bias products was generally less than 5 cm for GPS, whereas clock differences of greater than 10 cm were observed for BDS due to the large RT orbit error. Compared with PPP using global integer clock/bias products, the AR performance of PPP using the regional integer clock was obviously improved for four rover stations. For single GPS, the horizontal and vertical accuracies of ambiguity-fixed PPP solutions were improved by 56.2% and 45.3% on average, respectively, whereas improvements of 67.5% and 50.5% in the horizontal and vertical directions, respectively, were observed for the combined GPS/BDS situation. Based on a regional integer clock, the RMS error of a kinematic ambiguity-fixed PPP solution in the horizontal direction could reach 0.5 cm. In terms of initialisation time, the average time to first fix (TTFF) in combined GPS/BDS PPP was shortened from 18.2 min to 12.7 min. In view of the high AR performance realised with the use of regional integer clocks, this method can be applied to scenarios that require high positioning accuracy, such as deformation monitoring.     

GLONASS-based precise point positioning and performance analysis

Current precise point positioning (PPP) techniques are mainly based on GPS which has been extensively investigated. With the increase of available GLONASS satellites during its revitalization, GLONASS observations were increasingly integrated into GPS-based PPP. Now that GLONASS has reached its full constellation, there will be a wide interest in PPP systems based on only GLONASS since it provides a PPP implementation independent of GPS. An investigation of GLONASS-based PPP will also help the development of GPS and GLONASS combined PPP techniques for improved precision and reliability. This paper presents an observation model for GLONASS-based PPP in which the GLONASS hardware delay biases are addressed. In view of frequently changed frequency channel number (FCN) for GLONASS satellites, an algorithm has been developed to compute the FCN for GLONASS satellites using code and phase observations, which avoids the need to provide the GLONASS frequency channel information during data processing. The observation residuals from GLONASS-based PPP are analyzed and compared to those from GPS-based PPP. The performance of GLONASS-based PPP is assessed using data from 15 globally distributed stations.     

The update of BDS-2 TGD and its impact on positioning

Timing group delay (TGD) is an important parameter that affects the positioning performance of global navigation satellite systems (GNSS). The BeiDou navigation satellite system (BDS) broadcasts TGD corrections from B3I frequency to B1I and B2I frequencies, namely TGD₁ and TGD₂. On July 21, 2017, BDS updated TGD values with a maximum change of more than 4 ns. In this contribution, we explain the motivation for the BDS TGD update, which is due to the systematic bias between narrowly correlated and widely correlated pseudo-ranges in BDS monitoring receivers. To investigate the impact of the updated TGD, BDS signal-in-space range error (SISRE) and user positioning performance regarding single point positioning (SPP) and precise point positioning (PPP) are analyzed. Results show that after the update of TGD, the difference between the new TGD and multi-GNSS experiment (MGEX) differential code bias (DCB) decreases from 1.38 ns to 0.29 ns on TGD₁ and from 0.40 ns to 0.25 ns on TGD₂. With the contribution of more accurate TGD, the systematic bias of BDS radial SISRE no longer exists, and the overall BDS SISRE also reduces from 1.33 m to 0.87 m on B1I/B2I frequency, from 1.05 m to 0.89 m on B1I frequency, from 0.92 m to 0.91 m on B2I frequency, respectively, which proves the similar precision of BDS TGD and MGEX DCB. One week of statistical results from 28 globally distributed MGEX stations shows that the SPP performance improves on non-B3I frequencies after the TGD update, with a maximum improvement of more than 22% for the B1I/B2I or B1I/B3I combination. The new TGD mainly reduces SPP positioning bias in the East component. The updated TGD also slightly improves the PPP convergence performance for the B1I/B3I combination, but mostly contributes to a more accurate estimation of the receiver clock and ambiguities.     

BDS-2/BDS-3 uncalibrated phase delay estimation considering the intra-system bias

Intra-system biases (ISBs) between BDS-2 and BDS-3 are of critical importance when combining observations from the BDS-2 and BDS-3 systems, which is meaningful to fully take advantage of the BDS positioning capability. Meanwhile, ISBs should also be considered in the estimation of BDS uncalibrated phase delays (UPDs). In this research, we present a BDS-2/BDS-3 joint-processing scheme, as well as a method for estimating BDS UPDs. The characteristics of ISBs and the quality of BDS UPDs are analyzed based on 30-day data from 130 multi-GNSS experimental (MGEX) stations. Our results indicate that the ISBs are related to the type and version of the receiver. The ISBs can be regarded as constant across the course of a given day, and the mean standard deviation (STD) values of ISBs over one month for different types of receivers are generally within 0.2 m. Moreover, to assess the quality of UPD products, the residuals of the estimated UPDs and the utilization rates of the observation data are computed. The results show that the quality of BDS UPDs can be improved by correcting the satellite-induced pseudo-range variations, and by estimating the wide-lane (WL) UPD difference between BDS-2 and BDS-3. The average RMS values of the estimated residuals of WL UPD and narrow-lane (NL) UPD are 0.07 and 0.09 cycles, respectively; moreover, the utilization rate of the observation data of WL UPD and NL UPD can reach above 90 %. The performance of BDS precise point positioning (PPP) and PPP ambiguity resolution (PPP-AR) is analyzed in terms of positioning accuracy and convergence performance in both the static and kinematic modes. Compared with PPP ambiguity-float solutions, the positioning accuracy of PPP-AR is significantly improved, especially in the east direction. The impact of ISBs on PPP and PPP-AR is also analyzed, and the results indicate that ISBs can improve the convergence speed of float PPP, but can be disregarded in PPP-AR.     

Multi-GNSS contributions to differential code biases determination and regional ionospheric modeling in China

Due to the limited number and uneven distribution globally of Beidou Satellite System (BDS) stations, the contributions of BDS to global ionosphere modeling is still not significant. In order to give a more realistic evaluation of the ability for BDS in ionosphere monitoring and multi-GNSS contributions to the performance of Differential Code Biases (DCBs) determination and ionosphere modeling, we select 22 stations from Crustal Movement Observation Network of China (CMONOC) to assess the result of regional ionospheric model and DCBs estimates over China where the visible satellites and monitoring stations for BDS are comparable to those of GPS/GLONASS. Note that all the 22 stations can track the dual- and triple-frequency GPS, GLONASS, and BDS observations. In this study, seven solutions, i.e., GPS-only (G), GLONASS-only (R), BDS-only (C), GPS + BDS (GC), GPS + GLONASS (GR), GLONASS + BDS (RC), GPS + GLONASS + BDS (GRC), are used to test the regional ionosphere modeling over the experimental area. Moreover, the performances of them using single-frequency precise point positioning (SF-PPP) method are presented. The experimental results indicate that BDS has the same ionospheric monitoring capability as GPS and GLONASS. Meanwhile, multi-GNSS observations can significantly improve the accuracy of the regional ionospheric models compared with that of GPS-only or GLONASS-only or BDS-only, especially over the edge of the tested region which the accuracy of the model is improved by reducing the RMS of the maximum differences from 5–15 to 2–3 TECu. For satellite DCBs estimates of different systems, the accuracy of them can be improved significantly after combining different system observations, which is improved by reducing the STD of GPS satellite DCB from 0.243 to 0.213, 0.172, and 0.165 ns after adding R, C, and RC observations respectively, with an increment of about 12.3%, 29.4%, and 32.2%. The STD of GLONASS satellite DCB improved from 0.353 to 0.304, 0.271, and 0.243 ns after adding G, C, and GC observations, respectively. The STD of BDS satellite DCB reduced from 0.265 to 0.237, 0.237 and 0.229 ns with the addition of G, R and GR systems respectively, and increased by 10.6%, 10.4%, and 13.6%. From the experimental positioning result, it can be seen that the regional ionospheric models with multi-GNSS observations are better than that with a single satellite system model.     

Evaluation on real-time dynamic performance of BDS in PPP,RTK, and INS tightly aided modes

Since China’s BeiDou satellite navigation system (BDS) began to provide regional navigation service for Asia-Pacific region after 2012, more new generation BDS satellites have been launched to further expand BDS’s coverage to be global. In this contribution, precise positioning models based on BDS and the corresponding mathematical algorithms are presented in detail. Then, an evaluation on BDS’s real-time dynamic positioning and navigation performance is presented in Precise Point Positioning (PPP), Real-time Kinematic (RTK), Inertial Navigation System (INS) tightly aided PPP and RTK modes by processing a set of land-borne vehicle experiment data. Results indicate that BDS positioning Root Mean Square (RMS) in north, east, and vertical components are 2.0, 2.7, and 7.6?cm in RTK mode and 7.8, 14.7, and 24.8?cm in PPP mode, which are close to GPS positioning accuracy. Meanwhile, with the help of INS, about 38.8%, 67.5%, and 66.5% improvements can be obtained by using PPP/INS tight-integration mode. Such enhancements in RTK/INS tight-integration mode are 14.1%, 34.0%, and 41.9%. Moreover, the accuracy of velocimetry and attitude determination can be improved to be better than 1?cm/s and 0.1°, respectively. Besides, the continuity and reliability of BDS in both PPP and RTK modes can also be ameliorated significantly by INS during satellite signal missing periods.     

Estimation of QZSS differential code biases using QZSS/GPS combined observations from MGEX

As an important error source in Global Navigation Satellite System (GNSS) positioning and ionospheric modeling, the differential code biases (DCB) need to be estimated accurately, e.g., the regional Quasi-Zenith satellite system (QZSS). In this paper, the DCB of QZSS is estimated by adopting the global ionospheric modeling method based on QZSS/GPS combined observations from Multi-GNSS experiment (MGEX). The performance of QZSS satellite and receiver DCB is analyzed with observations from day of year (DOY) 275–364, 2018. Good agreement between our estimated QZSS satellite DCB and the products from DLR and CAS is obtained. The bias and root mean square (RMS) of DCB are mostly within ±0.3 ns. The day-to-day fluctuation of the DCB time series is less than 0.5 ns with about 96% of the cases for all satellites. However, the receiver DCB is a little less stable than satellite DCB, and their standard deviations (STDs) are within 1.9 ns. The result shows that the stability of the receiver DCBs is not significantly related to the types of receiver or antenna.     

精密进近阶段的多系统GNSS组合RAIM可用性算法及分析

给出了多系统全球卫星导航系统（GNSS）组合接收机自主完好性监测（ReceiverAutonomousIntegrityMonitoring,RAIM）可用性计算方法,在此基础上利用GPS、GLONASS实测数据与BDS、Galileo全星座仿真数据,分析了BDS、GPS、GLONASS和Galileo不同组合在精密进近阶段的RAIM可用性。通过试验分析发现,BDS的5颗地球同步轨道卫星和3颗倾斜地球同步轨道卫星对亚洲、非洲和欧洲大部分地区的RAIM可用性有很大的贡献。这些地区站星间几何观测结构得到改善,使得RAIM可用性相对于其他地区有很大幅度的提升。在亚太地区APV-I阶段单系统导航情况下,北斗导航系统RAIM可用性达到99.5%,高于其他三个导航系统。在精密进近阶段（APV-I、APV-II和CAT-I）,BDS与其他导航系统（GPS、GLONASS和Galileo）的组合导航可以满足全球大部分区域的RAIM可用性需求,大多可达到100%。     

Tightly combined GPS/Galileo RTK for short and long baselines: Model and performance analysis

To ensure the compatibility and interoperability with modernized GPS, Galileo satellites are capable of broadcasting navigation signals on carrier phase frequencies that overlap with GPS, i.e., GPS/Galileo L1-E1/L5-E5a. Moreover, the GPS/Galileo L2-E5b signals have different frequencies with wavelength differences smaller than 4.2?mm. Such overlapping and narrowly spaced signals between GPS and Galileo bring the opportunity to use the tightly combined double-differenced (DD) model for precise real-time kinematic (RTK) positioning, resulting in improved performance of ambiguity resolution and positioning with respect to the classical standard or loosely combined DD model. In this paper, we focus on the model and performance assessment of tightly combined GPS/Galileo L1-E1/L2-E5b/L5-E5a RTK for short and long baselines. We first investigate the tightly combined GPS/Galileo DD observational model for both short and long baselines with simultaneously considering the GPS/Galileo overlapping and non-overlapping frequencies. Particularly, we introduce a reparameterization approach to solve the rank deficiency that caused by the correlation between the DISB parameters and the DD ionospheric parameters for both overlapping and non-overlapping frequencies. Then we present performance assessment for the tightly combined GPS/Galileo RTK model with real-time estimation of the differential inter-system bias (DISB) parameters for short and long baselines in terms of ratio value, ambiguity dilution of precision (ADOP), ambiguity conditional number, decorrelation number, search count, empirical success rate, time-to-first-fix (TTFF), and positioning accuracy. Results from both static and kinematic experiments demonstrated that compared to the loosely combined model, the tightly combined model can deliver improved performance of ambiguity resolution and precise positioning with different satellite visibility. For the car-driven short baseline experiment with 10° elevation cut-off angle, the tightly combined model can not only significantly increase the ratio value by approximately 27.5% (from 16.0 to 20.4), but also reduce the ambiguity ADOP, the conditional number, and the search count in LAMBDA by approximately 22.2% (from 0.027 to 0.021 cycles), 14.9% (from 199.2 to 169.6), and 25.4% (from 150.1 to 112.0), respectively. Comparable decorrelation number, empirical success rate, and positioning accuracy are also obtained. For the car-driven long baseline experiment, it is also observed that the ambiguity resolution performance in terms of the ratio value, the decorrelation number, the condition number, and the search count are significantly improved by approximately 18.5% (from 2.7 to 3.2), 22.0% (from 0.186 to 0.227), 55.9% (from 937.6 to 413.7), and 10.3% (from 43.8 to 39.3), respectively. Moreover, comparable ADOP, empirical success rate, and positioning accuracy are obtained as well. Additionally, the TTFF can be reduced (from 54.1 to 51.8 epochs with 10° elevation cut-off angle) as well from the results of static experiments.     

Preliminary analysis of the quality and positioning performance of BDS-3 global interoperable signal B1C&B2a

BeiDou-3 Navigation Satellite System (BDS-3) satellites are equipped with the new generation GNSS signals B1C and B2a, which support the interoperability with GPS and Galileo systems. In this study, the pseudo-range multipath error and carrier phase observation noise of the BDS-3 B1C and B2a signals were evaluated based on zero baseline measurements from the day of year (DOY) 113 to 116, 2020. Further, the precision and performance of the single point positioning (SPP) and precise point positioning (PPP) are assessed at 9 stations. This assessment manifests that the standard deviations (STDs) of the pseudo-range multipath error are about 0.09 ~ 0.22 m, while STDs of the carrier phase observation noise are about 0.075 mm. For the single-frequency SPP, its positioning precision is about 2.03 ~ 4.85 m and 3.29 ~ 10.73 m at the 99.99% confidence level in horizontal and vertical directions, respectively, while the dual-frequency SPP precision is about 1.92 ~ 8.02 m and 4.81 ~ 12.77 m in horizontal and vertical directions, respectively. For the daily static PPP, the convergence time is about 20 ~ 30 min, while the daily positioning precision can reach 1.38 ~ 4.42 cm and -1.31 ~ 4.34 cm in horizontal and vertical directions, respectively. In general, the quality and the SPP and PPP performance of the BDS-3 B1C&B2a signals are comparable to the GPS and Galileo.     

Enhanced multipath mitigation method based on multi-resolution CNR model and adaptive statistical test strategy for real-time kinematic PPP

This study proposes an enhanced multipath mitigation method based on multi-resolution carrier-to-noise-ratio (CNR) model and adaptive statistical test strategy for real-time kinematic precise point positioning (PPP) applications. The multi-resolution CNR model is established with GPS observation data collected from DOY 152 to 181 of 2019 by 230 globally distributed IGS stations, which used to analyze the relevant factors affecting CNR. Statistical results indicate that the CNR is not only related to the satellite elevation, but also closely related to the receiver types and specific satellite. The maximum difference between different receivers can reach 20 dB for the same satellite at the same elevation. In addition, the performance of the CNR is also obviously different between each satellite, and the maximum difference between different satellites is about 10 dB for the same receiver at the same elevation. Hence, in terms of the method which is based on CNR information for multipath detection and mitigation, the independence of receiver types, satellite and frequency must be considered. With the above analysis, this study developed a multi-resolution CNR model based on different receiver types, different satellites and different elevation firstly. Then, combined with the adaptive statistical test strategy which is based on the difference of CNR between inter-frequency and the difference of CNR between adjacent epochs, the multipath can be detected effectively. For the epoch which affected by multipath, the down-weighted strategy based on CNR is adopted to mitigate the influence of multipath on positioning. Real-time kinematic PPP data are collected to assess the proposed method, and the results demonstrate that the proposed method can detect the multipath effectively, and the detection rate can reach 90.28%. Moreover, after adopting the mitigation strategy, the RMS bias of the east, north and up components are improved about 19.95%, 17.89% and 23.07% compared to the original results, respectively. It is worth noting that this proposed method is also suitable for other GNSS, such as GLONASS and BDS, but the corresponding CNR model must be established simultaneously.     

Correction model of BDS satellite-induced code bias and its impact on precise point positioning

There are code biases on the pseudo-range observations of the Beidou Navigation Satellite System (BDS) that range in size from several decimeters to larger than one meter. These biases can be divided into two categories, which are the code biases in the pseudo-range observations of Inclined Geo-Synchronous Orbit (IGSO) satellites and Medium Earth Orbit (MEO) satellites and the code biases in the pseudo-range observations of Geosynchronous Earth Orbit (GEO) satellites. In view of the code bias of the IGSO/MEO satellites, the code bias correction model is established using the weighted least square curve fitting method. After the correction, the code biases of the IGSO and MEO satellites are clearly mitigated. A methodology of correcting GEO code bias is proposed based on the empirical mode decomposition (EMD)-wavelet transform (WT) coupled model. The accuracies of the GEO multipath combination of the B1, B2 and B3 frequencies are improved by 39.9%, 17.9%, and 29.4%, respectively. Based on the corrections above, the ten days observations of three Multi-GNSS Experiment (MGEX) stations are processed. The results indicate that the convergence time of the precise point positioning (PPP) can be improved remarkably by applying a code bias. The mean convergence time can be improved by 14.67% after the IGSO/MEO code bias correction. By applying the GEO code bias, the mean convergence time can be further improved by 17.42%.     

The realization and convergence analysis of combined PPP based on raw observation

In order to speed up Precise Point Positioning (PPP)’s convergence, a combined PPP method with GPS and GLONASS which is based on using raw observations is proposed, and the positioning results and convergence time have been compared with that of single system. The ionospheric delays and receiver’s Differential Code Bias (DCB) corrections are estimated as unknown parameters in this method. The numerical results show that the combined PPP has not caused significant impacts on the final solutions, but it greatly improved Position Dilution of Precision (PDOP) and convergence speed and enhanced the reliability of the solution. Meanwhile, the convergence speed is greatly influenced by the receiver’s DCB, positioning results in horizontal which are better than 10 cm can be realized within 10 min. In addition, the ionosphere and DCB products can be provided with high precision.     

Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK

The overlapping carrier frequencies L1/E1, L5/E5a and B2/E5b from GPS/Galileo/BDS allow inter-system double-differencing of observations, which shows a clear advantage over differencing of the observations of each constellation independently. However, the inter-system biases destroy the integer nature of the inter-system double-differencing ambiguities. Two methods of direct rounding and parameter estimation are used to determine the ISB value. By analyzing data collected from Curtin University from 2015 to 2018, the phase and code inter-system bias (ISB) are related to the receiver type, firmware version and the selected overlapping frequency. Upgrade of receiver firmware version results in changes of ISB values. For example, the upgrade of Javad firmware in Dec, 15, 2017 causes the difference of 0.5 cycles ISB between BDS GEO and non-GEO satellites. By comparing the three dynamic models which include white noise process, random walk process, and random constant in the parameter estimation method, the ISB determined by the random constant model is consistent with the value obtained by the direct rounding method. After the calibration of ISBs, the performances of tightly combined positioning are assessed. The success rate of ambiguity resolution and accuracy of positioning for the tight combination (TC) are significantly improved in comparison with that for the loose combination (LC) over short baselines. For L5/E5a, on which only few satellites can be observed, the maximum increase in success rates of ambiguity resolution can reach 31.7%, i.e., from 54.9% of LC to larger than 86.6% of TC, and the positioning accuracies can even be increased by 0.13 m, i.e., from 0.208 of LC to 0.074 m of TC in East direction for the mix-receiver TRIMBLE NETR9-SEPT POLARX4 in 2018.     

观测资料和卫星钟差对PPP参数收敛影响分析

在GPS 单点定位中, 参数解算的收敛时间和收敛稳定性是重要的研究内容之一, 影响收敛时间和收敛稳定性的因素很多, 本文主要就观测资料的不同采样间隔、卫星钟差资料的不同采样间隔、不同的定位精度要求对精密单点定位中参数收敛时间的影响进行了深入的分析探讨, 以中国上海GPS综合应用网中的12个测站资料为例, 分析了采样间隔、定位精度要求与收敛时间的关系, 并对不同采样间隔的收敛时间进行了统计分析, 得出一些初步结论.