Search

To support real-time global navigation satellite systems (GNSS) precise applications, satellite clock corrections need to be precisely estimated at a high-rate update interval, which remains a challenge due to the rapid development of multi-GNSS constellations. In this study, we developed an undifferenced (UD) ambiguity resolution (AR) procedure to improve both the accuracy and computational efficiency for real-time multi-GNSS clock estimation realized by a square root information filter. In the proposed method, UD ambiguities are resolved after correcting the simultaneously estimated uncalibrated phase delays (UPD) and the fixed UD ambiguity parameters are eliminated immediately from the filter, so that the computational burden is significantly reduced. Moreover, based on the linear relationship between double-differenced (DD) and UD ambiguities, we investigated the difference between DD and UD AR in clock estimation. We found that the major reason why DD AR contributes little to the clock estimation while UD AR can speed up the convergence remarkably is that UD AR additionally provides a stable clock datum compared with DD AR. GNSS observations from about 100 globally distributed stations were processed with the proposed method to generate simulated real-time clocks and UPDs for GPS, Galileo, and BDS satellites over a one-month period. The results show that the percentage of wide-lane (WL) UPD residuals within ± 0.25 cycles and narrow-lane (NL) UPD residuals within ± 0.15 cycles are over 97.0% and 90.0%, respectively, which contributes to an ambiguity fixing rate of more than 90% for three systems. The mean daily standard deviation (STD) of the clocks of the UD-fixed solution with respect to Center for Orbit Determination in Europe 30 s final products is 0.021, 0.020, and 0.035 ns for GPS, Galileo, and BDS satellite, respectively, which is improved by 78.1%, 58.3%, and 79.8% compared to the float solution. Benefiting from the removal of fixed ambiguities, the average computation time per epoch was reduced from 3.88 to 1.05 s with a remarkable improvement of 72.9%. The quality of the satellite clock and UPD products was also evaluated by the performance of kinematic precise point positioning (PPP). The results show that fast and reliable multi-GNSS PPP-AR can be achieved with the derived UD-fixed clocks and UPDs, which outperforms that using DD-fixed clock and off-line UPD products with an average improvement of 7.9% and 19.9% in terms of convergence time and positioning accuracy, respectively. Furthermore, we demonstrated the effectiveness of the proposed UD AR method through a 7-day real-time clock estimation experiment. [ABSTRACT FROM AUTHOR]

Orbit prediction (OP) recently tends to be a very crucial step for supporting real-time GNSS orbit services due to the dynamic stability of navigation satellite orbits. The OP performance depends on the length of the predicted orbits and the accuracy of precise orbit determination (POD) as basis. Considering this, a new automatic processing engine is established for improving the multiple global navigation satellite systems (multi-GNSS) constellation OP performance. From the architecture-oriented high-performance parallel processing perspective, the multi-node and multi-core computer sources are fully exploited to implement the hourly update of the current multi-GNSS POD. For MEO-type satellites (e.g., Galileo satellites), the accuracy of predicted orbits is improved from 3.8 cm, 6.5 cm, and 12.3 cm to 3.5 cm, 4.3 cm, and 6.3 cm, in the radial, cross, and along directions, respectively, compared to the three-hour POD update. Despite the shortened OP length, the OP performance of regional navigation satellite system (RNSS) satellites is still limited due to their regional observability. The BDS-IGSO and QZSS-IGSO satellitesexhibit radial directional orbital errors of up to 36.9 cm and 28.9 cm, respectively. Therefore, an orbit fitting (OF) processing method with orbit reconstruction is implemented into the processing engine. By utilizing this method, the radial orbital errors for BDS-IGSO and QZSS-IGSO satellites can be reduced to 7.0 cm and 10.4 cm, respectively. The mean real-time positioning errors are thus reduced from 28.3 to 18.4 cm and from 24.4 to 18.2 cm in the horizontal and vertical components, respectively. [ABSTRACT FROM AUTHOR]

High-precision Global Navigation Satellite Systems (GNSS) orbits are critical for real-time clock estimation and precise positioning service; however, the prediction error grows gradually with the increasing prediction session. In this study, we present a new efficient precise orbit determination (POD) strategy referred to as the epoch-parallel processing to reduce the orbit update latency, in which a 24-h processing job is split into several sub-sessions that are processed in parallel and then stacked to solve and recover parameters subsequently. With a delicate handling of parameters crossing different sub-sessions, such as ambiguities, the method is rigorously equivalent to the one-session batch solution, but is much more efficient, halving the time-consuming roughly. Together with paralleling other procedures such as orbit integration and using open multi-processing (openMP), the multi-GNSS POD of 120 satellites using 90 stations can be fulfilled within 30 min. The lower update latency enables users to access orbits closer to the estimation part, that is, 30–60-min prediction with a 30-min update latency, which significantly improves the orbit quality. Compared to the hourly updated orbit, the averaged 1D RMS values of predicted orbit in terms of overlap for GPS, GLONASS, Galileo, and BDS MEO are improved by 39%, 35%, 41%, and 37%, respectively, and that of BDS GEO and IGSO satellites is improved by 47%. We also demonstrate that the boundary discontinuities of half-hourly orbit are within 2 cm for the GPS, GLONASS, and Galileo satellites, and for BDS the values are 2.6, 15.5, and 9.8 cm for MEO, GEO, and IGSO satellites, respectively. This method can also be implemented for any batch-based GNSS processing to improve the efficiency. [ABSTRACT FROM AUTHOR]

For the precise orbit determination (POD) of global navigation satellite systems (GNSS) constellation, it is very difficult to precisely model the solar radiation pressure (SRP) force acting on GNSS satellites. For GPS satellites, the ECOM model developed by the Center for Orbit Determination in Europe has been utilized by most of International GNSS Service (IGS) analysis centers. However, it should be adapted and optimized to the characteristics of satellites of each GNSS system or even individual satellites. It was extended to the ECOM2 model for GLONASS satellites and then for Galileo satellites by employing a box–wing model. Since November 2020, the third generation of the BeiDou satellite system (BDS-3) has been in its full operation and there are about 200 globally distributed IGS ground stations tracking BDS-3 signals, which creates a great potential to evaluate and optimize its SRP modeling. From the POD processing carried out in this study, we found significant fluctuations of up to 20 cm in overlapping orbit differences for satellites over eclipses in the radial direction and of about 20 and 50 cm in the cross and along directions for ECOM2 and ECOM models. Then, based on numerical analyses we demonstrate that the fourth- and sixth-order sine terms in the Sun direction can significantly reduce the overlapping orbit differences of ECOM. Therefore, an adapted SRP model by adding the fourth- and sixth-order sine periodical terms in the Sun direction to the ECOM model is presented. The adapted model is then validated for BDS-3 POD and orbit prediction. Results show that fluctuations in the amplitude of overlapping estimated orbits using ECOM models are reduced from 20 to < 10 cm in the radial-track component and satellite laser ranging residuals are reduced to half by the adapted SRP model. For the predicted BDS-3 satellite orbits, the RMS values over deep eclipses can be improved from about 7, 14 and 26 cm to about 3, 5 and 12 cm, in the radial, cross and along directions, respectively, compared to the ECOM model. [ABSTRACT FROM AUTHOR]

Global Navigation Satellite Systems (GNSS) Real-time Precise Point Positioning (RT-PPP) strongly relies on the precise satellite orbits and clocks, especially the latter requires high update rate, e.g., five seconds, due to its limited prediction accuracy. Such a high-rate update frequency is a burden to both the data analysis and communicating, and interruption in communicating is almost unavoidable. For the new constellations such as Galileo and BDS-3 with high-stability hydrogen clocks onboard, it is possible to predict satellite clocks to a longer arc. Moreover, recent developments in multi-GNSS ultra-rapid precise orbit determination enables the half-hourly update, ensuring the availability with a prediction arc of 30–60 min. We investigate multi-GNSS RT-PPP using half-hourly predicted products and demonstrate that a 3-D accuracy of 2.9 and 11.3 cm can be achieved for static and simulated kinematic solutions, respectively. We present the different clock prediction accuracies of different types of satellites and propose a satellite-specific weighting strategy in PPP, which exploits the benefits of the satellites of good performance. The method is based on the prediction accuracy of both, satellites and clocks, and shows an improvement of 15 to 60% compared to those without satellite-specific weighting or with simplified weighting strategies. We also demonstrate that Galileo satellites contribute the most in the quad-constellation solution, thanks to the highly stable satellite clocks. [ABSTRACT FROM AUTHOR]

