Evaluation of earth rotation parameters from modernized GNSS navigation messages

Modernized navigation messages of global navigation satellite systems like GPS CNAV include earth rotation parameters (ERPs), namely the pole coordinates and UT1-UTC (∆UT1) as well as their rates. Broadcast ERPs are primarily needed for space-borne GNSS applications that require transformations between earth-fixed and inertial reference frames like navigation in earth orbit as well as to the moon. Based on a global tracking network of 23 stations, broadcast ERP values are obtained for the global systems GPS and BeiDou as well as the regional QZSS and IRNSS. Subsequent data sets at daily intervals show polar motion discontinuities of 0.4 to 0.7 mas for GPS, QZSS, and IRNSS, whereas BDS is worse by a factor of about two. Discontinuities in ∆UT1 range from 0.17 to 0.45 ms. External comparison with the C04 series of the International Earth Rotation and Reference Systems Service results in polar motion RMS differences of 0.3 to 1.0 mas and ∆UT1 differences of about 0.13 ms for GPS, QZSS, and IRNSS. Due to less frequent update intervals, BDS performs worse by a factor of 2 – 4. In view of the current GNSS-based positioning errors at geostationary or even lunar distances, the accuracy of GPS, QZSS, and IRNSS broadcast ERPs is sufficient to support autonomous spacecraft navigation without the need for external data.


Introduction
The navigation messages of Global Navigation Satellite Systems (GNSSs) provide essential information such as satellite orbit and clock parameters needed for GNSS-based positioning. Modernized navigation messages like the civil navigation message (CNAV) include additional data such as earth rotation parameters (ERPs). The ERPs comprise the pole coordinates x p and y p , and the difference ∆UT1 between Universal Time 1 (UT1) and Universal Time Coordinated (UTC). The pole coordinates describe the location of the earth rotation axis on the earth's crust. ∆UT1 is based on the difference between earth rotation (UT1) and an atomic time scale (UTC).
The ERPs are part of the full set of earth orientation parameters (EOPs) which, in addition, include precession and nutation angles or celestial pole coordinates, and describe the overall transformation between the earth-centered earth-fixed (ECEF) and the earth-centered inertial (ECI) reference frame (Torge and Müller 2012;Capitaine 1986). Only the ERPs are included in the GNSS navigation message as nutation and precession can be described by analytical models with sufficient accuracy over long time scales. Therefore, only these parameters will be considered in the following.
For completeness, it should also be mentioned that the Interface Control Documents (ICDs) of the various GNSSs consistently use the term "EOP" when describing the format and usage of the corresponding parameters, even though the designation "ERP" would actually be more appropriate and consistent with established geodetic terminology. In fact, all currently available broadcast navigation messages are confined to polar motion and ∆UT1 data and do not provide additional parameters for the full ECEF/ECI transformation.
ERPs are of particular importance for space-borne GNSS applications as orbit computations and predictions are usually performed in an ECI frame whereas the GNSS satellite orbits are given in an ECEF frame. These applications include low earth orbit (LEO) satellites for, e.g., earth observation, navigation of geostationary satellites to their designated orbital slot as well as navigation to and around the moon. ERPs are also required for the numerical propagation of GNSS satellite orbits over extended periods of time. If the ERPs are included in the navigation message, such a propagation can be performed autonomously inside a GNSS receiver without the need for external data (Bar-Sever and Bertiger 2012; Montenbruck et al. 2020).
ERP data are currently transmitted by the United States Global Positioning System (GPS), the Chinese BeiDou Navigation Satellite System (BDS), the Japanese Quasi-Zenith Satellite System (QZSS), and the Indian Regional Navigation Satellite System (IRNSS). For our study, ERP data were decoded from raw navigation frames provided by 23 globally distributed receivers from two different manufacturers, see Fig. 1, ensuring continuous coverage. The stations are part of the Cooperative Network for GNSS Observations (CONGO, Montenbruck et al. 2011) operated by the German Aerospace Center (DLR) and the German Federal Agency for Cartography and Geodesy (BKG) as well as additional stations contributed by Deutsches GeoForschungsZentrum (GFZ, Ramatschi et al. 2019). For further processing, the extracted ERP data were stored in a draft version of the Receiver Independent Exchange format 4.00 (RINEX, IGS RWG and RTCM 2021). Our analysis covers the time period January 2020 until May 2021.
First, an overview of reference frame transformations and ERPs is given. The broadcast ERP data, including the transmission schemes of the different systems, are discussed in the following section. Finally, the quality of the broadcast ERPs is analyzed by misclosures of consecutive data sets as well as comparisons with an external reference series.

