Search

Global Navigation Satellite System (GNSS) network analysis benefits from precise point positioning with undifferenced (UD) or double differenced (DD) integer ambiguity resolution (IAR). Position time series with UD-IAR and a practical DD-IAR are compared for a mixture network consisting of 57 global and 82 regional stations with 4-year observations. The sparse global stations are more affected by different IAR methods than the regional stations of Crustal Movement Observation Network of China (CMONOC). The most significant differences between UD-IAR and the practical DD-IAR are observed in velocity estimations, with averages of 0.15 and 0.30 mm/yr for the horizontal and vertical components, respectively. UD-IAR can achieve smaller root mean square at 70–83 % of the daily solutions than the practical DD-IAR for global stations. Nevertheless, moderate differences of frame parameters and annual variations are introduced by different IAR methods. The IAR is also proved to be one of the potential origins for Global Positioning System (GPS) draconitic signals and white noise at high frequency. [ABSTRACT FROM AUTHOR]

With the capability for single-receiver AR (Ambiguity Resolution), real-time GNSS (Global Navigation Satellite System) satellite integer clock has been the research hotspot of GNSS community in the past decade. Robust ambiguity datum is one prerequisite for high-quality integer clock. In this paper, with undifferenced AR model derived, a robust ambiguity datum definition method is proposed to isolate the signal biases from the integer ambiguity. Firstly, based on the ambiguity graph, the ambiguity datum is selected with the minimum spanning tree algorithm. Secondly, the ambiguity datum is transferred to the subsequent epochs through signal bias connection. Moreover, additional datum ambiguities are introduced to accommodate the arising of receiver and satellite in real-time circumstance. With the combination of ambiguity graph and signal bias, the efficient definition and smooth transition of ambiguity datum are assured for the real-time processing. To validate the approach, simulated real-time integer clock estimation for GPS/Galileo/BDS-3 satellites are conducted with one-month (January 2022) data from about 100 global IGS (International GNSS Service) network stations. Experiments show that the proposed robust ambiguity datum can produce accurate real-time integer clock. The average WL (Wide-Lane) and NL (Narrow-Lane) fixed rates exceed 99 % and 95 %, respectively, for most of the GPS/Galileo satellites. While comparable WL fixed rate is achieved for the BDS-3 satellites, the NL fixed rate is about 5 % lower as a result of poor orbit quality and imperfect measurement model. With eclipsing satellites excluded, the standard deviations between real-time integer clock and post-processed reference clock reach 0.018 ns and 0.016 ns for the GPS and Galileo constellations, respectively, which are 28 % and 24 % smaller than those of corresponding float clocks. In contrast, the standard deviations of float and integer clocks are similar for the BDS-3 constellation. Further improvements on the background models, including unconservative force modeling, satellite attitude modeling and antenna calibration are required to achieve ultimate accuracy, especially for the BDS-3 satellites. [ABSTRACT FROM AUTHOR]

With the rapid progress of Global Navigation Satellite Systems (GNSS), integer ambiguity resolution has gradually become one of the key factors that affecting Precise Point Positioning (PPP) to achieve high-precision positioning results and fast convergence. Fortunately, the proposal of observable-specific code/phase biases (OSB) has greatly alleviated this problem tremendously. Different from obtaining OSB by dual-frequency linear combinations, this research proposed BeiDou Navigation Satellite System (BDS) OSB estimation method from raw observations. This method allows for more flexible estimation in a global network, without the need to compose the complex linear combinations when involving multi-frequency measurements processing. The BDS-3 observation data collected from 112 stations over one month (October 2022), containing B1C/B1I/B2a/B2b/B3I frequencies, were used for five-frequency phase OSBs estimation. Subsequently, static and kinematic PPP were executed utilizing two different ambiguity resolution strategies to corroborate the efficacy of this method. The results revealed that the BDS-3 five-frequency phase OSBs possessed good monthly stability, with standard deviation less than 0.3 cycle. High position accuracy was achieved in the static mode by fixing the ambiguity directly, with root mean square (RMS) of 3.0, 4.8 and 9.9 mm. Meanwhile, the ambiguity fixed rate was maintained at a high level around 96 % in multi-frequency tests. Simultaneously, the kinematic mode also presented the centimeter-level dynamic positioning results and the positioning accuracy improved by 55, 31 and 17 % in the East, North and Up components, respectively. Compared to the float solutions, the adoption of multi-frequency OSBs significantly accelerated the convergence speed by 55 % on average. [ABSTRACT FROM AUTHOR]

Multi-frequency integration can accelerate the ambiguity resolution (AR) and shorten the convergence time of precise point positioning (PPP). BDS-3 satellites can broadcast legacy B1I/B3I and new B1C/B2a/B2b/B2a + b signals, which facilitates the multi-frequency PPP. An uncalibrated phase delay (UPD) estimation method and a fast AR method are extended for BDS-3 six-frequency uncombined PPP. The characteristics of UPD estimates are analyzed in terms of stability, ambiguity residuals, and usage rate. The data sets on seven consecutive days from 224 globally distributed stations are selected for UPD estimation, and those from 21 stations are employed for PPP AR. The results indicate that the static convergence time of BDS-3 six-frequency PPP ambiguity-fixed solutions (with a threshold of 10 cm) can reach 8.2, 4.0, and 13.4 min in the east, north, and up directions, respectively, and the three-dimensional (3D) convergence time is shortened by 21.1 % and 26.5 % compared with the six-frequency ambiguity-float solutions and the dual-frequency ambiguity-fixed solutions, respectively. The corresponding kinematic convergence time can reach 13.0, 7.0, and 33.2 min in three directions, respectively, and the corresponding improvement rates are 21.0 % and 33.7 %, respectively. In addition, the convergence time of BDS-3 six-frequency PPP ambiguity-fixed solutions can be slightly shortened even compared with the five-frequency case. The static and kinematic positioning accuracy of BDS-3 six-frequency PPP ambiguity-fixed solutions can reach 4, 3, and 9 mm, and 18, 14, and 38 mm in three directions, respectively, which is slightly better than other cases. [ABSTRACT FROM AUTHOR]

• An improved SICB correction model is proposed based on elevation and azimuth angles. • The parameters are estimated for each BDS frequency with node-separation of 1°. • The threshold of elevation for the new model is determined based on the statistics. • Weight-average algorithm with SNR is used to obtain the optimal MP deviations. The satellite-induced code bias (SICB) resulting in code-phase divergences of nearly 1.0 m has been observed in the code measurements of BeiDou navigation satellite system (BDS). Since the SICB is proved to be elevation- and frequency-dependent, several SICB correction models have been proposed to remove the systematic bias, but most of them only consider elevation angles, overlooking the influence of azimuth angles, which would limit the accuracy of the correction models, especially for Medium Earth Orbit (MEO) satellites. To deal with this problem, we proposed an improved SICB piecewise correction model as a function of elevation and azimuth angles and estimated parameters for each BDS un-GEO (Geostationary Earth Orbit) satellite and each available frequency with node-separation of 1°. The epoch-difference method is employed to remove the effect of stable terms, such as phase-ambiguities, receiver- and satellite-dependent hardware delays. Moreover, to reduce the effect of code noise and multipath at low-elevation conditions, the threshold of elevation angle for the improved model is determined based on the statistics of Multipath (MP) combinations. Besides, weight-average algorithm is used to obtain the optimal MP deviations and improve the accuracy of proposed correction model in consideration of signal-to-noise (SNR). The experiment results with real GNSS data show that the SICB models are more suitable for BDS-2 MEO satellites than Inclined Geosynchronous Orbit (IGSO) satellites. The code measurements on B1 signal are more susceptible to the SICB than that on B3 signal. After applying the improved SICB correction models, the RMS errors of MP series are significantly reduced for BDS-2 un-GEO satellites. However, for BDS-3 satellites, only centimeter-level SICB deviations can be detected for all frequency-bands, and the RMS of MP variations have slight changes after applying the SICB corrections, indicating that the systematic bias in MP combinations is not obvious, especially for BDS-3 MEO satellites. [ABSTRACT FROM AUTHOR]