Introducing precise atmosphere information into precise point positioning enables rapid ambiguity resolution and introduces a significant accuracy improvement. However, it can only be implemented in regions with dense networks and stable communication links. For larger areas, e.g., an intercontinental level, there is a conflict between the accuracy of corrections and the amount of atmosphere information to be disseminated. We develop a hierarchical augmentation mode to combine the advantages of the fitting model and region interpolation model to relieve the communication burden. Relying on the fitting model with fewer coefficients applied over large areas as the essential information, the unmodeled errors are calculated at each reference station, and further correction information is optional compensation depending on the magnitude of the unmodeled residuals. We perform the proposed models on 103 EUREF Permanent Network stations with 200-km station spacing and 84 stations as the external validation. The ionosphere and troposphere fitting models have an average accuracy of about 4.2 and 1.3 cm, respectively, under meteorologically calm conditions. The unmodeled error transmission determined by the magnitude of residuals can be reduced by 61% and 96% for the ionospheric and tropospheric delays, respectively, with respect to the legacy interpolation mode. Further compensation implemented, i.e., unmodeled residuals, can achieve instantaneous convergence for 83.6% of all solutions, and the overall initialization time is within 1.0 min. Thus, the proposed hierarchical positioning mode satisfies real-time positioning convergence requirements and significantly reduces massive corrections in communication. [ABSTRACT FROM AUTHOR]

Since the orbit repeat cycle of Galileo satellites is about 10 days, existing spatial–temporal repeatability-based station multipath mitigation methods such as the modified sidereal filtering (MSF) and the multipath hemispherical map need more than one week of data in order to model the Galileo phase multipath correction. From this perspective, these methods are improper for most applications. We proposed a random walk multipath method (RWM) to resolve this problem. Since the elevation and azimuth angles between adjacent epochs do not change too much, the low-frequency phase multipath effects at adjacent epochs are similar. Multipath correction value can be estimated by the random walk model. Because the GPS satellite repeat cycle is short, GPS observations easily mitigate the multipath effect and have common coordinate parameters with Galileo. Multipath-reduced GPS signals can be a constant for separating coordinate parameters and multipath parameters. Galileo phase observation residuals of the latest day are used to calculate the variance for the random walk model. Experiment results show that compared to the traditional MSF model, the proposed RWM method can improve the Galileo residual reduction and positioning precision nearly without the need of any historical data. As for the practical real-time GNSS monitoring application in a severe multipath environment, the result shows that the new method can significantly reduce the Galileo multipath effect and subsequently yield a precise positioning solution. Moreover, the RWM method is invulnerable to the sampling rate and observation environment. [ABSTRACT FROM AUTHOR]

The computational efficiency is critical with the increasing number of GNSS satellites and ground stations since many unknown parameters must be estimated. Although only active parameters are kept in the normal equation in sequential least square estimation, the computational cost for parameter elimination is still a heavy burden. Therefore, it is necessary to optimize the procedure of parameter elimination to enhance the computational efficiency of GNSS network solutions. An efficient parallel algorithm is developed for accelerating parameter estimation based on modern multi-core processors. In the parallel algorithm, a multi-thread guided scheduling scheme, and cache memory traffic optimizations are implemented in parallelized sub-blocks for normal-equation-level operations. Compared with the traditional serial scheme, the computational time of parameter estimations can be reduced by a factor of three due to the new parallel algorithm using a six-core processor. Our results also confirm that the architecture of computers entirely limits the performance of the parallel algorithm. All the parallel optimizations are also investigated in detail according to the characteristics of CPU architecture. This gives a good reference to architecture-oriented parallel programming in the future development of GNSS software. The performance of the multi-thread parallel algorithm is expected to improve further with the upgrade of new multi-core coprocessors.s [ABSTRACT FROM AUTHOR]

Tropospheric delay modeling is challenging in high-precision Very Long Baseline Interferometry (VLBI) analysis due to the rapid water vapor variation and imperfect observation geometry, where observations from Global Navigation Satellite Systems (GNSS) co-locations can enhance the VLBI analysis. We investigate the impact of tropospheric ties in the VLBI and GNSS integrated processing during the CONT05–CONT17 campaigns, and present a method that automatically handles the systematic tropospheric tie biases. Applying tropospheric ties at VLBI–GNSS co-locations enhances the observation geometry and improves the solution reliability. The VLBI network is stabilized, with station coordinate repeatability improved by 12% horizontally and by 28% vertically, and the network scale improved by 32%. The Earth Orientation Parameters (EOP) improve by up to 20%. Both zenith delay and gradient ties contribute to the improvement of EOP, whereas the gradient ties contribute mainly to the improvement of length of day and celestial pole offsets. [ABSTRACT FROM AUTHOR]

The performance of precise point positioning (PPP) has been significantly improved thanks to the continuous improvements in satellite orbit, clock, and ambiguity resolution (AR) technologies, but the convergence speed remains a limiting factor in real-time PPP applications. To improve the PPP precision and convergence time, tropospheric delays from a regional network can be modeled to provide precise correction for users. We focus on the precise modeling of zenith wet delay (ZWD) over a wide area with large altitude variations for improving PPP-AR. By exploiting the water vapor exponential vertical decrease, we develop a modified optimal fitting coefficients (MOFC) model based on the traditional optimal fitting coefficients (OFC) model. The proposed MOFC model provides a precision better than 1.5 cm under sparse inter-station distances over the Europe region, with a significant improvement of 70% for high-altitude stations compared to the OFC model. The MOFC model with different densities of reference stations is further evaluated in GPS and Galileo kinematic PPP-AR solutions. Compared to the PPP-AR solutions without tropospheric delay augmentation, the positioning precision of those with 100-km inter-station spacing MOFC and OFC is improved by 25.7% and 17.8%, respectively, and the corresponding time to first fix (TTFF) is improved by 36.9% and 33.0% in the high-altitude areas. On the other hand, the OFC model only slightly improves the TTFF and positioning accuracy when using the 200 km inter-station spacing modeling and even degrades the positioning for high-altitude stations, whereas using the MOFC model, the PPP-AR solutions always improve. Moreover, the positioning precision improvement of MOFC compared with OFC is about 22.1%, 21.7%, and 25.7% for the Galileo-only, GPS-only, and GPS + Galileo PPP-AR solutions, respectively. [ABSTRACT FROM AUTHOR]

The article Integrated processing of ground- and space-based GPS observations: improving GPS satellite orbits observed with sparse ground networks, written by Wen Huang, Benjamin Männel, Pierre Sakic, Maorong Ge & Harald Schuh was originally published electronically on the publisher’s internet portal on [10 October 2020] without open access. [ABSTRACT FROM AUTHOR]

Wide-lane (WL) uncalibrated phase delay (UPD) is usually derived from Melbourne–Wübbena (MW) linear combination and is a prerequisite in Global Navigation Satellite Systems (GNSS) precise point positioning (PPP) ambiguity resolution (AR). MW is a linear combination of pseudorange and phase, and the accuracy is limited by the larger pseudorange noise which is about one hundred times of the carrier phase noise. However, there exist inconsistent pseudorange biases which may have detrimental effect on the WL UPD estimation, and further degrade user-side ambiguity fixing. Currently, only the large part of pseudorange biases, e.g., the differential code bias (DCB), are available and corrected in PPP-AR, while the receiver-type-dependent biases have not yet been considered. Ignoring such kind of bias, which could be up to 20 cm, will cause the ambiguity fixing failure, or even worse, the incorrect ambiguity fixing. In this study, we demonstrate the receiver-type-dependent WL UPD biases and investigate their temporal and spatial stability, and further propose the method to precisely estimate these biases and apply the corrections to improve the user-side PPP-AR. Using a large data set of 1560 GNSS stations during a 30-day period, we demonstrate that the WL UPD deviations among different types of receivers can reach ± 0.3 cycles. It is also shown that such kind of deviations can be calibrated with a precision of about 0.03 cycles for all Global Positioning System (GPS) satellites. On the user side, ignoring the receiver-dependent UPD deviation can cause significant positioning error up to 10 cm. By correcting the deviations, the positioning performance can be improved by up to 50%, and the fixing rate can also be improved by 10%. This study demonstrates that for the precise and reliable PPP-AR, the receiver-dependent UPD deviations cannot be ignored and have to be handled. [ABSTRACT FROM AUTHOR]

