Combined orbits and clocks from IGS second reprocessing

The Analysis Centers (ACs) of the International GNSS Service (IGS) have reprocessed a large global network of GPS tracking data from 1994.0 until 2014.0 or later. Each AC product time series was extended uniformly till early 2015 using their weekly operational IGS contributions so that the complete combined product set covers GPS weeks 730 through 1831. Three ACs also included GLONASS data from as early as 2002 but that was insufficient to permit combined GLONASS products. The reprocessed terrestrial frame combination procedures and results have been reported already, and those were incorporated into the ITRF2014 multi-technique global frame released in 2016. This paper describes the orbit and clock submissions and their multi-AC combinations and assessments. These were released to users in early 2017 in time for the adoption of IGS14 for generating the operational IGS products. While the reprocessing goal was to enable homogeneous modeling, consistent with the current operational procedures, to be applied retrospectively to the full history of observation data in order to achieve a more suitable reference for geophysical studies, that objective has only been partially achieved. Ongoing AC analysis changes and a lack of full participation limit the consistency and precision of the finished IG2 products. Quantitative internal measures indicate that the reprocessed orbits are somewhat less precise than current operational orbits or even the later orbits from the first IGS reprocessing campaign. That is even more apparent for the clocks where a lack of robust AC participation means that it was only possible to form combined 5-min clocks but not the 30-s satellite clocks published operationally. Therefore, retrospective precise point positioning solutions by users are not recommended using the orbits and clocks. Nevertheless, the orbits do support long-term stable user solutions when used with network processing with either double differencing or explicit clock estimation. Among the main benefits of the reprocessing effort is a more consistent long product set to analyze for sources of systematic error and accuracy. Work to do that is underway but the reprocessing experience already points to a number of ways future IGS performance and reprocessing campaigns can be improved.


Introduction
In early 2015, the International GNSS Service (IGS; Dow et al. 2009) Analysis Centers (ACs) completed a second reanalysis of Global Navigation Satellite System (GNSS) data collected for a global network of tracking stations. This second reprocessing, or repro2, updates the set of definitive IGS combination data products-station positions; satellite orbits and clocks; and Earth orientation parameters (EOPs)using the latest analysis models and methodologies. The data used for repro2 spanned ∼ 21 years, starting January 2, 1994, and continuing through February 14, 2015, or GPS Weeks B Jake Griffiths jake.griffiths@nrl.navy.mil 1 Naval Center for Space Technology, U.S. Naval Research Laboratory, Washington, DC, USA 730 through 1831. Participating analysis groups used data for the US Global Positioning System (GPS) and also the Russian GLONASS system in the case of three ACs. This paper focuses on results from the repro2 orbit and clock combinations, resulting in reprocessed IGS combined (IG2) orbit and clock products for GPS only. Please refer to Rebischung et al. (2016) for a description of the tracking network used and the associated reprocessed station positions and EOPs.

