An experimental combination of IGS repro3 campaign’s orbit products using a variance component estimation strategy

Over the past years, the International GNSS Service (IGS) has put efforts into reprocessing campaigns reanalyzing the full data collected by the IGS network since 1994. The goal is to provide a consistent set of orbits, station coordinates, and earth rotation parameters using state-of-the-art models. Different from the previous campaigns - namely: repro1 and repro2 - the repro3 includes not only GPS and GLONASS but also the Galileo constellation. The main repro3 objective is the contribution to the next realization of the International Terrestrial Reference Frame (ITRF2020). To achieve this goal, several Analysis Centers (AC) submitted their specific products, which are combined to provide the final solutions for each product type. In this contribution, we focus on the combination of the orbit products.We will present a consistent orbit solution based on a newly developed combination strategy where the weights are determined by a Least-Squares Variance Component Estimation (LSVCE). The orbits are combined in an iterative processing, first aligning all the products via a Helmert transformation, second defining which satellites will be used in the LSVCE, and finally normalizing the inverse of the variances as weights that are used to compute a weighted mean. Moreover, we will discuss the weight factors and their stability in the time evolution for each AC depending on the constellations. In addition, an external validation using a Satellite Laser Ranging (SLR) procedure will be shown for the combined solution.

An experimental combination of IGS repro3 campaign’s orbit products using a variance component estimation strategy

Over the past years, the International GNSS Service (IGS) has put efforts into reprocessing campaigns reanalyzing the full data collected by the IGS network since 1994. The goal is to provide a consistent set of orbits, station coordinates, and earth rotation parameters using state-of-the-art models. Different from the previous campaigns - namely: repro1 and repro2 - the repro3 includes not only GPS and GLONASS but also the Galileo constellation. The main repro3 objective is the contribution to the next realization of the International Terrestrial Reference Frame (ITRF2020). To achieve this goal, several Analysis Centers (AC) submitted their specific products, which are combined to provide the final solutions for each product type. In this contribution, we focus on the combination of the orbit products.We will present a consistent orbit solution based on a newly developed combination strategy where the weights are determined by a Least-Squares Variance Component Estimation (LSVCE). The orbits are combined in an iterative processing, first aligning all the products via a Helmert transformation, second defining which satellites will be used in the LSVCE, and finally normalizing the inverse of the variances as weights that are used to compute a weighted mean. Moreover, we will discuss the weight factors and their stability in the time evolution for each AC depending on the constellations. In addition, an external validation using a Satellite Laser Ranging (SLR) procedure will be shown for the combined solution.

Assessment and reduction of zenith path delay biases due to day boundary effect

The troposphere consists of dry air and water vapor, delaying the GNSS signal by about 2.4 m in the zenith direction. The water vapor only causes an error of about 0.2 m in distance measurement, but it is challenging to model and overcome. From 2003 the International GNSS Service (IGS) started to provide the new product of zenith path delay (ZPD) with an accuracy of 1.5-5 mm. However, we found an error in these products up to 30 mm at epochs between 2 days due to the day boundary effect and an average of 16mm RMS for nine days. Our research shows that for reducing the impact, the most critical factor is selecting the initial value for the ZPD, followed by satellite orbit/clock and, finally, the station coordinate values. By choosing an appropriate initial value for ZPD and employing a 3 days orbit/clock, the ZPD error due to the day boundary effect can be reduced to negligible. Meanwhile, the change in the station coordinate value in cm level does not impact the effect.

The cooperative IGS RT-GIMs: a reliable estimation of the global ionospheric electron content distribution in real time