The precise orbit determination (POD) of Global Navigation Satellite System (GNSS) satellites and low Earth orbiters (LEOs) are usually performed independently. It is a potential way to improve the GNSS orbits by integrating LEOs onboard observations into the processing, especially for the developing GNSS, e.g., Galileo with a sparse sensor station network and Beidou with a regional distributed operating network. In recent years, few studies combined the processing of ground- and space-based GNSS observations. The integrated POD of GPS satellites and seven LEOs, including GRACE-A/B, OSTM/Jason-2, Jason-3 and, Swarm-A/B/C, is discussed in this study. GPS code and phase observations obtained by onboard GPS receivers of LEOs and ground-based receivers of the International GNSS Service (IGS) tracking network are used together in one least-squares adjustment. The POD solutions of the integrated processing with different subsets of LEOs and ground stations are analyzed in detail. The derived GPS satellite orbits are validated by comparing with the official IGS products and internal comparison based on the differences of overlapping orbits and satellite positions at the day-boundary epoch. The differences between the GPS satellite orbits derived based on a 26-station network and the official IGS products decrease from 37.5 to 23.9 mm ( 34 % improvement) in 1D-mean RMS when adding seven LEOs. Both the number of the space-based observations and the LEO orbit geometry affect the GPS satellite orbits derived in the integrated processing. In this study, the latter one is proved to be more critical. By including three LEOs in three different orbital planes, the GPS satellite orbits improve more than from adding seven well-selected additional stations to the network. Experiments with a ten-station and regional network show an improvement of the GPS satellite orbits from about 25 cm to less than five centimeters in 1D-mean RMS after integrating the seven LEOs. [ABSTRACT FROM AUTHOR]

Along with the rapid development of GNSS, not only BeiDou, but also Galileo, and the newly launched GPS satellites can provide signals on three frequencies at present. To fully take advantage of the multi-frequency multi-system GNSS observations on precise point positioning (PPP) technology, this study aims to implement the triple-frequency ambiguity resolution (AR) for GPS, Galileo, and BeiDou-2 combined PPP using the raw observation model. The processing of inter-frequency clock bias (IFCB) estimation and correction in the context of triple-frequency PPP AR has been addressed, with which the triple-frequency uncalibrated phase delay (UPD) estimation is realized for real GPS observations for the first time. In addition, the GPS extra-wide-line UPD quality is significantly improved with the IFCB correction. Because of not being contaminated by the IFCB, the raw UPD estimation method is directly employed for Galileo which currently has 24 satellites in operation. An interesting phenomenon is found that all Galileo satellites except E24 have a zero extra-wide-lane UPD value. With the multi-GNSS observations provided by MGEX covering 15 days, the positioning solutions of GPS + Galileo + BeiDou triple-frequency PPP AR have been conducted and analyzed. The triple-frequency kinematic GNSS PPP AR can achieve an averaged 3D positioning error of 2.2 cm, and an averaged convergence time of 10.8 min. The average convergence time can be reduced by triple-frequency GNSS PPP AR by 15.6% compared with dual-frequency GNSS PPP AR, respectively. However, the additional third frequency has only a marginal contribution to positioning accuracy after convergence. [ABSTRACT FROM AUTHOR]

Multipath effect on carrier-phase observation is one of the bottlenecks for mm-level applications when using precise point positioning (PPP). Hence, we extract the multipath directly from raw carrier-phase residuals of GPS, GLONASS, Galileo, and BDS, by using PPP technique. Although the residuals for one frequency assimilate the errors from other frequencies, which is caused by error adjustment by the least squares estimator, the primary component of residuals is multipath. The results indicate that the residuals between frequencies have a significant linear negative correlation and synchronous time lag for each system. Besides BDS Geostationary Earth Orbit satellites, the residuals for other satellites can establish accurate mathematic relationship between the frequencies. For GLONASS, the residuals of R1 frequency recovered from R2 frequency with the mathematical relationship are better than 0.1 mm accuracy, which means the effect of inter-frequency bias can be neglected. These regularities double-reduce the complexity of data processing. Based on the multipath distribution, we propose a modified Multipath Hemispherical Map model (M-MHM), which constructs grids from residuals and is divided into three equal-elevation angle parts with an optimal resolution 0.2° × 0.2° × 1° from numerous experiments. In addition, the multipath manifests great consistency among satellites for GPS, GLONASS, and Galileo systems when elevation angles are higher than 15°, while is more satellite dependent for BDS. Although GPS L1 frequency is identical to Galileo E1, the model still has some systematic bias between GPS and Galileo. Compared with sidereal filtering and original MHM model, the M-MHM brings the highest improvement in both residual variance reduction and positioning accuracy. The positioning accuracy is on average 12% improvement compared to MHM and 29% improvement compared to SF. For four systems combined solutions with the M-MHM model, can reach an accuracy of 0.75, 0.55, and 2.08 cm in the east, north, and up components. [ABSTRACT FROM AUTHOR]

The correction of higher-order ionospheric (HOI) delays remaining in the dual-frequency ionosphere-free combined observations is suggested after the confirmation of its impact on precise Global Navigation Satellite System (GNSS) data processing. However, in the precise orbit determination (POD) for generating ultra-rapid orbits, the higher-order corrections are not always considered most likely because a RT ionospheric model needed for calculating the higher-order corrections is hardly available or the HOI impact is believed rather small compared to the accuracy of the predicted orbit. With the increasing requirement on the positioning performances from various applications, providing more accurate and reliable ultra-rapid orbits becomes an essential task of the real-time GNSS precise positioning services. In this contribution, the temporal–spatial characteristics of HOI effects on GNSS observables are investigated thoroughly using data collected from International GNSS Service (IGS) global ground stations and fluctuations of the higher-order delays up to several centimeters are detected during periods of high ionospheric activity. Hereafter, we evaluate the HOI effects on the multi-GNSS POD based on a network with globally distributed IGS stations. Results show that owing to the applied HOI corrections, the agreement of overlapping orbits can be improved significantly for all satellites and especially in radial direction. The three-dimensional RMS values of the overlapping differences are reduced from 1.6, 2.0, 4.6 and 1.7 to 1.0, 1.1, 3.4, and 1.5 cm for GPS, GALILEO, BDS, and GLONASS, respectively. Furthermore, the orbit improvement is also confirmed by the satellite laser ranging (SLR) observations over a 2-month time period where the STD of SLR residuals is reduced by HOI corrections from 6.4 to 5.3 cm for the BDS-IGSO satellites. [ABSTRACT FROM AUTHOR]

In multi-GNSS integration, fixing inter-system double-difference ambiguities to integers is still a challenge due to the existence of inter-system biases (ISB) when mixed types of GNSS receivers are used. It has been shown that when ISB is known, the inter-system ambiguities can be fixed and the reliability of ambiguity fixing can be improved significantly, especially under poor conditions when the number of observed satellites is small. In traditional methods, the intra-system ambiguity is fixed first; then, the ISB is estimated to ultimately fix the inter-system ambiguity. In our work, we use the particle filter-based method to estimate the ISB parameter and fix the inter-system ambiguities to integers at the same time. This method shows higher reliability and higher ambiguity fixing rate. Nevertheless, the existing particle filter approach for ISB parameter estimation is a one-dimensional algorithm. When satellites from three or more systems are observed, there are two or more ISB parameters. We extend the current one-dimensional particle filter approach to multi-dimensional case and estimate multi-ISB parameters in this study. We first present a multi-dimensional particle filter approach that can estimate multi-ISB parameters simultaneously. We also show that the RATIO values can be employed to judge the quality of multi-dimensional ISB values. Afterward, a two-dimensional particle filter approach is taken as an example to validate this approach. For example, in the experiment of GPS L5, Galileo E5a and QZSS L5 integration with 6 satellites using the IGS baseline SIN0-SIN1, only three ambiguities are resolved to integer when the ISBs are unknown. The integer ambiguity fixing rate is 41.0% with 53% of the ambiguity-fixed solutions having positioning errors larger than 3 cm. However, when our approach is adopted, the number of integer ambiguity parameters increases to five. The integer ambiguity fixing rate increases to 99.7% with 100% of ambiguity-fixed solutions having positioning errors smaller than 3 cm. [ABSTRACT FROM AUTHOR]

