3.1. PPP-B2b Orbit Accuracy Assessment
The PPP-B2b orbit correction sequences of the MEO satellites of BDS-3 decoded from the PPP-B2b message on 9 August 2021 (DoY 221) are shown in
Figure 2. The variations in orbit correction sequence for MEO satellites are within 0.1, 0.3, and 0.3 m in radial, along-track, and cross-track components, respectively. The centimeter-level corrections indicate that considerable accuracy could be achieved by the broadcast orbits of BDS-3 MEO satellites, which mainly benefits from the application of inter-satellite links (ISL) technology [
38].
Figure 2 also reveals the hourly discontinuity of the orbit correction series, which is related to the hourly update of the broadcast ephemeris [
39].
The real-time orbit errors in radial, along-track, and cross-track components for the PPP-B2b referenced to the final precise orbit from GFZ are presented in
Figure 3a, and those of broadcast ephemeris are presented in
Figure 3b. The WHU real-time orbit errors referenced to GFZ final products are plotted in
Figure 3c. The orbit errors for both real-time PPP-B2b orbit and broadcast orbit exhibit the same level of accuracy. The real-time orbit error series of PPP-B2b is smoother and more continuous than that of the broadcast orbit error. The orbit error of WHU real-time product is more continuous and smoother, and more consistent with the GFZ final orbit product. The satellite orbit errors of C41 and C42 of the PPP-B2b real-time orbits show some abnormal errors with 0.6 m in radial components, compared to the final precise orbit from GFZ. The significant difference in the orbits between the PPP-B2b and that of GFZ may relate to the PCO corrections released by IGS. To further analyze this orbit bias, the orbit errors in radial, along-track, and cross-track components that apply the PCO corrections provided by China Satellite Navigation Office (CSNO) are presented in
Figure 4. The abnormal orbit error of the satellites C41 and C42 no longer exist, but satellite orbit C44 appeared to have a significant bias, with about 1 m, on the radial component. Thus, the significant biases on satellite orbits are related to difference in the PCO corrections, which do not participate in real-time PPP solutions with PPP-B2b signals.
To further evaluate the accuracy of satellite orbits of PPP-B2b products, the Root Mean Square Error (RMSE) of PPP-B2b real-time orbits, broadcast orbits, and real-time orbits from WHU based on a 7-day dataset (6–12 August 2020) are presented in
Figure 5. The RMSE of different satellite orbits also indicates that the radial components of both PPP-B2b orbit and broadcast orbit have better performance than that of the along-track and cross-track components. The higher accuracy on the radial component is mainly related to the high-quality onboard hydrogen and rubidium clocks, which compensate for the systematic error in the radial component [
31]. The corresponding average RMSEs of the BDS-3 MEO satellites are listed in
Table 2. The average RMSE of PPP-B2b orbits is 8.5, 19.3, 14.0 cm in radial, along-track, and cross-track components, respectively. The difference between the PPP-B2b orbit and the broadcast orbit is only centimeter-level. This slight accuracy improvement may be related to the better continuity of the PPP-B2b orbit than that of the broadcast orbit. The WHU real-time orbits exhibit better accuracy than that of PPP-B2b orbits. The RMSE of WHU real-time orbits on radial, along-track, and cross-track components are 4.6, 8.1, and 6.1 cm, respectively. The better accuracy of the WHU real-time orbit is mainly benefited by the global GNSS network, while the PPP-B2b service only relies on the regional monitoring stations in China.
3.2. PPP-B2b Clock Offset Accuracy Assessment
The PPP-B2b clock offset correction sequences of BDS-3 MEO satellites on DoY 221 are shown in
Figure 6a. The PPP-B2b clock offset corrections of each satellite are within 3 m and contain systematic offsets. The details of clock offset correction parameters of satellites C24, C35, and C44 are presented in
Figure 6b. It is easy to find that the corrections of these three satellites fluctuate around constant values of about −1, 2, and 0 m, respectively. The corrections also exhibit discontinuity in hourly duration, which is also consistent with the update frequency of the broadcast ephemeris.
The error sequences of broadcast clock offset and the PPP-B2b clock offset referenced to GFZ precise clock offsets are presented in
Figure 7a,b, respectively. As shown in
Figure 7b, the continuity of the broadcast satellite clock offset is significantly improved by applying the PPP-B2b clock offset corrections. However, individual clock offset biases still exist in the clock offset error sequence. To effectively evaluate the accuracy of the PPP-B2b satellite clock offsets, the individual clock offset biases need to be corrected by subtracting the average clock offset, as described in re-edit Equations (9) and (10).
The series of re-edit clock offset differences of the broadcast satellite, and that of the PPP-B2b corrections, as well as the errors of WHU real-time clock offset are shown in
Figure 8. The re-edit clock offset bias of the PPP-B2b becomes smoother and more continuous. As shown in
Figure 8b, the clock offsets of PPP-B2b, however, clearly exhibit satellite-specified systematic errors. For example, the real-time PPP-B2b clock bias of satellite C35 reaches 5 ns, which results in more than 1.5 m systematic bias. The study revealed that these satellite-specified systematic errors are related to the pseudorange observations [
28]. If carrier phase observations are used for positioning, the satellite-specified systematic errors can be absorbed by ambiguities [
10,
28].
As the satellite-specified systematic errors existed in the BDS-3 real-time clock products, the STD values acquired from re-edited error sequences of BDS-3 MEO satellites over a 7-day period are summarized and presented in
Figure 9. The STD of PPP-B2b clock offsets is within 0.25 ns. The average STD of the PPP-B2b clock offsets is 0.124 ns, which is significantly improved compared to that of the broadcast clock offsets with the STD of 0.672 ns. The average STD of PPP-B2b clock offsets is also smaller than that of WHU real-time clock product, which is 0.360 ns.
3.5. Real-Time PPP Positioning Accuracy Assessment
The positioning results on DOY 221 of eight IGS MEGX stations in China and surrounding counties were used to assess the positioning accuracy of the PPP-B2b services, as shown in
Figure 14. The red dash indicates the service area of the PPP-B2b signal. The corresponding observation files were downloaded from the File Transfer Protocol Server (FTP) of GFZ. Both static and kinematic mode PPP were processed with B1I/B3I IF combinations. The processing strategies are shown in
Table 4 in detail. In this section, three stations, JFNG, MIZU, and SGOC are employed, to demonstrate the performance of the PPP-B2b service. The JFNG station is located in China, while the MIZU station and the SGOC station are situated in Japan and Sri Lanka, the boundary of the PPP-B2b service area, as shown in
Figure 14. The Position Dilution of Precision (PDOP) value, positioning accuracy, and convergence time of the three stations are analyzed in detail.
The visible satellites and available satellites with PPP-B2b corrections on DoY 221 are presented in
Figure 15. During the selected test period, at least nine satellites of BDS-3 are visible at these three stations. The number of satellites with available PPP-B2b corrections, however, is less than that of visible satellites. As shown in
Figure 15, the number of minimum available satellites is 4 at the SGOC station during the selected test period. The PPP-B2b service availability is further analyzed in these three stations. As shown in
Figure 16, the visible satellites and visible periods are marked as red lines, while the available satellites with PPP-B2b correction and the available duration are marked with purple lines. It is clearly illustrated that the visible satellites of each station are similar, but the number of PPP-B2b available satellites of each station is different. In JFNG station, the number of available satellites is similar to that of visible satellites, and the PPP-B2b availability rate is 84.3%, while the availability rates of MIZU and SGOC stations are 74.6% and 57.2%, respectively. The availability rates of the three stations decrease as the stations are closer to the boundary of the service area.
The PDOP is an important factor affecting the PPP positioning accuracy and convergence time [
17,
40,
41]. The Dilution of Precision (DOP) matrix can be expressed as:
where
is the design matrix of PPP solutions as expressed in Equation (16). The PDOP is expressed as:
where
,
, and
are the first three diagonal elements of the DOP matrix
. In our case, the PDOP is calculated according to the PPP-B2b available satellites.
The PDOP values of the JFNG, MIZU, and SGOC stations on DoY 221 are illustrated in
Figure 17. The PDOP values of the JFNG station fluctuate around 2 with a maximum value of 6. The PDOP values of MIZU and SGOC stations fluctuate more widely compared with those of the JFNG station. The maximum PDOP value of the MIZU station is 42 at epoch 7 h:12 m:00 s, while that of the SGOC station is 67 at epoch 10 h:18 m:00 s. To further analyze the fluctuation in PDOP values of these three stations, the sky images at 10 h:18 m:00 s are drawn in
Figure 18. At that epoch, the SGOC station has 11 visible satellites, including 5 satellites with PPP-B2b correction information, 4 of which are located on the north side of the SGOC station. Thus, the poor observation geometry results in the large value of PDOP.
The average number of visible satellites, the average number of PPP-B2b available satellites, and the corresponding availability rate, as well as the average PDOP values over 7 days for eight stations are listed in
Table 5. The average PPP-B2b correction availability decreases as distance from the center of the PPP-B2b service area. The maximum PPP-B2b correction availability rate is 88.76% at the JFNG station, and the minimum availability rate of that is 60.91% at the SGOC station. The variation in average PDOP values showed similar trends, ranging from 2.55 to 4.47.
To illustrate the positioning performance of the PPP-B2b service, the PPP with WHU real-time products (WHU PPP) is also carried out to compare with the PPP-B2b results (B2b PPP). The positioning errors of B2b PPP and WHU PPP on the east (E), north (N), and up (U) components for the static positioning mode compared to the reference positions are presented in
Figure 19 and
Figure 20, while the positioning errors in the kinematic positioning mode are presented in
Figure 21 and
Figure 22, respectively.
In
Figure 19 and
Figure 20, both B2b PPP and WHU PPP of the three stations achieve centimeter-level positioning accuracies after convergence. For B2b PPP, the positioning performance of the JFNG station, in terms of converge time and positioning accuracy, is better than that of the MIZU and the SGOC stations. For WHU PPP, the three stations present similar error levels, but the SGOC station needs longer convergence time on U components. The RMSE of E, N, and U components, as well as horizontal (H), and the three-dimensions (3D) for the eight stations on the static positioning mode based on the B2b PPP and WHU PPP are calculated and listed in
Table 6. For B2b PPP, the RMSE on the E, N, and U components of all stations are similar, achieving centimeter-level accuracy. The RMSE of the N component is smaller than that of the E and U components in all stations. The RMSE of JFNG station on the E, N, and U components is slightly better than those of other stations due to the station being within the service area. The 3D errors vary between 3.9 and 8.3 cm. The performance of WHU PPP presents a similar trend as B2b PPP in the analyzed area. The RMSE of 3D error of each station based on the WHU PPP, however, is larger than that of B2b PPP in all stations. For example, for the JFNG station, the RMSE of 3D error is 3.9 cm for B2b PPP, while that of WHU PPP is 7.2 cm.
In
Figure 21, the positioning errors of B2b PPP at the JFNG on E, N, and U components are within 0.3 m after convergence for the kinematic positioning mode. The horizontal components of the positions of the MIZU and the SGOC stations are similar to those of the JFNG. The U component error of the MIZU and the SGOC stations, however, is relatively large compared to that of the JFNG station. The maximum error of the U component is at SGOC station with 1.3 m, which is caused by poor observation geometrics and a lower number of available satellites with PPP-B2b corrections. The peak value can be easily observed on the SGOC station at epoch 10 h:18 m:00 s, 12 h:00 m:00 s, and 15 h:30 m:00 s, when the positioning error is larger than 0.9 m. The number of PPP-B2b available satellites at these epochs are less than 5, and the PDOP value is larger than 10, as shown in
Figure 15 and
Figure 17. In
Figure 22, the positioning errors of WHU PPP of the three stations are similar and all of them achieve decimeter-level accuracy. The positioning performance of the SGOC station with WHU PPP, compared to that of B2b PPP, is not affected by the number of PPP-B2b available satellites. The convergence time of WHU PPP, however, is slightly longer than that of B2b PPP.
The RMSE of E, N, and U components, horizontal (H), and 3-dimensions (3D) for eight stations on kinematic positioning mode regarding B2b PPP and WHU PPP are calculated and listed in
Table 7. The RMSE of the E, N, and U components of all stations achieve decimeter-level accuracy. The RMSE of the N component is smaller than that of the E and U components in all stations. For B2b PPP, the RMSE of the JFNG station on E, N, and U components are 7.3, 3.5, and 9.1 cm, respectively, while the errors of the three components at the SGOC station are 19.2, 9.2, and 21.3 cm, respectively, which are two times larger than that of JFNG station. Especially, for the RMSE of WHU PPP all eight stations present similar performance, achieving decimeter-level accuracy. For stations located on the service boundary, such as SGOC and LCK4, the RMSEs are smaller than those of B2b PPP. The RMSE values of different positioning strategies are summarized in
Figure 23.
To further evaluate the real-time performance of PPP-B2b services, the PPP with broadcast ephemeris is also carried out which is listed in
Table 8. For the static positioning mode, the RMSE of horizontal components are similar, with decimeter-level accuracy. For the kinematic positioning mode, the RMSE of horizontal component is around 1.5 m.
The convergence time of these stations is also used to analyze the positioning performance. The convergence time for the static positioning mode is defined as the positioning error less than 10 cm on the E and N components and less than 30 cm on the U component, and lasts for at least 10 min. For the kinematic positioning mode, the convergence time is defined as the error within 30 cm on the E and N components, within 60 cm on the U component, and lasts for at least 10 min. The daily convergence time over 7 days of the eight stations is summarized in
Table 9 with both B2b PPP and WHU PPP. In static positioning mode, for B2b PPP, the convergence time varies between 16.33 and 40.58 min except for the MIZU station. The convergence time of the MIZU station is 99.17 min due to the significant error on the U component at the beginning epochs, which is further caused by the considerable value of PDOP, as shown in
Figure 17. In kinematic positioning mode, B2b PPP, the convergence time varies between 11.08 and 45.66 min. The convergence time with B2b PPP is shorter than that of WHU PPP, which is mainly caused by the smaller STD on PPP-B2b clock offset.
From the results of the experiments, we found that the positioning accuracy of the selected stations with real-time PPP-B2b service can achieve centimeter-level accuracy in static positioning mode and decimeter-level accuracy in kinematic positioning mode. The overall positioning performance of PPP-B2b service within the test area is slightly better than that of WHU real-time products. The reason may be that the service areas of the PPP-B2b are consistent with those of its monitoring network, while the WHU real-time products are calculated based on the global GNSS network which is not consistent with the test area. The positioning performance of stations located at the boundary of the service area, however, is worse than that of the WHU real-time product. The PPP with broadcast ephemeris can achieve decimeter level accuracy in static positioning mode in the selected stations, which also corresponds to the smaller value of SISRE of broadcast ephemeris, which is 0.536 m.