A Decentralized Processing Schema for Efficient and Robust Real-time Multi-GNSS Satellite Clock Estimation

Real-time multi-GNSS precise point positioning (PPP) requires the support of high-rate satellite clock corrections. Due to the large number of ambiguity parameters, it is difficult to update clocks at high frequency in real-time for a large reference network. With the increasing number of satellites of multi-GNSS constellations and the number of stations, real-time high-rate clock estimation becomes a big challenge. In this contribution, we propose a decentralized clock estimation (DECE) strategy, in which both undifferenced (UD) and epoch-differenced (ED) mode are implemented but run separately in different computers, and their output clocks are combined in another process to generate a unique product. While redundant UD and/or ED processing lines can be run in offsite computers to improve the robustness, processing lines for different networks can also be included to improve the clock quality. The new strategy is realized based on the Position and Navigation Data Analyst (PANDA) software package and is experimentally validated with about 110 real-time stations for clock estimation by comparison of the estimated clocks and the PPP performance applying estimated clocks. The results of the real-time PPP experiment using 12 global stations show that with the greatly improved computational efficiency, 3.14 cm in horizontal and 5.51 cm in vertical can be achieved using the estimated DECE clock

Retrieving Precipitable Water Vapor From Shipborne Multi‐GNSS Observations

©2019. American Geophysical UnionPrecipitable water vapor (PWV) is an important parameter for climate research and a crucial factor to achieve high accuracy in satellite geodesy and satellite altimetry. Currently Global Navigation Satellite System (GNSS) PWV retrieval using static Precise Point Positioning is limited to ground stations. We demonstrated the PWV retrieval using kinematic Precise Point Positioning method with shipborne GNSS observations during a 20‐day experiment in 2016 in Fram Strait, the region of the Arctic Ocean between Greenland and Svalbard. The shipborne GNSS PWV shows an agreement of ~1.1 mm with numerical weather model data and radiosonde observations, and a root‐mean‐square of ~1.7 mm compared to Satellite with ARgos and ALtiKa PWV. An improvement of 10% is demonstrated with the multi‐GNSS compared to the Global Positioning System solution. The PWV retrieval was conducted under different sea state from calm water up to gale. Such shipborne GNSS PWV has the promising potential to improve numerical weather forecasts and satellite altimetry

Automatic Calibration of Process Noise Matrix and Measurement Noise Covariance for Multi-GNSS Precise Point Positioning

The Expectation-Maximization algorithm is adapted to the extended Kalman filter to multiple GNSS Precise Point Positioning (PPP), named EM-PPP. EM-PPP considers better the compatibility of multiple GNSS data processing and characteristics of receiver motion, targeting to calibrate the process noise matrix Qt and observation matrix Rt, having influence on PPP convergence time and precision, with other parameters. It is possibly a feasible way to estimate a large number of parameters to a certain extent for its simplicity and easy implementation. We also compare EM-algorithm with other methods like least-squares (co)variance component estimation (LS-VCE), maximum likelihood estimation (MLE), showing that EM-algorithm from restricted maximum likelihood (REML) will be identical to LS-VCE if certain weight matrix is chosen for LS-VCE. To assess the performance of the approach, daily observations from a network of 14 globally distributed International GNSS Service (IGS) multi-GNSS stations were processed using ionosphere-free combinations. The stations were assumed to be in kinematic motion with initial random walk noise of 1 mm every 30 s. The initial standard deviations for ionosphere-free code and carrier phase measurements are set to 3 m and 0.03 m, respectively, independent of the satellite elevation angle. It is shown that the calibrated Rt agrees well with observation residuals, reflecting effects of the accuracy of different satellite precise product and receiver-satellite geometry variations, and effectively resisting outliers. The calibrated Qt converges to its true value after about 50 iterations in our case. A kinematic test was also performed to derive 1 Hz GPS displacements, showing the RMSs and STDs w.r.t. real-time kinematic (RTK) are improved and the proper Qt is found out at the same time. According to our analysis despite the criticism that EM-PPP is very time-consuming because a large number of parameters are calculated and the first-order convergence of EM-algorithm, it is a numerically stable and simple approach to consider the temporal nature of state-space model of PPP, in particular when Qt and Rt are not known well, its performance without fixing ambiguities can even parallel to traditional PPP-RTK