Historical context
To date, two IGS reprocessing campaigns have occurred outside but in parallel with the IGS operational product streams (ultra-rapid, rapid, final) and were generally undertaken with the goal to advance Earth science research through updates to the International Terrestrial Reference Frame (ITRF).
Repro2 provided the IGS input  to ITRF2014 . It followed the successful first reprocessing campaign, or repro1, which provided the IGS input for ITRF2008 (Altamimi et al. 2011). The resulting reprocessed combination orbit and clock products plus the follow-on Finals operational products aim to disseminate the latest realization of the ITRF without a loss of fidelity. Therefore, the IGS reprocessing campaigns nominally aim to homogenize the full history of IGS combination data products in an internally consistent way, with ACs adopting the latest analysis models and methodologies available at the time of the campaign. In that approach, a posteriori quality assessments of the resulting reprocessed combination products can serve to inform future reprocessing campaigns and help to advance the state-of-the-art while providing quantitative measures of product accuracy. This is the ideal scenario; the realities are sometimes quite different. Often the ACs disagree on a common set of analysis models to be used or simply fail to apply timely changes, and that can lead to reduced internal precision and complicate interpretations of the combination results. This will become evident as repro2 orbit and clock combination results are discussed in the subsequent sections.
The idea that IGS reprocessing campaigns aim to homogenize the full history of IGS combination data products presumes that somehow heterogeneities appear in the established set of products at some point. It is not hard to see how those heterogeneities can get introduced. The evolution of the original IGS Finals operational data products (IGF) prior to repro1 is a good example. After the operational product streams were established in the early days of the IGS, the ACs naturally began adopting new analysis models and methods incrementally for their operational products. Mostly these newly adopted models led to improved AC products and therefore improved the IGS combination products going forward. But as analysis improvements compounded, the oldest products were less useful because they were less precise and inaccurate. Eventually, use of the earlier products in Earth science research began to limit studies of long-term geophysical processes. This was exacerbated in particular by periodic adoption of updated global reference frames and changes in calibration tables for ground and satellite antennas. It became clear that reprocessing the full history of data offered the potential to remove past heterogeneities, improve the IGS contribution to the ITRF and potentially advance Earth science research. For repro1, a set of common standard models and methods was established, and the operational ACs focused for several years on implementing them into their software. A full summary of those repro1 analysis standards is available online (http:// acc.igs.org/reprocess.html), but the main changes were due to: • a switch to absolute calibrations (Gendt, IGSMAIL-5272;Schmid et al. 2007) for receiver antennas and GLONASS and GPS satellite transmitter antennas; • adoption of the IGS05/igs05.atx framework (Ferland 2006, IGSMAIL-5447;Ferland and Piraszewski 2009), which is aligned to ITRF2005 (Altamimi et al. 2007); and • general implementation of the International Earth Rotation and Reference Systems Service (IERS) 2003 Conventions (McCarthy et al. 2003), including updated models for tropospheric propagation delay and for station displacements due to ocean tidal loading with whole-Earth center-of-mass corrections applied to SP3 orbits. Figure 1 shows WRMS time series statistics (smoothed) for each contributing AC in the repro1 orbit combinations. The trace of the IGF curve shows the evolution of the WRMS difference compared to IG1 and demonstrates historical improvements in the internal precision of the historical IGS operational products. The step decrease at GPS Week 1400 (November 05, 2006) is due to the switch to absolute antenna calibrations when the IGS05 reference frame was adopted for the operational product streams (Gendt, IGSMAIL-5438;Schmid et al. 2009). Other discrete analysis changes happened at other times and are less obvious in the WRMS statistics. For instance, the many IGS reference frames adopted over the years did not obviously impact the IGF WRMS statistics despite the fact that adopting a new frame is a discrete event. However, dramatic impacts from several of the frame changes are clearly seen in the Helmert rotation and translation parameters of some ACs (e.g., Fig. 2a,b) estimated in the repro1 orbit combination (Gendt et al. 2010).

Known remaining errors in IGS products
As Figs. 1 and 2 illustrate, repro1 marked an overall dramatic improvement in the orbit precision and long-term frame stability compared to the original IGS Finals operational products. As will be discussed in the following sections, the incremental improvements made for repro2 following repro1 are less dramatic because the overall errors in the IG1 products are already much smaller than they were in the original operational products. However, while the overall accuracy and precision of the IGS definitive products reached a remarkable level with repro1, significant deficiencies continue to affect the products at the few centimeters level and smaller (see Ray 2016 for a recent review). The largest of the errors include a combination of: unattributed subseasonal errors ; effects of background power-law noise in station coordinates on station velocities (Zhang et al. 1997;Santamaría-Gómez et al. 2011); effects of positional discontinuities on ground sta- orbits compared to IG1 orbits from repro1. The steady and dramatic decline in the IGF curve represents the accumulated effect of errors being removed in the operational products through the series of historical analysis changes. The IGF orbits were included for comparison only; they did not contribute to IG1. Black vertical lines along the top of the chart mark indicate when new terrestrial frames were adopted for the operational products tion velocities and frame stability (Williams 2003;Griffiths and Ray 2016); effects on EOPs and station time series residuals due to terrestrial frame misalignments to a longterm reference (Ray et al. 2017); subdaily EOP alias and draconitic errors in the satellite orbits Ray 2009, 2013) and those plus other harmonics in time series of ground station positions (Ray et al. 2008;Rebischung et al. 2016), which have been attributed to a combination of local near-field multipath, mismodeling of solar radiation pressure (SRP) and possibly other orbit-related errors; and various annual signals including unmodeled station displacements due to surface pressure loading, temporal changes in the low-degree geopotential coefficients and thermal expansion and flexure of ground antenna structures, among others. Some of these errors, like the background power-law noise in station coordinates, may be intrinsic to GPS. Others, like draconitic and tidal alias errors, can possibly be mitigated with continued efforts to update analysis models. Still other errors arise from inattentive operational and managerial diligence (e.g., excessive positional offsets caused by equipment changes, poor siting of tracking antennas and missing metadata needed to compute thermal variations).
In any case, most of these significant errors are expected to persist in the repro2 combination products and probably beyond.