Nowadays, the BDS-3 satellites transmit new signals B1c and B2a that are compatible with GPS and Galileo, which is meaningful to fully take advantage of the BDS positioning capability. However, further research and analysis are needed to investigate the precise point positioning ambiguity resolution (PPP-AR) performance of the new BDS-3 signal combination. In this manuscript, the characteristics and quality of differential code bias (DCB), inter frequency clock bias (IFCB) and uncalibrated phase delay (UPD) of BDS-3 new signals are analyzed based on 30 days observation data from 70 multi GNSS experiment (MGEX) stations. The results show that the mean value of IFCB between B1I/B3I and B1C/B2a is 1.18 cm, which indicates that the IFCB can be ignored when performing the clock errors reference conversion of BDS-3 satellites. Besides, the DCB of the new BDS-3 signals varies within 0.2 ns over 30 days, and with a difference of 0.08 ns compared to the DCB product provided by chinese academy of science (CAS). The percentage of the estimated residuals less than 0.25 cycles for new signals WL UPD and less than 0.15 cycles for new signals NL UPD are 99.76% and 97.19%, respectively. The above results demonstrate the reliability of the DCB and UPD products, which can be used to estimate the carrier phase observable specific signal bias (OSB) product for the new signals. The average positioning errors of the BDS-3 new signals static PPP-AR in the N, E and U directions is 0.50 cm, 0.43 cm and 1.62 cm respectively, with an improvement of 18.3%, 28.3% and 11.5% compared to the ambiguity float static PPP. The average positioning errors of the kinematic PPP-AR are 2.09 cm, 2.67 cm and 4.79 cm in three directions, respectively, with an improvement of 17.1%, 23.9% and 15.7% compared to the ambiguity float PPP. The convergence time of the static and kinematic PPP-AR are 26.36 min and 29.12 min, respectively. [ABSTRACT FROM AUTHOR]

Ambiguity resolution (AR) is critical for enhancing the orbit accuracy in precise orbit determination (POD) for low earth orbit (LEO) satellites. While orbit estimation using single-difference (SD) AR has been widely researched, the investigation of orbit estimation using undifferenced (UD) AR for LEO satellites is limited due to time-varying hardware biases at the LEO receiver end. To address this deficiency, we propose an improved UD AR method for LEO satellite orbit determination. The method employs the optimal integer datum ambiguity for a 1-day observation arc, and the random-walk clock model is utilized to transfer the integer ambiguity datum to the other epochs, which enables the arc-wise ambiguities to regain the integer property within the time frame. To validate the effectiveness of our improved method, numerical experiments are conducted. Moreover, we assess the performance of the random-walk clock model for a spaceborne ultra-stable oscillator (USO). The high frequency stability of the USO establishes the satisfactory requirements for this study. Both in the kinematic and dynamic modes, the AR success rates of our improved method are ∼ 2% higher than that of the current SD AR method, and the orbit results indicate a slight improvement. The 3-Dimensional (3D) root-mean-squares (RMS) of the orbit differences between JPL precise science orbits (PSO) and our orbits are reduced by up to 4% and 5% for kinematic and dynamic orbits, respectively. The subtle benefit is also proven by K-band ranging (KBR) system validations. In the case of the Gravity Recovery and Climate Experiment Follow-on (GFO) satellites, the results indicate that the effect of our improved UD AR can be equivalent to that of SD AR. Our improved UD AR method can be a good alternative for its superior ability to generate hardware delay at the LEO receiver end. [ABSTRACT FROM AUTHOR]

In addition to the legacy B1I and B3I signals, the BeiDou global navigation satellite system (Beidou-3) are capable of transmitting B1C/B2a/B2b/B2 several new navigation signals. The multi-frequency signals of Beidou-3 are expected to benefit precise point positioning with ambiguity resolution (PPP AR). Currently, most international GNSS service (IGS) analysis center adopt the Beidou-3 B1I/B3I ionosphere-free (IF) combination to estimate the satellite orbit and clock products, the Beidou-3 B1C/B2a/B2b/B2 phase biases should be aligned with satellite B1I/B3I clocks. Therefore, in this paper, the multi-frequency phase observable-specific phase biases (OSB) estimation method is developed for Beidou-3 system, which is compatible with satellite B1I/B3I legacy clock. Experimenting on the Beidou-3 multi-frequency observations over 30 days from the globally distributed stations, the characteristics of the estimated Beidou-3 multi-frequency phase OSB and the positioning performance of the IF kinematic PPP AR with different frequency combination are assessed. After the OSB are used in the raw observations, the results demonstrate that the wide-lane ambiguities residuals within ± 0.25 cycle are approximately 97%, 94%, 97% and 93% for Beidou-3 B1I/B3I, B1C/B2b, B1C/B3I and B1C/B2a combinations, respectively. When the OSB and satellite B1I/B3I legacy clock are used in different dual-frequency combination PPP AR tests, over 94% of the corresponding narrow-lane ambiguities are within ± 0.25 cycle. In addition, owing to the low noise of the B1C/B2a combination, the positioning performance of B1C/B2a PPP AR outperforms that of the other frequency combinations. More specifically, in case of B1I/B3I PPP AR, almost 22.7 min is required to successfully achieve 10 cm for east, north and up components, while only 18.8 min is needed for B1C/B2a PPP AR. The results show that such phase OSB, along with Beidou-3 satellite B1I/B3I legacy clock, can fully recover the PPP ambiguities based on all available frequency choices and observable combinations. [ABSTRACT FROM AUTHOR]

The current GLONASS system suffers from the inter-frequency code bias (IFCB) and the inter-frequency phase bias (IFPB) due to its employment of the frequency division multiple access (FDMA) technology. The IFCB affects the code observation's accuracy and further degrades the positioning accuracy, while the IFPB can be absorbed into the float ambiguity term. With increasing demand on precise point positioning (PPP) in mass market applications using low-cost single-frequency (SF) receivers, the impact of IFCB on SF-PPP with low-cost receivers need to be investigated. In this paper, we propose the use of the University of Calgary (UofC) model to estimate the IFCB in SF-PPP. The feasibility of the proposed approach is first verified using GNSS observations from the IGS stations. Comparing to the traditional SF-PPP model, the UofC model enables more stable IFCB estimations with 80.2% of the maximum changes less than 1.0 ns within one week. Then, the stability of IFCB in the low-cost SF receiver is analyzed with static data collected with u-blox M8T, and the STandard Deviation (STD) of IFCB for each satellite within 5 days is less than 1.5 ns when applying real-time orbit, clock and ionospheric corrections. Further, it's verified that the power off won't affect the IFCB's value significantly in the low-cost receiver. Finally, six different strategies are designed to evaluate the impact of the GLONASS IFCB on SF-PPP. Comparing to those strategies ignoring IFCB, the strategies considering IFCB in SF-PPP can improve the positioning accuracy in both static and kinematic tests. [ABSTRACT FROM AUTHOR]

Tight integration can enhance the model strength and positioning performance by considering the characteristic of differential inter-system bias (DISB), especially in obstructed environments. However, limited work emphasizes the comprehensive analysis of five-frequency DISBs between BDS-3 and other systems considering the receiver type, receiver configuration, and antenna type. In addition, the overlapping DISBs between BDS-3 and BDS-2 are also in great demand for further investigation since they are often regarded as one system. In this study, one DISB-float model is introduced to estimate the DISBs, and one DISB-fixed model and one DISB-free model are formulated to enhance the model strength of tight integration. Four dedicated datasets were collected to estimate the DISBs, which are also comprehensively analyzed considering the receiver type, receiver configuration, and antenna type. The results show that the DISBs between BDS-3 and other systems are rather stable over a certain period and are related to the receiver type and receiver configuration, whereas are not related to the antenna type. More interestingly, the B1I code DISB between BDS-3 and BDS-2 exhibits significant magnitude with a mean value of −1.44 m for the baseline composed of two different receivers. In this case, the B1I code DISB must be considered and the tight integration between BDS-3 and BDS-2 considering its calibration can improve the positioning performance. Besides, the tight integration of the DISB-fixed model can significantly improve the positioning accuracy between multiple GNSS. Compared to the loose integration, the improvement of 60.6 %, 56.6 %, and 61.2 % can be obtained in the E, N, and U directions, when only two satellites are available for each system. In real obstructed environments, the tight integration of the DISB-free model can also improve the positioning performance in terms of positioning availability and accuracy, as well as the ambiguity resolution performance. [ABSTRACT FROM AUTHOR]

• The ISB between the BDS-2 and BDS-3 is considered in UPD estimation. • The time-varying characteristics of ISB and the quality of BDS UPD are analyzed. • The satellite-induced bias of BDS-2 in UPD estimation is evaluated. • The impact of ISB is evaluated from the performance of BDS PPP and PPP-AR. Intra-system biases (ISBs) between BDS-2 and BDS-3 are of critical importance when combining observations from the BDS-2 and BDS-3 systems, which is meaningful to fully take advantage of the BDS positioning capability. Meanwhile, ISBs should also be considered in the estimation of BDS uncalibrated phase delays (UPDs). In this research, we present a BDS-2/BDS-3 joint-processing scheme, as well as a method for estimating BDS UPDs. The characteristics of ISBs and the quality of BDS UPDs are analyzed based on 30-day data from 130 multi-GNSS experimental (MGEX) stations. Our results indicate that the ISBs are related to the type and version of the receiver. The ISBs can be regarded as constant across the course of a given day, and the mean standard deviation (STD) values of ISBs over one month for different types of receivers are generally within 0.2 m. Moreover, to assess the quality of UPD products, the residuals of the estimated UPDs and the utilization rates of the observation data are computed. The results show that the quality of BDS UPDs can be improved by correcting the satellite-induced pseudo-range variations, and by estimating the wide-lane (WL) UPD difference between BDS-2 and BDS-3. The average RMS values of the estimated residuals of WL UPD and narrow-lane (NL) UPD are 0.07 and 0.09 cycles, respectively; moreover, the utilization rate of the observation data of WL UPD and NL UPD can reach above 90 %. The performance of BDS precise point positioning (PPP) and PPP ambiguity resolution (PPP-AR) is analyzed in terms of positioning accuracy and convergence performance in both the static and kinematic modes. Compared with PPP ambiguity-float solutions, the positioning accuracy of PPP-AR is significantly improved, especially in the east direction. The impact of ISBs on PPP and PPP-AR is also analyzed, and the results indicate that ISBs can improve the convergence speed of float PPP, but can be disregarded in PPP-AR. [ABSTRACT FROM AUTHOR]