Real-time satellite orbit and clock product is a key prerequisite for real-time precise positioning service based on precise point positioning (PPP). With the rapid development of the multiple global navigation satellite systems (Multi-GNSS), about 120 satellites will be processed for Multi-GNSS real-time clock estimation. Unfortunately, the computation is very time-consuming, especially for quality control since problematic observations are inevitable. Taking advantage of computer technology, sequential least square adjustment with an adapted online quality control procedure is developed to rapidly estimate Multi-GNSS real-time clocks, although various filtering estimators are commonly used now. A globally distributed network including 70 stations tracking mostly satellites of GPS, GLONASS, BDS, and Galileo is employed for experimental validation. The results show that the computation time per epoch is less than 3 s in average and can meet the 5 s update rate of the IGS real-time clock product. Compared to the GeoForschungsZentrum MGEX (GBM) final clock product, the averaged STD values of the estimated clocks of the four GNSS systems are 0.089 ns and 0.153 ns, respectively, for the clock solutions with and without the online quality control, which also confirms the importance of the quality control procedure. The Multi-GNSS kinematic PPP experiment using the estimated clocks with quality control shows that the positioning RMS is about 4 cm and generally 2 cm in vertical and horizontal components, respectively, and the corresponding convergence time is about 15 min. [ABSTRACT FROM AUTHOR]

To access the full capabilities of multi-frequency signals from the modernized GPS, GLONASS and newly deployed BDS, Galileo, the undifferenced and uncombined observable model in which the individual signal of each frequency is treated as independent observable has drawn increasing interest in GNSS community. The ionosphere delay is the major issue in the undifferenced and uncombined observable model. Though several ionosphere delay parameterization approaches have been promoted, we argue that the functional model with only deterministic characteristic may not follow the irregular spatial and temporal variations. On the contrary, when the ionosphere delay is estimated as random walk or even white noise with only stochastic characteristic, the ionosphere terms turn out to be non-estimable or not sensitive to their absolute value. In the authors' previous study, we have developed the deterministic plus stochastic ionosphere model, denoted as DESIGN, in which the deterministic part expressed with second-order polynomial is estimated as piece-wise constant over 5 min and the stochastic part is estimated as random walk with constrains derived based on statistics of 4 weeks data in 2010. In this contribution, we further model the deterministic part with Fourier series and update the variogram of the stochastic part accordingly based on two-year data collected by about 150 stations. From the statistic studies, it is concluded that the main frequency components are identical for different coefficients, different stations, as well as different ionosphere activity status, but with varying amplitude. Thus, in the Fourier series expression of the deterministic part, we fix the frequency and estimate the amplitude as daily constant unknowns. Concerning the stochastic component, the variation of variogram is both, geomagnetic latitude and ionosphere activity status dependent. Thus, we use the Gaussian function and Epstein function to model the variation of geomagnetic latitude and ionosphere activity status, respectively. Based on the undifferenced and uncombined observable model with ionosphere constrained with DESIGN, both dual-frequency and single-frequency PPP are carried out to demonstrate its efficiency with three-month data collected in 2010, 2014, and 2017 with different ionosphere activity status. The experimental results suggest that compared with ionosphere-free model and our previous method, the averaged 3D improvement of our new method is 17.8 and 7.6% for dual-frequency PPP, respectively. While for single-frequency PPP, the averaged 3D improvement is 37.0 and 14%, respectively. [ABSTRACT FROM AUTHOR]

During the past years, real-time precise point positioning has been proven to be an efficient tool in the applications of navigation, precise orbit determination of LEO as well as earthquake and tsunami early warning, etc. One of the most crucial issues of these applications is the high-precision real-time GNSS satellite clock. Though the performance and character of the GNSS onboard atomic frequency standard have been widely studied, the white noise model is still the most popular hypothesis that employed in the real-time GNSS satellite clock estimation. However, concerning the real-time applications, significant data discontinuity may arise either due to the fact that only regional stations involved, or the failure in the stations, satellites and network connections. These data discontinuity would result in an arbitrary clock jump between adjacent arcs when the clock offsets are modeled as white noise. In addition, it is also expected that the detection and identification of outliers would be benefited from the constrains of the satellite oscillator noise model. Thus in this contribution, based on the statistic analysis of almost 2-year multi-GNSS precise clock products, we developed the oscillator noise model for the satellites of GPS, GLONASS, BDS and Galileo according to the oscillator type as well as the block type. Then, the efficiency of this oscillator noise model in multi-GNSS satellite clock estimation is demonstrated with 2-months data for both regional and global networks in simultaneous real-time mode. For the regional network, the results suggest that compared with the traditional solution based on white noise model, the improvement is 44.4 and 12.1% on average for STD and RMS, respectively, and the improvement is mainly attributed to the efficiency of the oscillator noise model during the convergence period and the gross error resistance. Concerning the global experiment, since the stations guarantee the continuous tracking of the satellites with redundant observable, the improvement is not as evident as that of regional experiment for GPS, GLONASS and BDS. The STD of Galileo clock improves from 0.28 to 0.19 ns due to that, the satellites E14 and E18 still suffer significant data discontinuity during our experimental period. [ABSTRACT FROM AUTHOR]

Currently, with the rapid development of multi-constellation Global Navigation Satellite Systems (GNSS), the real-time positioning and navigation are undergoing dramatic changes with potential for a better performance. To provide more precise and reliable ultra-rapid orbits is critical for multi-GNSS real-time positioning, especially for the three merging constellations Beidou, Galileo and QZSS which are still under construction. In this contribution, we present a five-system precise orbit determination (POD) strategy to fully exploit the GPS + GLONASS + BDS + Galileo + QZSS observations from CDDIS + IGN + BKG archives for the realization of hourly five-constellation ultra-rapid orbit update. After adopting the optimized 2-day POD solution (updated every hour), the predicted orbit accuracy can be obviously improved for all the five satellite systems in comparison to the conventional 1-day POD solution (updated every 3 h). The orbit accuracy for the BDS IGSO satellites can be improved by about 80, 45 and 50% in the radial, cross and along directions, respectively, while the corresponding accuracy improvement for the BDS MEO satellites reaches about 50, 20 and 50% in the three directions, respectively. Furthermore, the multi-GNSS real-time precise point positioning (PPP) ambiguity resolution has been performed by using the improved precise satellite orbits. Numerous results indicate that combined GPS + BDS + GLONASS + Galileo (GCRE) kinematic PPP ambiguity resolution (AR) solutions can achieve the shortest time to first fix (TTFF) and highest positioning accuracy in all coordinate components. With the addition of the BDS, GLONASS and Galileo observations to the GPS-only processing, the GCRE PPP AR solution achieves the shortest average TTFF of 11 min with 7∘ cutoff elevation, while the TTFF of GPS-only, GR, GE and GC PPP AR solution is 28, 15, 20 and 17 min, respectively. As the cutoff elevation increases, the reliability and accuracy of GPS-only PPP AR solutions decrease dramatically, but there is no evident decrease for the accuracy of GCRE fixed solutions which can still achieve an accuracy of a few centimeters in the east and north components. [ABSTRACT FROM AUTHOR]

All BeiDou navigation satellite system (BDS) satellites are transmitting signals on three frequencies, which brings new opportunity and challenges for high-accuracy precise point positioning (PPP) with ambiguity resolution (AR). This paper proposes an effective uncalibrated phase delay (UPD) estimation and AR strategy which is based on a raw PPP model. First, triple-frequency raw PPP models are developed. The observation model and stochastic model are designed and extended to accommodate the third frequency. Then, the UPD is parameterized in raw frequency form while estimating with the high-precision and low-noise integer linear combination of float ambiguity which are derived by ambiguity decorrelation. Third, with UPD corrected, the LAMBDA method is used for resolving full or partial ambiguities which can be fixed. This method can be easily and flexibly extended for dual-, triple- or even more frequency. To verify the effectiveness and performance of triple-frequency PPP AR, tests with real BDS data from 90 stations lasting for 21 days were performed in static mode. Data were processed with three strategies: BDS triple-frequency ambiguity-float PPP, BDS triple-frequency PPP with dual-frequency (B1/B2) and three-frequency AR, respectively. Numerous experiment results showed that compared with the ambiguity-float solution, the performance in terms of convergence time and positioning biases can be significantly improved by AR. Among three groups of solutions, the triple-frequency PPP AR achieved the best performance. Compared with dual-frequency AR, additional the third frequency could apparently improve the position estimations during the initialization phase and under constraint environments when the dual-frequency PPP AR is limited by few satellite numbers. [ABSTRACT FROM AUTHOR]

