1. INTRODUCTION
Since 1985 China has been working on its own navigation satellite system BeiDou (BDS). The system is planned to be established in three steps or phases: demonstrational system, regional system and global system (China National Space Administration, 2003; China Satellite Navigation Office, 2013). On 27 December 2012, the regional system was officially announced to as providing positioning services over the Asia-Pacific region. Now the regional system is in full operation with its designed constellation comprising five Geostationary Earth Orbit (GEO) satellites, five Inclined Geosynchronous Earth Orbit (IGSO) satellites and four Medium Earth Orbit (MEO) satellites. The standard positioning service is free for all civilian users, while the precise positioning service and system-embedded wide area differential positioning are only for authorised users.
Besides the above-mentioned positioning services, precise relative positioning has already been investigated (Shi et al., Reference Shi, Zhao, Hu and Liu2013; Montenbruck et al., Reference Montenbruck, Hauschild, Steigenberger, Hugentobler, Teunissen and Nakamura2012) and precise orbit determination and clock estimation has also been demonstrated (Montenbruck et al., Reference Montenbruck, Hauschild, Steigenberger, Hugentobler, Teunissen and Nakamura2012; Ge et al., Reference Ge, Zhang, Jia, Song and Wickert2012b; Shi et al., Reference Shi, Zhao, Li, Tang, Hu, Lou, Zhang, Niu and Liu2012; Steigenberger et al., Reference Steigenberger, Hauschild, Montenbruck and Hugentobler2012; Zhao et al., Reference Zhao, Guo, Li, Qu, Hu, Shi and Liu2013; He et al., Reference He, Ge, Wang, Wickert and Schuh2013). Precise Point Positioning (Zumberge et al., Reference Zumberge, Heflin, Jefferson, Watkins and Webb1997) based on the products has also been undertaken in post-processing and simulated in real-time mode for validating the estimated orbit and clock products (Kouba and Heroux, Reference Kouba and Heroux2001).
Ge et al. (Reference Ge, Dousa, Li, Ramatschi, Nischan and Wickert2012a) carried out BDS PPP in both static and kimematic modes, but only for a single station in the middle of the tracking network. Shi et al. (Reference Shi, Zhao, Li, Tang, Hu, Lou, Zhang, Niu and Liu2012) demonstrated BDS PPP for a single station, which was not involved in orbit and clock estimation. The orbits and clocks were determined with zenith tropospheric delay and receiver clock estimated from GPS data, instead of using BDS as an independent system. In Steigenberger et al. (Reference Steigenberger, Hauschild, Montenbruck and Hugentobler2012) data from a station in Melbourne, Australia was utilized to test BDS PPP performance in static mode with and without GEOs. Although GEOs have a lower orbit accuracy compared to the IGSOs, they significantly improve the positioning accuracy because of the limited visibility. Zhao et al. (Reference Zhao, Guo, Li, Qu, Hu, Shi and Liu2013) undertook PPP for BDS and BDS+GPS data in both static and kinematic modes. However, again only one station, National Time Service Center (NTSC), China was used and it is several tens of km away from the tracking station Xi'An. In general, in most of the studies PPP is carried out for only a single station or a few stations which are not involved in the clock and orbit estimation because of the limited data availability. Real-time kinematic PPP is widely recognized as a mode of precise positioning service, due to its cost-efficiency, global coverage and flexibility (Gao et al., Reference Gao, Chen and Shen2003; Bisnath and Gao, Reference Bisnath and Gao2008). Its performance has been further enhanced through recent developments in PPP ambiguity-fixing (Ge et al., Reference Ge, Gendt, Rothacher, Shi and Liu2008; Laurichesse et al., Reference Laurichesse, Mercier, Berthias, Broca and Cerri2008, Collins et al., Reference Collins, Lahaye, Héroux and Bisnath2008) and regional augmentation (Li et al., Reference Li, Zhang and Ge2011; Ge et al., Reference Ge, Dousa, Li, Ramatschi, Nischan and Wickert2012a).
The International GNSS Service (IGS) established in 2007 its Real-Time Pilot Project Working Group (http://www.rtigs.net/) and recently announced after about six years of fruitful cooperation its operational real-time service. Thus, BDS PPP is of great interest, especially the capacity of its real-time PPP service. Theoretically, there should not be any significant difference between GPS PPP and BDS PPP thanks to their similar signal structure and analogous frequencies. However, the special BDS constellation with high GEO and IGSO and a few MEO satellites and tracking networks of the BDS regional system will bring some new issues in observation modelling and data processing.
In this contribution, a regional network comprising six stations with an inter-station distance of about 150 km and equipped with GPS- and BDS-capable dual-frequency receivers were employed for a PPP performance study. Precise orbits and clocks estimated from a global tracking network using the strategy by He et al. (Reference He, Ge, Wang, Wickert and Schuh2013) are employed in BDS PPP processing in static, kinematic and simulated real-time kinematic mode. The static results using various session lengths and kinematic trajectories are compared with the daily GPS PPP results using IGS precise orbits and clocks in order to assess the performance of BDS PPP in terms of accuracy and convergence time. After a short description of PPP principles and data processing of the PPP service system, the tracking network and precise products are explained. Then the software package is introduced. The data processing schemata are defined and results are analysed and assessed. Finally, conclusions are drawn accordingly.
2. EXPERIMENT DESIGN
As is well known, a PPP service system includes two parts: PPP server and client. At the server-end precise orbits and clocks and some additional corrections are generated using data from a reference network and transmitted to the client-end in real-time. At the client-end, most of the biases at the local receiver can be precisely removed by using the products provided by the server. Between the server and client a data link is required for transportation of correction information.
2.1. PPP Server
In this study, the BDS Experimental Test Service (BETS) network with BDS and GPS capability serves as the reference network and was established for the scientific purpose of a Positioning, Navigation and Time (PNT) service by the GNSS Research Center at Wuhan University, China. Since March 2011, 14 stations have been deployed in China and its neighbouring regions. Among them, nine stations are located inside China and five overseas. All stations are equipped with UR240 dual-frequency GPS/BeiDou receivers and UA240 antennas manufactured by the UNICORE Company in Beijing, China (http://www.unicorecomm.com/english/). A number of studies have been carried out based on data from this network (Ge et al., Reference Ge, Dousa, Li, Ramatschi, Nischan and Wickert2012a; Shi et al., Reference Shi, Zhao, Li, Tang, Hu, Lou, Zhang, Niu and Liu2012; Zhao et al., Reference Zhao, Guo, Li, Qu, Hu, Shi and Liu2013; He et al., Reference He, Ge, Wang, Wickert and Schuh2013).
Three-day data are used in a batch estimation to obtain a three-day solution, instead of combining three daily solutions on the level of normal equations (Steigenberger et al., Reference Steigenberger, Hugentobler, Montenbruck and Hauschild2011). In order to assess the quality of the estimated clocks and orbits, the differences over the overlapping time of two adjacent three-day solutions are utilized as usual. The orbit of the last day in a three-day solution is compared with that of the middle day of the next three-day solution. The RMS of the differences in along-track, cross-track and radial directions are taken as an orbit quality indicator. The statistical results for orbits and clocks are listed in Table 1.
Overlapping RMS of the estimated orbits and clocks for each individual satellites and the mean of each satellite type.
From the orbit RMS in Table 1, the along-track RMS is significantly larger than that of the other two directions as expected. GEOs have the largest RMS in the along-track direction of 114 cm, compared to 24 cm and 45 cm for IGSOs and MEOs, respectively. RMS in cross-track and radial are very similar for the three types of satellites, i.e., 10 cm and 6 cm, respectively, for GEOs, 15 cm and 7 cm, respectively, for IGSOs, and 13 cm and 12 cm, respectively, for MEOs.
2.2. PPP Client
Aside from the above mentioned network, a regional network comprising six different stations (Xu et al., Reference Xu, Xu, Ge, Xu, Zhu and Sui2013) equipped with the same UNICORE receivers and antenna were also employed as user stations for the tests of various precise point positioning modes. The inter-station distance is about 150 km on average. The distribution of the stations and their locations with respect to the BETS network are shown in Figure 1.
The regional BDS/GPS test network within the context of the global BDS reference network in China.
3. DATA PROCESSING
3.1. Software package
The PANDA (Positioning And Navigation Data Analyst) software is used for the data processing. The software was developed at the GNSS Research Center in Wuhan University (Liu and Ge, Reference Liu and Ge2003; Shi et al., Reference Shi, Zhao, Geng, Lou, Ge and Liu2008) and has recently been adapted for BDS data processing (Shi et al., Reference Shi, Zhao, Li, Tang, Hu, Lou, Zhang, Niu and Liu2012; Zhao et al., Reference Zhao, Guo, Li, Qu, Hu, Shi and Liu2013; He et al., Reference He, Ge, Wang, Wickert and Schuh2013). Similar to other scientific GNSS software packages, it includes the following basic modules: data pre-processing, orbit integration, parameter estimation, data editing and ambiguity-fixing. Two estimators have been developed: a least square estimator for post-mission processing and a square root information filter estimator for real-time processing.
3.2. Models
The undifferenced ionosphere-free code and phase observations are employed in the PPP processing. The a priori STD of B1 and B2 phase and range are 0·01 cycles and 1·0 m, respectively. Observations are down-weighted by factor of 2*sin(E) if their elevations are lower than 30° and are excluded if their elevations are below 7°. Satellite orbits and clocks and Earth rotation parameters are fixed to the products estimated from the BETS network at the server-end. The relativistic delay and station displacement caused by solid earth tides, pole tides and ocean tide loading are corrected following the instruction or using the models and subroutines recommended in IERS Convention 2003. The phase windup is calculated and corrected using the formulae in Wu et al. (Reference Wu, Wu, Hajj, Bertiger and Lichten1993).
Tropospheric delays are first corrected through ZTD computed using the Saastamoinen model and meteorological data from the Global Pressure and Temperature model (GPT) (Boehm et al., Reference Boehm, Heinkelmann and Schuh2007), as well as the Vienna Mapping Function (VMF) (Boehm et al., Reference Boehm, Werl and Schuh2006). The remaining ZTD are modeled by a random-walk stochastic process with a power density of 2·5 
$mm/\sqrt h $ (Bar-Sever and Kroger, Reference Bar-Sever and Kroger1998). Receiver clock error is estimated as white noise. Ambiguities are estimated as real-values and integer ambiguity resolution is not performed. Station coordinates are estimated for static mode as determined parameters with a loose constraint of about 100 m for each component and for kinematic mode as a first Gaussian-Markov process with 100 m constraint on initial coordinates and a power density of about 1 m/s.
Phase Centre Offset (PCO) and Phase Centre Variation (PCV) of BDS satellites and receivers are not currently available. Therefore, PCO and PCV corrections are neglected. Satellite attitude control mode has also not been announced. Moreover, the attitude control mode of the BDS satellites is assumed to be the same as that of the GPS Block IIR satellites (He et al., Reference He, Ge, Wang, Wickert and Schuh2013).
3.3. Processing Schemata
Three processing scenarios are defined for assessing the performance of the potential PPP services: static mode with various lengths of observations, kinematic positioning in post-mission mode, and simulated real-time kinematic positioning.
In static mode, data from the three selected days are divided into sessions of various lengths; estimated station coordinates are compared with the multi-day GPS solution; and their repeatability is also utilized as a quality indicator for their internal consistency. We used session lengths of 2 h, 3 h, 6 h, 12 h and 24 h.
Kinematic solutions with daily data are carried out in post-processing mode. The station coordinates are modelled by a first order Gaussian-Markov process. In the least squares adjustment, the state equation of the process is added as pseudo observation equations into the normal equation system. Station coordinates at the previous epochs are eliminated and they are recovered backwards after the parameters at the last epoch are solved. Therefore the estimated coordinates are equivalent to those from a total least squares adjustment using all of the observations – a smoothed solution. We tried to investigate the capability of the BDS real-time PPP service but in simulated real-time mode.
We used estimates of a three-day solution of precise orbit determination as the basis to predict the orbits of the day after. With the predicted orbits as known values, satellite clocks of the next day were estimated epoch-by-epoch with fixed station coordinates. Simulated real-time positioning using the estimated orbits and clocks was performed to validate the accuracy of BDS real-time positioning service.
As is well known, the convergence time of real-time kinematic PPP is a key indicator of its performance. In order to estimate the convergence time, we divided the data into two-hour sessions and each is processed in simulated kinematic mode. For each station there are 12 solutions every day and 216 solutions in total for six stations over three days. The estimated trajectories are compared with the known coordinates. The RMS of the coordinate differences at the same epoch number of all the solutions are computed to illustrate the relationship between length of observations and position accuracy.
4. RESULTS
Data from days 311–313 of 2012 of the six regional stations are selected as user stations in the three processing scenarios defined above. The estimated station coordinates are compared with those estimated from GPS data using IGS final products. The details are discussed in the following.
4.1. Static PPP
In the static test, PPP is performed with session lengths of two hours, three hours, six hours, 12 hours and 24 hours. Figure 2 shows the averaged RMS of the coordinate differences in east, north and up directions of the static solutions with various session lengths over the three selected days.
RMS of the positioning differences in east, north and up directions of the static PPP with various session lengths over the three selected days.
The accuracy of the PPP derived positions with two hours data is about 6 cm, 4 cm and 7 cm in east, north and up directions, respectively, compared to the GPS daily solution. The accuracy is improved along with the session lengths. If data of six hours or longer are involved in the processing, position accuracy in the east and north is about 1 cm, while there is a bias of about 3 cm in the up direction.
Figure 3 shows the time series of the differences for station HSSZ as an example. The three panels are for days 311, 312 and 313 from top to bottom, respectively. In each panel, all solutions of the 3 h, 6 h, 12 h and 24 h sessions are presented from left to right. The solution of the 3 h sessions is rather noisy and there is clearly a systematic bias in the vertical component. The bias in the vertical is most likely caused by the un-modelled PCO and PCV of the receiver and satellite antennas as well. Hence, vertical bias is present in all arc lengths.
Position differences of static PPP solutions with session length of 3 hours, 6 hours, 12 hours and 24 hours for days 311, 312 and 313 compared to the estimates using daily GPS data for station HSSZ.
4.2. Kinematic PPP
Kinematic PPP is also applied to the regional stations using the same orbit and clock products as for static PPP for the days 311–313 in 2012. The positions are estimated by means of the sequential least square adjustment with a loose constraint of 1 m to positions at two adjacent epochs. The epoch-wise station coordinates are eliminated at each epoch while processing forwards and are recovered backwards (Ge et al., Reference Ge, Gendt, Dick, Zhang and Rothacher2006).
The derived trajectories are compared to the daily PPP solution using GPS data. The bias and STD of the differences in east, north and up are shown in Figure 4 for the six stations. Both bias and STD are below 10 mm for east component, 10 mm and 15 mm for north component, respectively, about −40 mm and 50 mm for vertical component, respectively.
Bias (left) and STD (right) of the derived trajectories of the six regional stations in kinematic mode.
As an example, position differences of the trajectory for the six stations on day of year (DOY) 312 in 2012 are displayed in Figure 5 and the number of satellites used in each epoch is also shown. There is an obvious variation pattern in the vertical time series that differs from station to station. The reason is still under investigation and it could potentially be the multi-path effect or any inaccurate modelling.
Position differences of the kinematic PPP and number of satellites observed for the six regional stations on DOY 312 in 2012.
4.3. Simulated Real-Time PPP service
Figure 6 presents the position differences of the simulated real-time PPP service and the static daily GPS solution. The number of observed satellites is also plotted. Compared to the real-time PPP result in Figure 6 and the post-processing result in Figure 5, a long convergence time of about 1·5 hours is needed for real-time PPP in order to get positions of about 10 cm accuracy. Afterwards, the accuracy stays at the sub-dm-level.
Position differences of the simulated real-time and the static daily PPP for the six stations on DOY 311 in 2012.
According to the processing schema, for each epoch number, there are 216 estimated coordinates. Their differences with respect to the known coordinates reflect the position quality with the data up to this epoch. Therefore, the RMS of all the 216 coordinate differences can be a measure of the position accuracy. Figure 7 shows the relationship between the length of observations and the accuracy of the estimated positions. From Figure 7, we have similar accuracy for the east and up components although usually the east should be much better than the up component. After one hour convergence, the accuracy is about 8 cm, 20 cm and 20 cm for the north, east and up directions, respectively. With 90 minutes data, horizontal accuracy is better than 10 cm on average in the north and east components and about 15 cm for the vertical. The accuracy is still improved along with observing time and reaches to better than 10 cm for all components at the end of the two-hour session.
Convergence of station position along with the length of observations.
The bias and STD of the differences in simulated real-time PPP after 1·5 hours are shown on Figure 8 for the six stations. The bias is below 10 mm for east and north components, about −30 mm for the vertical component. The STD is about 20 mm for east and north components, about 60 mm for the vertical component. Therefore, with the BDS regional system, positioning accuracy of better than one decimetre is achievable in real-time PPP mode.
Bias (left) and STD (right) of the derived trajectories of the six stations in simulated real-time mode after 1·5 hours.
5. SUMMARY
The BDS navigation satellite system being established in China is moving steadily forward towards its final destination – an operational global navigation satellite system by 2020. The three-phase schedule enables the rapid emerging of the BDS system with its operational service over the Asia-pacific region now.
With the precise orbit and clock estimated from the BETS tracking network based on Chinese stations and several stations worldwide, PPP in static, kinematic and simulated real-time modes using BDS data only have been carried out and assessed.
In the static test, an accuracy of a few cm is attainable with about two-hour data compared to the GPS daily PPP solution. The accuracy in terms of RMS increases along with the observation time and reaches to better than 1 cm in the horizontal and 4 cm in the vertical with six hours data or longer. A clear bias changing along with observing time is found in the vertical direction, potentially being caused by the un-modelled PCO and PCV of the receiver and satellite antennas.
For kinematic post-positioning, the horizontal accuracy in STD is about 1 to 2 cm, but about 5 cm for vertical but with a bias of −4 cm which is about 6 cm in RMS. There are clearly long-term variations in the position time series, especially in the vertical direction, compared to the ground truth.
Simulated real-time PPP shows that after a convergence time of about 60 minutes, position accuracy reaches to 10 cm, 20 cm and 20 cm in north, east and up, respectively. On average, a convergence time of 90 minutes is needed for obtaining positions with an accuracy of about 10 cm. The accuracy of the converged simulated real-time kinematic positioning is 2 to 3 cm in the horizontal and 6 to 7 cm in the vertical.
Therefore, the operational BDS system can provide very precise positioning services similar to GPS and GLONASS, but only to China and its surroundings. Although the performance in terms of convergence time and accuracy is still slightly worse than GPS, these metrics can potentially be improved by applying accurate satellite and receiver PCO and PCV corrections and satellite attitude information.
ACKNOWLEDGEMENTS
The authors would like to thank the GNSS Research Center at Wuhan University and the Chinese Academy for Surveying and Mapping for providing the BDS data of the BETS network and the regional network for this study. This research is supported by the Liaoning Province Excellent Talent Program (LR201018), the PhD Programs Foundation of MOE of China (2012212111007), and the National Natural Science Foundation of China (41201454).
