Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

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.     

精密进近阶段的多系统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%。     

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.     

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.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

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 GFZ real-time GNSS precise positioning service system and its adaption for COMPASS

Motivated by the IGS real-time Pilot Project, GFZ has been developing its own real-time precise positioning service for various applications. An operational system at GFZ is now broadcasting real-time orbits, clocks, global ionospheric model, uncalibrated phase delays and regional atmospheric corrections for standard PPP, PPP with ambiguity fixing, single-frequency PPP and regional augmented PPP. To avoid developing various algorithms for different applications, we proposed a uniform algorithm and implemented it into our real-time software. In the new processing scheme, we employed un-differenced raw observations with atmospheric delays as parameters, which are properly constrained by real-time derived global ionospheric model or regional atmospheric corrections and by the empirical characteristics of the atmospheric delay variation in time and space. The positioning performance in terms of convergence time and ambiguity fixing depends mainly on the quality of the received atmospheric information and the spatial and temporal constraints. The un-differenced raw observation model can not only integrate PPP and NRTK into a seamless positioning service, but also syncretize these two techniques into a unique model and algorithm. Furthermore, it is suitable for both dual-frequency and sing-frequency receivers. Based on the real-time data streams from IGS, EUREF and SAPOS reference networks, we can provide services of global precise point positioning (PPP) with 5–10 cm accuracy, PPP with ambiguity-fixing of 2–5 cm accuracy, PPP using single-frequency receiver with accuracy of better than 50 cm and PPP with regional augmentation for instantaneous ambiguity resolution of 1–3 cm accuracy. We adapted the system for current COMPASS to provide PPP service. COMPASS observations from a regional network of nine stations are used for precise orbit determination and clock estimation in simulated real-time mode, the orbit and clock products are applied for real-time precise point positioning. The simulated real-time PPP service confirms that real-time positioning services of accuracy at dm-level and even cm-level is achievable with COMPASS only.     

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.     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

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.     

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.     

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.     

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.     

Real-time clock comparison and monitoring with multi-GNSS precise point positioning: GPS,GLONASS and Galileo

The precise point positioning (PPP) technique is widely used in time and frequency applications. Because of the real-time service (RTS) project of the International GNSS Service, we can use the PPP technique for real-time clock comparison and monitoring. As a participant in the RTS, the Centre National d’Etudes Spatiales (CNES) implements the PPPWIZARD (Precise Point Positioning with Integer and Zero-difference Ambiguity Resolution Demonstrator) project to validate carrier phase ambiguity resolution. Unlike the Integer-PPP (IPPP) of the CNES, fixing ambiguities in the post-processing mode, the PPPWIZARD operates in the real-time mode, which is also called real-time IPPP (RT-IPPP). This paper focuses on applying the RT-IPPP for real-time clock comparison and monitoring. We review the principle of real-time clock comparison and monitoring, and introduce the methodology of the RT-IPPP technique. The observations of GPS, GLONASS and Galileo were processed for the experiments. Five processing modes were provided in the experiment to analyze the benefits of ambiguity resolution and multi-GNSS. In the clock comparison experiment, the average reduction ratios of standard deviations with respect to the G PPP mode range from 9.7% to 35.0%. In the clock monitoring experiment, G PPP mode can detect clock jumps whose magnitudes are larger than 0.9 ns. The RT-IPPP technique with GRE PPP AR (G) mode allows for the detection of any clock jumps larger than 0.6 ns. For frequency monitoring, G PPP mode allows detection of frequency changes larger than 1.1 × 10⁻¹⁴. When the RT-IPPP technique is applied, monitoring with GRE PPP AR (G) mode can detect frequency changes larger than 6.1 × 10⁻¹⁵.     

Galileo real-time orbit determination with multi-frequency raw observations

Real-time GNSS-based applications require corresponding real-time orbit products. While traditional GNSS orbits are generated with the dual-frequency IF (Ionosphere-Free) model, the increase of multi-frequency signal satellites brings new challenges for the data processing. Therefore, real-time orbit determination with the multi-frequency UC (Uncombined) model is introduced in this study considering its flexibility. With the derived mathematical model conforming to IGS (International GNSS Service) dual-frequency clock definition and one-week triple-frequency Galileo observation data from 90 IGS network stations, the convergence and accuracy of real-time orbits is assessed and the characteristics of satellite IFCB (Inter-Frequency Clock Bias) are analyzed. Results indicate that the model differences, including dual-frequency IF model, dual-frequency UC model and triple-frequency UC model, contribute to only cm-level differences with CODE (Center for Orbit Determination in Europe) final orbits after a convergence time of around 12 h. The constellation-mean RMS (Root Mean Square) differences of the converged real-time orbits with the CODE final orbits reaches about 5.0 cm, 7.0 cm and 5.0 cm for the radial, tangential and normal directions. The convergence of satellite IFCB is much faster than that of satellite orbit, which reflects a loose correlation between these two parameters. While the Galileo satellite IFCB are temporally stable, the modeling of satellite IFCB may be unreliable when over constrained and becomes even more unstable with commonly encountered datum changes. In summary, real-time GNSS orbit determination with multi-frequency raw observations is feasible and extendable with proper treatment of IFCB.     

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.     

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

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

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.     

Initial accuracy and reliability of current BDS-3 precise positioning,velocity estimation,and time transfer (PVT)

The primary system of Chinese global BeiDou satellite system (BDS-3) was completed to provide global services on December 27, 2018; this was a key milestone in the development process for BDS in terms of its provision of global services. Therefore, this study analyzed the current performance of BDS-3, including its precise positioning, velocity estimation, and time transfer (PVT). The datasets were derived from international GNSS monitoring and assessment system (iGMAS) tracking networks and the two international time laboratories in collaboration with the International Bureau of Weights and Measures (BIPM). With respect to the positioning, the focus is on the real-time kinematic (RTK) positioning and precise point positioning (PPP) modes with static and kinematic scenarios. The results show that the mean available satellite number is 4.8 for current BDS-3 system at short baseline XIA1–XIA3. The RTK accuracy for three components is generally within cm level; the 3D mean accuracy is 8.9 mm for BDS-3 solutions. For the PPP scenarios, the convergence time is about 4 h for TP01 and BRCH stations in two scenarios. After the convergence, the horizontal positioning accuracy is better than cm level and the vertical accuracy nearly reaches the 1 dm level. With respect to kinematic scenarios, the accuracy stays at the cm level for horizontal components and dm level for the vertical component at two stations. In terms of velocity estimation, the horizontal accuracy stays at a sub-mm level, and the vertical accuracy is better than 2 mm/s in the BDS-3 scenario, even in the Arctic. In terms of time and frequency transfer, the noise level of BDS-3 time links can reach 0.096 ns for long-distances link NT01–TP02 and 0.016 ns for short-distance links TP01–TP02. Frequency stability reaches 5E–14 accuracy when the averaging time is within 10,000 s for NT01–TP02 and 1E–15 for TP01–TP02.