On the application of the raw-observation-based PPP to global ionosphere VTEC modeling: an advantage demonstration in the multi-frequency and multi-GNSS context

The ionospheric delay accounts for one of the major errors that the Global Navigation Satellite Systems (GNSS) suffer from. Hence, the global ionosphere Vertical Total Electron Content (VTEC) map has been an important atmospheric product within the International GNSS Service (IGS) since its early establishment. In this contribution, an enhanced method has been proposed for the modeling of the global VTECs, in which the enhancements include two aspects. Firstly, to cope with the rapid development of the newly established Galileo and BeiDou constellations in recent years, we extend the current dual-system (GPS/GLONASS) solution to a quad-system (GPS/GLONASS/Galileo/BeiDou) solution. More importantly, instead of using dual-frequency observations based on the Carrier-to-Code Leveling method, all available triple-frequency signals are utilized with a general raw-observation-based multi-frequency Precise Point Positioning model, which can process dual-, triple- or even arbitrary-frequency observations compatibly and flexibly. Benefiting from this, quad-system slant ionospheric delays can be retrieved based on multi-frequency observations in a more flexible, accurate and reliable way, which are finally used to establish global VTEC models with the spherical harmonic function. In this process, multi-GNSS Differential Code Biases (DCBs) are also estimated as by-products. More than 400 globally distributed stations from the IGS and the Multi-GNSS Experiment (MGEX) networks have been processed in both 2014 (with high solar activity) and 2018 (with low solar activity). Global VTECs have been compared with the IGS final products, and the over-ocean VTECs are validated with the results from the JASON altimeter. The mean RMS values of the VTEC differences are 1.84 (2014) and 1.23 (2018) TECUs with respect to the IGS final products. The standard deviations of the VTEC differences with respect to the JASON results are 4.71 (2014) and 2.82 (2018) TECUs, outperforming all the other products generated with the spherical harmonic function. Additionally, multi-GNSS satellite DCBs have also been validated with the existing products from the Center for Orbit Determination in Europe and MGEX. All the results prove that the proposed method can be used as an effective and accurate approach for global VTEC modeling and DCB estimation, especially in the future multi-frequency and multi-GNSS context.