Network-based triple-frequency carrier phase ambiguity resolution between reference stations using BDS data for long baselines

Network-based ambiguity resolution (AR) between reference stations is the prerequisite to realize a precise real-time kinematic positioning service. With the help of BDS triple-frequency signals, we can efficiently deal with the ionospheric delay and tropospheric delay, and achieve rapid and reliable AR. To overcome the inaccurate ionospheric delay estimated by the geometry-free three carrier ambiguity resolution (GF TCAR) technique, which leads to failure in the original ambiguity resolution, we propose an ionospheric-free (IF) TCAR method to resolve the ambiguity between the reference stations over long baselines. Taking full advantage of the known positions of the reference stations, the easily resolved extra-wide-lane (EWL) ambiguity, and the IF phase combinations, we can reliably fix the wide-lane (WL) ambiguity. A Kalman filter is applied to estimate precise IF ambiguities and the original ambiguity is resolved with the fixed WL ambiguity. A numerical analysis with triple-frequency BDS data from three long baselines of a CORS network is provided to compare the AR performance of GF TCAR with that of IF TCAR. The results show that both methods can reliably resolve the WL ambiguity with a remarkable correctly-fixed rate of higher than 99%, and the reliably-fixed rates of the IF TCAR slightly increase from 92.19, 94.67 and 94.61–98.26, 99.54 and 97.51% for the three baselines. Herein “correctly-fixed” and “reliably-fixed” mean the difference between the float ambiguity and the true one are less than ± 0.5 and ± 0.25 cycles, respectively. On the other hand, the AR performance of the original signals with the IF TCAR method is much better than that with the GF TCAR method attaining a 100% correctly-fixed rate, while the GF TCAR method can hardly fix the original ambiguity with the largest bias being as much as 4 cycles because of the amplified systematic bias.     

A modified mix-differenced approach for estimating multi-GNSS real-time satellite clock offsets

A modified mixed-differenced approach for estimating multi-GNSS real-time clock offsets is presented. This approach, as compared to the earlier presented mixed-differenced approach which uses epoch-differenced and undifferenced observations, further adds a satellite-differenced process. The proposed approach, based on real-time orbit products and a mix of epoch-differenced and satellite-differenced observations to estimate only satellite clock offsets and tropospheric zenith wet delays, has fewer estimated parameters than other approaches, and thus its implementing procedure is efficient and can be performed and extended easily. To obtain high accuracy, the approach involves three steps. First, the high-accuracy tropospheric zenith wet delay of each station is estimated using mixed-differenced carrier phase observations. Second, satellite clock offset changes between adjacent epochs are estimated using also mixed-differenced carrier phase observations. Third, the satellite clock offsets at the initial epoch are estimated using satellite-differenced pseudorange observations. Finally, the initial epoch clock results and clock offset changes are concatenated to obtain the clock results of the current epoch. To validate the real-time satellite clock results, multi-GNSS post-processing clock products from IGS ACs were selected for comparison. From the comparison, the standard deviations of the GPS, GLONASS, BeiDou and Galileo systems clock results are approximately 0.1–0.4 ns, except for the BeiDou GEO satellites. The root mean squares are about 0.4–2.3 ns, which are similar to those of other international real-time products. When the clock estimates were assessed based on a pseudo-kinematic PPP procedure, the positioning accuracies in the East, North and Up components reach 5.6, 5.5 and 7.6 cm, respectively, which meet the centimeter level and are comparable to the application of other products.     

On the accuracy of the GPS L2 observable for ionospheric monitoring

The introduction of the unencrypted global positioning system (GPS) L2 civil (L2C) signal has the potential to improve measurements made with the L2 frequency, an important observable in GPS-based ionospheric research and monitoring. Recent work has shown significant differences between the legacy L2P(Y) and L2C-derived total electron content rate of change index (ROTI). This difference is observed between L2P(Y) and L2C-derived ROTI with certain receiver models and between zero-baseline receiver pairs. We discuss the likely cause for these differences: L1-aided tracking used to track both the L2P(Y) and L2C signals. We also present L2C data that are confirmed to be from tracking independent of L1. Using the ionospheric-free linear combination, we show that the independently tracked carrier phase dynamics are significantly more accurate than the L1-aided observables. This result is confirmed by comparing the behavior of the L2C and L2P(Y) carrier phase observables upon a sudden antenna rotation.     

A new evil waveforms evaluating method for new BDS navigation signals