Reference frame transformations and earth rotation parameters
As a background for better understanding, this section describes the established models for transformations between celestial and terrestrial reference frames and highlights the role of ERPs in this transformation.

Reference frame transformations
The transformation from the ECEF to the ECI frame is conventionally composed of a sequence of three rotations (Petit and Luzum 2010): with the earth-centered inertial position vector r ECI and the earth-centered earth-fixed position vector r ECEF . The matrix Q describes essentially the combined effect of precession and nutation. It depends on the coordinates X and Y of the celestial intermediate pole (CIP) in the celestial intermediate reference system (CIRS), which are described by the IAU2006/IAU2000A model (Wallace and Capitaine 2006). The matrix R is based on the earth rotation angle (ERA), which, by definition, depends on UT1 in a strictly linear manner.
The transformation matrix for polar motion W is composed of rotation matrices R i around the axis i = 1, 2, 3 computed from the pole coordinates x p and y p and the locator s ′ of the terrestrial intermediate origin (TIO) The legacy reference frame transformation recommended in the IERS Conventions 1996(McCarthy 1996 is based on the nutation matrix N (IAU1980, Seidelmann 1982), the precession matrix P (IAU1976, Lieske et al. 1977) and Greenwich Apparent Sidereal Time (GAST) Other than the ERA, GAST depends on UT1 in a nonlinear manner and involves nutation-related periodic contributions (McCarthy 1996). Whereas this transformation is outdated for high precision applications, it still provides sufficient accuracy when using broadcast ERPs. Compared to the full IERS2010 models, it offers higher efficiency and lower computational effort, in particular, if the number of nutation model coefficients is further reduced (Montenbruck et al. 2020).

Earth rotation parameters
Variations in earth rotation are primarily driven by tidal forces of sun and moon as well as mass transport in the earth system. They are dominated by the typical tidal periods of 14 days and one month as well as variations at semiannual, annual, and Chandler periods. The Chandler period is about 435 days and caused by the free oscillation of the earth (Dehant and Mathews 2015 in the absence of leap second changes. ERPs are typically provided as daily estimates of offset and drift at noon or offset values at midnight. However, earth rotation is also affected by subdaily effects primarily introduced by ocean tides. These variations at diurnal and semidiurnal periods can reach up to 1 mas for polar motion and 0.1 ms for ΔUT1, see Fig. 2. The largest constituents for polar motion are the lunar diurnal tide K 1 with a period of 23.93 h and the principal lunar semidiurnal tide M 2 with a period of 12.42 h. For high-precision applications, a model for these subdaily effects has to be added to the daily ERP estimates, which, by convention, are mean values free of subdaily tide contributions. The model recommended by the International Earth Rotation and Reference Systems Service (IERS) can be found in the recent IERS conventions (Petit and Luzum 2010). The complete IERS2010 model consists of 41 diurnal and 30 semidiurnal terms. The superposition of different periods introduces beat periods, the most prominent ones with 14 days, half a year and a year. The beat period of about 14 days can clearly be seen in Fig. 2. If ERPs are estimated at daily intervals, the subdaily ERP variations have to be corrected in the observation model to avoid aliasing effects. As a consequence, the estimated ERPs are free of the subdaily effects, and they have to be added by the user after interpolating or propagating the tabulated ERPs to the desired epoch.
In addition to these subdaily variations, ΔUT1 and ΔLOD are affected by tidal deformations of the earth (4) caused by sun and moon with periods between 5 days and 18.6 years (zonal tides, Yoder et al. 1981). Whereas the linear ERP estimation model sufficiently approximates long-periodic variations, variations with periods up to roughly one month should be considered for best interpolation or extrapolation of tabular ERP data based on a daily epoch grid. Regularized versions of ΔUT1 and ΔLOD in which zonal tides with periods between 5 and 35 days as given in Table 8.1 of Petit and Luzum (2010) have been removed are designated as ΔUT1R and ΔLODR, respectively. A one-year time series of correction terms UT1 reg = ΔUT1 − ΔUT1R and LOD reg = ΔLOD − ΔLODR is given in Fig. 3. This plot is dominated by the largest terms of the 41 constituents, namely the lunisolar fortnightly tide M f with a period of 13.66 days and the lunar monthly tide M m with a period of 27.56 days.
In the determination of ∆UT1, the periodic part UT1 reg is removed from the observed ΔUT1 to get the smoother quantity ΔUT1R to improve the accuracy of linear inter-/ extrapolation. After the parameter estimation, UT1 reg is added back to obtain ΔUT1. The same procedure has to be applied when propagating ∆UT1 from the data reference epoch t ERP to an arbitrary epoch t. As illustrated in Fig. 4, ∆UT1 and ΔLOD given at t ERP are first reduced to ∆UT1R and ΔLODR. Then ∆UT1R is propagated to t with ΔLODR and UT1 reg at t is added back to obtain ∆UT1 at t. Finally, the subdaily corrections for polar motion and ΔUT1 have to be added to get x p , y p , and ∆UT1 at epoch t.
An overview about the characteristic errors of different commonly used ERP series is given in Table 1. Polar motion and ∆UT1 are given in units of milliarcseconds (mas) and milliseconds (ms) corresponding to about 3 and 46 cm on the earth's surface, respectively. Bulletin A (Wooden and Gambis 2004) is a multi-technique combined EOP series derived from GNSS, SLR, and VLBI observations as well as atmospheric angular momentum (AAM) data obtained from meteorological analysis and forecasts. It includes one year of predictions and is updated on a weekly basis.
IERS C04 (Bizouard et al. 2018) is also a combined multi-technique ERP series but has a higher latency than Bulletin A and does not include predictions. It is based on the combined ERP series of the International GNSS Service (IGS, Johnston et al. 2017), the International VLBI Service (IVS, Nothnagel et al. 2016), the International Laser Ranging Service (ILRS, Pearlman et al. 2019), and the International DORIS Service (IDS, Willis et al. 2015). In addition, input data of individual VLBI analysis centers (ACs) and the ERP data of the International Terrestrial Reference Frame (ITRF) combination are considered. Due to the latencies of the input series, IERS C04 is provided with a delay of one month after the last data point.
The IGS provides ERP series with different latencies obtained from a combination of the results contributed by several ACs. The IGS ultra-rapid ERP series (IGS 2021a) includes one day of observed ERP data as well as one day of predicted ERPs. As GNSSs are not able to determine ∆UT1, the ultra-rapid ∆UT1 values are derived from Bulletin A. The IGS final combined ERP series (IGS 2021b) has a latency of about two weeks and contains only values obtained from observations. The quality of the IGS final polar motion estimates as well as their error budget, is assessed in detail by Ray et al. (2017).

ERP data in GNSS navigation messages
ERPs are included in the GPS/QZSS CNAV and CNAV-2, BeiDou CNAV-1, CNAV-2 and CNAV-3, as well as IRNSS NAV messages, see Table 2. For GLONASS, transmission of ERPs is foreseen for the navigation messages of future L1, L2, and L3 code division multiple access (CDMA) signals, but no such data are broadcast at present. For Galileo, no ERP data are included in the current Open Service signals and no plans for future transmissions have been announced so far. Format descriptions and user guidelines are given in the ICDs of the individual systems, also listed in Table 2. Identical ERP parameters and formats are defined across the various navigation message types within each constellation. Furthermore, identical ERP values are provided in the different messages concurrently transmitted by an individual satellite. As such, no distinction is made between ERPs of individual messages   (Steigenberger et al. 2015). GPS CNAV is also transmitted by the latest generation of GPS III satellites launched since December 2018. It is emphasized, though, that the transmission of GPS CNAV is currently declared as pre-operational and "and should not be used for safety-of-life or other critical purposes until the government declares them operational." (https:// www. gps. gov/ syste ms/ gps/ moder nizat ion/ cnav/).
According to IS-GPS-200K (2019), GPS CNAV polar motion and polar motion rates already include the diurnal and semidiurnal tides, ∆UT1 and ∆ U T1 , in addition, include contributions of zonal tides. The same is stated for QZSS and IRNSS (see references in Table 2), whereas no information on this topic is given for BDS.
User algorithms as provided in the various ICDs specify that the pole coordinates and ∆UT1 at an epoch t in the vicinity of the reference epoch t ERP are computed as where both t and t ERP are referred to as the constellationspecific system time. None of the constellations foresees an explicit consideration of periodic corrections to this linear model, which naturally limits the quality of this model. Updated GPS ERP data sets are currently transmitted twice per day. The earliest transmit time t tm of the first set is between midnight and 03:00, but most records are transmitted before 02:00. For the second data set, the transmit time is between 22:15 and midnight. The reference epoch of the ERPs t ERP varies between 16:38-17:04 and 19:50-19:57 in a regular pattern and is typically 62 to 70 h ahead of the first transmit time. t ERP of the morning ERP values differs by roughly one day from the evening ERP values. However, the polar motion values of these two epochs agree within ± 1 µas when propagated to the same epoch with the corresponding polar motion rates (Eq. 5).
The ∆UT1 values agree within ± 50 ns. Furthermore, the polar motion rates and the ΔUT1 values of the morning and evening epochs are identical, also indicating that the ERP data of the evening transmission are based on prestored onboard data and do not reflect new information from an intermediate ground upload. Therefore, only the morning ERP values are considered in the following analysis.
For GPS, only ERP data after June 15, 2020, are considered in the present analysis due to inconsistent ERP data transmission schemes applied across the GPS constellation prior to that epoch. In early 2020, GPS III satellites already used a transmission scheme similar to the current one with t ERP 40 to 70 h ahead of transmit time. For Block IIR-M and IIF satellites, t ERP and t tm differed in general only by 10 to 30 min, but few records also reached differences of up to 10 h.
For completeness, we note that the specification of ERP data in GPS CNAV and CNAV-2 message was revised in mid 2020. Starting with IS-GPS-200L (2020), IS-GPS-705G (2020), and IS-GPS-800G (2020) ∆UT1 = UT1 − UTC has been replaced by ∆UT1 GPS = UT1 − GPST. In contrast to UTC, GPS time (GPST) is a continuous time scale and avoids potential problems in the vicinity of leap second changes. However, this change has not become effective yet, and the transmitted ERP data are based on the earlier specifications for the entire analysis period of our study.
According to NGA (2020), the National Geospatial-Intelligence Agency (NGA) provides ERPs to the GPS control segment that are based on data of the United States Naval Observatory (USNO). Based on the comparison with the broadcast ERPs, it can be concluded that the NGA data provide the current source of CNAV ERPs. The NGA parameters and hence the transmitted CNAV data are updated on a daily basis. The NGA predictions for polar motion and ∆UT1 at epoch t are described by a simple model with offset, drift, and harmonic terms: with and the polar motion index P: x p or y p , the ∆UT1 index U, the offset coefficient O, the drift coefficient D, the reference epoch for polar motion t p , the reference epoch for ∆UT1 t u , the semiannual period p 1 (182.625 d), the annual period p 2 (365.25 d), the Chandler period p 3 (435 d), the annual sine/ cosine polar motion terms S P,2 and C P,2 , the Chandler sine/ cosine polar motion terms S P,3 and C P,3 , the semiannual sine/ cosine ∆UT1 terms S U,1 and C U,1 , fixed to S U,1 = − 0.006 s and C U,1 = − 0.007 s, the annual sine/cosine ∆UT1 terms S U,2 and C U,2 , fixed to S U,2 = − 0.022 s and C U,2 = 0.012 s. For ∆UT1, the harmonic terms are fixed to the values given above, i.e., only offset and drift are estimated model parameters, whereas six independent parameters represent polar motion. The time interval of observations contributing to the estimation of this model is unknown. NGA publishes the model coefficients as well as tabulated ERP predictions based on Eq. 7/8 and includes zonal tides and subdaily terms at its ftp server (see data availability statement) in two different versions. The predicted ERPs cover seven days and 180 days, respectively. The latter interval might be driven by the requirement to maintain navigation function for 180 days without ground contact (AUTONAV, Menn and Bernstein 1994). As the ERPs published by NGA do not contain rates, these seem to be derived from the tabulated predictions. Based on inspection of the respective data, it may be inferred that a difference quotient covering the first three days of predicted NGA ERPs is used to obtain the CNAV ERP rates.
(7) P(t) = O P + D P ⋅ Δt p + 3 ∑ i=2 S P,i ⋅ sin P,i + C P,i ⋅ cos P,i (8) New BDS ERP records are transmitted once per day with updates at about 00:35 or 01:00. The immediately preceding midnight epoch is usually adopted as t ERP . The comparison with IGS or IERS ERP series suggests that the BDS broadcast ERPs are also derived from a harmonic prediction model as used by NGA. Unlike GPS, which uses new predictions every day, the BDS prediction model is updated irregularly at intervals of 4 -20 days. This resulted in larger errors and pronounced discontinuities of the transmitted ERPs at the model update epochs.
In contrast to the other systems, new QZSS ERP records are transmitted once every hour. The first transmission usually takes place at 3 min after the full hour. t ERP refers to the next full hour, i.e., 57 min after t tm . According to IS-QZSS-PNT-004 (2021), the QZSS ERPs have a fit interval of two hours and a validity interval of 12 h.
New IRNSS ERP records are transmitted once per day with updates between midnight and 0:03. t ERP always refers to the corresponding midnight epoch. Between November 20, 2020, and January 21, 2021, no updates of the underlying predictions were performed for polar motion data transmitted by IRNSS. Accordingly, the pole misclosures that will be introduced later are zero during that time period. Therefore, this time interval is excluded from the further analysis of the IRNSS ERPs. Unfortunately, no details about the origin of the ERP data and the parameterization of the ERP prediction model are provided for BDS, QZSS, and IRNSS.
In order to better understand the transmission scheme of the broadcast ERPs and to identify their zero age of data, the best agreement with the IERS C04 series (Bizouard et al. 2018) was identified by propagating the broadcast ERPs to different trial epochs t 0 . The resulting zero-age-of-data epochs provided in Table 3 are expressed by their offset dt from the midnight epoch preceding the epoch of first transmission of an ERP data set. In addition, it is evaluated if the application of the IERS2010 subdaily ERP model improves the RMS values or not. The broadcast ERP values are propagated from (9) t 0 = ⌊t tm ⌋ − dt t ERP to t 0 according to (5) and the subdaily ERPs are optionally removed at t 0 as the reference series is free of subdaily ERP variations. For BDS, only the first epoch of an updated ERP data set is used as the following epochs are strongly affected by prediction errors. These epochs have been determined by analyzing the discontinuities in polar motion and ∆UT1 time series w.r.t. a reference series as shown in Fig. 6b. As a consequence, only 11% of the available BDS ERP data are used for this analysis.
For GPS, significantly different values for dt are obtained for polar motion and ∆UT1. Whereas a shift of one day shows the lowest polar motion RMS values of about 0.15 mas, no shift is necessary for ∆UT1. Correction for subdaily polar motion improves the RMS values by a factor of more than two, whereas the improvement for ∆UT1 is only minor. Compared to GPS, the RMS values for BDS are larger by a factor of three for polar motion and a factor of two for ∆UT1. Like for GPS, a significant improvement could be achieved when applying the subdaily ERP model, thus demonstrating that the broadcast ERPs of both constellations include the subdaily tides at the zero age of data epoch, rather than those at the ERP epoch or the first transmit epoch.
The RMS values of QZSS are similar to those of GPS, but the best agreement is achieved without correcting for the subdaily ERP model. This is surprising as the spectra of the QZSS ERPs contain significant peaks at the typical subdaily periods. However, as the amplitudes of the subdaily terms included in the QZSS ERPs are smaller than the IERS model, the smallest RMS values are obtained without considering the subdaily ERP model. The latter is also true for IRNSS, but the RMS with respect to the reference series is worse by a factor of more than two for polar motion.
The ICDs of all considered systems state that subdaily corrections are included in the broadcast ERPs, but do not specify how exactly the subdaily tides are incorporated into the ERP values and rates at the ERP epoch. In fact, our analysis suggests that subdaily tides are applied in the ERP values at the zero age of data epoch for GPS and BDS, while IRNSS ERPs appear to be free of these corrections. In view of a linear propagation over several hours or even days, including periodic subdaily variations is questionable. For proper handling, the subdaily ERP model would have to be subtracted from the GPS and BDS ERPs at the zero age of data epoch and added back after propagation to the desired epoch with the ERP rates. However, as the errors of the broadcast ERPs at transmission time are dominated by the prediction errors of the ERPs, the error contribution of the subdaily ERPs is essentially negligible in the current context. The same is true for the linear propagation of ∆UT1, neglecting the reduction to UT1R. The latter introduces an RMS error of 10 µs with peak values of 50 µs.

Broadcast ERP quality
The precision of the broadcast ERPs is evaluated by misclosures of consecutive data sets, whereas comparisons with an external reference series assess the accuracy.

Internal consistency
The internal consistency of independent ERP estimates can be evaluated by ERP misclosures as illustrated in Fig. 5. The ERP values are evaluated at consecutive t tm epochs i and i + 1 with their rates, and the differences ∆ERP at these epochs allow to assess the internal consistency. t tm(i) is the epoch when a new set ERP(i) is transmitted the first time. The users apply these ERPs until transmission of the next data set ERP(i + 1) at t tm(i+1) illustrated by the shaded area in Fig. 5.
The RMS values of the ERP misclosures are given in Table 4. For GPS and IRNSS, the consistency is at roughly 0.5 mas for polar motion and 0.17 ms for ∆UT1. When looking at the hourly update rate of QZSS, the misclosures are smaller by a factor of ten for polar motion and seven for ∆UT1. However, when looking at daily intervals, polar   Page 9 of 12 50 motion is at a similar level as GPS and ∆UT1 is worse by a factor of more than two. Whereas the pattern of the misclosures is pretty random for GPS, QZSS, and IRNSS, the situation is different for BDS. Most of the misclosures are within ± 10 µas for polar motion (87%) and ± 10 µs for ∆UT1 (84%). The remaining misclosures reach up to 8 mas for polar motion and 3 ms for ∆UT1, thus causing the large misclosure values in Table 4. This behavior can be explained by the 4-20 day update interval of the underlying ERP forecast model for BDS, which results in very small discontinuities during transmission of the same model but large discontinuities when switching to the next parameter set.

Comparison with reference series
The IERS C04 ERP series consistent with ITRF2014 (Bizouard et al. 2018) is used as an external reference to evaluate the broadcast ERP quality. In the case of BeiDou, raw CNAV-1/2/3 ERP messages with void data, i.e., a null bit pattern, have been discarded in the data collection. Furthermore, occasional cases of erroneous data transmissions (from satellites C45/C46) resulting in grossly wrong ERP values have been identified and rejected in a coarse outlier screening. Unless otherwise noted, the comparison of C04 and broadcast ERPs was performed at first transmit time of the individual broadcast ERP records. At each epoch, subdaily ERP corrections were added to C04 ERP values according to IERS conventions, whereas no corrections were applied to the broadcast ERPs as specified in the constellation-specific ICDs (see Table 2). The ERP difference time series are shown in Fig. 6. None of the series shows a significant bias w.r.t. the C04 reference series, neither for polar motion nor for ∆UT1. The scatter of the GPS and QZSS ERP differences is quite random. GPS has maximum differences of ± 2 mas for polar motion and ± 0.5 ms for ∆UT. For BDS, steadily increasing differences with peak values of 7 mas and 3 ms are visible during time intervals of 4 -20 days due to less frequent ERP model updates already mentioned above. The anomalous ERP values of IRNSS in December 2020 and January 2021 are included in Fig. 6d. They reach a maximum of 10 mas for x p and ± 1.2 ms for ∆UT1 and are caused by missing updates of the IRNSS broadcast ERP model during that time period indicated by zero polar motion misclosures and ∆UT1 misclosures smaller than ± 10 µs.
Numerical values of the RMS differences of the broadcast ERPs w.r.t. the reference series are given in Table 5.
The RMS values are evaluated for three different epoch types: the transmission time and the next transmission time resulting in best and worst case values as well as for a time series with 15 min sampling to evaluate the average user performance. Due to the frequent updates, QZSS shows the smallest RMS values of about 0.4 mas for polar motion as well as the smallest differences for the three cases. The reason for the better agreement for a longer prediction time (t i+1 ) is unknown. GPS and IRNSS have a very similar performance that is worse by a factor of roughly two compared to QZSS. For ∆UT1, GPS, QZSS and IRNSS show the same level of RMS differences of 0.1 to 0.2 ms, but BDS is worse by a factor of three. The same comparisons were also done with the IGS final ERP series (IGS 2021b) but are not shown here as they are very similar to the C04 comparison differing by only up to 10 µas for polar motion and 20 µs for ∆UT1.
Compared to Table 3, the values in Table 4 are larger by a factor of up to eight, whereas the ∆UT1 differences are significantly smaller. Overall, the broadcast ERP performance experienced by GNSS users appears to be mainly driven by the time between the zero age of data epoch t 0 and the time of first transmission, t tm . This is particularly obvious for BDS, which uses less-frequent forecast model updates than other constellations. Furthermore, it may be noted that the forecast errors dominate the contributions of subdaily ERP variations and mask a possible neglect or ICD-related uncertainties in this correction.

Conclusions
Even though the availability of ERPs might be of interest for real-time precise point positioning (PPP) to properly model site displacement due to pole tides, their primary use is certainly in the field of space-borne navigation. Possible applications include navigation of satellites in low earth or geostationary orbit (Capuano et al. 2017), as well as navigation to the moon (Delépaut et al. 2020). Figure 7 illustrates the impact of ERP-induced errors on the transformation for ECEF to ECI for different orbit heights. Minimum and maximum errors evaluated at t tm from the analysis interval are given. With a single-frequency COTS GPS receiver connected to a signal simulator, Capuano et al. (2017) obtained a single epoch positioning error in GEO of about 17 m. As the ERP-induced position error at GEO is well below one meter for GPS, the broadcast ERP quality is certainly compatible with current performance needs in this application. The same is true for onboard orbit determination of the Chinese geostationary satellite TJS-2 presented by Wang et al. (2021). The GNSS receiver onboard TJS-2 is capable of tracking GPS, GLONASS and BDS-2 singlefrequency code observations. For the transformation to the inertial frame, predicted Bulletin-A ERPs are used. Wang et al. (2021) report a 3D position difference RMS value of 21.3 m between the onboard and a postprocessed solution using the same observation data but more sophisticated models as well as IGS final orbit and clock products. However, Capuano et al. (2017) could also demonstrate a positioning error of less than one meter for simulated Galileo dual-frequency observations. For this use case, the broadcast ERPs would significantly contribute to the position error in the inertial frame.
ERPs included in the broadcast messages of GPS, BDS, QZSS, and IRNSS allow for ECEF/ECI transformations and vice versa onboard GNSS-equipped spacecraft without the need for external data. For GPS, BDS, and QZSS, the current ERP accuracy appears sufficient in view of the overall positioning accuracy achievable in GNSS-based navigation. Due to the less frequent ERP model update intervals, the BeiDou ERPs show a reduced accuracy.
A conceptual problem in the use of broadcast ERPs results from the unclear, and most likely misleading, specifications for applications of subdaily tides in the ICDs of all GNSS presently providing ERP data as part of their navigation messages. While errors resulting from the associated uncertainties are generally smaller than the contribution of ERP forecast errors, it would be strongly advisable to transmit only mean ERP values and rates that do not include any subdaily tidal contributions. This would ensure consistency with other geodetic services, provide better transparency, and facilitate the handling of broadcast ERPs on the user side. With respect to ∆UT1, improved predictions could potentially be achieved when replacing UT1 with reduced UT1R. This is particularly true for GPS, where UT1 rates are based on a linear difference quotient approximation of predicted values over a 3-day interval and certainly affected by the zonal tides with periods of 5 days and up. However further analysis would be required to trade the possible accuracy gain against the increased complexity of evaluating the UT1 − UT1R model at the user end.
Finally, we recall that the European Galileo system does not, so far, provide ERPs as part of its INAV and FNAV navigation messages. This contrasts sharply with the benefit that Galileo tracking provides for space-borne users in view of its superior broadcast ephemeris accuracy (Montenbruck et al. 2022) and the upcoming introduction of a dedicated high-accuracy service (Fernandez-Hernandez et al. 2018). The addition of ERP data to the current navigation messages should therefore be investigated and is strongly encouraged.