The Real-Time Working Group (RTWG) of the International GNSS Service (IGS) is dedicated to providing high-quality data and high-accuracy products for Global Navigation Satellite System (GNSS) positioning, navigation, timing and Earth observations. As one part of real-time products, the IGS combined Real-Time Global Ionosphere Map (RT-GIM) has been generated by the real-time weighting of the RT-GIMs from IGS real-time ionosphere centers including the Chinese Academy of Sciences (CAS), Centre National d'Etudes Spatiales (CNES), Universitat Politècnica de Catalunya (UPC) and Wuhan University (WHU). The performance of global vertical total electron content (VTEC) representation in all of the RT-GIMs has been assessed by VTEC from Jason-3 altimeter for 3 months over oceans and dSTEC-GPS technique with 2 d observations over continental regions. According to the Jason-3 VTEC and dSTEC-GPS assessment, the real-time weighting technique is sensitive to the accuracy of RT-GIMs. Compared with the performance of post-processed rapid global ionosphere maps (GIMs) and IGS combined final GIM (igsg) during the testing period, the accuracy of UPC RT-GIM (after the improvement of the interpolation technique) and IGS combined RT-GIM (IRTG) is equivalent to the rapid GIMs and reaches around 2.7 and 3.0 TECU (TEC unit, 1016 el m−2) over oceans and continental regions, respectively. The accuracy of CAS RT-GIM and CNES RT-GIM is slightly worse than the rapid GIMs, while WHU RT-GIM requires a further upgrade to obtain similar performance. In addition, a strong response to the recent geomagnetic storms has been found in the global electron content (GEC) of IGS RT-GIMs (especially UPC RT-GIM and IGS combined RT-GIM). The IGS RT-GIMs turn out to be reliable sources of real-time global VTEC information and have great potential for real-time applications including range error correction for transionospheric radio signals, the monitoring of space weather, and detection of natural hazards on a global scale. All the IGS combined RT-GIMs generated and analyzed during the testing period are available at https://doi.org/10.5281/zenodo.5042622 (Liu et al., 2021b).

Feasibility analysis of GNSS-based ionospheric observation on a fast-moving train platform (GIFT)

The ionospheric effect plays a crucial role in the radio communications. For ionospheric observing and monitoring, the Global Navigation Satellite System (GNSS) has been widely utilized. The ionospheric condition can be characterized by the Total Electron Contents (TEC) and TEC Rate (TECR) calculated from the GNSS measurements. Currently, GNSS-based ionospheric observing and monitoring largely depend on a global fiducial network of GNSS receivers such as the International GNSS Service (IGS) network. We propose a new approach to observe the ionosphere by deploying a GNSS receiver on a Hong Kong Mass Transit Railway (MTR) train. We assessed the TECR derived from the MTR-based GNSS receiver by comparing it with the TECR derived from a static GNSS receiver. The results show that the Root-Mean-Squares (RMS) errors of the TECR derived from the MTR-based GNSS receiver is consistently approximately 23% higher than that derived from the static GNSS receiver. Despite the increased error, the findings suggest that the GNSS observation on a fast-moving platform is a feasible approach to observe the ionosphere over a large region in a rapid and cost-effective way.

Status of CAS global ionospheric maps after the maximum of solar cycle 24

As a new Ionosphere Associate Analysis Center (IAAC) of the International GNSS Service (IGS), Chinese Academy of Sciences (CAS) started the routine computation of the real-time, rapid, and final Global Ionospheric Maps (GIMs) in 2015. The method for the generation of CAS rapid and final GIMs and recent updates are presented in the paper. The quality of CAS post-processed GIMs is assessed during 2015–2018 after the maximum of solar cycle 24. To perform an independent and fair assessment, Jason-2/3 Vertical Total Electron Contents (VTEC) are first used as the references over the ocean. GPS differential Slant TECs (dSTEC) generated from 55 Multi-GNSS Experimental (MGEX) stations of the IGS are also employed, which provides a complementing way to evaluate the ability of electron content models to reproduce the spatial and temporal gradients in the ionosphere. During the test period, Jet Propulsion Laboratory (JPL) GIMs present significantly positive deviations compared to the Jason VTEC and GPS dSTEC. Technical University of Catalonia (UPC) rapid GIM UQRG exhibits the best performance in both Jason VTEC and GPS dSTEC analysis. The CAS GIMs show comparable performance with the results of the first four IAACs of the IGS. As expected, the poor performance of all GIMs is in equatorial regions and the high latitudes of the southern hemisphere. The consideration of generating multi-layer or three-dimensional ionospheric maps is emphasized to mitigate the inadequacy of ionospheric single-layer assumption in the presence of pronounced latitudinal gradients. The use of ionospheric observations from the new GNSS constellations and other space- or ground-based observation techniques is also suggested in the generation of future GIMs, given the sparse GPS/GLONASS stations in the southern hemisphere.

