Abstract
The International GNSS Service (IGS) Analysis Center Coordinator initiated in 2019 an experimental multi-GNSS orbit combination service by adapting the current combination software that has been used for many years for IGS GPS and GLONASS combinations. The multi-GNSS orbits are based on individual products generated by IGS and multi-GNSS Pilot Project analysis centers. However, the combinations are not yet considered to be the final products at this time. The goal of this research is to provide a quality assessment of the very first IGS experimental multi-GNSS combined orbits based on Satellite Laser Ranging (SLR) observations and the mean position errors from the orbit combinations. The errors available in the combined orbit files provide information about the consistency between orbits from different analysis centers, whereas SLR provides independent orbit validation results even for those satellites which are considered only by one analysis center, and thus, the quality of the combination is not provided in the orbit files. We found that the BeiDou-3 satellites manufactured by China Academy of Space Technology and Shanghai Engineering Center for Microsatellites are characterized by opposite SLR residual dependencies with respect to the position of the sun which means that the orbit models for BeiDou-3 need further improvement. Smallest SLR residuals are obtained for Galileo, GLONASS-K1, and GLONASS-M+ . However, the latter is characterized by a bias of + 29 mm. The mean standard deviations of SLR residuals are 23, 29, 87, 51, 40, and 72 mm for Galileo, GLONASS, BeiDou GEO, BeiDou IGSO, BeiDou MEO, and QZSS, respectively. The mean orbit combination errors in the radial direction are three times lower than those from SLR residuals in the case of MEO satellites and vary between 8 and 14 mm, whereas the orbit errors are four times lower than SLR residuals in the case of GEO and IGSO and equal to 11–21 mm.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
Introduction
The experimental multi-GNSS orbit combination service was initiated in 2019 by the International GNSS Service (IGS, Johnston et al. 2017) Analysis Centre Coordinator (ACC). The current combination software, used for many years for the IGS GPS and GLONASS combinations (Beutler et al. 1995; Kouba and Mireault, 1998; Weber and Springer 2001), has been employed for the generation of multi-GNSS orbits. The combined orbit products are not yet considered to be the final IGS products at this time, but instead, they can be used for the comparison of individual analysis centers contributing to IGS and the multi-GNSS Pilot Project (MGEX, Montenbruck et al. 2017). The combined products include a plethora of different satellite systems and generations; GPS: Block IIA, IIR, IIR-M, IIF, III; GLONASS: M, M+ , K1; Galileo: In-Orbit-Validation (IOV), Full Operational Capability (FOC), and FOC in eccentric orbits; BeiDou-2: Geostationary (GEO), Inclined Geosynchronous (IGSO), and Medium Earth Orbiters (MEOs); BeiDou-3 and BeiDou-3S: IGSO and MEO, manufactured by China Academy of Space Technology (CAST) and Shanghai Engineering Center for Microsatellites (SECM, Zhao et al. 2018); QZSS: GEO and eccentric IGSO.
All new GNSS satellites, except for GPS, are equipped with laser retroreflector arrays (LRA) for Satellite Laser Ranging (SLR). The International Laser Ranging Service (ILRS, Pearlman et al. 2019a, b) initiated in 2015 a series of intensive GNSS tracking campaigns, which resulted in a substantial increase in collected SLR observations to GNSS (Bury et al. 2019a), which was possible thanks to the optimized and enhanced tracking strategy at ILRS stations. SLR observations to GNSS can be used for the orbit modeling validation (Springer et al. 1999; Arnold et al. 2015; Kazmierski et al. 2018), precise orbit determination of GNSS satellites (Bury et al. 2019c), co-location of GNSS and SLR techniques in space (Thaller et al. 2011; Bruni et al. 2018), and deriving global geodetic parameters (Sośnica et al. 2018a, 2019).
The goal of this study is to assess the quality of the IGS experimental multi-GNSS combined orbits based on SLR observations and on the mean position errors from the orbit combinations. SLR provides independent orbit validation results even for those satellites, which are considered only by one analysis center, and the standard deviation values (STDs) of the combination are not provided. We also analyze STD values provided in the combined SP3 orbit files, which reflect the consistency level between orbits from different analysis centers.
Methodology
The experimental multi-GNSS products are made available on a weekly basis, with a delay of about 20 days, which is similar to the final IGS combination delay (Griffiths and Ray 2009; Johnston et al. 2017; Griffiths 2019). The combination includes final and rapid MGEX submissions, as well as final IGS submissions, all combined together because the purpose of the experimental combination is to include as many submissions as possible. Only the MGEX submission is employed in the case of analysis centers providing both the final IGS and MGEX submissions, whereas IGS products are used then only for comparison purposes. The same weighted average approach developed by Beutler et al. (1995) is used for generating combined orbits. Orbits are combined as weighted averages computed over all centers. Each center and satellite are given a position weight computed from the center’s absolute deviation to the unweighted average orbit. Correlations between components of the position vector and positions referring to different epochs, that is, the time correlations, are neglected. Under these simplifications, the combination procedure for each coordinate reads as follows:
First, a simple mean of the orbit solutions by different analysis centers is calculated:
where \( x_{k}\) is the orbit solution of the kth analysis center for each \(X, Y, Z\) orbital component at the epoch \(t\), and \(n\) is the total number of analysis centers providing orbit solution for a satellite.
For each orbit solution, seven-parameter Helmert transformation between the individual orbit and the mean orbit is estimated by minimizing the mean absolute deviation (MAD) of the position components of the two orbits as a robust function. This L1-norm minimization is performed using a bracketing and bisection algorithm for finding the root of the derivative of the robust function, as described by Press et al. (2007). Each orbit solution is transformed using the above estimated Helmert parameters, and the weights for the analysis centers are calculated based on the deviations of each orbit and its transformed version:
where \(d^{3D}_{k} \left( t \right) = \left| {x_{k}^{T} \left( t \right) - x_{k} \left( t \right)} \right|\), and MAD is the mean absolute deviation of the position components of the orbit solution \(x_{k}\) from the positions of the transformed orbit \(x_{k}^{T}\) for a particular day. By using these weights for the analysis center data, the weighted mean of the orbits is calculated:
Another set of the seven-parameter Helmert transformation is estimated between each orbit solution and the above weighted mean using the same L1-norm minimization as described for the previous step. Each orbit is then transformed using the estimated transformation parameters. Finally, the combined orbit is calculated as the weighted mean of the transformed orbit positions \(x_{k}^{T}\), with the weights as calculated before:
In the case of very large RMS values for a satellite in the combined orbit exceeding triple STD, that satellite is removed from the orbit solution, and in case of a very large RMS for an orbit solution, or very large transformation parameters estimated for an orbit solution, that orbit solution is completely removed from the combination. In such cases, the above algorithm is repeated until none of the above issues persist in the combination.
Position components \(x_{k} \) represent the individual orbit solutions in the earth-fixed frame. Individual orbit components (X, Y, Z) could be considered independently. However, in this approach, the weights are calculated based on differences of full position 3D vectors. According to Beutler et al. (1995), the weighted mean of the orbits satisfies the equation of motion provided that the weights are constant.
The weighted average software was originally developed by Timon Springer and Gerhard Beutler in 1993 and used by IGS for the combination of GPS and later for GLONASS orbits (Beutler et al. 1995; Weber and Springer 2001; Kouba 2009). The IGS ACC initiated a web service providing weighted RMS (WRMS) of the individual orbit solutions with respect to the combined orbits and estimated transformation parameters between individual solutions and the IGS combined orbits. The service also generates plots showing the time series of transformation parameters, see: https://acc.igs.org/mgex_experimental.html. The mean WRMS of individual contributors in the combination is 10, 25, 14, 52, and 46 mm for GPS, GLONASS, Galileo, BeiDou, and QZSS, respectively. Transformations to the IGS reference frame are not yet applied due to the variety of the products used for the combinations. Not all analysis centers provide earth rotation parameters and products in SINEX files (Mansur et al. 2020). The errors resulting from this neglect are much smaller than the errors in the orbit determination process. For most of the centers, the rotation errors do not exceed 0.1 mas, which corresponds to 3 mm on the earth surface and 12 mm at GNSS heights. The weighting scheme used in this experimental combination is the same as the one used in the current operational combination. An improved weighting that considers the differences between different orbit types is being developed in an upgraded combination software, which aims at a fully operational multi-GNSS combination. So far, only GNSS orbits are combined, which satisfies the needs of users processing double-differenced GNSS observations, whereas a proper GNSS clock combination is still pending for undifferenced solutions. In this research, the analysis period spans from April 29 to September 29, 2019.
GPS
Combined GPS orbits are generated on the basis of twelve analysis centers, six of which emerge from the MGEX solutions with the suffix -M, and another six from IGS final, see Table 1. Two IGS solutions from CODE and GFZ and the combined IGS rapid (IGR) and final (IGS) solutions are included for statistical reasons in the combination reports but not included in the combination solutions due to their redundancy. The current agreement between experimental IGS combination and IGR or IGS products is at the level of 3 and 4 mm, respectively, whereas most of the individual centers have an RMS of about 10 mm when compared to the combined orbits. So far, only two GPS satellites from Block IIA, SVN 35 and 36, have been equipped with LRA for SLR. These satellites were deactivated in 2013 and 2014 and recently are occasionally being reactivated again providing very sparse SLR data, which is why the SLR analysis for GPS is omitted here. However, it is planned that future GPS III satellites, launched after 2026, will carry LRAs for SLR tracking.
GLONASS
The combined GLONASS solutions are based on six MGEX centers, two IGS centers (denoted with a suffix ‘-X’), and two GLONASS-only contributors: IAC, and the SLR-only solution provided by MCC. The MCC contribution and IGS GLONASS final products (IGL) are not considered for the combination, but only for comparison purposes. The GLONASS constellation is mainly composed of GLONASS-M satellites broadcasting navigation signal on two frequencies L1 and L2, three M+ and two GLONASS-K1 with the additional L3 signal (Montenbruck et al. 2015, 2017). GLONASS-M and M+ are equipped with rectangular LRAs with 112 corner cubes, whereas K1s are equipped with a ring retroreflector surrounding the microwave transmitter antennas consisting of 123 corner cubes. Despite a large area of the ring retroreflector, the precision of the SLR measurements should be better for single-photon detectors because the probability of the reflection from each corner cube is the same. Therefore, when an SLR station collects hundreds or thousands of single full-rate reflections and generates one normal point based on 300 s of observations, the mean SLR observation corresponds to the centroid of the LRA onboard the GLONASS-K1 (Sośnica et al. 2015; Rodríguez et al. 2019). The future GLONASS-K2 will be equipped with ring LRAs with the number of corner cubes reduced to 36. GLONASS-K R26 (SVN 801) was the first, purely experimental satellite of the new type that was launched in 2011. This satellite cannot be tracked by all stations and is not considered by all analysis centers. The official status of SVN 801 is ‘flight test.’ GLONASS-K R09 (SVN 802), launched in 2014, is a newer K1-type satellite that is fully operational and transmits a signal on a channel that is accessible to all receivers.
Two GLONASS-M R02 (SVN 747) and R15 (SVN 857), one M + R05 (SVN 856), and K1 R09 (SVN 802) are intensively being tracked by the ILRS stations with the highest priority since July 2019, whereas the tracking of the remaining constellation depends on the capabilities and time availability of SLR stations. Before that, the ILRS recommended tracking of R02, R09, R11, and R14 (SVNs 747, 802, 853, and 732) and with a lower priority also R03, R18, R21, and R17 (SVNs 744, 854, 855, and 851).
Galileo
The Galileo constellation consists of three active IOV satellites, one IOV transmitting signal only on one frequency and thus not included in the combination, 19 active FOC satellites, and two FOC satellites launched into eccentric orbits that are used for geodesy and general relativity studies (Steigenberger and Montenbruck 2017; Hadas et al. 2019; Delva et al. 2015). Galileo IOV: E12 (SVN 102), and FOC: E01 (SVN 210), E08 (SVN 208), as well as FOC eccentric: E14 (SVN 202), were selected by the ILRS for intensive tracking. Before July 2019, the ILRS recommended tracking E12, E14, E09, E01, E36, E13, E15, and E04 (SVNs 102, 202, 209, 210, 219, 220, 221, and 213) with decreasing priorities. Galileo IOV satellites are equipped with rectangular retroreflectors with 84 corner cubes of a diameter of 33 mm, whereas FOCs are equipped with 60 corner cubes with an optical diameter of just 28.2 mm. Therefore, SLR tracking of FOC satellites is much more challenging than the tracking of IOV for SLR stations. However, the signature effect and the RMS of SLR normal points for FOC are reduced, especially for those SLR stations that use multi-photon detectors (Sośnica et al. 2018b).
BeiDou
In October 2019, the BeiDou constellation consisted of 34 active satellites and nine satellites not included in the operational mode (4 experimental BeiDou-3S and 5 BeiDou-3). BeiDou-2 includes five GEO, seven IGSO, and three MEO satellites, whereas BeiDou-3 comprises 18 MEO and one IGSO active spacecraft manufactured by SECM and CAST (Lv et al. 2020). The characteristics of the Chinese constellation can be followed at the MGEX Web site (https://mgex.igs.org/IGS_MGEX_Status_BDS.php). BeiDou-2 IGSO and MEO combined orbit solutions are based on five centers (see Table 1). The BeiDou-2 GEO orbits are based on GFM and WUM, because other centers do not consider GEO satellites or, as in the case of TUM, provide incomplete coverage of BeiDou-2 GEO satellites with the differences of mean orbits exceeding several meters. BeiDou-3 MEO orbits are generated only by WUM. The ILRS selected for intensive tracking four BeiDou-3: 3-M2 (C20, SVN 202), 3-M3 (C21, SVN 203), 3-M9 (C29, SVN 207), and 3-M10 (C30, SVN 208). Before July 2019, also BeiDou-2: C08, C10, C11, C13 (SVNs 008, 010, 012, 017), were scheduled for SLR tracking. BeiDou-2 MEO and BeiDou-3 SECM are equipped with LRA consisting of 42 corner cubes, whereas GEO and IGSO with 90 corner cubes, all of which have corner cube diameters of 33 mm. Two BeiDou-3 CAST from the ILRS priority list are equipped with retroreflectors consisting of 38 corner cubes.
QZSS
QZSS combined orbits are based on five MGEX centers; however, COM does not include solutions for the QZSS GEO satellite (J07, SVN 003), whereas JAM provides solutions only for the very first QZS-1 IGSO spacecraft (J01, SVN 001). All QZSS satellites are equipped with LRAs of 56-corner cubes of 41 mm in diameter and are scheduled in the ILRS priority list for tracking.
Internal evaluation of combined orbits
Figure 1 shows the number of GNSS satellites included in the combined IGS products. GPS, Galileo, and GLONASS are on a stable level of 32, 24, and 23 active spacecrafts. The number of BeiDou satellites in the combined products varies between 10 and 33 because BeiDou-3 data are currently being processed only by one center, which is subject to outages in case of product delays or processing errors. The number of considered QZSS satellites varies between 1 and 4.
The combined orbits in SP3 format include not only the precise positions and clocks of the satellites but also information on the quality of the combination, i.e., STD of weighted position residuals from the combination, for those satellite orbits which are generated by more than one analysis center. However, the STD from SP3 files should be interpreted as a function of the consistency between orbits from different centers rather than the absolute orbit accuracy. Most of the analysis centers employ similar orbit models, e.g., the Empirical CODE Orbit Model (ECOM, Beutler et al. 1994), or an extended version with twice-per-revolution parameters included, the so-called ECOM2 (Arnold et al. 2015). Some of the centers employ, in addition, a priori box-wing models (Springer et al. 2019). Therefore, the consistency level of the orbit solutions employing the same models may be at a high level, whereas all solutions may be affected by the same systematic errors. The systematic errors in orbit positions can thus better be assessed using independent observational techniques, such as SLR.
Figure 2 shows the orbital STD values from the SP3 files for 3D positions and the radial component, which has a fundamental meaning for positioning due to the extensive impact on the signal-in-space range error (SISRE, Montenbruck et al. 2018). The STD values in radial \(m_{r}\) are calculated using the variance–covariance propagation law as:
where the radial direction is calculated as \({\varvec{e}}_{r} = \frac{{{\varvec{r}}_{{\varvec{s}}} }}{{r_{s} }}\), with \({\varvec{r}}_{{\varvec{s}}}\) being the satellite position vector from SP3, and \(C_{S}\) is the diagonal variance matrix containing the squares of STD of the X, Y, and Z position components from SP3.
The radial precision for GPS satellites is at the level of 5 mm for all satellites, except for G04 with STD of 11 mm. G04 was occupied by the new GPS III (SVN74) satellite between January and July 2019 and again after October 2019. Apparently, the orbit models and antenna offsets used at different IGS analysis centers might be inconsistent or improper for the new GPS III satellite.
The radial precision for Galileo satellites is 8 mm, with no prominent difference between IOV, FOC, and FOC in eccentric orbits. For GLONASS, the radial precision is 9 mm for all satellites and 11 mm for R01 (SVN 730), R19 (SVN 720), and the experimental GLONASS-K1 R26 (801). The smallest radial STD of 7.6 mm is for R09 (SVN 802)—the second experimental GLONASS-K1 satellite.
For BeiDou-2, two groups can be distinguished: GEO with the mean radial STD of 21 mm, and IGSO and MEO with mean STD of 14 mm. The radial STD of QZSS IGSO is equal to 15 mm and thus corresponds well to other IGSO satellites. The 3D mean STD values are 37, 24, 15, 9, 26, and 16 mm for BeiDou GEO, BeiDou IGSO&MEO, Galileo, GPS, QZSS, and GLONASS, respectively.
Validation of IGS experimental combined multi-GNSS orbits using SLR
The GLONASS system is supported by the Russian network of SLR stations established in Asia, Europe, Africa, and South America, which results in a substantial number of SLR observations to selected GLONASS satellites (Fig. 3). The Galileo constellation is equally tracked by the European SLR stations, whereas other SLR stations track the Galileo satellites included in the ILRS priority list. GEO satellites are very challenging targets for laser ranging and have limited visibilities from SLR stations. Thus, the number of SLR observations to C01 and J07 is very low. One BeiDou-2 MEO (C11), three IGSO (C08, C10, and C13), and four BeiDou-3 MEO (CAST C20, C21; SECM C29, C30) were tracked by SLR in the period considered in this research.
The precision of the SLR measurements to GNSS strongly depends on the laser pulse width, timer and detector type used at the ILRS station, the number of collected photons, LRA size and corner cube arrangement, and data screening procedures used to generate SLR normal points. For satellites observed in zenith, the precision is 4–8 mm and decreases to 20–30 mm for satellites at low elevation angles due to signature effect of flat LRAs (https://ilrs.cddis.eosdis.nasa.gov/missions/satellite_missions/current_missions/ga02_stadata.html). The overall precision of SLR normal points to GNSS is between 10 and 20 mm for most of the stations and 4 mm for Matera due to a different screening procedure employed.
The SLR residuals are differences between SLR measurements and the theoretical distances to satellites based on IGS orbit positions. The SLR measurements are corrected by tropospheric delay models (Mendes and Pavlis 2004), effects emerging from general relativity, satellite eccentricities, and LRA offset corrections provided by mission operators. The processing strategy is the same as typically employed at the ILRS Associated Analysis Center at UPWr. (Zajdel et al. 2017; Otsubo et al. 2019) using the modified version of Bernese GNSS Software (Dach et al. 2015). The Bernese GNSS Software v.5.2 has the full capability of processing GPS and GLONASS observations, whereas Galileo processing is possible with some limitations. The modified version allows for processing Galileo, BeiDou, and QZSS data with a large number of GNSS satellites as used in the service www.govus.pl. The ILRS version of the International Terrestrial Reference Frame, SLRF2014, the ILRS eccentricity and data handling files, and SLR observations have been obtained from the ILRS Data Centers (Noll et al. 2019).
Figure 4 shows the results of the SLR validation for individual satellites. The smallest STD of residuals is 20, 24, and 24 mm for Galileo FOC eccentric, IOV, and FOC, respectively. For GLONASS, the STD of SLR residuals is 29, 27, 22, and 31 mm for GLONASS-M, M + , K1 R09, and the experimental K1 R26. For BeiDou-2, the STD is 87, 51, and 29 mm for GEO, IGSO, and MEO satellites, whereas for BeiDou-3, STD is 41 and 37 mm for CAST and SECM, respectively. QZSS has the STD of 59, 74, and 81 mm for QZS-1, other QZSS IGSO, and GEO, respectively.
The mean SLR offset is below 2 mm for Galileo FOC, FOC eccentric, and GLONASS-K1 R09. These satellites are equipped with either small-size LRAs or LRAs with the ring-shape corner cube arrangement so that they minimize the SLR signature effect, which is a systematic effect emerging from the laser reflection from multiple corner cubes. The large positive offsets of SLR residuals of 26 and 29 mm is visible for M + R05, R12, R21, and M R15, (SVNs 856, 858, 855, and 857) despite the very low STD of SLR residuals for the latter satellite (Fig. 4). This can be associated with the wrong value of LRA or antenna offset for the M+ satellites. GLONASS-M R15 (SVN 857) was launched recently in November 2018 and shows a similar SLR offset of + 31 mm to that of all other M+ satellites, despite it still officially belongs to the previous M-generation (Steigenberger et al. 2019). We assume that the construction of the R15 must be similar to M+ , and hence, in the statistics, R15 is considered together with the M+ satellites. One Galileo satellite, E22 (SVN 204), was affected by the failure of both hydrogen masers. Thus, the satellite was deactivated in February 2018. However, E22 was activated for a 2-week test period in June 2019, which resulted in a large negative offset of SLR residuals of − 40 mm.
BeiDou-3 satellites have small mean offsets between − 8 and 0 mm, the offset of BeiDou-2 MEO is − 9 mm, GLONASS-M satellites have a mean offset of − 10 mm, and the values for Galileo IOV are − 11 mm. The remaining satellites have larger offsets, which should be associated with the inferior quality of the orbit determination for GEO and IGSO.
From the comparison between SLR residuals from Fig. 4 and the mean radial precision from Fig. 2, one may conclude that for MEO satellites, the mean orbit STD values in the radial direction are three times lower than those from SLR residuals. In the case of GEO and IGSO, the orbit errors are four times lower than SLR residuals. The SLR residuals provide an independent validation tool for GNSS orbits, especially for the radial component. However, they are affected as well by SLR-specific systematic errors, such as systematic range and time biases, the blue-sky effect, and the satellite signature effect (Arnold et al. 2019). The SLR signature effect is a dependency between the real reflection point and the real distance to a satellite LRA centroid due to different laser incidence angles at the laser retroreflector array. For SLR stations operating in the multi-photon regime, the strongest registered reflection is always associated with the nearest edge of the flat retroreflector array when inclined with respect to the SLR station. The registered distance is shorter as the mean reflection point does not reach the optical centroid of the retroreflector array. For single-photon SLR stations operating with low-energy laser pulses, the probability of the reflection from each corner cube is the same, and after collecting several thousands of reflections, the mean observation corresponds to the optical centroid of the SLR retroreflector (Sośnica et al. 2015, 2018b). The blue-sky effect is related to the weather dependency of SLR observations. SLR observations are typically collected under blue skies when the atmospheric pressure deforms the earth crust downwards because of the atmospheric pressure loading and generates a systematic effect that is detectable in SLR products (Sośnica et al. 2013; Bury et al. 2019b).
The SLR residuals typically overestimate the total error of the radial orbit component due to SLR-specific systematic errors and LRA signature effect, whereas the radial consistency factor from SP3 files underestimates the total orbit error due to some systematic errors that are common in all IGS orbit solutions. The latter errors include the solar radiation pressure modeling errors, satellite attitude modeling errors, satellite and receiver antenna offsets and variations for individual frequencies, as well as errors in signal propagation models, such as tropospheric delay and higher-order ionospheric delay modeling errors.
Figure 5 shows the SLR residuals plotted as a function of the sun elevation angle above the orbital plane (absolute value of β) and the satellite argument of latitude with respect to the latitude of the sun (Δu). Galileo IOV and FOC show a pattern with SLR residuals shifted toward negative values for Δu close to 180°. This pattern was very large for the reduced ECOM and could substantially be reduced when introducing the extended ECOM2 or even almost eliminated when using the a priori box-wing model based on Galileo metadata with ECOM (Bury et al. 2019a). The pattern is due to the cuboid shape of Galileo satellites with the ratio between the X- and Z-bus area of 1.3:3.0, where the latter carries the transmitter antennas. Galileo IOVs show large negative residuals of about − 50 mm for maximum |β| angles when the sun is almost perpendicular to the orbital plane. The highest maximum |β| angles of 76° occur currently only for the Galileo C-plane (Sośnica et al. 2018b). Fortunately, no prominent issues with the orbit determination of eclipsing satellites, when |β|< 12.3°, are visible for Galileo, despite the asymmetrically distributed radiators on the satellite bus.
GLONASS-M and M+ are characterized by positive SLR residuals for Δu close to 180°, which is opposite to the situation of Galileo satellites. This can be explained by the fact that GLONASS-M and M+ are cylindrical with a small surface of the Z-bus side compared to the surface of the cylinder side. Neglecting the estimation of twice-per-revolution accelerations in the satellite–sun direction typically results in systematic effects similar to those observed for M and M+ (Arnold et al. 2015). SLR residuals to GLONASS-K1 are similar to those satellites that have well-performing orbit models due to K1′s regular cuboid shape (Fig. 5). SLR residuals to GLONASS-M+ are systematically shifted toward positive residuals almost for all β angles.
BeiDou-2 IGSO shows very prominent systematic errors. These range from − 160 mm when Δu is close to 0° to + 100 mm for Δu close to 180° (Fig. 4). The pattern is similar to that of GLONASS-M satellites, which could be explained by a much smaller area of the Z-bus surface when compared to the X-bus surface and mismodeling of the higher-order solar radiation pressure perturbing terms. The pattern for BeiDou-2 MEO is similar to that of GLONASS-M satellites, whereas BeiDou GEO and QZSS satellites show large-scale SLR residuals with no evident systematic patterns (Fig. 6). The dominant problem of orbit determination for BeiDou GEO and QZSS satellites is threefold: (1) not employing the normal orbit mode by all analysis centers, (2) satellite orbit instability in the geostationary region due to the resonance with earth rotation and gravity field, and (3) no change of the satellite observation geometry for ground-borne receivers.
Interestingly, BeiDou-3 CAST and SECM are characterized by opposite patterns of SLR residuals (Fig. 5). CAST satellites have patterns similar to those of elongated GLONASS-M satellites, whereas SECM satellites have patterns similar to those of Galileo satellites suggesting that the Z-bus area with transmitter antennas is much larger than the X-bus surface area and vice versa for CAST. The results agree with the BeiDou-3 metadata released in December 2019, disclosing that the effective surface area for SECM is 2.59 and 1.25 m2 for the Z and X panels, respectively, whereas for CAST, it is 2.18 and 2.86 m2 for the Z and X panels. Figure 6 shows the dependencies of the SLR residuals as a function of the satellite elongation angle with respect to the position of the sun with, again, an opposite pattern for SECM and CAST satellites. The systematic patterns of SLR residuals can greatly be reduced when using the a priori box-wing model based on BeiDou-3 metadata and when estimating phase center correction models (Yan et al. 2019). We may thus conclude that the orbit modeling of all GEO, IGSO, and BeiDou-3 satellites should be improved to eliminate the systematic effects, whereas BeiDou-2 MEO has today acceptable orbit processing strategies for high-accuracy geodetic applications.
Conclusions
The IGS ACC established in 2019 an experimental multi-GNSS orbit combination service by adapting the well-established IGS combination procedures. The multi-GNSS orbits are based on individual products generated by IGS and MGEX analysis centers. The number of included GPS, Galileo, and GLONASS is at a stable level of 32, 24, and 23 satellites, respectively. The number of BeiDou satellites in the combined products varies between 10 and 33 because BeiDou-3 satellites are currently being processed only by one center, whereas the number of considered QZSS satellites varies between 1 and 4.
From the orbit quality assessment based on WRMS of the orbit combination, the 3D mean errors are 9, 15, 16, 24, 37, and 26 mm for GPS, Galileo, GLONASS, BeiDou IGSO and MEO, BeiDou GEO, and QZSS, respectively (Table 2). The mean STD values in the radial direction are 5, 8, 9, 14, 21, and 15 mm for the satellites in the same order. STD of SLR residuals is 23, 29, 40, 51, 87, and 72 mm for Galileo, GLONASS, BeiDou MEO, IGSO, GEO, and QZSS, respectively. In all tests, Galileo turned out to be the system with the highest quality of orbit products out of all new GNSS, which confirms that Galileo is fully suitable for the future realizations of International Terrestrial Reference Systems, such as the International Terrestrial Reference Frame ITRF2020.
We found that GPS III G04 (SVN 74) has combination STD twice as large than other GPS satellites. The GLONASS-M+ satellites (SVNs 855, 856, 858) and R15 (SVN 857) have a large positive SLR bias of + 29 mm. GLONASS-K1 R09 (SVN 802) has the smallest STD of SLR residuals out of all GLONASS satellites, which can be explained by a regular cuboid shape of the satellite bus and a new type of SLR retroreflector arranged in the form of a ring.
Galileo IOVs show large negative SLR residuals when the sun is near perpendicular to the orbital plane. BeiDou-3 SECM and CAST satellites have opposite patterns of SLR residuals when analyzed as a function of the sun elongation angle and the satellite latitude, which means that the construction and surface properties of BeiDou-3 MEO satellites from different manufacturers may be different. The orbit modeling of QZSS, BeiDou-2 GEO, IGSO, BeiDou-3, and GPS III needs enhancement to meet the demand for high-accuracy GNSS products. Interestingly, the smallest bias of just − 1 mm and the STD of SLR residuals of 18 mm were obtained for E14 (SVN 202)—the Galileo satellite that had never been intended to fly in an eccentric orbit. Galileo E12 (SVN 102) also has a very low STD value of SLR residuals of 17 mm. However, the bias for E12 is equal to − 7 mm.
The current combination procedure was optimized to fulfill the needs of the IGS repro3 for the ITRF2020 combination, in which double difference GPS–GLONASS–Galileo solutions will be employed. The combination procedure will be improved in near future by considering a proper clock combination, aligning the orbits with the IGS reference frame, and improving the robustness of the combination for BeiDou-3 as well as GEO and IGSO orbits from the BeiDou and QZSS constellations, which are currently missing in many single-day solutions. All these changes are indispensable for a transition from experimental to the fully operational multi-GNSS combination.
Data availability
IGS products, including experimental orbit combinations, are available at https://acc.igs.org/. Information on LRA offsets for GNSS satellites is available at the ILRS Web site: https://ilrs.cddis.eosdis.nasa.gov/missions/satellite_missions/current_missions/index.html and by MGEX: https://mgex.igs.org/. Galileo metadata are available at https://www.gsc-europa.eu/support-to-developers/galileo-satellite-metadata, whereas BeiDou-2 metadata are available at https://en.beidou.gov.cn/WHATSNEWS/201912/t20191209_19641.html and for BeiDou-3 at https://www.beidou.gov.cn/yw/gfgg/201912/t20191209_19613.html. The LRA offsets for BeiDou CAST and SECM in this study are consistent with values from Yang et al. (2019).
References
Arnold D, Meindl M, Beutler G, Dach R, Schaer S, Lutz S, Prange L, Sośnica K, Mervart L, Jäggi A (2015) CODE’s new solar radiation pressure model for GNSS orbit determination. J Geod 89(8):775–791. https://doi.org/10.1007/s00190-015-0814-4
Arnold D, Montenbruck O, Hackel S, Sośnica K (2019) Satellite laser ranging to low earth orbiters: orbit and network validation. J Geod 93(11):2315–2334. https://doi.org/10.1007/s00190-018-1140-4
Beutler G, Brockmann E, Gurtner W, Hugentobler U, Mervart L, Rothacher M, Verdun A (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(6):367–386
Beutler G, Kouba J, Springer T (1995) Combining the orbits of the IGS analysis centers. Bull Géodésique 69(4):200–222. https://doi.org/10.1007/BF00806733
Bruni S, Rebischung P, Zerbini S, Altamimi Z, Errico M, Santi E (2018) Assessment of the possible contribution of space ties on-board GNSS satellites to the terrestrial reference frame. J Geod 92(4):383–399. https://doi.org/10.1007/s00190-017-1069-z
Bury G, Zajdel R, Sośnica K (2019a) Accounting for perturbing forces acting on Galileo using a box-wing model. GPS Solut 23(3):74. https://doi.org/10.1007/s10291-019-0860-0
Bury G, Sośnica K, Zajdel R (2019b) Impact of the atmospheric non-tidal pressure loading on global geodetic parameters based on satellite laser ranging to GNSS. IEEE Trans Geosci Remote Sens 57(6):3574–3590. https://doi.org/10.1109/TGRS.2018.2885845
Bury G, Sośnica K, Zajdel R (2019c) Multi-GNSS orbit determination using satellite laser ranging. J Geod 93(12):2447–2463. https://doi.org/10.1007/s00190-018-1143-1
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
Delva P, Puchades N, Schönemann E, Dilssner F, Courde C, Bertone S, Gonzalez F, Hees A, Le Poncin-Lafitte Ch, Meynadier F (2015) A gravitational redshift test using eccentric Galileo satellites. Phys Rev Lett 121(23):231101. https://doi.org/10.1103/PhysRevLett.121.231101
Duan B, Hugentobler U, Selmke I (2019) The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites. Adv Space Res 63(5):1803–1812. https://doi.org/10.1016/j.asr.2018.11.007
Griffiths J, Ray JR (2009) On the precision and accuracy of IGS orbits. J Geod 83(3–4):277–287. https://doi.org/10.1007/s00190-008-0237-6
Griffiths J (2019) Combined orbits and clocks from IGS second reprocessing. J Geod 93(2):177–195. https://doi.org/10.1007/s00190-018-1149-8
Guo F, Li X, Zhang X, Wang J (2017a) Assessment of precise orbit and clock products for Galileo, BeiDou, and QZSS from IGS multi-GNSS Experiment (MGEX). GPS Solut 21(1):279–290. https://doi.org/10.1007/s10291-016-0523-3
Guo J, Xu X, Zhao Q, Liu J (2016) Precise orbit determination for quad-constellation satellites at Wuhan university: strategy, result validation and comparison. J Geod 90(2):143–159. https://doi.org/10.1007/s00190-015-0862-9
Guo J, Chen G, Zhao Q, Liu J, Liu X (2017b) Comparison of solar radiation pressure models for BDS IGSO and MEO satellites with emphasis on improving orbit quality. GPS Solut 21(2):511–522. https://doi.org/10.1007/s10291-016-0540-2
Hadas T, Kazmierski K, Sośnica K (2019) Performance of Galileo-only dual-frequency absolute positioning using the fully serviceable Galileo constellation. GPS Solut 23(4):108. https://doi.org/10.1007/s10291-019-0900-9
Johnston G, Riddell A, Hausler G (2017) The International GNSS Service. In: Teunissen PJG, Montenbruck O (eds) Springer handbook of global navigation satellite systems, 1st edn. Springer, Cham, pp 967–982. https://doi.org/10.1007/978-3-319-42928-1
Katsigianni G, Loyer S, Perosanz F, Mercier F, Zajdel R, Sośnica K (2019) Improving Galileo orbit determination using zero-difference ambiguity fixing in a multi-GNSS processing. Adv Space Res 63(9):2952–2963. https://doi.org/10.1016/j.asr.2018.08.035
Kazmierski K, Sośnica K, Hadas T (2018) Quality assessment of multi-GNSS orbits and clocks for real-time precise point positioning. GPS Solut 22:11. https://doi.org/10.1007/s10291-017-0678-6
Kouba J (2009) A guide to using international GNSS service (IGS) products. Geodetic Survey Division Natural Resources, Canada, Ottawa, 6, 34, p 2015
Kouba J, Mireault Y (1998) IGS orbit, clock and EOP combined products: an update. In: Brunner FK (ed) Advances in positioning and reference frames. International association of geodesy symposia, vol 118. Springer, Berlin.https://doi.org/10.1007/978-3-662-03714-0_39
Li X, Ge M, Dai X, Ren X, Fritsche M, Wickert J, Schuh H (2015) Accuracy and reliability of multi-GNSS real-time precise positioning: GPS, GLONASS, BeiDou, and Galileo. J Geod 89(6):607–635. https://doi.org/10.1007/s00190-015-0802-8
Lv Y, Geng T, Zhao Q, Xie X, Zhou R (2020) Initial assessment of BDS-3 preliminary system signal-in-space range error. GPS Solut 24:16. https://doi.org/10.1007/s10291-019-0928-x
Mansur G, Sakic P, Männel B, Schuh H (2020) Multi-constellation GNSS orbit combination based on MGEX products. Adv Geosci 50:57–64. https://doi.org/10.5194/adgeo-50-57-2020
Mendes V, Pavlis E (2004) High-accuracy zenith delay prediction at optical wavelengths. Geophys Res Lett. https://doi.org/10.1029/2004GL020308
Montenbruck O, Schmid R, Mercier F, Steigenberger P, Noll C, Fatkulin R, Ganeshan AS (2015) GNSS satellite geometry and attitude models. Adv Space Res 56(6):1015–1029. https://doi.org/10.1016/j.asr.2015.06.019
Montenbruck O, Steigenberger P, Prange L, Deng Z, Zhao Q, Perosanz F, Schmid R (2017) The Multi-GNSS experiment (MGEX) of the international GNSS service (IGS)–achievements, prospects and challenges. Adv Space Res 59(7):1671–1697. https://doi.org/10.1016/j.asr.2017.01.011
Montenbruck O, Steigenberger P, Hauschild A (2018) Multi-GNSS signal-in-space range error assessment–Methodology and results. Adv Space Res 61(12):3020–3038. https://doi.org/10.1016/j.asr.2018.03.041
Noll C, Ricklefs R, Horvath J, Mueller H, Schwatke C, Torrence M (2019) Information resources supporting scientific research for the international laser ranging service. J Geod 93(11):2211–2225. https://doi.org/10.1007/s00190-018-1207-2
Otsubo T, Müller H, Pavlis EC, Torrence MH, Thaller D, Glotov VD, Wilkinson MJ (2019) Rapid response quality control service for the laser ranging tracking network. J Geod 93(11):2335–2344. https://doi.org/10.1007/s00190-018-1197-0
Pearlman MR, Noll CE, Pavlis EC, Lemoine FG, Combrink L, Degnan JJ, Schreibe U (2019a) The ILRS: approaching 20 years and planning for the future. J Geod 93(11):2161–2180. https://doi.org/10.1007/s00190-019-01241-1
Pearlman M, Arnold D, Davis M, Barlier F, Biancale R, Vasiliev V, Sośnica K, Bloßfeld M (2019b) Laser geodetic satellites: a high-accuracy scientific tool. J Geod 93(11):2181–2194. https://doi.org/10.1007/s00190-019-01228-y
Prange L, Orliac E, Dach R, Arnold D, Beutler G, Schaer S, Jäggi A (2017) CODE’s five-system orbit and clock solution the challenges of multi-GNSS data analysis. J Geod 91(4):345–360. https://doi.org/10.1007/s00190-016-0968-8
Press WH, Teukolsky SA, Vetterling WT, Flannery BP (2007) Numerical recipes 3rd edition: the art of scientific computing. Cambridge University Press, Cambridge
Rodríguez J, Appleby G, Otsubo T (2019) Upgraded modelling for the determination of centre of mass corrections of geodetic SLR satellites: impact on key parameters of the terrestrial reference frame. J Geod 93(12):2553–2568. https://doi.org/10.1007/s00190-019-01315-0
Sośnica K, Thaller D, Dach R, Jäggi A, Beutler G (2013) Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results. J Geod 87(8):751–769. https://doi.org/10.1007/s00190-013-0644-1
Sośnica K, Thaller D, Dach R, Steigenberger P, Beutler G, Arnold D, Jäggi 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, Bury G, Zajdel R (2018a) Contribution of Multi-GNSS constellation to SLR-derived terrestrial reference frame. Geophys Res Lett 45(5):2339–2348. https://doi.org/10.1002/2017GL076850
Sośnica K, Prange L, Kaźmierski K, Bury G, Drożdżewski M, Zajdel R, Hadas T (2018b) 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
Sośnica K, Bury G, Zajdel R, Strugarek D, Drożdżewski M, Kazmierski K (2019) Estimating global geodetic parameters using SLR observations to Galileo, GLONASS, BeiDou, GPS, and QZSS. Earth Planets Space 71(1):20. https://doi.org/10.1186/s40623-019-1000-3
Springer TA, Beutler G, Rothacher M (1999) Improving the orbit estimates of GPS satellites. J Geod 73(3):147–157. https://doi.org/10.1007/s001900050230
Springer T, Dach R, Sibthorpe A (2019) GNSS orbit RPR model considerations. In: Proceedings of international GNSS service: Analysis Centre Workshop, Potsdam 2019. https://s3-ap-southeast-2.amazonaws.com/igs-acc-web/igs-acc-website/workshop2019/Orbit_PosPaper.pdf
Steigenberger P, Montenbruck O (2017) Galileo status: orbits, clocks, and positioning. GPS Solut 21(2):319–331. https://doi.org/10.1007/s10291-016-0566-5
Steigenberger P, Thoelert S, Montenbruck O (2019) GPS and GLONASS satellite transmit power: update for IGS repro3. Technical note, DLR/GSOC TN 19-01, 11 July 2019. https://acc.igs.org/repro3/TX_Power_20190711.pdf
Steigenberger P, Montenbruck O (2019) Consistency of MGEX orbit and clock products. Engineering. https://doi.org/10.1016/j.eng.2019.12.005
Thaller D, Dach R, Seitz M, Beutler G, Mareyen M, Richter B (2011) Combination of GNSS and SLR observations using satellite co-locations. J Geod 85(5):257–272. https://doi.org/10.1007/s00190-010-0433-z
Uhlemann M, Gendt G, Ramatschi M, Deng Z (2016) GFZ global multi-GNSS network and data processing results. In: Rizos C, Willis P (eds) IAG 150 years. International association of geodesy symposia. Springer, Berlin, pp 673–679. https://doi.org/10.1007/1345_2015_120
Villiger A, Dach R (eds) (2019) International GNSS service technical report 2018 (IGS Annual Report). IGS Central Bureau and University of Bern, Bern Open Publishing. https://doi.org/10.7892/boris.130408
Weber R, Springer T (2001) The international GLONASS experiment: products, progress and prospects. J Geod 75(11):559–568. https://doi.org/10.1007/s001900100199
Weiss J, Steigenberger P, Springer T (2017) Orbit and clock product generation. In: Teunissen PJG, Montenbruck O (eds) Springer handbook of global navigation satellite systems. Springer, Cham, pp 983–1010. https://doi.org/10.1007/978-3-319-42928-1_34
Yan X, Huang G, Zhang Q, Wang L, Qin Z, Xie S (2019) Estimation of the antenna phase center correction model for the BeiDou-3 MEO satellites. Remote Sens 11(23):2850. https://doi.org/10.3390/rs11232850
Yang H, Xu T, Nie W, Gao F, Guan M (2019) Precise orbit determination of BDS-2 and BDS-3 using SLR. Remote Sens 11(23):2735. https://doi.org/10.3390/rs11232735
Zajdel R, Sośnica K, Bury G (2017) A new online service for the validation of multi-GNSS orbits using SLR. Remote Sens 9(10):1049. https://doi.org/10.3390/rs9101049
Zajdel R, Sośnica K, Dach R, Bury G, Prange L, Jäggi A (2019) Network effects and handling of the geocenter motion in multi-GNSS processing. J Geophys Res Solid Earth 124(6):5970–5989. https://doi.org/10.1029/2019JB017443
Zhao Q, Wang C, Guo J, Wang B, Liu J (2018) Precise orbit and clock determination for BeiDou-3 experimental satellites with yaw attitude analysis. GPS Solut 22(1):4. https://doi.org/10.1007/s10291-017-0673-y
Acknowledgements
We would like to thank IGS and MGEX for providing multi-GNSS data, orbit and clock products, and ILRS for providing the SLR observations and state-of-the-art data processing standards. This work was funded by the National Science Centre, Poland (NCN), grants UMO-2018/29/N/ST10/00289 and UMO‐2018/29/B/ST10/00382. We also thank the Wroclaw Center of Networking and Supercomputing computational grant using MATLAB Software License No. 101979 (https://www.wcss.wroc.pl). We would like to thank Editor-in-Chief and two anonymous reviews for their valuable comments.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sośnica, K., Zajdel, R., Bury, G. et al. Quality assessment of experimental IGS multi-GNSS combined orbits. GPS Solut 24, 54 (2020). https://doi.org/10.1007/s10291-020-0965-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10291-020-0965-5