With the advent of new global navigation satellite systems (GNSSs) and new signals, GNSS users will rely more on them to obtain higher-accuracy positioning. Evil waveform monitoring and assessment are of great importance for GNSS to achieve its positioning, velocity, and timing service with high accuracy. However, the advent of new navigation signals introduces the necessity to extend the traditional analyzing techniques already accepted for binary phase-shift keying modulation to new techniques. First, the well-known second-order step thread model adopted by the International Civil Aviation Organization is introduced. Then the extended new general thread models are developed for the new binary offset carrier modulated signals. However, no research has been done on navigation signal waveform symmetry yet. Simulation results showed that, waveform asymmetry may also cause tracking errors, range biases, and position errors in GNSS receivers. It is thus imperative that the asymmetry be quantified to enable the design of appropriate error budgets and mitigation strategies for various application fields. A novel evil waveform analysis method, called waveform rising and falling edge symmetry (WRaFES) method, is proposed. Based on this WRaFES method, the correlation metrics are provided to detect asymmetric correlation peaks distorted by received signal asymmetry. Then the statistical properties of the proposed methods are analyzed, and a proper deformation detection threshold is calculated. Finally, both simulation results and experimentally measured results of Beidou navigation satellite system (BDS) M1-S B1Cd signal are given, which show the effectiveness and robustness of the proposed thread models.     

Position-domain integrity risk-based ambiguity validation for the integer bootstrap estimator

Integrity monitoring for ambiguity resolution is of significance for utilizing the high-precision carrier phase differential positioning for safety–critical navigational applications. The integer bootstrap estimator can provide an analytical probability density function, which enables the precise evaluation of the integrity risk for ambiguity validation. In order to monitor the effect of unknown ambiguity bias on the integer bootstrap estimator, the position-domain integrity risk of the integer bootstrapped baseline is evaluated under the complete failure modes by using the worst-case protection principle. Furthermore, a partial ambiguity resolution method is developed in order to satisfy the predefined integrity risk requirement. Static and kinematic experiments are carried out to test the proposed method by comparing with the traditional ratio test method and the protection level-based method. The static experimental result has shown that the proposed method can achieve a significant global availability improvement by 51% at most. The kinematic result reveals that the proposed method obtains the best balance between the positioning accuracy and the continuity performance.     

Beidou satellite maneuver thrust force estimation for precise orbit determination

Beidou satellites, especially geostationary earth orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites, need to be frequently maneuvered to keep them in position due to various perturbations. The satellite ephemerides are not available during such maneuver periods. Precise estimation of thrust forces acting on satellites would provide continuous ephemerides during maneuver periods and could significantly improve orbit accuracy immediately after the maneuver. This would increase satellite usability for both real-time and post-processing applications. Using 1 year of observations from the Multi-GNSS Experiment network (MGEX), we estimate the precise maneuver periods for all Beidou satellites and the thrust forces. On average, GEO and IGSO satellites in the Beidou constellation are maneuvered 12 and 2 times, respectively, each year. For GEO satellites, the maneuvers are mainly in-plane, while out-of-plane maneuvers are observed for IGSO satellites and a small number of GEO satellites. In most cases, the Beidou satellite maneuver periods last 15–25 min, but can be as much as 2 h for the few out-of-plane maneuvers of GEO satellites. The thrust forces acting on Beidou satellites are normally in the order of 0.1–0.7 mm/s². This can cause changes in velocity of GEO/IGSO satellites in the order of several decimeters per second. In the extreme cases of GEO out-of-plane maneuvers, very large cross-track velocity changes are observed, namely 28 m/s, induced by 5.4 mm/s² thrust forces. Also, we demonstrate that by applying the estimated thrust forces in orbit integration, the orbit errors can be estimated at decimeter level in along- and cross-track directions during normal maneuver periods, and 1–2 m in all the orbital directions for the enormous GEO out-of-plane maneuver.     

A new multipath mitigation method based on adaptive thresholding wavelet denoising and double reference shift strategy

Multipath disturbance is one of the major error sources in high-accuracy positioning for global navigation satellite system (GNSS). Although various methods based on software and hardware strategies have been developed to mitigate this error, they are still limited by different kinds of factors and the effect is not ideal. After analyzing the existing methods, a new single-difference sidereal filtering method, based on adaptive thresholding wavelet denoising and double reference shift strategy (ATDR), is proposed to mitigate multipath effects for static short-baseline GNSS applications. The key idea of the proposed method is the use of both the adaptive thresholding wavelet denoising to extract an accurate multipath correction model from the reference Day and the double reference shift strategy to mitigate multipath for subsequent Day 2 more accurately and efficiently. By applying the introduced adaptive thresholding method, the average improvement rate of the RMS values of the single-difference residuals can reach about 15.82% compared with the constant thresholding method. Moreover, after applying the proposed ATDR method, the 3D positioning precision is improved by about 37.73% for the single epoch mode with 30 s data sampling rate and about 31.22% for the continuous mode with 1 s high sampling rate compared with the original results. Even compared with the constant thresholding single orbital reference (CTSR) method, the improvement percentage is about 33.94% in single epoch mode and about 25.40% in continuous mode for 3D positioning precision, respectively. In conclusion, the results of the two experiments indicate that the proposed ATDR method performs much better than the CTSR method in mitigating multipath for different sampling rates and different processing modes in the measurement domain for GNSS static short-baseline postprocessing applications.     