A Comparative Study on the Solar Radiation Pressure Modeling in GPS Precise Orbit Determination

For Global Positioning System (GPS) precise orbit determination (POD), the solar radiation pressure (SRP) is the dominant nongravitational perturbation force. Among the current SRP models, the ECOM and box-wing models are widely used in the International GNSS Service (IGS) community. However, the performance of different models varies over different GPS satellites. In this study, we investigate the performances of different SRP models, including the box-wing and adjustable box-wing as a priori models, and ECOM1 and ECOM2 as parameterization models, in the GPS POD solution from 2017 to 2019. Moreover, we pay special attention to the handling of the shadow factor in the SRP modeling for eclipsing satellites, which is critical to achieve high-precision POD solutions but has not yet been fully investigated. We demonstrate that, as an a priori SRP model, the adjustable box-wing has better performance than the box-wing model by up to 5 mm in the orbit day boundary discontinuity (DBD) statistics, with the largest improvement observed on the BLOCK IIR satellites using the ECOM1 as a parameterization SRP model. The box-wing model shows an insignificant orbit improvement serving as the a priori SRP model. For the eclipsing satellites, the three-dimensional (3D) root mean square (RMS) values of orbit DBD are improved when the shadow factor is applied only in the D direction (pointing toward to Sun) than that in the three directions (D, Y, and B) in the satellite frame. Different SRP models have comparable performance in terms of the Earth rotation parameter (ERP) agreement with the IERS EOP 14C04 product, whereas the magnitude of the length of day (LoD) annual signal is reduced when the shadow factor is applied in the D direction than in the three directions. This study clarifies how the shadow factor should be applied in the GPS POD solution and demonstrates that the a priori adjustable box-wing model combined with ECOM1 is more suitable for high-precision GPS POD solutions, which is useful for the further GNSS data analysis.

IONORING: Real-Time Monitoring of the Total Electron Content over Italy

IONORING (IONOspheric RING) is a tool capable to provide the real-time monitoring and modeling of the ionospheric Total Electron Content (TEC) over Italy, in the latitudinal and longitudinal ranges of 35°N-48°N and 5°E-20°E, respectively. IONORING exploits the Global Navigation Satellite System (GNSS) data acquired by the RING (Rete Integrata Nazionale GNSS) network, managed by the Istituto Nazionale di Geofisica e Vulcanologia (INGV). The system provides TEC real-time maps with a very fine spatial resolution (0.1° latitude x 0.1° longitude), with a refresh time of 10 min and a typical latency below the minute. The TEC estimated at the ionospheric piercing points from about 40 RING stations, equally distributed over the Italian territory, are interpolated using locally (weighted) regression scatter plot smoothing (LOWESS). The validation is performed by comparing the IONORING TEC maps (in real-time) with independent products: i) the Global Ionospheric Maps (GIM) - final product- provided by the International GNSS Service (IGS), and ii) the European TEC maps from the Royal Observatory of Belgium. The validation results are satisfactory in terms of Root Mean Square Error (RMSE) between 2 and 3 TECu for both comparisons. The potential of IONORING in depicting the TEC daily and seasonal variations is analyzed over 3 years, from May 2017 to April 2020, as well as its capability to account for the effect of the disturbed geospace on the ionosphere at mid-latitudes. The IONORING response to the X9.3 flare event of September 2017 highlights a sudden TEC increase over Italy of about 20%, with a small, expected dependence on the latitude, i.e., on the distance from the subsolar point. Subsequent large regional TEC various were observed in response to related follow-on geomagnetic storms. This storm is also used as a case event to demonstrate the potential of IONORING in improving the accuracy of the GNSS Single Point Positioning. By processing data in kinematic mode and by using the Klobuchar as the model to provide the ionospheric correction, the resulting Horizontal Positioning Error is 4.3 m, lowering to, 3.84 m when GIM maps are used. If IONORING maps are used as the reference ionosphere, the error is as low as 2.5 m. Real-times application and services in which IONORING is currently integrated are also described in the conclusive remarks.