The main objective of this research was to examine the feasibility of Multi-GNSS precise point positioning (PPP) in precision agriculture (PA) through a series of experiments with different working modes (i.e. stationary and moving) under different observation conditions (e.g. open sky, with buildings or with canopy). For the stationary test carried out in open space in the UK, the positioning accuracy achieved was 13.9 mm in one dimension by a PPP approach, and the repeatability of positioning results was improved from 19.0 to 6.0 mm by using Multi-GNSS with respect to GPS only. For the moving test carried out in similar location in the UK, almost the same performance was achieved by GPS-only and by Multi-GNSS PPP. However, for a moving experiment carried out in China with obstruction conditions, Multi-GNSS improved the accuracy of baseline length from 126.0 to 35.0 mm and the repeatability from 110.0 mm to 49.0 mm, The results suggested that the addition of the BeiDou, Galileo and GLONASS systems to the standard GPS-only processing improved the positioning repeatability, while a positioning accuracy was achieved at about 20 mm level in the horizontal direction with an improvement against the GPS-only PPP results. In space-constrained and harsh environments (e.g. farms surrounded with dense trees), the availability and reliability of precise positioning decreased dramatically for the GPS-only PPP results, but limited impacts were observed for Multi-GNSS PPP. In addition, compared to real time kinematic (RTK) GNSS, which is currently most commonly used for high precision PA applications, similar accuracy has been achieved by PPP. In contrast to RTK GNSS, PPP can provide high accuracy positioning with higher flexibility and potentially lower capital and running costs. Hence, PPP might be a great opportunity for agriculture to meet the high accuracy requirements of PA in the near future. [ABSTRACT FROM AUTHOR]

Relative positioning using multi-GNSS (global navigation satellite systems) can improve accuracy, reliability, and availability compared to the use of a single constellation system. Intra-system double-difference (DD) ambiguities (ISDDAs) refer to the DD ambiguities between satellites of a single constellation system and can be fixed to an integer to derive the precise fixed solution. Inter-system ambiguities, which denote the DD ambiguities between different constellation systems, can also be fixed to integers on overlapping frequencies, once the inter-system bias (ISB) is removed. Compared with fixing ISDDAs, fixing both integer intra- and inter-system DD ambiguities (IIDDAs) means an increase of positioning precision through an integration of multiple GNSS constellations. Previously, researchers have studied IIDDA fixing with systems of the same frequencies, but not with systems of different frequencies. Integer IIDDAs can be determined from single-difference (SD) ambiguities, even if the frequencies of multi-GNSS signals used in the positioning are different. In this study, we investigated IIDDA fixing for multi-GNSS signals of narrowly spaced frequencies. First, the inter-system DD models of multi-GNSS signals of different frequencies are introduced, and the strategy for compensating for ISB is presented. The ISB is decomposed into three parts: 1) a float approximate ISB number that can be considered equal to the ISB of code pseudorange observations and thus can be estimated through single point positioning (SPP); 2) a number that is a multiple of the GNSS signal wavelength; and 3) a fractional ISB part, with a magnitude smaller than a single wavelength. Then, the relationship between intra- and inter-system DD ambiguity RATIO values and ISB was investigated by integrating GPS L1 and GLONASS L1 signals. In our numerical analyses with short baselines, the ISB parameter and IIDDA were successfully fixed, even if the number of observed satellites in each system was small. [ABSTRACT FROM AUTHOR]

Global navigation satellite systems (GNSS) are acting as an indispensable tool for geodetic research and global monitoring of the Earth, and they have been rapidly developed over the past few years with abundant GNSS networks, modern constellations, and significant improvement in mathematic models of data processing. However, due to the increasing number of satellites and stations, the computational efficiency becomes a key issue and it could hamper the further development of GNSS applications. In this contribution, this problem is overcome from the aspects of both dense linear algebra algorithms and GNSS processing strategy. First, in order to fully explore the power of modern microprocessors, the square root information filter solution based on the blocked QR factorization employing as many matrix-matrix operations as possible is introduced. In addition, the algorithm complexity of GNSS data processing is further decreased by centralizing the carrier-phase observations and ambiguity parameters, as well as performing the real-time ambiguity resolution and elimination. Based on the QR factorization of the simulated matrix, we can conclude that compared to unblocked QR factorization, the blocked QR factorization can greatly improve processing efficiency with a magnitude of nearly two orders on a personal computer with four 3.30 GHz cores. Then, with 82 globally distributed stations, the processing efficiency is further validated in multi-GNSS (GPS/BDS/Galileo) satellite clock estimation. The results suggest that it will take about 31.38 s per epoch for the unblocked method. While, without any loss of accuracy, it only takes 0.50 and 0.31 s for our new algorithm per epoch for float and fixed clock solutions, respectively. [ABSTRACT FROM AUTHOR]

In recent years, the precise orbit determination (POD) of the regional Chinese BeiDou Navigation Satellite System (BDS) has been a hot spot because of its special constellation consisting of five geostationary earth orbit (GEO) satellites and five inclined geosynchronous satellite orbit (IGSO) satellites besides four medium earth orbit (MEO) satellites since the end of 2012. GEO and IGSO satellites play an important role in regional BDS applications. However, this brings a great challenge to the POD, especially for the GEO satellites due to their geostationary orbiting. Though a number of studies have been carried out to improve the POD performance of GEO satellites, the result is still much worse than that of IGSO and MEO, particularly in the along-track direction. The major reason is that the geostationary characteristic of a GEO satellite results in a bad geometry with respect to the ground tracking network. In order to improve the tracking geometry of the GEO satellites, a possible strategy is to mount global navigation satellite system (GNSS) receivers on MEO satellites to collect the signals from GEO/IGSO GNSS satellites so as that these observations can be used to improve GEO/IGSO POD. We extended our POD software package to simulate all the related observations and to assimilate the MEO-onboard GNSS observations in orbit determination. Based on GPS and BDS constellations, simulated studies are undertaken for various tracking scenarios. The impact of the onboard GNSS observations is investigated carefully and presented in detail. The results show that MEO-onboard observations can significantly improve the orbit precision of GEO satellites from metres to decimetres, especially in the along-track direction. The POD results of IGSO satellites also benefit from the MEO-onboard data and the precision can be improved by more than 50% in 3D direction. [ABSTRACT FROM AUTHOR]

Single-frequency precise point positioning (SF-PPP) is a potential precise positioning technique due to the advantages of the high accuracy in positioning after convergence and the low cost in operation. However, there are still challenges limiting its applications at present, such as the long convergence time, the low reliability, and the poor satellite availability and continuity in kinematic applications. In recent years, the achievements in the dual-frequency PPP have confirmed that its performance can be significantly enhanced by employing the slant ionospheric delay and receiver differential code bias (DCB) constraint model, and the multi-constellation Global Navigation Satellite Systems (GNSS) data. Accordingly, we introduce the slant ionospheric delay and receiver DCB constraint model, and the multi-GNSS data in SF-PPP modular together. In order to further overcome the drawbacks of SF-PPP in terms of reliability, continuity, and accuracy in the signal easily blocking environments, the inertial measurements are also adopted in this paper. Finally, we form a new approach to tightly integrate the multi-GNSS single-frequency observations and inertial measurements together to ameliorate the performance of the ionospheric delay and receiver DCB-constrained SF-PPP. In such model, the inter-system bias between each two GNSS systems, the inter-frequency bias between each two GLONASS frequencies, the hardware errors of the inertial sensors, the slant ionospheric delays of each user-satellite pair, and the receiver DCB are estimated together with other parameters in a unique Kalman filter. To demonstrate its performance, the multi-GNSS and low-cost inertial data from a land-borne experiment are analyzed. The results indicate that visible positioning improvements in terms of accuracy, continuity, and reliability can be achieved in both open-sky and complex conditions while using the proposed model in this study compared to the conventional GPS SF-PPP. [ABSTRACT FROM AUTHOR]