IGS repro2 analysis models
A full list of analysis standards used for repro2 is summarized online (http://acc.igs.org/reprocess2.html), but the main changes from repro1 are: • a switch from weekly to daily terrestrial frame integrations to facilitate the study of station displacements at higher temporal resolution (Griffiths 2012, IGSMAIL-6613;Griffiths and Choi 2013;Rebischung et al. 2013); • inclusion of GLONASS data by some (three) but not all ACs; • the implementation of the IGb08/igs08.atx reference frame and calibration framework (Rebischung et al. 2012, IGSMAIL-6354, IGSMAIL-6663); • general implementation of the IERS 2010 Conventions (Petit et al. 2010), and of particular note the addition of the ocean pole tide displacement and the change to the cubic plus linear conventional model for the mean pole motion; • higher-order (at least 2nd order) ionospheric and updated tropospheric models for propagation delays; • implementation of new attitude models for eclipsing satellites (Kouba 2009;Dilssner 2010;Dilssner et al. 2011) by some but not all ACs; and (b) repro1 Y-axis orbital frame translations (a) repro1 Z-axis orbital frame rotations As with repro1, ACs were asked not to apply model corrections for the load displacements caused by large-scale non-tidal atmosphere, ocean and hydrological surface fluid motions. The primary seasonal components of the load effect on station positions and EOPs, as well as all other annual and semiannual signals, were removed empirically by explicit fitting in the subsequent long-term stacking process used by Altamimi et al. (2016) to form ITRF2014. While the above repro2 standards were generally implemented, some AC software changes departed from the recommended standards (see Table 2 for explanation of AC abbreviations). The main known departures are: • COD • assumes nominal attitude during eclipses for GPS and GLONASS (applies to operational products only since broadcast clocks provided in CF2 SP3 files) • ocean pole tide not applied (displacements preapplied for SINEX combination) • EMR • assumes nominal attitude during eclipses for GPS Block IIF satellites; other GPS satellites use yawsteering model • ESA • assumes nominal attitude during eclipses for GPS • no modeling of Earth albedo and antenna thrust • ocean pole tide not applied (displacements preapplied for SINEX combination) • 2nd-order ionospheric correction not applied

• GFZ
• site displacements due to atmospheric S1/S2 tidal loading applied without applying associated centerof-mass offset corrections to SP3 file • GRG • site displacements and SP3 center-of-mass corrections due to atmospheric S1/S2 tidal loading applied • JPL • 30-hour data spans and orbital arcs • MIT • ocean pole tide not applied (displacements preapplied for SINEX combination) • SRP parameters constrained between days over 9 d window • applied non-tidal atmospheric pressure loading in processing; effects removed from SINEX • no modeling of orbit perturbations due to ocean tidal geopotential variations Some of the departures listed are small and have limited implications for the combination products. Others, like the application of day-boundary constraints on SRP parameters by MIT, have unknown effects. Still others have potentially dramatic impacts, depending on weighting of the offending AC product in the combination. For instance, radial orbit errors caused by neglecting the effects of ocean tidal geopotential variations can reach 20 centimeters on average over the long term for a geodetic satellite in low earth orbit (Petit et al. 2010). For GPS, the impact is < 0.1 mm on average over the long term, but subdaily perturbations can exceed 42 centimeters 3D RSS. So, the effects are of course dampened for GPS satellites due to their altitude, but are still significant for high-accuracy applications. Also, the lack of consistent yaw modeling among clock ACs has a large impact on the satellite clock combination. As a practical matter, mixing input clock solutions derived with differing yaw models complicates usage of the IGS combination clock product because the user is unable to match their own selected model with that of the IGS clocks. At least one AC made changes to their operational products since their repro2 submissions. ESA made several orbit modeling changes, iterating on the most suitable box-wing model for GPS satellites (Springer 2017). A change was made to their operational products at GPS Week 1892, dropping altogether the box-wing model for GPS Block IIF. This put the ESA orbits into strong alignment with JPL until about GPS Week 1938 (http://acc.igs.org/igsacc_final.html). Consequently, and also because the National Geodetic Survey (NGS) and Scripps Institute of Oceanography (SIO) ACs did not contribute to the repro2 effort, the consistency between the IG2 products and the follow-on operational Finals is not optimal.

Repro2 AC contributions and usage
Nine analysis groups submitted solutions for repro2 (Table 1): seven IGS Finals operational ACs and two centers from the IGS tide gauge benchmark monitoring working group (TIGA; Schöne et al. 2009), which primarily serve to densify the tracking network with GNSS stations that are co-located with tide gauges, tying the tide gauge measurements directly to ITRF2014. However, one operational AC and both TIGA centers were ultimately excluded from the orbit and clock combinations for reasons described below. As Table 1 and Fig. 3 show, each contributed solution spans a different segment of time. They start as early as January 2, 1994, and extend to at least the end of 2013. Generally, those end dates correspond to when an AC completed their repro2 software changes in their IGS Finals operational products. For the TIGA groups, the time spans were determined by their internal group-specific requirements. AC operational products were used to extend AC solutions uniformly through GPS Week 1831. The approach outlined in Table 1 is consistent with Table 1 of Rebischung et al. (2016) for the IG2 station products.
AC repro2 orbits and clocks were used to form the combination products in a way that closely matches what was done in forming the IGS combined SINEX files . This was done to maintain consistency between the orbital and terrestrial reference frames expressed by the  Fig. 3 Graphical representation of the repro2 orbit and clock inputs described in Table 1 joint set of IG2 products (Table 2). One exception is the use of GFZ orbits and clocks instead of GTZ. This was done for purely conventional reasons because the GFZ products contribute to the operational product stream and the GFZ orbits and clocks are nominally identical to those of GTZ.
The tracking network used for GTZ was simply expanded to include TIGA stations in the GTZ submission. This decision has no measurable impact on the combination results. The other exception to the IG2 SINEX selections is in the handling of MIT contributions to the orbit and clock combi-  nations. Biases exceeding 0.50 ppb in their orbit scale and many clock jumps and biases exceeding 500 ps should preclude their use entirely, but excluding their orbits over the full repro2 span introduces a 5-mm WRMS disparity in IG2 orbits compared to IG1. So MIT orbits were included after GPS Week 1249 in order to minimize discrepancies with the IG1/F products in the earlier years 1 . The MIT clocks were also problematic due to large biases and frequent jumps. They were largely excluded, though they were used to fill-in gaps in the combination for the earliest years (GPS Weeks 730-779 and 1111-1115) because there are too few good AC clock solutions to form a robust and precise IG2 CLK product. Moreover, no 30-sec satellite clock combinations were formed due to an insufficient number (i.e., < 3) of useable (e.g., detectable jumps and removable biases) AC solutions. 1 A test was performed to quantify the impact of including/excluding MIT orbits from the combination. Including MIT orbits after GPS Wk 1249 did not degrade the combination with respect to the IG1/IGF. Apparently, there are sufficient AC inputs to largely suppress the ocean tidal geopotential errors contributed by the MIT orbits, but not the large orbit scale issues prior to GPS Wk 1249. The issues around including/excluding MIT were discussed with the IGS Analysis Center Coordinator priori to finalizing the IG2 combinations; no concerns were raised with the proposal to include MIT after GPS Week 1249.

IGorbit and clock combinations
The combination models used for generating the IG2 orbits and clocks are identical to those used by the Analysis Center Coordinator (ACC; acc.igs.org) in generating the IGS operational Final products. The underlying orbit combination model has been the same since the early days of the IGS (Beutler et al. 1995;Kouba et al. 1995;Kouba andMireault 1996, 1997), though a change to the combination pre-processing steps was introduced in early 2000 in order to co-align the IGS orbital and terrestrial reference frames (Kouba et al. 1998;Ferland et al. 2000, IGSMAIL-2740Ferland et al. 2000;Springer 2000, IGSMAIL-2750 As with the IG1 reprocessed products, the full suite of repro2 orbit and clock combination products (see Table 3) are available the IGS Global Data Center at NASA (i.e., ftp://cddis.gsfc.nasa.gov/gps/products/wwww/repro2 or ftp: //cddis.gsfc.nasa.gov/gps/products/repro2). The ig2ww ww7.erp and ig215p01.erp files are from the ACC orbit and clock combination system. They are derived for the purpose of monitoring the official ERP products generated in the SINEX combination (Ferland et al. 2000;Rebischung et al. 2016). The main differences between the ACC and SINEX ERPs are the use of a posteriori AC orbit residual statistics as a priori weights to the AC ERPs and the absence of station covariances when they are combined in the ACC system. The definitive EOP products are provided by the Reference Frame Combination Center at the Institut National de l'Information Géographique et Forestière (IGN) in France and are the ones to be used with the orbits.

Orbit combination statistics
In addition to the resulting combined ephemerides, the orbit combination generates post-fit residuals for each AC after removing seven Helmert transformation parameters (rotation, translation, scale) for each daily input orbit relative to the daily combined orbit. The time series of daily WRMS of the input orbit residuals computed from the repro2 combination results are shown in Fig. 4 (smoothed). The WRMS in Fig. 4 are akin to those in Fig. 1, but for the repro2 inputs. With exception of UL2 and GR2, the WRMS agreement of the input solutions and IG2 trends approximately exponentially over the repro2 time span. The WRMS agreement between IG1/F and IG2 improves from ∼ 28 mm in the earliest years converging to ∼ 5 mm at GPS Week 1250. Their mutual agreement improves further to about ∼ 2.5 mm over a time span of about 27 weeks, from GPS Week 1773 (December 29, 2013) until about GPS Week 1800 (July 6, 2014). Recall that this is the period when COD, EMR and ESA operational products were introduced to the repro2 combinations (Table 1). These are three of the more highly weighted inputs in the operational Finals combinations (http://acc.igs. org/igsacc_final.html). The remaining 2.5-mm WRMS difference between IG2 and IG1/F after GPS Week 1800 is attributed mostly to the absence of a repro2 solution from the NGS AC (the missing SIO AC solution having little impact). This result suggests that the absence of an NGS repro2 solution introduces a background difference in IG2 compared to IG1/F at the ∼ 2.5 mm WRMS level.

Translations
Helmert orbital frame translations estimated in the repro2 orbit combinations are plotted in Fig. 5 for each geocentric coordinate. All ACs have nearly annual motions in the Ztranslation with various phases, but approximately similar amplitudes, except for GRG whose annual motions some-times reach 20 mm for unexplained reasons. Meanwhile, the equatorial translations are nearly featureless for all ACs.

Rotations
Helmert rotations about each geocentric axis coordinate from the repro2 orbit combinations are plotted in Fig. 6. Large variations sometimes occur about all three axes for most ACs prior to GPS Week 1150(January 20, 2002. The most strik-ing result, however, is the large (up to 100 µas) Rx and Ry offsets for IG1/F prior to GPS Week 1702, which was caused by the ACC frame alignment software bug described by Griffiths (2012) and mentioned earlier.

Scale
Helmert scale offsets estimated in the series of repro2 orbit combinations are plotted in Fig. 7 , but the MIT products were excluded during this period so they do not impact the IG2 orbit scale. The EMR and JPL orbital scales track each other closely because they both use GIPSY-OASIS for their data analysis. ESA and IG1/F vary similarly for GPS Weeks 1050-1550, but then afterward drift independently. The average scale offset of ∼ 0.55 ppb (1.46 cm at GPS) between IG2 and IG1/F before Week ∼ 1700 is smaller than the ∼ 0.72 ppb (Ray 2012) offset expected from adopting Earth albedo and antenna thrust models in the repro2 analyses. The missing 0.17 ppb part of the IG2 scale change is attributed primarily to the fact that ESA did not adopt these two orbit modeling changes for repro2. The IG1/F scale converges toward IG2 starting at GPS Weeks 1702 (August 19, 2012), which, as mentioned above, coincides with the switch to daily products and the software fix for applying AC SINEX rotations to the AC orbits. The faster convergence of IG2 and IG1/F after GPS Week 1740 (August 21, 2011) is likely due to operational ACs adopting the repro2 Earth albedo and antenna thrust models in their contributions to the operational IGS Finals products. With the exception of ESA, these model changes were completed by about GPS Week 1773 (December 29, 2013). Figure 8 shows weekly averages for the RMS statistics of individual repro2 AC clock solutions with respect to the combined 5-min clock product. The RMS statistics include separate biases for each satellite and station clock, which are computed and removed before generating the final IG2 clocks, when possible. The step in the combination process that computes the biases requires three reliable (few jumps, distinguishable biases) input clock products. When the number of reliable clock products falls below three, the ability to reliably determine the clock biases fails and a large number of AC clock rejections can occur (see for example Fig. 9). When the step for determining AC biases fails, the combination becomes unstable and the combined clock product gets contaminated with AC biases, which appear as spikes in the AC RMS curves like those near GPS Week 1523 in the JP2 series. The step decrease in the RMS curve for IGF at GPS Week 1631 (April 17, 2011) corresponds to when the IGS switched from the IGS05 to IGS08 framework (including station coordinates and antenna calibrations) in the IGS Finals operational product stream. The step increase in the RMS curve for JP2 at GPS Week 1773 (December 29, 2013) corresponds to when the CODE operational clocks (and orbits) were incorporated into the repro2 combination. Recall from Fig. 3 that CODE did not contribute repro2 clock estimates.

Clock combination statistics
It was mentioned earlier that the clock combination requires at least three good input AC solutions to form a robust product. Another complicating factor negatively impacting the combined clocks is the fact that ESA assumes nominal attitude through satellite eclipse season. Other ACs adopted a yaw-steering model, but not necessarily the same one. During these periods the resulting clock inconsistencies cause automatic rejections to increase in the combination, which can sometimes become unstable and unreliable if insufficient numbers of usable clock ACs remain.

Quality assessment of the IG2 orbits and clocks
Determining the quality of the repro2 products provides insights into their utility for accessing the IGS reference frame for high-accuracy positioning, navigation and timing applications, as well as indicating the progress of AC modeling efforts. While there is no method to validate IGS orbits, several approaches have been developed for detecting systematic errors. Analysis of positional discontinuities at the midpoint epoch between successive daily SP3 files (Griffiths 2009a, b; Griffiths and Ray 2009) has been extremely valuable for detecting limiting errors in the IERS Conventions model for subdaily EOP variations , among other things. A refined version of that 2013 analysis is currently underway for the IG2 orbits and is the subject of a manuscript in preparation. Satellite laser ranging (SLR) has also been used to great effect (e.g., Urschl et al. 2005;Ziebart et al. 2007;Sośnica et al. 2015Sośnica et al. , 2017). An anonymous reviewer of this paper used IG2 orbits to repro- cess International Laser Ranging Service (ILRS; ilrs.org) SLR data for the two GPS satellites (SNV35 and SVN36) equipped with a laser retroreflector array. They apparently found that the RMS difference between the optical and microwave ranges is reduced by 1.5 cm compared to those when using IG1/IGF orbits. That result is consistent with the main IG2-IG1/IGF orbit scale offset discussed earlier (Fig. 7), which was an expected result given the repro2 orbit model changes for Earth albedo and satellite antenna thrust-ing. The reviewer also found that the standard deviation of SLR residuals decreased by 1-2 mm, which may be due to improved rotational stability of the orbital frame after correcting the combination software bug discussed earlier.
One must be cautious, however, not to overstate SLR's role in validating orbit accuracy. Validation is a very high standard. With ranging precisions approaching 1 mm, SLR is highly sensitive to variations in the orbit radial directionthe most accurately modeled orbit component due to Kepler's  (or May 1, 2000, andAugust 14, 2005). May 1, 2000, has some significance because that is the date when the US Air Force turned off Selective Availability (SA). It is unlikely that turning off SA is responsible for the spike in AC clock rejections. Rather, better inter-AC clock agreement in the pre-SA-off period is relative to some background noise level. In the pre-SA-off era, when there is a tremendously large noise background, it is likely that the clock combo algorithm simply cannot detect when an AC clock value is bad or not, except when they are egregious. On the other hand, when SA gets turned off, the background satellite clock noise drops by a very large amount and the algorithm is able to better detect discrepant AC values 3rd law-but SLR is relatively insensitive to along-track and cross-track errors 2 . Moreover, the SLR technique continues 2 The sensitivity decreases in proportion to the square root of the range. For GPS SVN35, a 1-mm SLR error is consistent with an orthogonal to suffer other limitations affecting its accuracy, as illustrated by long-standing issues with range biases (e.g., Appleby Footnote 2 continued position offset of 199.5 m to 202.3 m, depending on its position in the slightly elliptical orbit.  (Petit et al. 2010), which for GPS can exceed 2.5 cm in the along-track and cross-track directions . For all of these reasons, the next two subsections are focused on results from a longarc orbit analysis, following , and a precise point positioning (PPP) analysis to determine the utility of IG2 orbit and clock products for long-term reprocessing by IGS users.

Long-arc analysis
Estimates of orbit precision were derived from a long-arc analysis of each weekly set of IG2 and IG1/F orbits, where seven days of combined SP3c files are fit for 15 parameters of the extended CODE orbit model (ECOM; Springer et al. 1998)-a six-parameter orbit state vector (position and velocity) and the nine terms of an empirical harmonic SRP forcing model expressed in a satellite body-fixed frame. Then RMS, WRMS and median statistics are computed for the orbit residuals on a daily basis. These statistics are reported in the companion *.sum files noted in Table 3 of Section 2. This procedure is identical to that used for the operational Finals summaries. Table 4 shows the ensemble averages and standard deviations for the daily median statistics from the long-arc analysis for IG1/F and IG2 orbits. These ensemble statistics are computed over the full repro2 time span (GPS Weeks 0730 thru 1831, or 1994.003 thru 2015.121) and are taken to represent the internal long-term precision of the IG2 and IG1/F orbits. The question is whether the precision of IG2 is significantly different from that of IG1/F. A z-test comparing the average median from IG2 and IG1/F indicates that the precision of the two solution series is indeed significantly different at the 99.9% confidence interval, and that IG2 is less precise than IG1/F. The root of the square differences for the values in Table 4 over the full time span is computed to approximate the magnitude in the loss of precision for IG2: 2 As mentioned earlier, the orbit combination results (Fig. 4) indicate that about half of the WRMS discrepancy between IG2 and IG1/F can be attributed to not having a repro2 contribution from the NGS AC. It then follows that half of the 7.31-mm loss in precision for IG2 is attributable to the absence of an NGS repro2 solution. That is, not having a repro2 solution from NGS introduces a significant discrepancy with respect to IG1/F and reduces the overall precision of the IG2 orbits. In that case, the remaining half of the 7.31 mm RDS medi are likely due to other errors in other highly weighted AC orbits (i.e., EMR, ESA, JPL) that do not appear to exist in their operational solutions. For instance, errors committed in ES2 solutions associated with using a boxwing a priori solar radiation pressure model (Springer et al. 2014) that was later found to be erroneous (Springer 2017). Of course, there is a question as to whether ECOM is appropriate for the IG2 long-arc analysis given inter-AC differences for SRP modeling. That is, perhaps ECOM (Beutler et al. 1994) is simply more consistent with the IG1/F orbits because a so-called reduced-parameter version (i.e., 6-parameter state vector and five terms of the full nine term truncated Fourier series) of ECOM is the most common among the ACs, and so the differences in the long-arc statistics largely reflect model inconsistency. We have no way to answer this question directly. However, the large number of degrees of freedom in the full ECOM should be adequate to represent the IG2 orbits for this analysis. The AC-specific long-arc results in Table 5, where no bias in favor of ACs using ECOM is seen, seem to support this assertion. Therefore, the author believes it is the IG2 orbits that are less precise even for the second half of the repro2 time span. This loss in Table 6 Average and standard deviations for post-fit residuals after Helmert transformation comparison of a PPP network of coordinates to the IG2 daily terrestrial frames, which are expressed in IGb08. Statistics are given over the full repro2 time span and for a more recent subset in order to gauge their utility over the last decade of repro2. The smaller biases and scatter in the station position residuals indicate higher precision in the latter span of the products, but the IG2 products are still not as precise as IG1/F

PPP analysis
A PPP analysis was performed to assess the utility of the IG2 orbit and clock products as a means for accessing the IGS reference frame in a long-term PPP reprocessing by generating coordinate sets consistent with IGS08 and transforming them to IGS14 via the IERS transformation parameter values . The overall analysis approach follows the one used operationally by the IGS ACC. The Bernese v5.2 software  was employed to estimate daily positional coordinates for 163 core IGS14 stations (Rebischung et al. 2016, IGSMAIL-7399) separately using IG1/F and IG2 orbit and clock products. IERS 2010 Conventions are generally implemented, including GMF/GPT2 (Böhm et al. 2006;Lagler et al. 2013) troposphere modeling with daily tropospheric delay gradients following Chen and Herring (1997). Plus, • RINEX files for each IGS14 core station were downloaded from NASA/CDDIS (cddis.gsfc.nasa.gov); • daily SP3 satellite files were fit using IG2 official EOPs and extended CODE (6+9) orbit model in order to provide satellite positions at the sampling interval of the data files; • IG2 clocks and station data were preprocessed to detect and eliminate clock jumps and cycle slips; • code and phase observational data were reduced with iterative post-fit cleaning; outliers are phase residuals exceeding 25 mm and code residuals exceeding 2.5 m; • floating-point phase ambiguity parameters were estimated; and • daily coordinate sets were compared to IG2 combined daily SINEX solutions, which were aligned to IGb08 by Rebischung et al. (2016), resulting in daily sets of PPP-derived terrestrial frame parameters and station coordinate residuals.
Compared to IG1/F, the IG2 orbits and clocks deliver longterm global coordinate sets that are noisier and maybe slightly more biased (Table 6). The added noise in the long-term PPP results is caused by a high number of station data rejections in the PPP processing prior to GPS Week 1336. This period corresponds to when there were large numbers of satellite clock rejections in the clock combination ( Fig. 9), implying that the AC input clocks were of poor quality that degrades the PPP results. The time series of PPP-derived Helmert parameters for the daily networks are shown in Figs. 10 and 11 with ensemble statistics summarized in Table 7. Like the position residuals, the frame parameters based on IG2 orbits and clocks are also clearly more scattered than for IG1/F, especially before GPS Week 1336. This is evident in the standard deviations shown in Table 7. In Fig. 10, the convergence of the IG2 and IG1/F Tz series after adopting IGS08 in the IGS operational products is consistent with the Tz offset between IGS05 and IGS08. The elimination of spurious Rx and Ry rotations for IG2 compared to IG1/F when the s/w bug fix was implemented at GPS Week 1702 in the IGF operational series (Fig. 11) is consistent with the orbit combination results (Fig. 6a,b).

Summary
IG2 GPS orbits and clocks have been derived from the contributions of 9 analysis groups (6 included with nonzero weights) that reprocessed more than 21 years of ground tracking data. All IG2 orbit and clock product files (ig2wwwwd.sp3, ig2wwwwd.clk), along with the weekly summary reports (ig2wwww7.sum), have been adopted by the IGS and are publicly available. These are in addition to the IG2 terrestrial reference frame products, which formed the IGS contribution to ITRF2014 ). The overall inter-AC orbit agreement is at about the ∼ 4 mm WRMS level (1D) with outlier ACs at exceedingly large levels, reaching several centimeters (GRG) to decimeters (ULR). The agreement between IG2 and IG1/F persists at the ∼ 5 mm WRMS level for most of the repro2 time span, until the most recent years when operational AC Table 7 Average and standard deviations of Helmert transformation parameters for the same PPP-based coordinate sets as used for the results in Table 6, and also covering two different time spans: full repro2 and since 2005-08-14. The large IG2 rotations and translations are due to frequently high numbers of station rejections. The stations were rejected during the iterative outlier rejection steps because of large biases in the IG2 satellite clocks, for instance, prior to 2005-08-14. In the more recent years, the IG2 products have smaller long-term frame errors than IG1/F, but the frame less stable at short time intervals as indicated by the larger standard deviations in the IG2 frame parameters since 2005-08-14. The instability arises also because of poorer IG2 clocks compared to IG1/F clocks products were used to extend their submissions uniformly through GPS Week 1831 (February 14, 2015). For those later weeks, the IG2-IG1/F WRMS disagreement drops to 2.5 mm. Missing from this mix, however, is a reprocessed contribution from NGS (as well as SIO), and that is probably the biggest source for the background 2.5-mm WRMS difference between IG2 and IG1/F orbits. Large rotational offsets of the orbital frame prior to GPS Week 1702 have been eliminated due to a software bug that was fixed in the operational Finals products and inherited by the repro2 combinations. Absence of Earth albedo and satellite antenna thrust models by ESA dampens an expected 0.72 ppb scale change due to these effects. The scale shift realized in the IG2 orbits was ∼ 0.55 ppb, leaving a ∼ 0.17 ppb residual error when IG2 orbits are used and Earth albedo and antenna thrust models are adopted by the user.
A long-arc analysis was performed as a measure of IG2 orbit precision. The root square difference of long-arc post-fit residuals for the IG2 and IG1/F orbits indicates an increased uncertainty of ∼ 7.3 mm (1D RMS) for IG2. However, while less precise than IG1/F, the IG2 orbits offer at least two advantages in a long-term reprocessing: they are internally more self-consistent over the full history, and the orbital frame is much better aligned to the ITRF. The ∼ 0.17 ppb orbit scale error present should be mitigated in applications using network processing with either double differencing or explicit clock estimation. The same is probably true for mitigating the effects of IG2 lower precision orbits.
The satellite clocks are severely limited by large residual biases and incompatible satellite attitude models adopted by the ACs, and therefore the IG2 clocks should not be used for long-term PPP reprocessing.
Recommendations for the next IGS reprocessing are based on results from this study, along with results from other published work as cited below for the sake of completeness, and include: • Need improved orbit modeling and AC participation • Forego clock submission/combination if number and quality of submissions are insufficient for robust combined product, including 30-s satellite clocks ideally • Need improved models for subdaily variations in Earth orientation due to ocean tides ( As discussed in the introduction of this paper, other improvements at GNSS tracking stations could benefit IGS products overall, for example by mitigating multipath effects (especially from near-field reflectors), reducing unnecessary discontinuities due to equipment changes and enabling thermal expansion corrections to be computed by collecting the relevant metadata, including but not limited to monument dimensions, materials properties, descriptions of cabling lengths above and below ground.