In this contribution, an adaptive autoregressive model is proposed and developed to predict global ionospheric vertical total electron content maps (VTEC). Specifically, the spherical harmonic (SH) coefficients are predicted based on the autoregressive model, and the order of the autoregressive model is determined adaptively using the F-test method. To test our method, final CODE and IGS global ionospheric map (GIM) products, as well as altimeter TEC data during low and mid-to-high solar activity period collected by JASON, are used to evaluate the precision of our forecasting products. Results indicate that the predicted products derived from the model proposed in this paper have good consistency with the final GIMs in low solar activity, where the annual mean of the root-mean-square value is approximately 1.5 TECU. However, the performance of predicted vertical TEC in periods of mid-to-high solar activity has less accuracy than that during low solar activity periods, especially in the equatorial ionization anomaly region and the Southern Hemisphere. Additionally, in comparison with forecasting products, the final IGS GIMs have the best consistency with altimeter TEC data. Future work is needed to investigate the performance of forecasting products using the proposed method in an operational environment, rather than using the SH coefficients from the final CODE products, to understand the real-time applicability of the method.

In May 2016, the availability of GNSS raw measurements on smart devices was announced by Google with the release of Android 7. It means that developers can access carrier-phase and pseudorange measurements and decode navigation messages for the first time from mass-market Android-devices. In this paper, an improved Hatch filter algorithm, i.e., Three-Thresholds and Single-Difference Hatch filter (TT-SD Hatch filter), is proposed for sub-meter single point positioning with raw GNSS measurements on Android devices without any augmentation correction input, where the carrier-phase smoothed pseudorange window width adaptively varies according to the three-threshold detection for ionospheric cumulative errors, cycle slips and outliers. In the mean time, it can also eliminate the inconsistency of receiver clock bias between pseudorange and carrier-phase by inter-satellite difference. To eliminate the effects of frequent smoothing window resets, we combine TT-SD Hatch filter and Kalman filter for both time update and measurement update. The feasibility of the improved TT-SD Hatch filter method is then verified using static and kinematic experiments with a Nexus 9 Android tablet. The result of the static experiment demonstrates that the position RMS of TT-SD Hatch filter is about 0.6 and 0.8 m in the horizontal and vertical components, respectively. It is about 2 and 1.6 m less than the GNSS chipset solutions, and about 10 and 10 m less than the classical Hatch filter solution, respectively. Moreover, the TT-SD Hatch filter can accurately detect the cycle slips and outliers, and reset the smoothed window in time. It thus avoids the smoothing failure of Hatch filter when a large cycle-slip or an outlier occurs in the observations. Meanwhile, with the aid of the Kalman filter, TT-SD Hatch filter can keep continuously positioning at the sub-meter level. The result of the kinematic experiment demonstrates that the TT-SD Hatch filter solution can converge after a few minutes, and the 2D error is about 0.9 m, which is about 64%, 89%, and 92% smaller than that of the chipset solution, the traditional Hatch filter solution and standard single point solution, respectively. Finally, the TT-SD Hatch filter solution can recover a continuous driving track in this kinematic test.

No abstract

Rapid precise point positioning ambiguity resolution (PPP-AR) is of great importance to improving precise positioning efficiency. There is an expectation that Galileo multi-frequency (three or more frequencies) data processing will offer a promising way to accelerate PPP-AR. However, the performance of different combination observables out of raw Galileo multi-frequency data is still unclear, and the adverse impacts of missing receiver antenna phase center corrections have not been quantified in detail. We therefore studied uncombined Galileo PPP-AR by contrasting three typical triple-frequency combinations, which are E1/E5a/E5b, E1/E5a/E6, and E1/E5/E6 signals, using 30 days of data from 15 stations across Australia. We carried out triple-frequency PPP-AR by separately applying the official GPS receiver antenna phase centers, as currently employed in most relevant literatures, as well as the pilot Galileo receiver antenna phase centers preliminarily measured by the International GNSS Service. We found that, compared to dual-frequency (E1/E5a) PPP-AR, triple-frequency PPP-AR based on E1/E5a/E5b signals shortened the convergence time by only 7.6%, while those based on E1/E5a/E6 and E1/E5/E6 increased unexpectedly the convergence time by 17.6% and 12.7%, respectively, if the GPS receiver antenna corrections were presumed for Galileo signals. However, after using the pilot Galileo phase center corrections, triple-frequency PPP-AR based on E1/E5a/E5b, E1/E5a/E6, and E1/E5/E6 signals could speed up the convergence on average by about 16.2%, 30.3%, and 17.7%, respectively. Therefore, we demonstrate the critical impact of correct Galileo receiver antenna phase centers on multi-frequency PPP-AR convergences. Moreover, the triple-frequency signal combination E1/E5a/E6 is advantageous over others in achieving rapid triple-frequency Galileo PPP-AR.

Observing subdaily surface deformations is important to the interpretation of rapidly developing transient events. However, it is not known whether Global Navigation Satellite System (GNSS) is able to identify millimeter-level transient displacements over various subdaily timescales. We studied non-tidal ocean loading (NTOL) using 18 GNSS stations along the southern North Sea for November-December 2013 and compared 3-h GPS/GLONASS displacements with NTOL predictions. It was found that they overall agreed well with a mean correlation coefficient of 0.6 and their vertical differences had an RMS of 5.7 mm, but a 10-mm subsidence prediction for December 5th could only be marginally detected. Hence, the spatial coherence among the loading signatures at the 18 stations was harnessed to improve subdaily GNSS, and then the predicted displacements of 5-10-mm over the subdaily timescales could be discriminated successfully. We envision that adding Galileo/BeiDou signals to GPS/GLONASS can further improve the resolution of subdaily GNSS, which can also enhance the spatial coherence of transient signals captured by regional GNSS stations. GENG ET AL.