Precise positioning with the current Chinese BeiDou Navigation Satellite System is proven to be of comparable accuracy to the Global Positioning System, which is at centimeter level for the horizontal components and sub-decimeter level for the vertical component. But the BeiDou precise point positioning (PPP) shows its limitation in requiring a relatively long convergence time. In this study, we develop a numerical weather model (NWM) augmented PPP processing algorithm to improve BeiDou precise positioning. Tropospheric delay parameters, i.e., zenith delays, mapping functions, and horizontal delay gradients, derived from short-range forecasts from the Global Forecast System of the National Centers for Environmental Prediction (NCEP) are applied into BeiDou real-time PPP. Observational data from stations that are capable of tracking the BeiDou constellation from the International GNSS Service (IGS) Multi-GNSS Experiments network are processed, with the introduced NWM-augmented PPP and the standard PPP processing. The accuracy of tropospheric delays derived from NCEP is assessed against with the IGS final tropospheric delay products. The positioning results show that an improvement in convergence time up to 60.0 and 66.7% for the east and vertical components, respectively, can be achieved with the NWM-augmented PPP solution compared to the standard PPP solutions, while only slight improvement in the solution convergence can be found for the north component. A positioning accuracy of 5.7 and 5.9 cm for the east component is achieved with the standard PPP that estimates gradients and the one that estimates no gradients, respectively, in comparison to 3.5 cm of the NWM-augmented PPP, showing an improvement of 38.6 and 40.1%. Compared to the accuracy of 3.7 and 4.1 cm for the north component derived from the two standard PPP solutions, the one of the NWM-augmented PPP solution is improved to 2.0 cm, by about 45.9 and 51.2%. The positioning accuracy for the up component improves from 11.4 and 13.2 cm with the two standard PPP solutions to 8.0 cm with the NWM-augmented PPP solution, an improvement of 29.8 and 39.4%, respectively. [ABSTRACT FROM AUTHOR]

Although double-differenced (DD) observations between satellites from different systems can be used in multi-GNSS relative positioning, the inter-system DD ambiguities cannot be fixed to integer because of the existence of the inter-system bias (ISB). Obviously, they can also be fixed as integer along with intra-system DD ambiguities if the associated ISBs are well known. It is critical to fix such inter-system DD ambiguities especially when only a few satellites of each system are observed. In most of the existing approaches, the ISB is derived from the fractional part of the inter-system ambiguities after the intra-system DD ambiguities are successfully fixed. In this case, it usually needs observations over long times depending on the number of observed satellites from each system. We present a new method by means of particle filter to estimate ISBs in real time without any a priori information based on the fact that the accuracy of a given ISB value can be qualified by the related fixing RATIO. In this particle filter-based method, the ISB parameter is represented by a set of samples, i.e., particles, and the weight of each sample is determined by the designed likelihood function related to the corresponding RATIO, so that the true bias value can be estimated successfully. Experimental validations with the IGS multi-GNSS experiment data show that this method can be carried out epoch by epoch to provide precise ISB in real time. Although there are only one, two, or at most three Galileo satellites observed, the successfully fixing rate increases from 75.5% for GPS only to 81.2%. In the experiment with five GPS satellites and one Galileo satellites, the first successfully fixing time is reduced to half of that without fixing the inter-system DD ambiguities. [ABSTRACT FROM AUTHOR]

Precise point positioning (PPP) using the Global Positioning System (GPS) is widely recognized as an efficient approach for providing precise positioning services. However, its accuracy and reliability could be significantly degraded by unexpected observation discontinuities and unfavorable tracking geometry which are unavoidable, especially in severe environments such as city canyons. Therefore, in the last decades inertial navigation system (INS) has been integrated to overcome such drawbacks. Recently, multi-Global Navigation Satellite Systems (GNSS) were applied to enhance the PPP performance by appropriate usage of the increased number of satellites. We present a new approach to tightly integrate the multi-GNSS PPP and INS together in the observation level. The inter-system bias and inter-frequency bias of multi-GNSS and the hardware errors of INS sensors are estimated to improve the position accuracy and to shorten the convergence time of PPP. In order to demonstrate the impact of multi-GNSS observations and INS data on the derived position, velocity, attitude, and the convergence time of PPP, the new approach is validated through an experimental test with a set of land vehicle data. The results show that the position accuracy can be improved by multi-GNSS and INS significantly, but very little improvement in velocity and attitude is achieved. The position root-mean-square improves from 23.3, 19.8, and 14.9 cm of the GPS PPP/INS tightly coupled integration (TCI) solution to 7.9, 3.3, and 5.1 cm of multi-GNSS PPP/INS TCI in north, east, and up components, respectively. Furthermore, GNSS outages are simulated and their effect on the performance of multi-GNSS PPP/INS TCI is investigated to demonstrate the contribution of the multi-GNSS PPP/INS TCI during GNSS outages. In addition, the convergence test also shows that both multi-GNSS and INS can improve the PPP convergence performance noticeably. [ABSTRACT FROM AUTHOR]

The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5-2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing. [ABSTRACT FROM AUTHOR]

The frequency division multiple access adopted in present GLONASS introduces inter-frequency bias (IFB) at the receiver-end both in code and phase observables, which makes GLONASS ambiguity resolution rather difficult or even not available, especially for long baselines up to several thousand kilometers. This is one of the major reasons that GLONASS could hardly reach the orbit precision of GPS, both in terms of consistency among individual International GNSS Service (IGS) analysis centers and discontinuity at the overlapping day boundaries. Based on the fact that the GLONASS phase IFB is similar on L1 and L2 bands in unit of length and is a linear function of the frequency number, several approaches have been developed to estimate and calibrate the IFB for integer ambiguity resolution. However, they are only for short and medium baselines. In this study, a new ambiguity resolution approach is developed for GLONASS global networks. In the approach, the phase ambiguities in the ionosphere-free linear combination are directly transformed with a wavelength of about 5.3 cm, according to the special frequency relationship of GLONASS L1 and L2 signals. After such transformation, the phase IFB rate can be estimated and corrected precisely and then the corresponding double-differenced ambiguities can be directly fixed to integers even for baselines up to several thousand kilometers. To evaluate this approach, experimental validations using one-month data of a global network with 140 IGS stations was carried out for GLONASS precise orbit determination. The results show that the GLONASS double-difference ambiguity resolution for long baselines could be achieved with an average fixing-rate of 91.4 %. Applying the fixed ambiguities as constraints, the GLONASS orbit overlapping RMS at the day boundaries could be reduced by 37.2 % in ideal cases and with an averaged reduction of about 21.4 %, which is comparable with that by the GPS ambiguity resolution. The orbit improvement is also confirmed by the better agreement with the independent satellite laser ranging observations. [ABSTRACT FROM AUTHOR]

GLONASS could hardly reach the positioning performance of GPS, especially for fast and real-time precise positioning. One of the reasons is the phase inter-frequency bias (IFB) at the receiver end prevents its integer ambiguity resolution. A number of studies were carried out to achieve the integer ambiguity resolution for GLONASS. Based on some of the revealed IFB characteristics, for instance IFB is a linear function of the received carrier frequency and L1 and L2 have the same IFB in unit of length, most of recent methods recommend estimating the IFB rate together with ambiguities. However, since the two sets of parameters are highly correlated, as demonstrated in previous studies, observations over several hours up to 1 day are needed even with simultaneous GPS observations to obtain a reasonable solution. Obviously, these approaches cannot be applied for real-time positioning. Actually, it can be demonstrated that GLONASS ambiguity resolution should also be available even for a single epoch if the IFB rate is precisely known. In addition, the closer the IFB rate value is to its true value, the larger the fixing RATIO will be. Based on this fact, in this paper, a new approach is developed to estimate the IFB rate by means of particle filtering with the likelihood function derived from RATIO. This approach is evaluated with several sets of experimental data. For both static and kinematic cases, the results show that IFB rates could be estimated precisely just with GLONASS data of a few epochs depending on the baseline length. The time cost with a normal PC can be controlled around 1 s and can be further reduced. With the estimated IFB rate, integer ambiguity resolution is available immediately and as a consequence, the positioning accuracy is improved significantly to the level of GPS fixed solution. Thus the new approach enables real-time precise applications of GLONASS. [ABSTRACT FROM AUTHOR]