Global navigation satellite system (GNSS)-based attitude determination has been widely adopted in a wide variety of terrestrial, sea, air, and space applications. Recently, the emergence of commercial multi-GNSS common-clock receivers has brought new opportunities for high-precision GNSS-based attitude determination with single-differenced (SD) model. However, the key requirement of using this approach is the accurate estimation of the troublesome line bias (LB) in real-time. In this contribution, we propose a particle filter-based real-time phase LB estimation approach that apply to SD model with single-system single-frequency observations from common-clock receiver. We first analyzed the relationship between the integer ambiguity ratio value and the phase LB. It is proved that the accuracy of a given phase LB value can be qualified by the related ambiguity resolution ratio value, and the normalized ratio value can therefore be used to represent the likelihood function of observations. Then, we presented the particle filter-based real-time phase LB estimation procedure, and assessed its performance using GPS L1/BDS B1I observations from two datasets collected with different types of common-clock receivers in terms of the accuracy and convergence time of phase LB estimation, the computation load, and the positioning and attitude determination accuracy with respect to the double-differenced (DD) model. Experimental results demonstrated that the phase LB could be accurately estimated with short convergence time (generally within 15 epochs). Moreover, compared with the classical DD approach, the particle filter-based SD approach delivers comparable positioning root-mean-square (RMS) errors in the North and East components but significantly smaller RMS errors in the Up component. Accordingly, the achievable yaw accuracy is comparable whereas the pitch accuracy is remarkably improved. The improvements of positioning accuracy in the Up component and pitch accuracy are approximately 35.7 % to 63.7 %, and 33.3 % to 63.1 %, respectively. Additionally, the single-epoch computation time with our particle filter-based SD approach is generally 0.08 s, which is obviously larger than the DD approach but could still meet the requirements of real-time applications below 10 Hz sampling. [ABSTRACT FROM AUTHOR]

• Attitude and inertial parameter are estimated without 3D model. • Ill-posed point cloud is detected and avoided. • Combination of MEKF and pose graph reduces error of estimate. Determining the attitude and inertial parameters of a noncooperative target is essential in an on-orbit servicing mission. Various methods based on machine vision have been proposed, but most of them require the 3D model of the target. This paper proposes a model-free method through sequentially registering point clouds captured by a depth camera. Our main contributions are the avoidance of the ambiguity in registration, and the combination of the multiplicative extended Kalman filter and the pose graph optimization to reduce the effect of measurement noise and drift error. A hardware experiment was performed to generate the sequence of point clouds of a three-axis free-floating target and validate our method. The result shows that the proposed method outperforms existing methods and effectively identifies the inertial parameters, including the normalized principal moments of inertia and the orientation of principal axes. [ABSTRACT FROM AUTHOR]

Global Navigation Satellite System (GNSS) precise positioning can be significantly affected by severe multipath effects and outliers in harsh environments, and highly relies on quality control strategies. Previous studies mainly focus on the posterior residuals to check and exclude the outliers in GNSS observations, limited work emphasizes the combined quality control method considering both the prior and posterior knowledge simultaneously. This paper proposed a real-time combined quality control method to process the multipath effects and outliers in harsh environments simultaneously. Specifically, in the prior stage, a modified multipath processing strategy is proposed for both phase and code observations, then a modified detection, identification, and adaptation (DIA) method considering the maximum times of data snooping is studied in the posterior stage. Two dedicated experiments in real harsh environments were carried out to evaluate the performance of the proposed combined quality control method. For the static experiment, the proposed method exhibits smaller positioning errors, the best positioning accuracy, and the highest availability in this study. Specifically, the proposed method exhibits an improved percentage of 55.4 %, 56.3 %, and 59.7 % for positioning accuracy compared to those without the quality control method in the E, N, and U directions, respectively. Besides, the proposed method can further improve the performance of ambiguity resolution with an improved percentage of 32.2 %. For the kinematic experiment, the three-dimensional positioning accuracy of the proposed method is 0.577 m, which exhibits a 40.0 % improvement compared to those without the quality control method. Also, the proposed method exhibits better performance under relatively strong multipath effects. In this sense, the proposed real-time combined quality control method is highly appreciated in terms of positioning availability, accuracy, and ambiguity resolution for GNSS precise positioning, especially in harsh environments. [ABSTRACT FROM AUTHOR]

The full operational capability of the Chinese BeiDou Satellite Navigation System (BDS) has injected additional energy into the GNSS field. Along with the visibility of more BDS satellites, as well as the precise products generated by analysis centers, ambguity float and fixed Precise Point Positioning (PPP) performance with the complete BDS constellation are evaluated in this study. Inter-system biases (ISBs) between BDS-2 and BDS-3 are first evaluated, and Fractional Cycle Bias (FCB) with and without consideration of ISBs are assessed. Results show that when considering ISBs, the standard deviations of ambiguity residuals are smaller, thus the ISB parameter should be taken into consideration in BDS PPP. BDS-3 FCBs are more stable than those of BDS-2, at the same time, the ambiguity residuals are smaller for BDS-3 satellites. In terms of PPP performance, BDS-3 outperforms BDS-2 in positioning accuracy, convergence time and time to first fix. For the test data, the convergence time of BDS-2 static solution is 48.8 min, while for BDS-2 + 3 it is 19.7 min, representing an improvement of 59.6% when BDS-3 satellites are included. The positioning accuracy of static BDS-2, BDS-3 and BDS-2 + 3 PPP can reach the same millimeter-level after a long time of convergence, while kinematic solution with BDS-2 + 3 has the highest accuracy compared to BDS-2 and BDS-3, which could reach 1.7, 1.3 and 4.3 cm in East, North and Up components, respectively, in fixed solution. The performance of BDS PPP with the whole constellation is comparable with that of GPS and Galileo, and the PPP performance is promising with receiver updates for MGEX stations in the near future. [ABSTRACT FROM AUTHOR]

Multi-GNSS multi-frequency observations improve the accuracy, availability and reliability of the single-site RTK positioning but also simultaneously raise the communication bandwidth needed for the transmission of differential corrections. A RTK positioning algorithm is proposed that transmit the phase-only differential corrections to achieve the ambiguity-fixed RTK positioning, significantly reducing the communication bandwidth requirements. The essence of the presented method is that the code differential positioning solution implied in the conventional RTK positioning procedure can be replaced by the single point positioning solution of the rover receiver to realize RTK positioning. For typical application scenarios, the transmitted bits in a single epoch can be reduced by more than 25% when the corrections are sent with the RTCM protocol. The experiments by the GPS dual-frequency data over the 10 km and 180 km baselines shown that the performance of the presented RTK algorithm are comparable to that of the conventional RTK positioning, which verifies its feasibility and effectiveness. [ABSTRACT FROM AUTHOR]

With the development of commercial services in the European Union, the advanced designs of Galileo E6 signals can be freely offered to Global Navigation Satellite System (GNSS) users. However, only a few studies have analyzed the benefit of Galileo E6 signals, and their application in the real-time instantaneous precise point positioning (PPP) with ambiguity resolution (AR) has not been investigated. In this study, based on the triple-frequency ionosphere-free combination, the positioning performance of instantaneous decimeter-level PPP AR was assessed with real-time orbit, clock and phase biases provided by the Centre National d'Études Spatiales Navigation Team (CNES/NAV). Experimenting on the GPS/Galileo observations from the 42 globally distributed stations, the wide-lane (WL) ambiguity residuals within ± 0.25 cycle for GPS L1/L2/L5, Galileo E1/E5a/E5b and E1/E5a/E6 frequency combinations were 92.54%, 90.72% and 96.48% with single-epoch data, respectively. When the WL ambiguities are fixed instantaneously, owing to the low combination noise of Galileo E1/E5a/E6, the positioning performance of instantaneous PPP AR with GPS/Galileo (E1/E5a/E6 used) outperforms that of GPS/Galileo (E1/E5a/E5b used). If the position dilution of precision is <1.5, the Root Mean Square (RMS) of GPS/Galileo (E1/E5a/E6 used) for east, north and up components are 0.13 m, 0.13 m and 0.57 m, respectively, while those of GPS/Galileo (E1/E5a/E5b used) are 0.17 m, 0.19 m and 0.62 m, respectively. Moreover, the WL ambiguity fixing rate of GPS/Galileo (E1/E5a/E6 used) is higher than that of GPS/Galileo (E1/E5a/E5b used), which shows an improvement of approximately 9% in the relatively low satellite elevation angles (20°< elevation angle < 40°). The results indicate that the Galileo E6 signal can bring a significant improvement in real-time instantaneous decimeter-level PPP AR. [ABSTRACT FROM AUTHOR]