Mitigation of Unmodeled Error to Improve the Accuracy of Multi-GNSS PPP for Crustal Deformation Monitoring

High-rate multi-constellation global navigation satellite system (GNSS) precise point positioning (PPP) has been recognized as an efficient and reliable technique for large earthquake monitoring. However, the displacements derived from PPP are often overwhelmed by the centimeter-level noise, therefore they are usually unable to detect slight deformations which could provide new findings for geophysics. In this paper, Global Positioning System (GPS), GLObalnaya NAvigatsionnaya Sputnikovaya Sistema (GLONASS), and BeiDou navigation satellite system (BDS) data collected during the 2017 Mw 6.5 Jiuzhaigou earthquake were used to further exploit the capability of BDS-only and multi-GNSS PPP in deformation monitoring by applying sidereal filtering (SF) in the observation domain. The equation that unifies the residuals for the uncombined and undifferenced (UCUD) PPP solution on different frequencies was derived, which could greatly reduce the complexity of data processing. An unanticipated long-term periodic error term of up to ± 3 cm was found in the phase residuals associated with BDS satellites in geostationary Earth orbit (GEO), which is not due to multipath originated from the ground but is in fact satellite dependent. The period of this error is mainly longer than 2000 s and cannot be alleviated by using multi-GNSS. Compared with solutions without sidereal filtering, the application of the SF approach dramatically improves the positioning precision with respect to the weekly averaged positioning solution, by 75.2%, 42.8%, and 56.7% to 2.00, 2.23, and 5.58 cm in the case of BDS-only PPP in the east, north, and up components, respectively, and 71.2%, 27.7%, and 37.9% to 1.25, 0.81, and 3.79 cm in the case of GPS/GLONASS/BDS combined PPP, respectively. The GPS/GLONASS/BDS combined solutions augmented by the SF successfully suppress the GNSS noise, which contributes to the detection of the true seismic signal and is beneficial to the pre- and post-seismic signal analysis

Modeling wide-area tropospheric delay corrections for fast PPP ambiguity resolution

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

A square root information filter for multi-GNSS real-time precise clock estimation

Real-time satellite orbit and clock estimations are the prerequisite for Global Navigation Satellite System (GNSS) real-time precise positioning services. To meet the high-rate update requirement of satellite clock corrections, the computational efficiency is a key factor and a challenge due to the rapid development of multi-GNSS constellations. The Square Root Information Filter (SRIF) is widely used in real-time GNSS data processing thanks to its high numerical stability and computational efficiency. In real-time clock estimation, the outlier detection and elimination are critical to guarantee the precision and stability of the product but could be time-consuming. In this study, we developed a new quality control procedure including the three standard steps: i.e., detection, identification, and adaption, for real-time data processing of huge GNSS networks. Effort is made to improve the computational efficiency by optimizing the algorithm to provide only the essential information required in the processing, so that it can be applied in real-time and high-rate estimation of satellite clocks. The processing procedure is implemented in the PANDA (Positioning and Navigation Data Analyst) software package and evaluated in the operational generation of real-time GNSS orbit and clock products. We demonstrated that the new algorithm can efficiently eliminate outliers, and a clock precision of 0.06 ns, 0.24 ns, 0.06 ns, and 0.11 ns can be achieved for the GPS, GLONASS, Galileo, and BDS-2 IGSO/MEO satellites, respectively. The computation time per epoch is about 2 to 3 s depending on the number of existing outliers. Overall, the algorithm can satisfy the IGS real-time clock estimation in terms of both the computational efficiency and product quality

Improving the vertical modeling of tropospheric delay