Precise knowledge and consistent modeling of the yaw-attitude of GNSS satellites are essential for high-precision data processing and applications. As the exact attitude control mechanism for the satellites of the BeiDou Satellite Navigation System (BDS) is not yet released, the reverse kinematic precise point positioning (PPP) method was applied in our study. However, we confirm that the recent precise orbit determination (POD) processing for GPS satellites could not provide suitable products for estimating BDS attitude using the reverse PPP because of the special attitude control switching between the nominal and the orbit-normal mode. In our study, we propose a modified processing schema for studying the attitude behavior of the BDS satellites. In this approach, the observations of the satellites during and after attitude switch are excluded in the POD processing, so that the estimates, which are needed in the reverse PPP, are not contaminated by the inaccurate initial attitude mode. The modified process is validated by experimental data sets and the attitude yaw-angles of the BDS IGSO and MEO satellites are estimated with an accuracy of better than $$9^{\circ }$$ . Furthermore, the results confirm that the switch is executed when the Sun elevation is about $$4^{\circ }$$ and the actual orientation is very close to its target one. Based on the estimated yaw-angles, a preliminary attitude switch model was established and reintroduced into the POD, yielding to a substantial improvement in the orbit overlap RMS. [ABSTRACT FROM AUTHOR]

The rapid development of the Chinese BeiDou Navigation Satellite System (BDS) brings a promising prospect for the real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV), which is of great benefit for supporting the time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a real-time ZTD/PWV processing method based on Global Positioning System (GPS) and BDS observations. The performance of ZTD and PWV derived from BDS observations using real-time precise point positioning (PPP) technique is carefully investigated. The contribution of combining BDS and GPS for ZTD/PWV retrieving is evaluated as well. GPS and BDS observations of a half-year period for 40 globally distributed stations from the International GNSS Service Multi-GNSS Experiment and BeiDou Experiment Tracking Network are processed. The results show that the real-time BDS-only ZTD series agree well with the GPS-only ZTD series in general: the RMS values are about 11-16 mm (about 2-3 mm in PWV). Furthermore, the real-time ZTD derived from GPS-only, BDS-only, and GPS/BDS combined solutions are compared with those derived from the Very Long Baseline Interferometry. The comparisons show that the BDS can contribute to real-time meteorological applications, slightly less accurately than GPS. More accurate and reliable water vapor estimates, about 1.3-1.8 mm in PWV, can be obtained if the BDS observations are combined with the GPS observations in the real-time PPP data processing. The PWV comparisons with radiosondes further confirm the performance of BDS-derived real-time PWV and the benefit of adding BDS to standard GPS processing. [ABSTRACT FROM AUTHOR]

In this contribution, we present a GPS+GLONASS+BeiDou+Galileo four-system model to fully exploit the observations of all these four navigation satellite systems for real-time precise orbit determination, clock estimation and positioning. A rigorous multi-GNSS analysis is performed to achieve the best possible consistency by processing the observations from different GNSS together in one common parameter estimation procedure. Meanwhile, an efficient multi-GNSS real-time precise positioning service system is designed and demonstrated by using the multi-GNSS Experiment, BeiDou Experimental Tracking Network, and International GNSS Service networks including stations all over the world. The statistical analysis of the 6-h predicted orbits show that the radial and cross root mean square (RMS) values are smaller than 10 cm for BeiDou and Galileo, and smaller than 5 cm for both GLONASS and GPS satellites, respectively. The RMS values of the clock differences between real-time and batch-processed solutions for GPS satellites are about 0.10 ns, while the RMS values for BeiDou, Galileo and GLONASS are 0.13, 0.13 and 0.14 ns, respectively. The addition of the BeiDou, Galileo and GLONASS systems to the standard GPS-only processing, reduces the convergence time almost by 70 %, while the positioning accuracy is improved by about 25 %. Some outliers in the GPS-only solutions vanish when multi-GNSS observations are processed simultaneous. The availability and reliability of GPS precise positioning decrease dramatically as the elevation cutoff increases. However, the accuracy of multi-GNSS precise point positioning (PPP) is hardly decreased and few centimeter are still achievable in the horizontal components even with 40 $$^{\circ }$$ elevation cutoff. At 30 $$^{\circ }$$ and 40 $$^{\circ }$$ elevation cutoffs, the availability rates of GPS-only solution drop significantly to only around 70 and 40 %, respectively. However, multi-GNSS PPP can provide precise position estimates continuously (availability rate is more than 99.5 %) even up to 40 $$^{\circ }$$ elevation cutoff (e.g., in urban canyons). [ABSTRACT FROM AUTHOR]

Although the computational burden of global navigation satellite systems (GNSS) data processing is nowadays already a big challenge, especially for huge networks, integrated processing of denser networks with data of multi-GNSS and multi-frequency is desired in the expectation of more accurate and reliable products. Based on the concept of carrier range, in this study, the precise point positioning with integer ambiguity resolution is engaged to obtain the integer ambiguities for converting carrier phases to carrier ranges. With such carrier ranges and pseudo-ranges, rigorous integrated processing is realized computational efficiently for the orbit and clock estimation using massive networks. The strategy is validated in terms of computational efficiency and product quality using data of the IGS network with about 460 stations. The experimental validation shows that the computation time of the new strategy increases gradually with the number of stations. It takes about 14 min for precise orbit and clock determination with 460 stations, while the current strategy needs about 82 min. The overlapping orbit RMS is reduced from 27.6 mm with 100 stations to 24.8 mm using the proposed strategy, and the RMS could be further reduced to 23.2 mm by including all 460 stations. Therefore, the new strategy could be applied to massive networks of multi-GNSS and multi-frequency receivers and possibly to achieve GNSS data products of higher quality. [ABSTRACT FROM AUTHOR]

The complementary advantages of GPS and seismic measurements are well recognized in seismotectonic monitoring studies. Therefore, integrated processing of the two data streams has been proposed recently in an attempt to obtain accurate and reliable information of surface displacements associated with earthquakes. A hitherto still critical issue in the integrated processing is real-time detection and precise estimation of the transient baseline error in the seismic records. Here, we report on a new approach by introducing the seismic acceleration corrected by baseline errors into the state equation system. The correction is performed and regularly updated in short epochs (with increments which may be as short as seconds), so that station position, velocity, and acceleration can be constrained very tightly and baseline error can be estimated as a random-walk process. With the adapted state equation system, our study highlights the use of a new approach developed for integrated processing of GPS and seismic data by means of sequential least-squares adjustment. The efficiency of our approach is demonstrated and validated using simulated, experimental, and real datasets. The latter were collected at collocated GPS and seismic stations around the 4 April 2010, E1 Mayor-Cucapah earthquake (Mw, 7.2). The results have shown that baseline errors of the strong-motion sensors are corrected precisely and high-precision seismic displacements are real-timely obtained by the new approach. [ABSTRACT FROM AUTHOR]

An increasing number of GNSS reference stations are installed around the world to provide real-time precise positioning services. In most of the current services, a full network solution is required for the precise determination of biases. Such a network solution is time consuming and difficult to achieve for very large regions such as Europe or China. Therefore, we developed a multi-layer processing scheme for precise point positioning (PPP) regional augmentation to avoid processing large networks. Furthermore, we use L1 and L2 raw observations and estimate atmospheric delays, which were properly constrained to the atmospheric corrections derived from the reference stations. Therefore, inaccurate representation of atmospheric delays due to temporal and/or spatial atmospheric fluctuations in the processing can be compensated. The proposed scheme of PPP regional augmentation was implemented into the operational real-time PPP service system at GFZ for validation. The real-time orbit and clock corrections, the uncalibrated phase delays, and regional augmentation corrections are generated by this system. The augmentation corrections from the regional network are investigated and the positioning performance in terms of positioning accuracy and time for fixed solution is demonstrated in real-time. Our results indicate that a reliable fixing is possible after 5 s on average. The positioning accuracy is about 12, 10, and 25 mm in east, north, and vertical direction, respectively. [ABSTRACT FROM AUTHOR]