The integrated navigation of single-frequency global navigation satellite system (GNSS) and inertial navigation system (INS) has broad application prospects with the development of multi-constellation GNSS. However, GNSS signals are susceptible to interference from a complex environment, which results in low-quality observations. The combination of frequencies cannot be applied for a single-frequency signal, so it is difficult to fix the ambiguity and process cycle slip. Therefore, the method of ambiguity resolution and cycle slip processing with INS aiding is proposed in this paper. When full ambiguities cannot be completely fixed, a subset of ambiguity is determined based on the elevation angle and posteriori residuals. The satellite is firstly eliminated one by one if the absolute value of posteriori residual is over the set threshold. And if the ambiguities still cannot be fixed, the satellite with the minimum elevation angle is eliminated each one. The short-term and high-accuracy position predicted by the INS was used as the constraint equation of the cycle slip detection. The unit weight mean square error of the cycle slip detection equation was used to identify the epoch where the cycle slip occurs, and the residual calculated by robust estimation was used to detect and repair the cycle slip. A set of data collected by a vehicle was used to analyze the proposed method. The experimental results showed that the proposed partial ambiguity resolution can obtain the optimal ambiguity fixing rate and positioning performance compared with the full ambiguity resolution and partial ambiguity resolution based on the elevation angle. The proposed algorithm can effectively detect and accurately repair the cycle slip for the simulated small and large cycle slips, and the cycle slips can still be detected when the global positioning system (GPS), global navigation satellite system of Russia (GLONASS), and BeiDou navigation satellite system (BDS) were outage for 10 s, 5 s, and 10 s without sufficient visible satellites. [ABSTRACT FROM AUTHOR]

The ionospheric delay is one of the major error sources in GNSS precise data processing. In most scenarios, only the first-order ionosphere delay was fully considered, which contributes to about 99% of the total ionosphere delay. Recent research reveals that the higher-order ionospheric (HOI) delay can also introduce up to 2 cm biases for the GPS L1 observations. This bias is generally not considered in the underlying function model and consequently, it introduces unintentional biases, which may degrade the reliability of the ambiguity estimation in long-baseline RTK. This study examined the impact of the HOI on the observations and the ambiguity estimation procedure. We considered the HOI bias in the underlying model to evaluate the unintentional HOI bias on the success rate (SR) of ambiguity estimation. The numerical results showed that the double-differenced HOI effect reaches about 7 mm in the ionospheric active period, and can affect the ambiguity estimation success rate by up to 30% in the GNSS data processing and the HOI impact is more significant when a new satellite arose. The presence of the HOI biases may change about 0.3% ambiguity fixing results. In addition, the HOI bias can also introduce several millimeters biases in positioning solutions. [ABSTRACT FROM AUTHOR]

AR (Ambiguity Resolution) is one of the critical issues in high-accuracy GNSS (Global Navigation Satellite System) applications. Commonly, the double-difference AR technique is employed in GNSS network data processing. With the demand of high-accuracy, low-latency and high-integrity GNSS orbit, the orbit estimation in sequential mode becomes a candidate of that in batch mode. Unlike orbit estimation with double-difference AR in either batch or sequential mode, researches on real-time orbit estimation with undifferenced AR in sequential mode is still limited. In this study, real-time orbit determination with epoch-independent undifferenced AR in sequential mode is assessed. With one-week observation data from about 100 global distributed IGS (International GNSS Service) network stations, simulated real-time Galileo orbit determination are conducted with GFZ post-processed solution as reference for validation. Numerical results indicate that, first, the 3D error RMS (Root Mean Square) of fixed orbit is about 30% smaller than that of float orbit; second, the Galileo orbit error RMS reach 3.66, 3.60 and 3.09 cm for the radial, tangential and normal directions respectively for the fixed solution; third, small jumps exist in the real-time orbit and may be attributed to the temporal variations of UPD (Uncalibrated Phase Delay) and wrongly fixed ambiguities. In summary, the orbit quality is comparable to that of IGS real-time solutions, and further improvements may be achieved with improved AR constraint strategy and background orbit model. [ABSTRACT FROM AUTHOR]

In road transport, continuous high-accuracy positioning is required in real time. To ensure the proper functioning and safety of vehicular applications, integrity monitoring (IM) is needed to protect from the positioning errors under a certain alert limit (AL) with a pre-defined probability of misleading information (MI). In this study, a detailed threat model is developed for real-time kinematic (RTK) positioning application of short baselines. The model distinguishes between ambiguity-float and -fixed scenarios, and considers the influences of phase and code multipath as well as between-receiver atmospheric residuals. With the float ambiguities temporally constrained, the bias contribution that propagates with time-updated ambiguities was studied analytically for the horizontal protection level (HPL) in IM. Based on real data from both static and kinematic experiments, HPL was computed along the direction of the semi-major axis of the horizontal error ellipse. In ambiguity-float and -fixed cases, the HPL was mostly several meters and decimetres, respectively. It was found that time-propagated biases play a dominant role in the ambiguity-float HPL, and among them, phase and code multipath had in general the largest contributions. For ambiguity-fixed case, the phase multipath was found to play a dominant role in the HPL. This shows the importance of considering the biases in the RTK IM for both the ambiguity-float and -fixed scenarios. Given a horizontal alert limit (HAL) of 5 m, the availabilities of ambiguity-float solutions were low, i.e., below 50% for the static roof tests and below 5% for the kinematic road tests. For the ambiguity-fixed scenario, with HAL at 0.5 m, integrity availability was nearly 100% for the static roof tests and above 85% for the kinematic road tests. [ABSTRACT FROM AUTHOR]

This article proposes a new method for uncalibrated phase delay (UPD) estimation to improve the accuracy of precise point positioning (PPP), which uses only observation station data. This means that the station used to generate the UPDs is the same station to which they are applied. First, dual-frequency observation equations based on a raw PPP model are developed. Then, the UPDs are calculated from integer linear combinations of float ambiguities. Third, with the UPD corrections, the least-squares ambiguity decorrelation adjustment (LAMBDA) method is utilized to obtain the integer ambiguities. Since only observation station data are used for UPD estimation, the partial ambiguity resolution (PAR) method is adopted to increase the possibility of finding a subset of integer ambiguities. The UPD estimation and ambiguity resolution are performed in each epoch. To obtain the correct integer ambiguity, the ratio test and success rate (bootstrapping) are used to evaluate the estimated integer ambiguity. Finally, by treating the integer ambiguities as constants, fixed solutions can be obtained. Quality control is also applied throughout the entire data processing procedure to obtain high quality float and fixed solutions. Data from 22 stations of the International Global Navigation Satellite System (GNSS) Service (IGS) in East Asia on day of year (DOY) 206, 2017, are used to verify the feasibility of this method. The experimental results show that compared with the float solution, the proposed method can significantly improve the accuracy in the east, north and up directions by 24%, 21% and 18% for static PPP and 36%, 18% and 34% for dynamic PPP, respectively. However, the accuracy of the proposed method is still lower than that of the fixed solutions obtained by the PRIDE-PPPAR software, in which the fractional cycle bias is computed based on reference network data. These findings sufficiently show that the proposed method can offer better solution accuracy than the float solution. However, the quality of the UPDs estimated only from observation station data is not as good as that of the estimates obtained based on reference network data. [ABSTRACT FROM AUTHOR]

