Status of IGS Reprocessing Activities at GFZ

Based on a large network of continuously operated GNSS tracking stations the International GNSS Service (IGS) has a valuable contribution for the realization of the International Terrestrial Reference System (ITRS). In order to contribute to its next realization, the IGS is preparing for a new reprocessing of the GNSS data from 1994 to 2020 including GPS, GLONASS, and – for the ﬁrst time – Galileo. A ﬁrst test campaign including single- and multi-system solutions for 2017 and 2018 was performed to derive consistent transmitter phase center corrections for all systems. Preliminary results of the test solutions derived at GFZ show well determined orbits with overlaps of 28 mm for GPS, 67mm for GLONASS, and 40 mm for Galileo and an overall RMS of satellite laser ranging residuals for Galileo of 58 mm. Using multi-GNSS antenna calibrations (including also E5a and E5b calibrations) horizontal coordinate differences are almost zero between a GPS+GLONASS and a Galileo-only solutions. Due to the mixture of estimated (GPS, GLONASS) and measured (Galileo) transmitter phase center offsets a scale difference of 1:16 ˙ 0:27 ppb is found between both solutions which agrees nicely to results derived by other analysis centers.


Introduction
To provide the best possible GNSS solution for the realization of the International Terrestrial Reference System, the Analysis Centers (ACs) of the International GNSS Service (IGS, Johnston et al. 2017) are preparing for a full reprocessing of GNSS data from 1994 to 2020. Like the previous efforts (repro1 and repro2) the upcoming reprocessing will provide a fully consistent set of orbits, station coordinates and Earth rotation parameters derived with the best and most consistent models available. It is well B. Männel ( ) · A. Brandt · M. Bradke · P. Sakic · A. Brack · T. Nischan Helmholtz Centre Potsdam -German Research Centre for Geosciences, Potsdam, Germany e-mail: benjamin.maennel@gfz-potsdam.de known that in terms of reference frame parameters the most critical issues for GNSS are, firstly, the transmitter phase center offsets (which are highly correlated with the terrestrial scale, e.g., Zhu et al. 2003) and, secondly, the modeling of the solar radiation pressure on the orbits (main reason for draconitic period in geocenter results, e.g., Meindl et al. 2013). While trying to reduce or to solve both issues several additional topics have to be considered like the 13.63/13.66 day signal in GNSS time series (see for example Ray et al. 2013) or remaining modeling inconsistencies compared to other space geodetic techniques. Compared to the last reprocessing, new satellite systems like Galileo and BeiDou became almost fully operational. As their signals were tracked by an increasing number of IGS stations during the past years, the set of considered systems has to be redefined from GPS+GLONASS in repro2 to an up-to-date multi-GNSS solution. During the IGS Analysis Workshop 2019 held in Potsdam, Germany, the IGS ACs agreed to strive for an combined GPS (G), GLONASS (R), and Galileo (E) solution in the upcoming repro3. However, to avoid systematic distortions, so far missing, receiver antenna corrections for the Galileo signals E5a, E5b, and E6 and the GPS L5 frequency as well as consistent transmitter phase center offsets (PCOs) are required (see e.g. Schmid et al. 2016). Whereas the first issue was solved for many antenna types used in the IGS as Geo++ provided robot-based calibrations for these signals it was agreed to solve the second issue by setting up a test campaign. This campaign includes multi-and singlesystem solutions (if possible GRE, GR, G, R, and E) for 2017 and 2018 which will be used to derive phase center offsets for GPS and GLONASS based on the Galileo offsets which are known thanks to published chamber calibrations (GSA 2017). As these Galileo PCOs are measured -and not estimated from observations itself -they are independent of the terrestrial scale which has to be fixed to the ITRF scale otherwise (Schmid et al. 2007). Therefore, an independent GNSS scale will become available if the final repro3 could be performed with this consistent set of re-estimated and calibrated PCOs. It was also agreed to run a second test campaign using the final repro3 setup including the station selection as benchmark test before starting the processing tasks.
This paper summarizes the current reprocessing status at GFZ during the first test campaign and highlights preliminary outcomes. Section 2 describes investigations regarding station selection and testing some models. Initial results based on the different GFZ solutions in the first test campaign will be presented in Sect. 3. Section 4 provides an outlook to the upcoming tasks.

Data and Processing
This section discusses the processed data, the station network, and the selected models for the test campaign but also for the final reprocessing.

Data Selection
The station selection process is based on the pre-selection and station classification which was provided to the Analysis Centres by Paul Rebischung via the IGS AC mailing list (IGS-ACS-1235). 1 According to software and time capabilities we will process stations listed in Categories 1 (revised set of IGS14 core stations), 2 (stations with local ties to other techniques), 3 (redundant local tie stations), and 4 (remaining IGS14 stations) as well as IGS stations operated by GFZ placed in lower categories. In order to assess this selection we re-imported the whole data set into our archive with dedicated checks for formal correctness and consistency with the provided site logs. In addition, we processed the GPS observations from all selected IGS stations using the EPOS.P8 software in PPP mode to identify the stations temporal behavior. 2 The processing was done based on orbit and clock products derived within a GFZ internal reprocessing effort which we carried out in 2018 to derive consistent products in the IGS14 frame. Figure 1 shows the derived height coordinate repeatability for each station with the symbol identifying the station's category. Overall, 325 stations are contained while the average repeatability is 7:5˙1:9 mm (median is 7.2 mm). The GFZ contribution to the first test campaign was, however, processed using the station selection used in GFZ's operational processing (IGS final line, 220 stations) supported by additional 30 IGS stations selected to achieve a basic coverage for Galileo also in early 2017. However, the number of selected Galileo sites was still rather low for an independent Galileo-only solution. For January 1st, 2017 the number of stations was 222/137/68 for GPS, GLONASS, and Galileo while rising to 210/145/124 for Dec. 31st, 2018 due to ongoing station upgrades within the IGS tracking network. With respect to the number of Galileo satellites it might be interesting to know that the number of satellites increased from 13 to 22 within the same time span. 3 Table 1 provides the processing strategy applied for the test campaign. The same parametrization can be expected for the final reprocessing. However, the orbit parametrization might be modified for some satellite blocks as tests regarding the optimal setup are performed currently. In general, the settings follow the IERS Conventions 2010 (Petit and Luzum 2010) and the repro2 standards. Using the EPOS.P8 software GFZ provided four solutions (GRE, GR, G, and E).