Standard network RTK has been widely used since it was proposed in the mid-1990s. Rovers can obtain high-precision estimates of position by resolving double-differenced (DD) ambiguities. The focus of this study is a new undifferenced network RTK method, abbreviated as URTK hereafter, based on undifferenced (UD) observation corrections whose single-differenced (SD) ambiguities between satellites can be resolved in several seconds. The tools for studying the real-time realization of the new method are our developments of logical schemes that have the capability for the real-time modeling of a reference network and the instantaneous resolution of SD ionosphere-free (IF) ambiguities at a single station. This research demonstrates the validity of modeling regional UD-unmodeled errors on the ground and examines the maximum differences when compared to modeling the errors using ionospheric pierce points (IPP). With data collected at 48 stations from a CORS network in Shanxi Province (SXCORS) in China through May 21, 2010, the efficiency of the presented real-time strategies is validated using IGS final products in a postprocessing mode. The results verify that more than 83 % of SD wide-lane (WL) ambiguity can be fixed with 5 s of observation data, and the average resolution time of all the WL tests is 4.96 s. More than 80 % of SD L1 ambiguity can be fixed within 5 s, and the average resolution time is only 6.66 s. Rovers could gain rapidly centimeter-level absolute positioning service, comparable to standard network RTK. In addition, the URTK method transforms the fixed DD-ambiguities of the reference network into UD-ambiguities, and it does not need to set the base station and base satellite. Since the UD-corrections are modeled for each common visible satellite, it breaks down the connections between stations and satellites of the DD-corrections in the current network RTK. The UD-corrections can be broadcast by the base station and automatically selected and optimized by a rover during the real-time kinematic processing, thus avoiding ambiguity in reinitialization due to the change of reference, so it should be very flexible and useful for a wide range of applications. [ABSTRACT FROM AUTHOR]

In order to improve the performance of precise point positioning (PPP), this paper presents a new data processing scheme to shorten the convergence time and the observation time required for a reliable ambiguity-fixing. In the new scheme, L1 and L2 raw observations are used and the slant ionospheric delays are treated as unknown parameters. The empirical spatial and temporal constraints and the ionospheric delays derived from a real-time available ionospheric model are all considered as pseudo-observations into the estimation for strengthening the solution. Furthermore, we develop a real-time computational procedure for generating uncalibrated phase delays (UPDs) on L1 and L2 frequencies. The PPP solution is first carried out on all reference stations based on the proposed scheme, undifferenced float ambiguities on L1 and L2 frequencies can be directly obtained from the new scheme. The L1 and L2 UPDs are then generated and broadcasted to users in real-time. This data product and also the performance of the new PPP scheme are evaluated. Our results indicate that the new processing scheme considering ionospheric characteristics can reduce the convergence time by about 30 % for float kinematic solutions. The observation time for a reliable ambiguity-fixing is shortened by 25 % compared to that of the traditional ambiguity-fixed kinematic solution. When the new method is used for static reference stations, the observation time for ambiguity-fixing is about 10 min in static mode and only 5 min if the coordinates are fixed to well-known values. [ABSTRACT FROM AUTHOR]

Ambiguity resolution dedicated to a single global positioning system (GPS) station can improve the accuracy of precise point positioning. In this process, the estimation accuracy of the narrow-lane fractional-cycle biases (FCBs), which destroy the integer nature of undifferenced ambiguities, is crucial to the ambiguity-fixed positioning accuracy. In this study, we hence propose the improved narrow-lane FCBs derived from an ambiguity-fixed GPS network solution, rather than the original (i.e. previously proposed) FCBs derived from an ambiguity-float network solution. The improved FCBs outperform the original FCBs by ensuring that the resulting ambiguity-fixed daily positions coincide in nature with the state-of-the-art positions generated by the International GNSS Service (IGS). To verify this improvement, 1 year of GPS measurements from about 350 globally distributed stations were processed. We find that the original FCBs differ more from the improved FCBs when fewer stations are involved in the FCB estimation, especially when the number of stations is less than 20. Moreover, when comparing the ambiguity-fixed daily positions with the IGS weekly positions for 248 stations through a Helmert transformation, for the East component, we find that on 359 days of the year the daily RMS of the transformed residuals based on the improved FCBs is smaller by up to 0.8 mm than those based on the original FCBs, and the mean RMS over the year falls evidently from 2.6 to 2.2 mm. Meanwhile, when using the improved rather than the original FCBs, the RMS of the transformed residuals for the East component of 239 stations (i.e. 96.4% of all 248 stations) is clearly reduced by up to 1.6 mm, especially for stations located within a sparse GPS network. Therefore, we suggest that narrow-lane FCBs should be determined with ambiguity-fixed, rather than ambiguity-float, GPS network solutions. [ABSTRACT FROM AUTHOR]

Realtime satellite clock corrections are usually estimated using undifferenced phase and range observations from a global network. Because a large number of ambiguity parameters must be estimated, the computation is time-consuming. Consequently, only a sparse global network of limited number of stations is processed by most IGS Realtime Analysis Centers with an update rate of 5 s. In addition, it is very desirable to build the capability to simultaneously estimate clock corrections for multi-GNSS constellations. Although the estimation can be sped up by epoch-differenced observations that eliminate ambiguities, the derived clocks can contain a satellite-specific bias that diminishes the contribution of range observations. We introduce a computationally efficient approach for realtime clock estimation. Both the epoch-differenced phase and undifferenced range observations are used together to estimate the epoch-differenced satellite clocks and the initial clock bias for each satellite and receiver. The biased clock corrections accumulated from the estimated epoch-differenced clocks are then aligned with the estimated clock biases and provided as the final clock corrections to users. The algorithm is incorporated into the EPOS-RT software developed at GFZ (GeoForschungsZentrum) and experimentally validated with the IGS global network. The comparison with the GFZ rapid products shows that the accuracy of the clock estimation with the new approach is comparable with that of the undifferenced approach, whereas the computation time is reduced to one-tenth. As a result, estimation of high-rate satellite clocks from a large reference network and tracking satellites of multi-GNSS constellations becomes achievable. [ABSTRACT FROM AUTHOR]

Integer ambiguity resolution at a single receiver can be achieved if the fractional-cycle biases are separated from the ambiguity estimates in precise point positioning (PPP). Despite the improved positioning accuracy by such integer resolution, the convergence to an ambiguity-fixed solution normally requires a few tens of minutes. Even worse, these convergences can repeatedly occur on the occasion of loss of tracking locks for many satellites if an open sky-view is not constantly available, consequently totally destroying the practicability of real-time PPP. In this study, in case of such re-convergences, we develop a method in which ionospheric delays are precisely predicted to significantly accelerate the integer ambiguity resolution. The effectiveness of this method consists in two aspects: first, wide-lane ambiguities can be rapidly resolved using the ionosphere-corrected wide-lane measurements, instead of the noisy Melbourne-Wübbena combination measurements; second, narrow-lane ambiguity resolution can be accelerated under the tight constraints derived from the ionosphere-corrected unambiguous wide-lane measurements. In the test at 90 static stations suffering from simulated total loss of tracking locks, 93.3 and 95.0% of re-convergences to wide-lane and narrow-lane ambiguity resolutions can be achieved within five epochs of 1-Hz measurements, respectively, even though the time latency for the predicted ionospheric delays is up to 180 s. In the test at a mobile van moving in a GPS-adverse environment where satellite number significantly decreases and cycle slips frequently occur, only when the predicted ionospheric delays are applied can the rate of ambiguity-fixed epochs be dramatically improved from 7.7 to 93.6% of all epochs. Therefore, this method can potentially relieve the unrealistic requirement of a continuous open sky-view by most PPP applications and improve the practicability of real-time PPP. [ABSTRACT FROM AUTHOR]

Global navigation satellite system (GNSS) precise point positioning (PPP) requires continuous carrier-phase observations to achieve a solution of high precision. Precisely correcting cycle slips caused by signal interruptions is crucial for recovering the data continuity. Most of the existing approaches usually employ only data of one epoch after the interruption for real-time cycle slip processing. In this study, we propose to introduce and estimate cycle slip parameters together with standard PPP parameters, such as position, ionospheric delay, and ambiguities in the case that possible cycle slips are detected, using a Kalman-filter-based procedure with the undifferenced and uncombined PPP model. The integer search strategy is used to fix cycle slips. To reduce the probability of wrong integer fixing, a strict integer validation threshold is suggested. As a result, it is not easy to fix all cycle slips with only one epoch of data. Our approach can be easily extended to use multi-epoch observations to enhance the cycle slip estimation. Once the cycle slips are correctly determined, continuous PPP can be achieved instantaneously. This new approach is tested and validated with three groups of experiments using GPS and GLONASS stations operated by the International GNSS Service from DOY 1–10, 2017, and a real vehicle kinematic data. Numerous experimental results showed that the proposed method can correctly fix the cycle slips for more than 99.5% of epochs suffering from re-convergence. On average, this method takes observation information from about 1.5–2.5 epochs to fix cycle slips and realize rapid re-convergence. Consequently, positioning performance is significantly improved. [ABSTRACT FROM AUTHOR]