The integer ambiguity resolution (AR) of carrier phase is significant for Global Navigation Satellite System (GNSS) precise positioning. However, in kinematic case, single-epoch AR methods based on alone GNSS are usually not reliable due to the instable pseudorange accuracy. Moreover, the computation of classical AR method Least Squares Ambiguity Decorrelation Adjustment (LAMBDA) is large. Thus, the inertial measurement unit (IMU) is introduced, a new inertial-aided AR method that directly rounds the float ambiguity of BeiDou triple-frequency combined observations, which is characterized by long wavelength, low carrier-phase noise and ionospheric delay, is proposed. The mathematical model of the new method is derived first. Then the impacts of the carrier-phase noise, ionospheric delay and inertial navigation system (INS) position error on the AR success ratio of combined observation are analyzed through probabilistic approach. Based on above investigation, the combinations (0, −1, 1), (1, 4, −5) and (4, −2, −3) are selected to resolve the original ambiguity. A vehicular integrated navigation test is performed to demonstrate the proposed method. The results show that the average AR success ratios of the three selected combinations, whose float ambiguity errors are 0.041, 0.146, 0.279 cycle respectively, are above 97.25% without regard to low-elevation C05. With respect to positioning accuracy based on our AR method when compared with IE software, the east, north, up error RMS of position are 0.042, 0.024, 0.069 m, respectively. In terms of the AR recover after the BeiDou signals outage, as long as 62 s BeiDou signal complete outage, all the ambiguities of all satellites could be re-fixed immediately. Besides, during the 90 s signals partial outage, the AR is not influenced by the position error, since the float ambiguity errors are all below half-cycle. The research of this contribution demonstrates the effectiveness of the proposed new method, which indicates it is applicable to kinematic positioning, even in BDS degraded and denied environments. [ABSTRACT FROM AUTHOR]

Obtaining reliable GNSS uncalibrated phase delay (UPD) or integer clock products is the key to achieving ambiguity-fixed solutions for real-time (RT) precise point positioning (PPP) users. However, due to the influence of RT orbit errors, the quality of UPD/integer clock products estimated with a globally distributed GNSS network is difficult to ensure, thereby affecting the ambiguity resolution (AR) performance of RT-PPP. In this contribution, by fully utilising the consistency of orbital errors in regional GNSS network coverage areas, a method is proposed for estimating regional integer clock products to compensate for RT orbit errors. Based on Centre National d'Études Spatiales (CNES) RT precise products, the regional GPS/BDS integer clock was estimated with a CORS network in the west of China. Results showed that the difference between the estimated regional and CNES global integer clock/bias products was generally less than 5 cm for GPS, whereas clock differences of greater than 10 cm were observed for BDS due to the large RT orbit error. Compared with PPP using global integer clock/bias products, the AR performance of PPP using the regional integer clock was obviously improved for four rover stations. For single GPS, the horizontal and vertical accuracies of ambiguity-fixed PPP solutions were improved by 56.2% and 45.3% on average, respectively, whereas improvements of 67.5% and 50.5% in the horizontal and vertical directions, respectively, were observed for the combined GPS/BDS situation. Based on a regional integer clock, the RMS error of a kinematic ambiguity-fixed PPP solution in the horizontal direction could reach 0.5 cm. In terms of initialisation time, the average time to first fix (TTFF) in combined GPS/BDS PPP was shortened from 18.2 min to 12.7 min. In view of the high AR performance realised with the use of regional integer clocks, this method can be applied to scenarios that require high positioning accuracy, such as deformation monitoring. [ABSTRACT FROM AUTHOR]

The main challenge in real-time precise point positioning (PPP) is that the data outages or large time lags in receiving precise orbit and clock corrections greatly degrade the continuity and real-time performance of PPP positioning. To solve this problem, instead of directly predicting orbit and clock corrections in previous researches, this paper presents an alternative approach of generating combined corrections including orbit error, satellite clock and receiver-related error with broadcast ephemeris. Using ambiguities and satellite fractional-cycle biases (FCBs) of previous epoch and the short-term predicted tropospheric delay through linear extrapolation model (LEM), combined corrections at current epoch are retrieved and weighted with multiple reference stations, and further broadcast to user for continuous enhanced positioning during outages of orbit and clock corrections. To validate the proposed method, two reference station network with different inter-station distance from National Geodetic Survey (NGS) network are used for experiments with six different time lags (i.e., 5 s, 10 s, 15 s, 30 s, 45 s and 60 s), and one set of data collected by unmanned aerial vehicle (UAV) is also used. The performance of LEM is investigated, and the troposphere prediction accuracy of low elevation (e.g., 10–20degrees) satellites has been improved by 44.1% to 79.0%. The average accuracy of combined corrections before and after LEM is used is improved by 12.5% to 77.3%. Without LEM, an accuracy of 2–3 cm can be maintained only in case of small time lags, while the accuracies with LEM are all better than 2 cm in case of different time lags. The performance of simulated kinematic PPP at user end is assessed in terms of positioning accuracy and epoch fix rate. In case of different time lags, after LEM is used, the average accuracy in horizontal direction is better than 3 cm, and the accuracy in up direction is better than 5 cm. At the same time, the epoch fix rate has also increased to varying degrees. The results of the UAV data show that in real kinematic environment, the proposed method can still maintain a positioning accuracy of several centimeters in case of 20 s time lag. [ABSTRACT FROM AUTHOR]

The quality and availability of Uncalibrated Phase Delay (UPD) solutions are crucial to the Precise Point Positioning (PPP) service, and the long-term temporal variability and its contributing factors should be better understood. In this paper, we comprehensively investigate the long-term time-varying characteristics of each UPD product respectively generated by a global and regional network and their interoperable application in PPP-AR (ambiguity resolution), the sampling of the WL and NL UPDs are daily and 30 s, respectively. Firstly, in terms of our 30 day Wide-Lane (WL) UPD products of 31 satellites, the Standard Deviation (STD) of each satellite WL UPDs ranges from 0.04 to 0.06 cycles, indicating that the long-term prediction accuracy of satellite WL UPD is sufficient for fixing Wide-Lane ambiguities. Secondly, when a satellite in eclipsing the discontinulity may corrupt the determination of Narrow-Lane (NL) UPD in form of offset, as a result of lacking or poor satellite attitude dynamic modeling. When the influence of discontinuity is removed, the STD of our estimated satellite NL UPDs is less than 0.05 cycles. Thirdly, the STD of our estimated receiver WL UPDs is mainly below 0.2 cycles, which implies that its stability is one order poorer that of the satellite. In addition, if they are used for stations in and around the network covered region, the stability of the UPD products from the CMONOC (Crustal Movement Observation Network of China) is better than that from a global network, benefit from the fact that all the CMONOC stations are equipped with the same receiver type. Finally, the PPP-AR results show that a rate of 82.9% for stations with a WL-ambiguity-fixed rate of over 90% while 69.5% for stations with an NL-ambiguity-fixed rate of over 80% can be achieved when using UPD from the global network, which is worse than that of using UPD from the CMONOC (85.7% for stations with a WL-ambiguity-fixed rate of over 90% while 75% for stations with an NL-ambiguity-fixed rate of over 80%). The results of the experiment on the UPD interoperable application in PPP show that the global network UPD products can provide a fast AR at any single station, and the convergence time is well below 25 min. Particularly, when the location of a station is in and around the regional network, our results show that the PPP results obtained using regional UPDs enable the consistent use of global UPDs. When the location of a station is far away from the regional network, using the regional UPDs can not achieve PPP-AR. Finally, the WL UPDs of the previous day is used for forecasting to estimate the NL UPDs, the stability analysis results of NL UPDs solution and positioning results are demonstrate the validity of forecasted UPD products. [ABSTRACT FROM AUTHOR]

The integer characteristics of the undifferenced ambiguity can be recovered after the correction of fractional cycle biases (FCB), which can significantly shorten the initialization time of precise point positioning (PPP). In the former work, the ambiguity resolution model of the client end needs to be consistent with the FCB products server. However, the client may face the problem that inconsistent FCB products are provided and cannot be used to fix the ambiguity. We propose a new method for client of undifferenced and uncombined PPP ambiguity resolution based on FCB product of ionosphere-free combination PPP model. The hardware delay in server of ionosphere-free combination FCB (IF FCB) and undifferenced and uncombined FCB (L1L2 FCB) are deduced in detail. Constructing the narrow-lane ambiguity by using the raw ambiguity and the wide-lane ambiguity, the FCB product of the ionosphere-free combination PPP model is used to fix the ambiguity step by step. A globally distributed network consisting of about 118 MGEX stations is selected as the server end. Two types of FCB products, i.e., IF FCB and L1L2 FCB, are generated. 50 additional MGEX stations in global are selected as client end for performance assessment. The experimental results show that the new method for client of UU PPP model based on the IF FCB products has a comparable performance with the traditional method based on L1L2 FCB products, which improves the utilization of FCB products. [ABSTRACT FROM AUTHOR]

