Abstract
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.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 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 \(\sim 21\) years, starting January 2, 1994, and continuing through February 14, 2015, or GPS Weeks 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.
1.1 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 (Rebischung et al. 2016) to ITRF2014 (Altamimi et al. 2016). 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.
Time series of AC weighted RMS (WRMS) statistics from the repro1 orbit combination (Gendt et al. 2010). PD1 is for the Potsdam/Dresden group, which contributed as an extra solution. SI1 is for the solution contributed by Scripps Institution of Oceanography, an operational IGS AC. The WRMS statistics are computed using the post-fit residuals of the combination for each contributed solution. The curve for IGF represents the WRMS difference for the operational IGS Finals 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
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).
1.2 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 (Ray et al. 2013); 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 station 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 long-term reference (Ray et al. 2017); subdaily EOP alias and draconitic errors in the satellite orbits (Griffiths and 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.
In addition to WRMS statistics for post-fit residuals in Fig. 1, the orbit combination estimates AC orbital frame transformation parameters with respect to the combined orbit. a Time series of Helmert orbit frame rotations about the Z-axis for AC repro1 contributions and IGF, and b time series of Helmert origin translations in the Y-direction for the same AC solutions (Griffiths et al. 2009). The two plots in (a, b) illustrate historical frame errors removed with the adoption of IGS05 and reprocessing the full history of GPS data for the IG1 data products
1.3 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
-
modeling of Earth radiation pressure (Rodriguez-Solano et al. 2011a, b) including also the thrust acting on the spacecraft caused by signal transmission along the satellite antenna boresight (Rodriguez-Solano et al. 2012; acc.igs.org/orbits/thrust-power.txt).
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 pre-applied for SINEX combination)
-
-
EMR
-
assumes nominal attitude during eclipses for GPS Block IIF satellites; other GPS satellites use yaw-steering model
-
-
ESA
-
assumes nominal attitude during eclipses for GPS
-
no modeling of Earth albedo and antenna thrust
-
ocean pole tide not applied (displacements pre-applied for SINEX combination)
-
2nd-order ionospheric correction not applied
-
-
GFZ
-
site displacements due to atmospheric S1/S2 tidal loading applied without applying associated center-of-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 pre-applied 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
-
Graphical representation of the repro2 orbit and clock inputs described in Table 1
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\,\hbox {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.
1.4 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 (Rebischung et al. 2016). This was done to maintain consistency between the orbital and terrestrial reference frames expressed by the 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 combinations. 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 yearsFootnote 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.
WRMS statistics for repro2 AC post-fit residuals from the orbit combination. The IG1 and IGF orbits do not contribute to the combined IG2 orbits; they were included for comparison only. While UL2 orbits were also included for comparison only, the UL2 residuals exceeded 110 mm and do not fit in the scale of this plot
2 IG2 orbit 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 and Mireault 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-2740; Ferland et al. 2000; Springer 2000, IGSMAIL-2750). The change meant applying small Helmert rotations, determined from the IGS SINEX combination, to the AC Final orbits prior to combining them into the IGS Final orbits. During the switch to daily SINEX integrations in 2012, however, a bug in the combination software was discovered. The sets of SINEX TRF rotational offsets used for co-aligning the AC orbits were being applied with incorrect signs (Griffiths 2012). Correcting the software bug starting with GPS Week 1702 (August 19, 2012) eliminated spurious rotations of the IGS Final orbital frame, which reached \(100\,{\upmu }\)as about the Y-axis. The combination model and strategy used for the clocks are the one introduced by Kouba and Springer (2001).
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 ig2wwww7.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.
2.1 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 \(\sim 28\,\hbox {mm}\) in the earliest years converging to \(\sim 5\,\hbox {mm}\) at GPS Week 1250. Their mutual agreement improves further to about \(\sim 2.5\,\hbox {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 \(\sim 2.5\,\hbox {mm}\) WRMS level.
2.2 Orbital frame
2.2.1 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 Z-translation with various phases, but approximately similar amplitudes, except for GRG whose annual motions sometimes reach 20 mm for unexplained reasons. Meanwhile, the equatorial translations are nearly featureless for all ACs.
![](http://media.springernature.com/lw610/springer-static/image/art%3A10.1007%2Fs00190-018-1149-8/MediaObjects/190_2018_1149_Fig5b_HTML.gif)
2.2.2 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 striking result, however, is the large (up to \(100\,{\upmu }\)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.
Same as in Fig. 5, but for rotations about the geocentric a X-axis, b Y-axis and c Z-axis
![](http://media.springernature.com/lw610/springer-static/image/art%3A10.1007%2Fs00190-018-1149-8/MediaObjects/190_2018_1149_Fig6b_HTML.gif)
2.2.3 Scale
Helmert scale offsets estimated in the series of repro2 orbit combinations are plotted in Fig. 7. MIT orbits have large spurious offsets prior to GPS Week 1250 (December 21, 2003), 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 \(\sim 0.55\) ppb (1.46 cm at GPS) between IG2 and IG1/F before Week \(\sim 1700\) is smaller than the \(\sim \,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).
2.3 Clock combination statistics
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.
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.
Number of AC GPS satellite clock rejections per week in the combination. The rejections are determined per epoch of the clock sampling, so \(\sim 500\) rejections compared to clock estimates for 32 satellites sampled at 5-minute intervals, for instance. The same rejection criteria were used for the entire combination period. The ACs contributing to the clock combinations are listed, with the earliest years beginning with EM2, GF2, JP2 and MI2. MI2 was given zero weight when ES2 was added at GPS Week 780, except for when MIT was used to fill-in gaps between GPS Week 1111 and 1115. The CODE operational clocks (COF) were added at GPS Week 1772. The number of rejections is especially large for the period between GPS Weeks 1060 and 1335 (or May 1, 2000, and August 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
3 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 (Griffiths and Ray 2013), 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. 2015, 2017). An anonymous reviewer of this paper used IG2 orbits to reprocess 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 thrusting. 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 direction—the most accurately modeled orbit component due to Kepler’s 3rd law—but SLR is relatively insensitive to along-track and cross-track errorsFootnote 2. Moreover, the SLR technique continues to suffer other limitations affecting its accuracy, as illustrated by long-standing issues with range biases (e.g., Appleby et al. 2016). This is compounded by the fact that SLR to GPS is limited to two old satellites with relatively infrequent observations that ended in 2013.347 (SVN35) and 2015.835 (SVN36). All of this ignores inter-technique modeling errors inherent in the IERS Conventions (Petit et al. 2010), which for GPS can exceed 2.5 cm in the along-track and cross-track directions (Griffiths and Ray 2013). For all of these reasons, the next two subsections are focused on results from a long-arc orbit analysis, following Griffiths and Ray (2009), and a precise point positioning (PPP) analysis to determine the utility of IG2 orbit and clock products for long-term reprocessing by IGS users.
3.1 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:
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 \(\hbox {RDS}_{\mathrm{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 box-wing 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 precision impacts PPP, but probably not results where double differencing is used because the errors should be largely common mode and minimized or eliminated by the differencing technique.
3.2 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 (Altamimi et al. 2016). The overall analysis approach follows the one used operationally by the IGS ACC. The Bernese v5.2 software (Dach et al. 2015) 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 long-term 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).
PPP-derived frame translations for IG2-based (gray) and IG1/F-based (purple) solutions compared to IG2 SINEX frames that are aligned to IGS08 (Rebischung et al. 2016). The solid curves are from a 7d boxcar average. The large scatter appearing as noise in the earlier weeks of the gray series reflects instability of the PPP-based IG2 frame caused by large numbers of station rejections (due to satellite clock issues). The blue vertical line marks the switch from the IGS05 frame to IGS08 in the IGF operational stream. The magenta vertical line marks the timing of the s/w bug fix in IGF
Same as Fig. 10 but for rotation and scale
4 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 (Rebischung et al. 2016).
The overall inter-AC orbit agreement is at about the \(\sim 4\,\hbox {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 \(\sim 5\,\hbox {mm}\) WRMS level for most of the repro2 time span, until the most recent years when operational AC 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 \(\sim 0.55\) ppb, leaving a \(\sim 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 \(\sim 7.3\,\hbox {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 \(\sim 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 (Griffiths and Ray 2013)
-
Full implementation of IERS Conventions by all ACs (with changes below)
-
Modify IERS mean pole model to older linear form to agree with Wahr et al. (2015)
-
Adopt model for seasonal variations of the low-degree geopotential terms, preferably in agreement with International Laser Ranging Service (ILRS) and International Doris Service (IDS) (e.g., Melachroinos et al. 2014)
-
Enforce a common improved satellite attitude model for all ACs to improve inter-AC clock consistency, especially during satellite Earth eclipses
-
Require full consistency with IGS and updated IERS Conventions by all ACs
-
Establish well-defined analysis standards and reject non-compliant AC solutions
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.
Notes
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.
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 position offset of 199.5 m to 202.3 m, depending on its position in the slightly elliptical orbit.
References
Altamimi Z, Collilieux X, Legrand J, Garayt B, Boucher C (2007) ITRF2005: a new release of the International Terrestrial Reference Frame based on time series of station positions and Earth Orientation Parameters. J Geophys Res Solid Earth 112:B09401. https://doi.org/10.1029/2007JB004949
Altamimi A, Collilieux X, Métivier L (2011) ITRF2008: an improved of the International Terrestrial Reference Frame. J Geod 85(8):457–473. https://doi.org/10.1007/s00190-011-0444-4
Altamimi A, Rebischung P, Métivier L, Collilieux X (2016) ITRF2014: a new release of the International Terrestrial Reference Frame modeling nonlinear station motions. J Geophys Res 121(8):2169–9356. https://doi.org/10.1002/2016JB013098
Appleby G, Rodriguez J, Altamimi Z (2016) Assessment of the accuracy of global geodetic satellite laser ranging observations and estimated impact on ITRF scale: Estimation of systematic errors in LAGEOS observations 1993–2014. J Geod 90(12):1371–1388. https://doi.org/10.1007/s00190-016-0929-2
Beutler G, Brockmann E, Gurtner W, Hugentobler U, Mervart L, Rothacher M (1994) Extended orbit modeling techniques at the CODE processing center of the International GPS Service for geodynamics (IGS): theory and initial results. Manuscr Geod 19:367–386
Beutler G, Kouba J, Springer T (1995) Combining the orbits of the IGS Analysis Centers. Bull Geod 69:200–222
Böhm J, Niell A, Tregoning P, Schuh H (2006) Global Mapping Function (GMF): a new empirical mapping function based on numerical weather model data. Geophys Res Lett. https://doi.org/10.1029/2005GL025546
Chen G, Herring T (1997) Effects of atmospheric azimuthal asymmetry on the analysis of space geodetic data. J Geophys Res 102(B9):20489–20502. https://doi.org/10.1029/97JB01739
Dach R, Lutz S, Walser P, Fridez P (eds) (2015) Bernese GNSS software version 5.2. user manual. Astronomical Institute, University of Bern, Bern Open Publishing. https://doi.org/10.7892/boris.72297, ISBN: 978-3-906813-05-9
Dilssner F (2010) GPS IIF-1 satellite antenna phase center and attitude modeling. Inside GNSS 5(6):59–64
Dilssner F, Springer T, Gienger G, Dow J (2011) The GLONASS-M satellite yaw-attitude model. Adv Space Res 47(1):160–171. https://doi.org/10.1016/j.asr.2010.09.007
Dow JM, Neilan RE, Rizos C (2009) The International GNSS Service in a changing landscape of Global Navigation Satellite Systems. J Geod 83(3–4):191–198. https://doi.org/10.1007/s00190-008-0300-3
Ferland R, Piraszewski M (2009) The IGS-combined station coordinates, earth rotation parameters and apparent geocenter. J Geod 83:385–392. https://doi.org/10.1007/s00190-008-0295-9
Ferland R, Kouba J, Hutchison D (2000) Analysis methodology and recent results of the IGS network combination. Earth Planets Space 52(11):953–957. https://doi.org/10.1186/BF03352311
Gendt G, Griffiths J, Nischan Th, Ray J (2010) IGS reprocessing—summary of orbit/clock combination and first quality assessment. IGS 2010 Workshop, Newcastle upon Tyne, UK, June 2010 http://acc.igs.org/repro1/repro1_IGSW10.pdf
Griffiths J (2009a) IGS reprocessed GPS satellite orbits: preliminary results for period 2005–2007. Unpublished http://acc.igs.org/orb-repro1_old.html. Accessed 19 May 2017
Griffiths J (2009b) IGS reprocessed GPS satellite orbits: preliminary results for 3 equal periods from Aug 1999 to Dec 2007. Unpublished http://acc.igs.org/orbrepro1.html. Accessed 19 May 2017
Griffiths J (2012) A review of the IGS Final orbit combination procedures. Unpublished report, http://acc.igs.org/orbits/acc_report_final_rotations.pdf
Griffiths J, Choi K (2013) Analysis center coordinator report. In: Dach R, Jean Y (eds) International GNSS Service technical report 2012. University of Bern, Bern, pp 21–34. https://doi.org/10.7892/boris.80303
Griffiths J, Ray J (2009) On the precision and accuracy of IGS orbits. J Geod 83:277–287. https://doi.org/10.1007/s00190-008-0237-6
Griffiths J, Ray J (2013) Sub-daily alias and draconitic errors in the IGS orbits. GPS Solut 17(3):413–422. https://doi.org/10.1007/s10291-012-0289-1
Griffiths J, Ray J (2016) Impacts of GNSS position offsets on global frame stability. Geophys J Int 204:480–487. https://doi.org/10.1093/gji/ggv455
Griffiths J, Gendt G, Nischan T, Ray J (2009) Assessment of the orbits from the 1st IGS reprocessing campaign. EOS Trans AGU 90(52), Fall Meet. Suppl., Abstract G13A-04. http://acc.igs.org/repro1/repro1-orbits_agu09.pdf
Kouba J (2009) A simplified yaw-attitude model for eclipsing GPS satellites. GPS Solut 13(1):1–12. https://doi.org/10.1007/s10291-008-0092-1
Kouba J, Mireault Y (1996) IGS analysis coordinator report. International GPS Service for Geodynamics (IGS) 1995 Annual Report, pp 45–76
Kouba J, Mireault Y (1997) Analysis coordinator report. International GPS Service for Geodynamics (IGS) 1996 Annual Report, pp 55–100
Kouba J, Springer T (2001) New IGS station and satellite clock combination. GPS Solut 4(4):31–36. https://doi.org/10.1007/PL00012863
Kouba J, Mireault Y, Lahaye F (1995) IGS orbit/clock combination and evaluation, Appendix I of the analysis coordinator report. International GPS Service for Geodynamics (IGS) 1994 Annual Report, pp 70–94
Kouba J, Ray J, Watkins MM (1998) IGS reference frame realization. In: 1998 IGS analysis center workshop proceedings. European Space Operations Centre, Darmstadt
Lagler K, Schingelegger M, Böhm J, Krásná H, Nilsson T (2013) GPT2: Empirical slant delay model for radio space geodetic techniques. Geophys Res Lett 40(6):1069–1073. https://doi.org/10.1002/grl.50288
McCarthy DD, Petit G, Conventions IERS (2003) IERS technical note 32. Verlag des Bundesamts für Kartographie und Geodäsie, Frankfurt am Main
Melachroinos S, Lemoine FG, Chinn DS, Zelensky NP, Nicholas JB, Beckley BD (2014) The effect of seasonal and long-period geopotential variations on the GPS orbits. GPS Solut 18(4):497–507. https://doi.org/10.1007/s10291-013-0346-4
Petit G, Luzum B, Conventions IERS (2010) IERS technical note 36. Verlag des Bundesamts für Kartographie und Geodäsie, Frankfurt am Main
Ray J (2012) Personal communication
Ray J (2016) Precision, accuracy, and consistency of GNSS products. In: Grafarend EW (ed) Encyclopedia of geodesy. Springer, Cham, pp 1–5. https://doi.org/10.1007/978-3-319-02370-0_97-1
Ray J, Altamimi Z, Collilieux X, van Dam T (2008) Anomalous harmonics in the spectra of GPS position estimates. GPS Solut 12(1):55–64. https://doi.org/10.1007/s10291-007-0067-7
Ray J, Griffiths J, Collilieux X, Rebischung P (2013) Subseasonal GNSS errors in IGS products. Geophys Res Lett 40:5854–5860. https://doi.org/10.1002/2013GL058160
Ray J, Rebischung P, Griffiths J (2017) IGS polar motion measurement accuracy. J Geod Geodyn. https://doi.org/10.1016/j.geog.2017.01.008
Rebischung P, Griffiths J, Ray J, Schmid R, Collilieux X, Garayt B (2012) IGS08: the IGS realization of ITRF2008. GPS Solut 16(4):483–494. https://doi.org/10.1007/s10291-011-0248-2
Rebischung P, Garayt B, Collilieux X, Altamimi Z (2013) IGS reference frame working group coordinator report. In: Dach R, Jean Y (eds) International GNSS Service technical report 2012. University of Bern, Bern, pp 171–178. https://doi.org/10.7892/boris.80303
Rebischung P, Altamimi A, Ray J, Garayt B (2016) The IGS contribution to ITRF2014. J Geod 90(7):611–630. https://doi.org/10.1007/s00190-016-0897-6
Rodriguez-Solano CJ, Hugentobler U, Steigenberger P (2011a) Earth radiation pressure model for GNSS satellites. In: Abstract EGU2011-9824 presented at EGU general assembly, Vienna
Rodriguez-Solano CJ, Hugentobler U, Steigenberger P, Lutz S (2011b) Impact of earth radiation pressure on GPS position estimates. J Geod 86(5):309–317. https://doi.org/10.1007/s00190-011-0517-4
Rodriguez-Solano CJ, Hugentobler U, Steigenberger P (2012) Adjustable box-wing model for solar radiation pressure impacting GPS satellites. Adv Space Res 7(1):1113–1128. https://doi.org/10.1016/j.asr.2012.01.016
Santamaría-Gómez A, Bouin M-N, Collilieux X, Wöppelmann G (2011) Correlated errors in GPS position time series: implications for velocity estimates. J Geophys Res. https://doi.org/10.1029/2010JB007701
Schmid R, Steigenberger P, Gendt G, Ge M, Rothacher M (2007) Generation of a consistent absolute phase center correction model for GPS receiver and satellite antennas. J Geod 81:781–798. https://doi.org/10.1007/s00190-007-0148-y
Schmid R, Steigenberger P, Dach R, Schmitz M, Dilssner F, Hugentobler U (2009) Interactions of the IGS reprocessing and the IGS antenna phase center model. EOS Transactions AGU, Fall Meeting Supplement, Abstract G13A-01, Vol. 90, Nr. 52, AGU, 2009. http://acc.igs.org/repro1/sat-ant-pcv_agu09.pdf
Schöne T, Schön N, Thaller D (2009) IGS Tide Gauge Benchmark Monitoring Pilot Project (TIGA): scientific benefits. J Geod 83:249–261. https://doi.org/10.1007/s00190-008-0269-y (with more info online at http://kb.igs.org/hc/en-us/sections/200423353-Technical-Information)
Sośnica K, Thaller D, Dach R, Steigenberer P, Beutler G, Arnold D, Jaggi A (2015) Satellite laser ranging to GPS and GLONASS. J Geod 89(7):725–743. https://doi.org/10.1007/s00190-015-0810-8
Sośnica K, Prange L, Kaźmierski K, Bury G, Drożdżewski M, Zajdel R, Hadas T (2017) Validation of Galileo orbits using SLR with a focus on satellites launched into incorrect orbital planes. J Geod 92(2):131–148. https://doi.org/10.1007/s00190-017-1050-x
Springer (2017) Personal communication
Springer TA, Beutler G, Rothacher M (1998) A new solar radiation pressure model for the GPS satellites. In: IGS workshop proceedings, ESOC, Darmstadt, Germany, 9–11 February 1998. www.iapg.bgu.tum.de/mediadb/10128/10129/ORBMOD98.ps
Springer T, Otten M, Flohrer C, Pereira F, Gini F, Enderle W (2014) GNSS satellite orbit modeling at ESOC. 2014 IGS Analysis Center Workshop Presentation, Pasadena, California
Urschl C, Gurtner W, Hugentobler U, Schaer S, Beutler G (2005) Validation of GNSS orbits using SLR observations. Adv Space Res 36(3):412–417. https://doi.org/10.1016/j.asr.2005.03.021
Wahr J, Nerem RS, Bettadpur SV (2015) The pole tide and its effect on GRACE time-variable gravity measurements: implications for estimates of surface mass variations. J Geophys Res 102(6):4597–4615. https://doi.org/10.1002/2015JB011986
Williams SDP (2003) Offsets in Global Positioning System time series. J Geophys Res. https://doi.org/10.1029/2002JB002156
Ziebart M, Sibthorpe A, Cross P, Bar-Sever Y, Haines B (2007) Cracking the GPS-SLR orbit anomaly. In: Proceedings of ION-GNSS-2007, pp 2033–2038
Zhang J, Bock Y, Johnson H, Fang P, Williams S, Genrich J, Wdowinski S, Behr J (1997) Southern California Permanent GPS geodetic array: error analysis of daily position estimates and site velocities. J Geophys Res 102(B8):18035–18055. https://doi.org/10.1029/97JB01380
Acknowledgements
Support for this work was provided by the NASA ROSES-2016 program under NRA number NNH16AC75I. Thanks to all participating Analysis Centers for their efforts in contributing to the second IGS reprocessing campaign. Thanks to Jim Ray and Paul Rebischung for their advice in forming the repro2 orbit and clock combinations. They also reviewed an early draft of this paper, and their comments/suggestions improved it significantly. Thanks to Nacho Romero for a third early review. Thanks to three anonymous reviewers for their feedback, which helped to improve this paper. Thanks to Kevin Choi for helping to verify IGS software performance on NRL computers.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Griffiths, J. Combined orbits and clocks from IGS second reprocessing. J Geod 93, 177–195 (2019). https://doi.org/10.1007/s00190-018-1149-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00190-018-1149-8