Accurate ionospheric delay model for real-time GPS-based positioning of LEO satellites using horizontal VTEC gradient estimation

Ionospheric delays compensation is a mandatory step for precise absolute and relative positioning of Low Earth Orbit Satellites (LEO) by GPS measurements. The most frequently used ionosphere model for real-time GPS-based navigation in LEO is an isotropic model proposed by Lear, which uses the Vertical Total Electron Content (VTEC) above the receiver and a mapping function for TEC evaluation along a given ray path. Based on significant assessed results available for ground-based GPS receivers, we propose the use of a different model relying on the thin shell assumption and a bilinear horizontal variation of the VTEC as a function of latitude and longitude in the shell. It is expected that this model is capable of better describing horizontal gradients in the ionosphere, thus improving ionospheric delay estimation, especially in intense ionospheric conditions. This model is referred to as Linear Thin Shell (LTS). LTS performance in estimating undifferenced and double-differenced ionospheric delays is checked by comparing measured and predicted delays computed using flight data from the GRACE mission. Results show that the LTS always outperforms the isotropic model, especially in case of high solar activity. Moreover, the LTS model provides a higher performance uniformity over a wide range of ionospheric delays, thus ensuring good performance in different conditions. The results obtained demonstrate that the LTS model improves the ionosphere delays estimation accuracy by 20 and 40% for undifferenced and double-differenced delays, respectively. This suggests the LTS model can effectively contribute to improving precision in LEO positioning applications.     

Carrier phase bias estimation of geometry-free linear combination of GNSS signals for ionospheric TEC modeling

The ionosphere can be modeled and studied using multi-frequency GNSS signals and their geometry-free linear combination. Therefore, a number of GNSS-derived ionospheric models have been developed and applied in a broad range of applications. However, due to the complexity of estimating the carrier phase ambiguities, most of these models are based on low-accuracy carrier phase smoothed pseudorange data. This, in turn, critically limits their accuracy and applicability. Therefore, we present a new methodology of estimating the phase bias of the scaled L1 and L2 carrier phase difference which is a function of the ambiguities, the ionospheric delay, and hardware delays. This methodology is suitable for ionospheric modeling at regional and continental scales. In addition, we present its evaluation under varying ionospheric conditions. The test results show that the carrier phase bias of geometry-free linear combination can be estimated with a very high accuracy, which consequently allows for calculating ionospheric TEC with the uncertainty lower than 1.0 TECU. This high accuracy makes the resulting ionosphere model suitable for improving GNSS positioning for high-precision applications in geosciences.     

A new GPS/BDS tropospheric delay resolution approach for monitoring deformation in super high-rise buildings

A new approach for deformation monitoring of super high-rise building using GPS/BDS technology is proposed for the case when prior coordinates are known and the baseline is short but has a large height difference. The approach is based on the ambiguity function method (AFM). Considering that the double-differenced (DD) troposphere delay residual error cannot be ignored, the relative zenith tropospheric delay (RZTD) parameter is introduced into the original AFM equation. Thus, the RZTD and 3D coordinate parameters are together obtained through the modified AFM (MAFM). Due to the low computational efficiency of conventional AFM, an improved particle swarm optimization (IPSO) algorithm is used to search the four optimal parameters X/Y/Z/RZTD and replaces the grid search method. In this study, GPS/BDS deformation monitoring data for buildings with approximately 290 m height difference were used to verify the feasibility of the proposed MAFM. Numerical results show a single-epoch average computation time of approximately 0.3 s, which meets the requirements of near-real-time dynamic monitoring. The average accuracy of the GPS single-epoch RZTD solution is better than 1 cm, the combined GPS/BDS MAFM performance outperforms the GPS-only system, and using multi-epoch observations can further improve the accuracy of the RZTD solution. After RZTD correction, GPS/BDS monitoring precision can be improved, particularly the height dimension, whose precision is improved by approximately 6 cm.