A single-receiver integer ambiguity resolution-enabled precise point positioning (PPP-RTK) user experiences a long convergence time when the rather weak single-constellation dual-frequency ionosphere-float model is used. Nowadays, the rapid development of Global Navigation Satellite Systems (GNSS) provides a multitude of available satellites and frequencies that can serve in improving the user's model strength and, therefore, its ambiguity resolution and positioning capabilities. In this study, we provide insight into and analyze the global impact of a multi-GNSS (GPS, Galileo, BeiDou-3) multi-frequency integration on the expected ambiguity resolution and positioning performance of the ionosphere-float uncombined PPP-RTK user model, and demonstrate whether it is the increased number of satellites or frequencies, or a combination thereof, that speeds up ambiguity-resolved positioning. Moreover, we explore the capabilities of both full (FAR) and partial (PAR) ambiguity resolution, considering the full ambiguity information content with the LAMBDA method, and investigate whether PAR is an efficient solution to the multi-dimensional ambiguity case. The performance of our solutions is assessed in terms of the ambiguity success rate (ASR), the number of epochs (TTFA) to achieve both an ASR criterion and a horizontal positioning precision better than 10 cm, as well as the gain in precision improvement. Based on multi-system multi-frequency simulated data from nine globally distributed stations and a large number of kinematic solutions over a day, we found that the increase in number of frequencies enhances the ambiguity resolution performance, with PAR achieving a TTFA reduction of 70% when five instead of two Galileo frequencies are used, while the ambiguity-float solution is only slightly improved. Further, our numerical results demonstrated that the increase in number of satellites leads to an improvement in both the positioning and ambiguity resolution performance, due to the improved geometry strength. It is shown that the GPS+Galileo+BeiDou solutions can achieve a TTFA of 6.5 and 4.5 min (at 90%) on a global scale when two and three frequencies are used, respectively, without any a priori information on the ionospheric delays. Finally, we analyzed the sensitivity of the PPP-RTK user's performance to changes in the precision of the measurements. [ABSTRACT FROM AUTHOR]

Knowledge of differential inter-systems phase/code biases are critical to integrate observations from several different Global Navigation Satellite Systems. Many studies on tightly-combined model have mostly focused on the overlapping frequencies and short baselines, and fewer have focused on non-overlapping frequencies and medium-long baselines. Considering differential inter-systems code biases and inter-frequency code biases, we propose a tightly combined real-time kinematic model by a reparameterization process for non-overlapping frequencies among GPS/GLONASS/BDS systems. Compared with the traditional double-difference ambiguity parameter, we use single-difference ambiguity as the parameter to be estimated, which is not affected by reference satellite switching and wavelength limitation. The paper first analyzes the stability of differential inter-systems code biases and inter-frequency code biases with several zero and ultra-short baselines. Then, the paper tests the positioning performance under occlusion environment. The experiment shows that the tightly-combined model is more robust than the loosely-combined model and the average time to first fix is shortened by approximately 22.4%. At last, the paper compares the positioning performance of tightly-combined model, loosely-combined model, and tightly-combined model after calibration under the medium and long baselines. The experiment shows that the tightly-combined model after calibration is better than the tightly-combined model and loosely-combined model in the ambiguity resolution rate, initialization speed and positioning accuracy, respectively. Compared with intra-system loosely-combined model, the positioning accuracy of tightly-combined model after calibration is improved by 42.1%, 25.5%, and 32.3% in the three-dimensional directions, respectively. [ABSTRACT FROM AUTHOR]

The real-time precise point positioning (RT PPP) with ambiguity resolution (AR) technique has attracted increasing attention due to its high accuracy and low cost. With the availability of RT precise orbits, clocks, and multi-frequency phase bias products for multi-GNSS constellations, which are provided freely by the Centre National d'Etudes Spatiales (CNES), it is possible now to investigate RT PPP AR based on undifferenced and uncombined PPP model. However, only a few studies described these phase bias model tersely, especially in the multi-frequency and multi-system PPP AR. Moreover, the impact on inter-frequency clock bias (IFCB) for triple-frequency GPS/Galileo RT PPP by using the phase bias products have not been investigated. In this paper, the transformation between multi-frequency phase bias model and integer recovery clock (IRC) model was formulated where IFCBs are well considered. Meanwhile, the convergence time, positioning accuracy, narrow-lane ambiguity N1 fixing rate of GPS/Galileo PPP AR were analyzed by using a modified version of the Precise Point Positioning with Integer and Zero-difference Ambiguity Resolution Demonstrator (PPP-WIZARD). Experiments with GPS/Galileo observation from 50 stations over 32 days were performed in static and kinematic modes. Data were processed with four strategies: dual-frequency float and fixed PPP, triple-frequency float and fixed PPP. Results showed that the effect of the IFCB of GPS Block IIF could be mitigated by using the phase bias products from CNES. Among four solutions, the triple-frequency PPP AR achieved the fastest convergence time of 10.43 and 12.55 min in static and kinematic modes, respectively. The average RMS values are 1.2, 1.5, 3.1 cm in static mode while 2.0, 2.5, 3.8 cm in kinematic mode with the 3 h observation session by triple-frequency PPP AR in the east, north and up directions, respectively. Compared with dual-frequency PPP AR, triple-frequency PPP AR could contribute to improving the performance of convergence time and positioning accuracy, especially when few satellites can be observed. [ABSTRACT FROM AUTHOR]

Precise point positioning with ambiguity resolution (PPP-AR) is a useful tool for high-precision geodetic and geophysical applications, while phase bias products are the prerequisite to implement PPP-AR. Wuhan University has been providing the final (the best operationally post-processing solution based) phase clock/bias products with a latency of two weeks since March of 2019, while a dedicated open-source software package PRIDE PPP-AR is released to leverage these products for high-precision positioning. In order to satisfy some both time and precision critical applications, such as rapid earthquake response, Wuhan University also released rapid (with comparable quality but with much shorter delivery latency) phase clock/bias products with a latency of less than 24 h and updated PRIDE PPP-AR in July 2019. We first introduce the phase clock/bias generation and validation schemes and the maintenance of routine products provision. Then, with 14 days (July 2 to July 15 in 2019) of GPS data collected from 146 globally distributed IGS (International GNSS Service) stations, we evaluated the positioning performance of the rapid products with respect to their final counterparts. It is found that positioning precision of PPP-AR using rapid products is comparable to that using final products, especially in kinematic positioning mode. When rapid products are used, the RMS of PPP-AR in static mode with respect to IGS weekly solutions can reach 1.7 mm, 1.8 mm and 5.5 mm in the east, north and up components, respectively. Furthermore, the RMS of epoch-wise positions with respect to daily solutions for the east, north and up components are 0.51 cm, 0.57 cm and 1.51 cm for PPP-AR with rapid products in kinematic mode. It demonstrates that the rapid phase clock/bias products can sufficiently meet the precision requirement of most geodetic and geophysical applications yet with much shorter time delay. Finally, we study the July 6th M7.1 2019 Ridgecrest, California earthquake using the rapid phase clock/bias products and demonstrate their comparable performance against the final products. [ABSTRACT FROM AUTHOR]

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

The Ambiguity Dilution of Precision (ADOP) is a well-known scalar measure that can be used to infer the strength of the Global Navigation Satellite System (GNSS) model of the carrier phase ambiguities involved in precise relative GNSS positioning. Odijk and Teunissen (2008a) derived closed-form expressions for single-baseline GNSS models that allow verifying the factors affecting the ambiguity resolution as well as the probability of its correct resolution as integer values. However, this weighted-ionosphere ADOP closed form assumes that the standard deviation of the ionosphere delay is dependent on a function associated with the baseline length. This means that a baseline of the same length at different locations around the world and in different moments in time would have the same standard deviation, which in practice is not true. This becomes even worse in regions like Brazil, where the ionospheric anomalies are more intense and frequent, especially in periods of high solar activity. In this work, a new method for calculating the ionospheric delay standard deviation that considers the ionospheric scintillation S4 index is proposed to improve the closed-form ADOP performance. Experiments with a baseline of 280 m located in Presidente Prudente, Brazil (magnetic latitude of around −13°), were carried out. The results showed that in periods of both weak and strong ionospheric scintillation, the introduction of S4 into the ionospheric delay standard deviation was beneficial for ADOP. The average improvement of the ADOP closed form was of around 72% in the analyzed period of weak scintillation and of 35% in the period of strong scintillation. In addition, the results showed that in periods of strong scintillation, relative positioning accuracy is around 100 times worse, up to 30 m, when compared to periods of weak ionospheric scintillation. [ABSTRACT FROM AUTHOR]

Precise orbit determination was investigated for a satellite constellation comprised of two different missions, the CHAllenging Minisatellite Payload (CHAMP) satellite and the Gravity Recovery And Climate Experiment (GRACE) twin satellites. The orbital planes of these two missions aligned closely during March to May 2005, allowing precise baseline determinations between the associated three satellites based on their onboard BlackJack Global Positioning System (GPS) receivers. The GRACE-A/B satellites fly in tandem formation with a baseline of around 220 km, whereas the baselines between CHAMP and the GRACE tandem vary from about 110 to 7500 km during 24-h orbital arcs centered around the points of closest approaches. A number of factors had to be dealt with for orbit determinations, including the cross-talk between the CHAMP GPS main navigation and occultation antennas, the different levels of non-gravitational accelerations, and the rapidly changing geometry that complicates the fixing of integer ambiguities for the GPS carrier-phase observations. Quality assessments of the orbit solutions were based on comparisons with Satellite Laser Ranging (SLR) observations, best orbit solutions had a precision of typically 1.7–2.3 cm. Consistency checks between reduced-dynamic and kinematic orbit solutions were done. For the GRACE baselines, the reduced-dynamic/kinematic baseline consistency was typically better than 1 cm, with an ambiguity fixing success rate of around 94 %. The agreement with the K/Ka-Band Radar Ranging (KBR) measurements was about 0.6 mm. For the CHAMP/GRACE pairs, the reduced-dynamic/kinematic baseline consistency varied from 0.5 to 2.5 cm, where better consistency was obtained for shorter arcs. [ABSTRACT FROM AUTHOR]