Updating Models for Repro 3
According to the discussions between the IGS ACs and the IERS, several models have to be updated for the final reprocessing compared to Table 1. For computing the rotational deformation (pole tide) the linear mean pole will be used as adopted by the IERS in 2018. Regarding the gravity field, a static gravity field up to degree and order 12 was used whereas a time-variable gravity field should be used in the reprocessing. The ocean tides and ocean loading model will be updated to a more recent FES2014b model (Carrere et al. 2016). In order to consider high-frequency variations in Earth orientation parameters (EOP) it was agreed to use the model provided by Desai and Sibois (2016) instead of the model provided in the IERS Conventions.

InitialResults
Within this section initial results derived within the test campaign will be discussed with respect to the derived orbits and stations coordinates. Figure 2 shows the orbit misclosures (orbit overlap error) for all satellites processed in   the full solution (GRE). The overlaps are computed for 2 h each while estimating transformation parameters between the two orbit solutions. An averaged RMS of 28 mm is achieved for the GPS satellites. As shown in Table 2 the RMS is independent of the solution type (GRE, GR, or G) for GPS with differences of 1 mm. It is obvious, that the GPS orbits are not downgraded in terms of overlap errors by adding other systems to the solutions. For GLONASS a larger mean RMS of 67 mm is observed while five satellites exceed 90 mm overlap error (see Table 3). As shown in this table, Dach et al. (2019) reported similarly large overlaps for two of the three satellites (R715 and R719 are not contained in Table 1 of Dach et al. (2019)). Without these satellites the remaining average is 57.8 mm. Galileo orbit overlaps are in general larger compared to GPS mis-closures. Without large variations between the satellites an averaged RMS of 40 mm is determined for Galileo which is also achieved by the Galileo-only solution.
In order to further asses the satellite orbits a validation based on satellite laser ranging (SLR) was performed. While fixing the SLR telescope positions to their ITRF2014 coordi-nates (Altamimi et al. 2016) and estimating no other parameters, the derived residuals (i.e., the differences between observed and computed ranges) provide insights into the absolute orbit accuracy. The number of SLR observations varies between 350 normal points (E217 and E218) and more than 17,000 (R802, R853). On average 5,230 normal points are collected per satellite within the 2 years. During the processing 3.6 % and 2.1 % of the observations where excluded as outliers for GLONASS and Galileo, respectively. Figure 3 shows the derived statistics. The determined biases, i.e., mean values over all residuals, reach up to 25 mm for GLONASS with some large variations between the satellites. For R856 with only 742 observations a larger bias of 44 mm is determined. In general, remaining biases indicate systematics contained in the orbits (or applied sensor offsets). In the current solution positive biases are visible for most of the GLONASS satellites which needs further investigations. The biases for Galileo are small (few millimeters) but also mostly positive for the whole constellation. The larger bias of 13 mm for E218 might be related to the low number of 350 SLR observations available for this satellite launched in December 2017. The standard deviation of all SLR residuals reaches 91:5˙13:3 mm for GLONASS and 58:3˙8:7 mm for Galileo.
With respect to the different solutions a comparison of the station coordinates is very important. In theory, the estimated station coordinate should be independent of the processed GNSS. However, it was shown for example by Villiger et al. (2019) that one has to expect significant coordinate differences between system-specific solutions. These differences are mostly related to the considered antenna corrections for transmitter and receiver. As stated earlier, transmitter PCOs are either estimated in a global adjustment or -as for Galileo -chamber calibrated. In addition, robot-calibrations of receiving antennas were not available for several GNSS signals so far. Figure 4 shows the mean differences (and standard deviation) in North, East, and height coordinates derived in a combined (GR) and a Galileo-only solution using the provided multi-GNSS antenna corrections. It has to

Summary and Outlook
The presented repro3 test campaign was performed at GFZ as contribution to the re-determination of transmitter phase center offsets for GPS and GLONASS in preparation for the final reprocessing. The current setup for this reprocessing is a three-system solution (GPS, GLONASS, Galileo) while some models like the time-variable gravity field are still in discussion between the IGS ACs and the IERS. The preliminary results show acceptable overlap errors for the individual satellites with on average 28 mm, 67 mm and 40 mm for GPS, GLONASS, and Galileo. An SLR orbit validation revealed also good orbit quality without significant biases for the Galileo satellites. Comparing the derived station coordinates a good agreement between the GR-and the Galileo-only solution can be found in the horizontal components with a scale difference of around 1.1 ppb. According to the reprocessing schedule, a second test campaign will be performed after final decisions on the models. In addition, a few open questions have to be addressed like the reason for large orbit overlap errors for some GLONASS satellites or the coordinate difference for some stations equipped with JAVRINGANT_G3T and JAVRINGANT_G5T antennas.
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons. org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence and indicate if changes were made. The images or other third party material in this chapter are included in the chapter's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.