Estimation and Analysis of the Observable-Specific Code Biases Estimated Using Multi-GNSS Observations and Global Ionospheric Maps

Observable-specific bias (OSB) parameterization allows observation biases belonging to various signal types to be flexibly addressed in the estimation of ionosphere and global navigation satellite system (GNSS) clock products. In this contribution, multi-GNSS OSBs are generated by two different methods. With regard to the first method, geometry-free (GF) linear combinations of the pseudorange and carrier-phase observations of a global multi-GNSS receiver network are formed for the extraction of OSB observables, and global ionospheric maps (GIMs) are employed to correct ionospheric path delays. Concerning the second method, satellite and receiver OSBs are converted directly from external differential code bias (DCB) products. Two assumptions are employed in the two methods to distinguish satellite- and receiver-specific OSB parameters. The first assumption is a zero-mean condition for each satellite OSB type and GNSS signal. The second assumption involves ionosphere-free (IF) linear combination signal constraints for satellites and receivers between two signals, which are compatible with the International GNSS Service (IGS) clock product. Agreement between the multi-GNSS satellite OSBs estimated by the two methods and those from the Chinese Academy of Sciences (CAS) is shown at levels of 0.15 ns and 0.1 ns, respectively. The results from observations spanning 6 months show that the multi-GNSS OSB estimates for signals in the same frequency bands may have very similar code bias characteristics, and the receiver OSB estimates present larger standard deviations (STDs) than the satellite OSB estimates. Additionally, the variations in the receiver OSB estimates are shown to be related to the types of receivers and antennas and the firmware version. The results also indicate that the root mean square (RMS) of the differences between the OSBs estimated based on the CAS- and German Aerospace Center (DLR)-provided DCB products are 0.32 ns for the global positioning system (GPS), 0.45 ns for the BeiDou navigation satellite system (BDS), 0.39 ns for GLONASS and 0.22 ns for Galileo.

Impact of GPS/BDS Satellite Attitude Quaternions on Precise Point Positioning with Ambiguity Resolution

Precise point positioning with ambiguity resolution (PPP-AR) based on multiple global navigation satellite system (multi-GNSS) constellations is an important high-precision positioning tool. However, some unmodeled satellite and receiver biases (such as errors in satellite attitude) make it difficult to fix carrier-phase ambiguities. In order to fix ambiguities of eclipsing satellites, accurate integer clock and satellite attitude products (i.e., attitude quaternion) have been provided by the International GNSS Service (IGS). Nevertheless, the quality of these products and their positioning performance in multi-GNSS PPP-AR have not been investigated yet. Using the PRIDE PPP-AR II software associated with the corresponding rapid satellite orbit, integer clock and attitude quaternion products of Wuhan University (WUM), we carried out GPS/BDS PPP-AR using 30 days of data in an eclipsing season of 2020. We found that about 75% of GPS, 60% of BDS-2 and 57% of BDS-3 narrow-lane ambiguity residuals after integer clock corrections fall within ±0.1 cycles in the case of using nominal attitudes. However, when using attitude quaternions, these percentages will rise to 80% for GPS, 70% for BDS-2 and 60% for BDS-3. GPS/BDS daily kinematic PPP-AR after integer clock and nominal attitude corrections can usually achieve a positioning precision of about 10, 10 and 30 mm for the east, north and up components, respectively. In contrast, the counterparts are 8, 8 and 20 mm when using attitude quaternions. Compared with the case of using attitude quaternions only at the network end for the integer clock estimation, using attitude quaternions only at the user end shows a pronounced improvement of 15% in the east component and less than 10% in the north and up components. Therefore, we suggest PPP users apply integer clock and satellite attitude quaternion products to realize more efficient ambiguity fixing, especially in satellite eclipsing seasons.