Precise point positioning (PPP) usually takes about 30 min to obtain centimetre-level accuracy, which greatly limits its application. To address the drawbacks of convergence speed and positioning accuracy, we develop a PPP model with integrated GPS and BDS observations. Based on the method, stations with global coverage are selected to estimate the fractional cycle bias (FCB) of GPS and BDS. The short-term and long-term time series of wide-lane (WL) FCB, and the single day change of narrow-lane (NL) FCB are analysed. It is found that the range of GPS and BDS non-GEO (IGSO and MEO) WL FCB is stable at up to a 30-day-time frame. At times frame of up to 60 days, the stability is reduced a lot. Whether for short-term or long-term, the changes in the BDS GEO WL FCB are large. Moreover, BDS FCB sometimes undergoes a sudden jump. Besides, 17 and 10 stations were used respectively to investigate the convergence speed and positioning errors with six strategies: BDS ambiguity-float PPP (Bfloat), GPS ambiguity-float PPP (Gfloat), BDS/GPS ambiguity-float PPP (BGfloat), BDS ambiguity-fixed PPP (Bfix), GPS ambiguity-fixed (Gfix), and BDS/GPS ambiguity-fixed (BGfix). The average convergence speed of the ambiguity-fixed solution is greatly improved compared with the ambiguity-float solution. In terms of the average convergence time, the Bfloat is the longest and the BGfix is the shortest among these six strategies. Whether for ambiguity-float PPP or ambiguity-fixed PPP, the convergence reduction time in three directions for the combined system is the largest compared with the single BDS. The average RMS value of the Bfix in three directions (easting (E), northing (N), and up (U)) are 2.0 cm, 1.5 cm, and 5.9 cm respectively, while those of the Gfix are 0.8 cm, 0.5 cm, and 1.7 cm. Compared with single system, the BDS/GPS combined ambiguity-fixed system (BGfix) has the fastest convergence speed and the highest accuracy, with average RMS as 0.7 cm, 0.5 cm, and 1.9 cm for the E, N, U components, respectively. [ABSTRACT FROM AUTHOR]

Abstract Ambiguity fixing to integer numbers of the phase measurements has been proven to ameliorate the accuracy of GNSS data processing. Until now there are two strategies to achieve ambiguity resolution; by forming double differences and by using undifferenced phase measurements. In this article, the potentiality of using the latter method for fixing the phase ambiguities for the Galileo system combined with GPS measurements in a Multi-GNSS solution is proved. The integer property of the Galileo phase clocks is demonstrated. GPS and Galileo phase fixed orbit and "integer" clock products are compared to the float solutions. Both orbit overlaps and orbit validation using satellite laser ranging (SLR) validation methods showed that there is an improvement mainly in the normal and the along track direction. [ABSTRACT FROM AUTHOR]

Abstract The performance of Real-Time (RT) Precise Point Positioning (PPP) ambiguity resolution can be further improved by incorporating multiple Global Navigation Satellite System (GNSS) observations. RT PPP with ambiguity resolution based on raw GNSS observations, aided by the satellite clock offsets and GPS un-calibrated phase delays (UPDs), is realized to evaluate the effect of GNSS observations. Compared with previous studies, Galileo observations are integrated in this analysis and the effect of new BeiDou/Galileo systems in RT PPP ambiguity resolution is analyzed for global users. Comparing with the GPS-only and GPS/GLONASS solution, the time to first fix (TTFF) is 18.8% and 7.9% improved for those users in the Asia-Pacific regions, while it is 9.6% and 3.8% improved for other users by adding the BeiDou/Galileo observations. The improvement of the fixing rate and positioning accuracy is similar for users in different regions. The fixing rate is 2.1% and 0.5% improved when adding BeiDou/Galileo observations on top of GPS and GPS/GLONASS respectively. The accuracy of the north, east and up components are 0.4, 0.0 and 0.5 mm improved when adding on top of GPS, while 0.2, 0.1 and 0.4 mm when adding on top of GPS/GLONASS. By utilizing observations of four systems, an accuracy of 9.0, 9.1 and 28.3 mm is achieved for the fixed kinematic PPP solution after a mean initialization time of 797.9 s. [ABSTRACT FROM AUTHOR]

Rapid initialization with ambiguity fixed resolution plays a key role in real-time Global Navigation Satellite Systems (GNSS) displacement monitoring of bridges. In this study, we propose a rapid initialization method with triple-frequency (TF) BeiDou Navigation Satellite System (BDS) and Global Positioning System (GPS) observations. The key step of this method is to form extra-wide-lane (EWL) ambiguity with TF BDS and GPS carrier-phase to reduce the ambiguity search space. In order to improve the speed of ambiguity resolution (AR) further, the strategies of fixing ambiguities in groups sequentially from long wavelength to the short and partial AR (PAR) with elevation and signal-to-noise ratio (SNR) criterions were proposed. Based on the real monitoring data on the Baishazhou Yangtze River Bridge, the AR efficiency was tested. The results illustrate that over 90% EWL ambiguities can be fixed by only single-epoch with TF BDS + GPS. After EWL ambiguity is fixed, the wide-lane (WL) ambiguity success rate could improve from 20% with dual-frequency (DF) BDS + GPS data to 70% with TF BDS + GPS data, 50–80% and 80–90% for the three stations respectively. Since the area of ambiguity search space is reduced by TF observations, the ambiguity search time decrased significantly compared with the DF BDS + GPS. In the Time-to-first-fix (TTFF) test, over 85% of epochs can complete the initialization process by only one epoch with one hour of 10 Hz TF BDS + GPS data. However, DF BDS + GPS has only 25%. As for the initialization time longer than two epochs, the counts of TTFFs of TF BDS + GPS are significantly less than that of DF BDS + GPS. Meanwhile, the integrated BDS and GPS can further enhance the geometry of satellites, which will be better for the precision improvement than the increase of signal frequencies. [ABSTRACT FROM AUTHOR]

The key problem of dynamic single epoch GPS positioning is the rapid decomposition of ambiguity. In some special dynamic positioning, the coordinates of the GPS receiver placed on the rover station often have the constraint condition, if this useful information is taken into account, it will help to determine the ambiguity accurately. For this reason, the straight motion trajectory of GPS receiver is used as the constraint condition for single epoch ambiguity resolution, in this algorithm, the wide lane ambiguities are calculated by using the MW combination, and the corresponding alternative combinations of wide lane ambiguity are established, which are determined according to the minimum residual sum of squares criterion firstly, it is used to lay the foundation for the subsequent L 1 and L 2 frequency ambiguity solution. And then the straight trajectory constraint condition and the phase observation equation are used to solve the integer ambiguity N 1 . When the constraint information of the trajectory is unknown, the space linear equation is established as the constraint condition according to the coordinate solution by wide lane relative positioning to solve the unknown problem. By designing the corresponding experiments, the results of this algorithm are proved to be feasible by using different experimental schemes, such as without constraint condition, the constraint condition with known prior information, introducing dynamic constraint condition and so on, to solve the ambiguity, and to compare and analyze the effect of linear constraint condition on single epoch ambiguity resolution. [ABSTRACT FROM AUTHOR]

Correct ambiguity resolution (AR) is the key for high-precision positioning. In order to reliably resolve ambiguity, an improved single-epoch segmented ambiguity algorithm is proposed in this paper. First, the double-difference wide-lane ambiguities resolved by the Melbourne-Wubbena (MW) method are divided into two categories based on signal-to-noise ratio (SNR). One is the primary ambiguity group that can be resolved easily, and the other is the secondary ambiguity group that is difficult to resolve. Second, the primary ambiguity group is divided into different levels in a (4,1,1…1) pattern based on SNR, and the different levels’ ambiguities are resolved successively, and the method being that the integer ambiguity is selected from two candidate values built by rounding up and down. Third, the resolved primary ambiguity group is used as a constraint to solve the secondary ambiguity group. Finally, the wide-lane observation equation and the basic observation equation are combined to solve the ambiguity of the L1 carrier phase, and the resolved ambiguity can be obtained by directly rounding. The candidate values of the primary ambiguity group is limited to two cycles through this step, and the (4,1,1…1) pattern can reduce ambiguity search times. Since an appropriate selection enables the correct resolution of primary wide-lane ambiguities, the selection of the threshold for SNR is the key to group wide-lane ambiguity, and we use the data from Novatel receiver as an example to discuss how to confirm the SNR threshold in this paper. Static short baseline data and kinematic short baseline that were less than 20 km were collected to validate the proposed algorithm. These results were compared with the commercial software IE8.60. These experiments show that the successful rate reaches 100%, 99.78% and 99.42%, respectively, for the new single-epoch segmented ambiguity method in the three data sets, and the accuracy is similar to that of IE8.60 for the combined double difference. [ABSTRACT FROM AUTHOR]