Accurate tropospheric delays from Numerical Weather Models (NWM) are an important input to space geodetic techniques, especially for precise real-time Global Navigation Satellite Systems, which are indispensable to earthquake and tsunami early warning systems as well as weather forecasting. The NWM-based tropospheric delays are currently provided either site-specific with a limited spatial coverage, or on two-dimensional grids close to the Earth surface, which cannot be used for high altitudes. We introduce a new method of representing NWM-derived tropospheric zenith hydrostatic and wet delays. A large volume of NWM-derived data is parameterized with surface values and additional two or three coefficients for their vertical scaling to heights up to 14 km. A precision of 1–2 mm is achieved for reconstructing delays to the NWM-determined delays at any altitudes. The method can efficiently deliver NWM-derived tropospheric delays to a broader community of space geodetic techniques.DFG, 434617780, SFB 1464: Relativistische und quanten-basierte Geodäsie (TerraQ

Validating HY-2A CMR precipitable water vapor using ground-based and shipborne GNSS observations

The calibration microwave radiometer (CMR) on board the Haiyang-2A (HY-2A) satellite provides wet tropospheric delay correction for altimetry data, which can also contribute to the understanding of climate system and weather processes. The ground-based global navigation satellite system (GNSS) provides precise precipitable water vapor (PWV) with high temporal resolution and could be used for calibration and monitoring of the CMR data, and shipborne GNSS provides accurate PWV over open oceans, which can be directly compared with uncontaminated CMR data. In this study, the HY-2A CMR water vapor product is validated using ground-based GNSS observations of 100 International GNSS Service (IGS) stations along the global coastline and 56 d shipborne GNSS observations over the Indian Ocean. The processing strategy for GNSS data and CMR data is discussed in detail. Special efforts were made in the quality control and reconstruction of contaminated CMR data. The validation result shows that HY-2A CMR PWV agrees well with ground-based GNSS PWV with 2.67 mm as the root mean square (rms) within 100 km. Geographically, the rms is 1.12 mm in the polar region and 2.78 mm elsewhere. The PWV agreement between HY-2A and shipborne GNSS shows a significant correlation with the distance between the ship and the satellite footprint, with an rms of 1.57 mm for the distance threshold of 100 km. Ground-based GNSS and shipborne GNSS agree with HY-2A CMR well

Observation of the Ωc0\Omega_{c}^{0} Charmed Baryon at CLEO

The CLEO experiment at the CESR collider has used 13.7 fb$^{-1}$ of data to search for the production of the $\Omega_c^0$ (css-ground state) in $e^{+}e^{-}$ collisions at $\sqrt{s} \simeq 10.6$ {\rm GeV}. The modes used to study the $\Omega_c^0$ are $\Omega^- \pi^+$, $\Omega^- \pi^+ \pi^0$, $\Xi^- K^- pi^+ \pi^+$, $\Xi^0 K^- pi^+$, and $\Omega^- \pi^+ \pi^- \pi^+$. We observe a signal of 40.4$\pm$9.0(stat) events at a mass of 2694.6$\pm$2.6(stat)$\pm$1.9(syst) {\rm MeV/$c^2$}, for all modes combined.Comment: 10 pages postscript, also available through http://w4.lns.cornell.edu/public/CLN

Observation of B→ϕKB\to \phi K and B→ϕK∗B\to \phi K^{*}

We have studied two-body charmless hadronic decays of $B$ mesons into the final states phi K and phi K^*. Using 9.7 million $B\bar{B}$ pairs collected with the CLEO II detector, we observe the decays B- -> phi K- and B0 -> phi K*0 with the following branching fractions: BR(B- -> phi K-)=(5.5 +2.1-1.8 +- 0.6) x 10^{-6} and BR(B0 -> phi K*0)=(11.5 +4.5-3.7 +1.8-1.7) x 10^{-6}. We also see evidence for the decays B0 -> phi K0 and B- -> phi K*-. However, since the statistical significance is not overwhelming for these modes we determine upper limits of <12.3 x 10^{-6} and <22.5 x 10^{-6} (90% C.L.) respectively.Comment: 9 pages postscript, also available through http://w4.lns.cornell.edu/public/CLN