Integer ambiguity fixing with uncalibrated phase delay (UPD) products can significantly shorten the initialization time and improve the accuracy of precise point positioning (PPP). Since the tracking arcs of satellites and the behavior of atmospheric biases can be very different for the reference networks with different scales, the qualities of corresponding UPD products may be also various. The purpose of this paper is to comparatively investigate the influence of different scales of reference station networks on UPD estimation and user ambiguity resolution. Three reference station networks with global, wide-area and local scales are used to compute the UPD products and analyze their impact on the PPP-AR. The time-to-first-fix, the unfix rate and the incorrect fix rate of PPP-AR are analyzed. Moreover, in order to further shorten the convergence time for obtaining precise positioning, a modified partial ambiguity resolution (PAR) and corresponding validation strategy are presented. In this PAR method, the ambiguity subset is determined by removing the ambiguity one by one in the order of ascending elevations. Besides, for static positioning mode, a coordinate validation strategy is employed to enhance the reliability of the fixed coordinate. The experiment results show that UPD products computed by smaller station network are more accurate and lead to a better coordinate solution; the PAR method used in this paper can shorten the convergence time and the coordinate validation strategy can improve the availability of high precision positioning. [ABSTRACT FROM AUTHOR]

Precision Orbit Determination (POD) is a prerequisite for the success of many Low Earth Orbiting (LEO) satellite missions. With high-quality, dual-frequency Global Positioning System (GPS) receivers, typically precisions of the order of a few cm are possible for single-satellite POD, and of a few mm for relative POD of formation flying spacecraft with baselines up to hundreds of km. To achieve the best precision, the use of Phase Center Variation (PCV) maps is indispensable. For LEO GPS receivers, often a-priori PCV maps are obtained by a pre-launch ground campaign, which is not able to represent the real space-borne environment of satellites. Therefore, in-flight calibration of the GPS antenna is more widely conducted. This paper shows that a further improvement is possible by including the so-called Code Residual Variation (CRV) maps in absolute/undifferenced and relative/Double-differenced (DD) POD schemes. Orbit solutions are produced for the GRACE satellite formation for a four months test period (August-November, 2014), demonstrating enhanced orbit precision after first using the in-flight PCV maps and a further improvement after including the CRV maps. The application of antenna maps leads to a better consistency with independent Satellite Laser Ranging (SLR) and K-band Ranging (KBR) low-low Satellite-to-Satellite Tracking (ll-SST) observations. The inclusion of the CRV maps results also in a much better consistency between reduced-dynamic and kinematic orbit solutions for especially the cross-track direction. The improvements are largest for GRACE-B, where a cross-talk between the GPS main antenna and the occultation antenna yields higher systematic observation residuals. For high-precision relative POD which necessitates DD carrier-phase ambiguity fixing, in principle frequency-dependent PCV maps would be required. To this aim, use is made of an Extended Kalman Filter (EKF) that is capable of optimizing relative spacecraft dynamics and iteratively fixing the DD carrier-phase ambiguities. It is found that PCV maps significantly improve the baseline solution. CRV maps slightly enhance the baseline precision, more significantly they lead to a much better initialization of the ambiguity fixing. The GRACE single-satellite orbit solutions compare to within a few cm 3-dimensionally with state-of-the-art external orbit solutions and SLR observations, whereas for the baseline a consistency of better than 0.7 mm with KBR observations is achieved. [ABSTRACT FROM AUTHOR]

The mitigation of ionospheric delay is still of crucial interest in GNSS positioning, especially in precise solutions such as instantaneous RTK positioning. Thus, several effective algorithms and functional models were developed, and also numerous investigations of ionospheric correction properties in RTK positioning have been performed so far. One of the most highly effective approaches in precise relative positioning is the application of the ionosphere-weighted model with network-derived corrections. This contribution investigates the impact of the accuracy of the network ionospheric corrections on time-to-fix in RTK-OTF positioning. Also, an attempt has been made to estimate the desirable accuracy of the network ionospheric corrections, allowing for reliable instantaneous ambiguity resolution. The experiment is based on a multi-baseline GPS RTK positioning supported with network-derived ionospheric corrections for medium length baselines. The results show that in such scenario, the double-differenced ionospheric correction residuals should not exceed ∼1/3 of the L1 wavelength for successful single-epoch ambiguity resolution. [ABSTRACT FROM AUTHOR]

The Extended Kalman filter is widely used for its robustness and simple implementation. Parameters estimated for solving dynamical systems usually require certain time to converge and need to be smoothed by a dedicated algorithms. The purpose of our study was to implement smoothing algorithms for processing both code and carrier phase observations with Precise Point Positioning method. We implemented and used the well known Rauch–Tung–Striebel smoother (RTS). It has been found out that the RTS suffer from significant numerical instability in smoothed state covariance matrix determination. We improved the processing with algorithms based on Singular Value Decomposition, which was more robust. Observations from many permanent stations have been processed with final orbits and clocks provided by the International GNSS service (IGS), and the smoothing improved stability and precision in every cases. Moreover, (re)convergence of the parameters were always successfully eliminated. [ABSTRACT FROM AUTHOR]

Abstract: Ambiguity resolution (AR) is a critical step for successful attitude determination using carrier phase measurements of a satellite navigation system such as Beidou. This paper proposes an improved method for AR in support of Beidou attitude determination based on the concept of a “virtual baseline”. In the traditional long-short baseline method, the short baseline is limited to a length less than half of the carrier wave length of the Beidou signals. In the proposed method, a virtual short baseline is formed by differencing two collinear baselines. The AR equations for virtual short and long baselines are derived and the factors impacting the AR accuracy are analysed. Numerical simulation studies were carried out to evaluate the performance of the proposed AR method. The simulation results confirmed that the proposed method is an improvement over the traditional approach -- not only is it easier to deploy collinear antennas but also it keeps the capability of epoch-by-epoch AR, which makes it immune to cycle slips and there is no need for initialisation of ambiguity searching. [Copyright &y& Elsevier]

Abstract: In this paper, Global Positioning System-based (GPS) Orbit Determination (OD) for the KOrea-Multi-Purpose-SATellite (KOMPSAT)-2 using single- and double-differenced methods is studied. The requirement of KOMPSAT-2 orbit accuracy is to allow 1m positioning error to generate 1-m panchromatic images. KOMPSAT-2 OD is computed using real on-board GPS data. However, the local time of the KOMPSAT-2 GPS receiver is not synchronized with the zero fractional seconds of the GPS time internally, and it continuously drifts according to the pseudorange epochs. In order to resolve this problem, an OD based on single-differenced GPS data from the KOMPSAT-2 uses the tagged time of the GPS receiver, and the accuracy of the OD result is assessed using the overlapping orbit solution between two adjacent days. The clock error of the GPS satellites in the KOMPSAT-2 single-differenced method is corrected using International GNSS Service (IGS) clock information at 5-min intervals. KOMPSAT-2 OD using both double- and single-differenced methods satisfies the requirement of 1-m accuracy in overlapping three dimensional orbit solutions. The results of the SAC-C OD compared with JPL’s POE (Precise Orbit Ephemeris) are also illustrated to demonstrate the implementation of the single- and double-differenced methods using a satellite that has independent orbit information available for validation. [Copyright &y& Elsevier]

Abstract: Next-generation Very Long Baseline Interferometry (VLBI) system designs are aiming at 1mm global position accuracy. In order to achieve this, it is not only necessary to deploy improved VLBI systems, but also to develop analysis strategies that take full advantage of the observations taken. Since the new systems are expected to incorporate four independent radio frequency bands, it should be feasible to resolve phase ambiguities directly from post-correlation data, providing roughly an order of magnitude improvement in precision of the delay observable. As the unknown ambiguities are of integer nature, it is discussed here how they the can be resolved analytically using algorithms which have been developed for Global Navigation Satellite System (GNSS) applications. Furthermore, it will be shown that ionosphere contribution and source structure effects, so-called core-shifts, can be solved simultaneously with the delay, which is the main geodetic observable for follow-on analysis. In order to verify the proposed algorithm, simulated observations were created using parameters from actual design studies. It is shown that, even in the case of low signal-to-noise ratio observations, reliable phase ambiguity resolution can be achieved and it is discussed how the integer ambiguity recovery depends on the number of observations and signal-to-noise ratio. [Copyright &y& Elsevier]