New IGS Station and Satellite Clock Combination

As part of the International GNSS Service (IGS) multi-GNSS Pilot Project (MGEX), precise orbit and clock products for multi-GNSS constellations have been submitted by several analysis centers (ACs) since 2012. Based on the 30 s satellite clocks from 6 MGEX ACs, the multi-GNSS clocks are combined and the consistency of the AC clocks is assessed in this study. Usually, a linear transformation between the combined and AC solution is used for clock combination, and the clock residuals of the AC solution w.r.t the combined solution are used to determine the weights of the AC. However, any inconsistent satellite-specific bias or nonlinear variations in the clocks induced by the AC’s processing strategy can contaminate the linear transformation as well as the determination of the weight. In this study, the analysis center and satellite-specific bias (ASB) of the MGEX AC clock solutions is first identified and estimated by using observations from globally distributed stations. Moreover, the clock solutions with nonlinear variations induced by the reference clock or nonconstant intersystem bias (ISB) are corrected by aligning the clocks to the selected reference solution before clock combination. With the correction of the ASB, the root-mean-square of the clock residuals decreases significantly and reaches 14-26, 37-91, 33-48 and 12-44 ps for GPS, GLONASS, BDS-2 and Galileo, respectively. In general, the consistency of AC solutions w.r.t the combination reaches 8-16, 27-58, 13-27 and 9-36 ps for GPS, GLONASS, BDS-2 and Galileo in terms of the standard deviation. Finally, the individual AC and combined orbit as well as clock solutions of different constellations are assessed by precise point positioning, and the combined multi-GNSS solutions show competitive performance with the best AC solution in terms of both the positioning accuracy and stability of the reference frame parameters. However, inconsistent scale parameters of both the AC and combined solutions are identified and require more investigation.Graphical

The World Geodetic System (WGS) was conceived as a practical geodetic reference system that maintains consistency with the best scientific terrestrial reference system at the time but also retains some stability. It has evolved from its beginnings as WGS 60 to its present manifestation as WGS 84. A significantly improved Earth Gravitational Model and global geoid were released in 1996, and the realization of the WGS 84 reference frame is consistent with the International Terrestrial Reference Frame (ITRF) 1994 at the 5-cm level. Evaluations of WGS 84 relative to ITRF94 have been made by comparing International GPS Service (IGS) GPS orbits with the National Imagery and Mapping Agency’s precise orbits and by computing WGS 84 positions for a globally distributed set of IGS stations with known ITRF94 coordinates. The results indicate that there is no practical difference between ITRF94 and WGS 84 for mapping, charting, navigation and many survey applications.

Ground station receivers or networks’ location uses different geodetic datums. The geodetic datum defines size and shape of the Earth model and origin of orientation of coordinate systems used to map the Earth. They are modified to converge to the International Terrestrial Reference Frame (ITRF). The Thailand geodetic reference frame has been maintained by the Royal Thai Survey Department (RTSD); currently, it is based on GNSS permanent networks. Precise orbit and clock solutions, earth orientation parameters and tropospheric gradients are applied as fixing parameters to estimate station positions using a network approach where all station coordinates are simultaneously adjusted; defined as a zero-order network. They are determined with respected to ITRF with centimetre precisions and accuracies. The shape of the Earth varies over time due to oceanic tides and plate movements. Changing are more gradually when earthquakes and volcanic eruptions happen. ITRF and other geodetic datums are changing continually; therefore, frame transformations are required. The latest adjusted coordinates are based on ITRF2005, 2008 and 2014 frame and velocities are set at the epoch 2008.83, 2013.81 and 2021.93 respectively. Constraints are cast in form of transformation parameters. The grid shift is then generated to transform from global coordinate systems to local geocentric reference frames. GNSS measurements obtained from another set of ground stations are used as check coordinates. Transformed coordinates are at 1 centimetre-level of accuracy both horizontally and vertically at 95 percent confidence level. Coordinate transformations can be applied to map Thailand based on one map on one datum policy.

<p>The ITRF2020 is the upcoming realization of the International Terrestrial Reference Frame. As the successor of the ITRF2014, it is based on an inter-technique combination of all four space-geodetic techniques, i.e., VLBI, GNSS, SLR and DORIS, and it is based on contributions from different institutions around the world. In this context, the Combination Centre of the International VLBI Service for Geodesy and Astrometry (IVS) – operated by the Federal Agency for Cartography and Geodesy (BKG, Germany) and the Deutsches Geodätisches Forschungsinstitut (DGFI-TUM, Germany) – generates the VLBI intra-technique combination for ITRF2020 utilizing the individual contributions of multiple IVS Analysis Centres (AC).</p><p>For the contribution to the ITRF2020 solution, sessions containing 24h VLBI observations from 1979 until the end of 2020 are reprocessed by 11 ACs and submitted to the IVS Combination Centre. All individual sessions include datum-free normal equations containing station coordinates and source positions as well as full sets of Earth Orientation Parameters (EOP) in the required SINEX format. For ensuring consistently combined solutions, time series of EOP and station coordinates, as well as a VLBI-only Terrestrial Reference Frame (VTRF), have been investigated.</p><p>This contribution focuses on detailed investigations concerning the session-dependency of the scale and the impact of the individual AC contributions to the combination. Thereby significant differences of the scale estimates from the different session types are investigated. Also, the evaluation of the ACs’ contributions to the combined solution will be presented.</p>

This paper evaluates the effect of the accuracy of reference radio sources on the daily estimates of station positions, nutation angle offsets, and the estimated site coordinates determined by very long baseline interferometry (VLBI), which are used for the realization of the international terrestrial reference frame (ITRF). Five global VLBI solutions, based on VLBI data collected between 1979 and 2006, are compared. The reference solution comprises all observed radio sources, which are treated as global parameters. Four other solutions, comprising different sub-sets of radio sources, were computed. The daily station positions for all VLBI sites and the corrections to the nutation offset angles were estimated for these five solutions. The solution statistics are mainly affected by the positional instabilities of reference radio sources, whereas the instabilities of geodetic and astrometric time-series are caused by an insufficient number of observed reference radio sources. A mean offset of the three positional components (Up, North, East) between any two solutions was calculated for each VLBI site. From a comparison of the geodetic results, no significant discrepancies between the respective geodetic solutions for all VLBI sites in the Northern Hemisphere were found. In contrast, the Southern Hemisphere sites were more sensitive to the selected set of reference radio sources. The largest estimated mean offset of the vertical component between two solutions for the Australian VLBI site at Hobart was 4 mm. In the worst case (if a weak VLBI network observed a limited number of reference radio sources) the daily offsets of the estimated height component at Hobart exceeded 100 mm. The exclusion of the extended radio sources from the list of reference sources improved the solution statistics and made the geodetic and astrometric time-series more consistent. The problem with the large Hobart height component offset is magnified by a comparatively small number of observations due to the low slewing rate of the VLBI dish (1°/ s). Unless a minimum of 200 scans are performed per 24-h VLBI experiment, the daily vertical positions at Hobart do not achieve 10 mm accuracy. Improving the slew rate at Hobart and/or having an increased number of new sites in the Southern Hemisphere is essential for further improvement of geodetic VLBI results for Southern Hemisphere sites.

Along with other space geodetics, the Global Navigation Satellite Systems (GNSS) technique has a profound impact on the realization of terrestrial reference frames and the provision of geodetic and geophysical parameters such as Earth Rotation Parameters (ERP) and geocenter (GCC) coordinates. In GNSS data processing ambiguity resolution (AR) is critical to achieve high-precision estimates. Due to the uncalibrated phase delay (UPD), currently most International GNSS Service (IGS) Analysis Centers (ACs) adopt the double-differenced (DD) ambiguity resolution strategy, which eliminates UPD by between-satellite and between-station differencing. Undifferenced (UD) ambiguity resolution, which recovers the integer characteristics of the UD ambiguities by directly estimating the UPD, was originally developed for Precise Point Positioning and nowadays is also applied to network processing, for example for Precise Orbit Determination (POD). In this study, we conducted a comparative study on the UD and DD ambiguity resolution for GNSS POD data processing. We demonstrated that the UD strategy significantly improves the solution compared to the DD solutions. The GPS orbit precision is improved by 20%, 12%, and 5% in the along, cross, and radial component, respectively. Moreover, the ERP are significantly improved, and the agreement with the IGS final products is improved by 28% and 23% for the offsets of x-pole and y-pole, respectively, and by 23%, and 17% for the corresponding rates. Additionally, the repeatability of station coordinates and GCC coordinates is also slightly improved. More important, we confirmed that the two strategies are theoretically equivalent and the outperformance of UD strategy could be explained by its strong robustness to wrong-fixings and its capability to fix the most UD ambiguities, while wrong-fixings and missing fixings are hardly avoidable in the DD strategy due to the requirement of forming DD-ambiguities.

The BeiDou-3 Navigation Satellite System (BDS-3) generates navigation messages from a regional ground tracking network and inter-satellite links, relying on predicted station coordinates and Earth rotation parameters (ERPs). These procedures, however, may lead to orientation errors in the broadcast ephemeris, termed constellation rotation errors. Traditionally, these errors are monitored through Helmert transformation analysis comparing broadcast and post-processed orbits, a method that highlights discrepancies between the global navigation satellite system (GNSS) terrestrial reference frame (TRF) and the International TRF (ITRF) but suffers from a two-week evaluation delay. This study proposes a near real-time approach for estimating constellation rotation errors using high-precision satellite laser ranging (SLR) data, representing the difference between GNSS TRF and ITRF through the International laser ranging service TRF. The impact of SLR data quantity and distribution on rotation estimation accuracy is investigated, and the effectiveness of this method on BDS-3 orbits with ERP-induced orientation deficiencies is validated. Results show that the root mean square of broadcast orbit errors and the orbit user ranging error reduce from 0.8/0.7/0.2 m to 0.1/0.1/0.1 m, respectively. This demonstrates the method’s capability for real-time monitoring and correcting orientation errors within BDS-3 navigation messages.

Evaluation of BDS-2 real-time orbit and clock corrections from four IGS analysis centers

Following the first reprocessing campaign performed by the International GNSS Service (IGS) in 2008, a second reprocessing campaign (repro2) was finalized in 2015. Nine different Analysis Centers (ACs) reanalyzed the history of GNSS data collected by a global tracking network back to 1994 using the latest available models and methodology, and provided daily terrestrial frame solutions among other products. Daily combinations of the AC terrestrial frame solutions provided the IGS input to the next release of the International Terrestrial Reference Frame (ITRF2014). From weighted root mean squares values of the residuals of the daily repro2 combinations, the overall inter-AC level of agreement is assessed to be 1.5 mm for the horizontal components and 4 mm for the vertical component of station positions, 25–40 $$\upmu $$ as for pole coordinates, 140–200 $$\upmu $$ as/day for pole rates, 8–20 $$\upmu $$ s/day for calibrated length-of-day estimates, 4 mm for the X and Y components of geocenter motion, 8 mm for its Z component and 0.5 mm for the terrestrial scale. On the long term, the origins (resp. scales) of the AC terrestrial frames show relative offsets and rates within $$\pm $$ 3 mm and $$\pm $$ 0.3 mm/year (resp. $$\pm $$ 0.5 mm and $$\pm $$ 0.05 mm/year). The combination residuals also present AC-specific features, some of which are explained by known analysis specifics, while others remain under investigation.

The contribution of the International VLBI Service for Geodesy and Astrometry (IVS) to the ITRF2005 (International Terrestrial Reference Frame 2005) has been computed by the IVS Analysis Coordinator’s office at the Geodetic Institute of the University of Bonn, Germany. For this purpose the IVS Analysis Centres (ACs) provided datum-free normal equation matrices in Solution INdependent EXchange (SINEX) format for each 24 h observing session to be combined on a session-by-session basis by a stacking procedure. In this process, common sets of parameters, transformed to identical reference epochs and a prioris, and especially representative relative weights have been taken into account for each session. In order to assess the quality of the combined IVS files, Earth orientation parameters (EOPs) and scaling factors have been derived from the combined normal equation matrices. The agreement of the EOPs of the combined normal equation matrices with those of the individual ACs in terms of weighted root mean square (WRMS) is in the range of 50–60 μas for the two polar motion components and about 3 μs for UT1−UTC. External comparisons with International GNSS Serive (IGS) polar motion components is at the level of 130–170 μas and 21 μs/day for length of day (LOD). The scale of the terrestrial reference frame realized through the IVS SINEX files agrees with ITRF2000 at the level of 0.2 ppb.

Abstract. Solar cycle variations of the amplitudes of the 27-day solar rotation period reflected in the geomagnetic activity index Ap, solar radio flux F10.7cm and critical frequency foF2 for mid-latitude ionosonde station Moscow from the maximum of sunspot cycle 18 to the maximum of cycle 23 are examined. The analysis shows that there are distinct enhancements of the 27-day amplitudes for foF2 and Ap in the late declining phase of each solar cycle while the amplitudes for F10.7cm decrease gradually, and the foF2 and Ap amplitude peaks are much larger for even-numbered solar cycles than for the odd ones. Additionally, we found the same even-high and odd-low pattern of foF2 for other mid-latitude ionosonde stations in Northern and Southern Hemispheres. This property suggests that there exists a 22-year cycle in the F2-layer variability coupled with the 22-year cycle in the 27-day recurrence of geomagnetic activity. Key words. Ionosphere (mid-latitude ionosphere; ionosphere- magnetosphere interactions) – Magnetospheric physics (solar wind-magnetosphere interactions)

This study investigates the contribution of different CH4 sources to the seasonal cycle of δ13C during 2000–2012 by using the TM5 atmospheric transport model, including spatially varying information on isotopic signatures. The TM5 model is able to produce the background seasonality of δ13C, but the discrepancies compared to the observations arise from incomplete representation of the emissions and their source-specific signatures. Seasonal cycles of δ13C are found to be an inverse of CH4 cycles in general, but the anti-correlations between CH4 and δ13C are imperfect and experience a large variation (p=−0.35 to −0.91) north of 30° S. We found that wetland emissions are an important driver in the δ13C seasonal cycle in the Northern Hemisphere and Tropics, and in the Southern Hemisphere Tropics, emissions from fires contribute to the enrichment of δ13C in July–October. The comparisons to the observations from 18 stations globally showed that the seasonal cycle of EFMM emissions in the EDGAR v5.0 inventory is more realistic than in v4.3.2. At northern stations (north of 55° N), modeled δ13C amplitudes are generally smaller by 12–68%, mainly because the model could not reproduce the strong depletion in autumn. This indicates that the CH4 emission magnitude and seasonal cycle of wetlands may need to be revised. In addition, results from stations in northern latitudes (19–40° N) indicate that the proportion of biogenic to fossil-based emissions may need to be revised, such that a larger portion of fossil-based emissions is needed during summer.

The present communication contains results of a detailed study of the horizontal drifts in the ionosphere recorded at selected northern and southern hemisphere stations. Certain important general conclusions about the differences in the behavior of drifts at the northern and southern hemisphere stations and the variation of the drift parameters with latitude are presented.

LAser RElativity Satellite 2 (LARES-2) successfully joined the constellation of geodetic satellites tracked by Satellite Laser Ranging (SLR) stations on July 13, 2022. LARES-2 has a spherical shape and a very favorable area-to-mass ratio that minimizes the non-gravitational orbit perturbations. Due to very small retroreflectors, the spread of center-of-mass corrections for different detectors installed at SLR sites is much smaller than for LAGEOS satellites. LARES-2 orbits at a similar height as LAGEOS-1, however, with a complementary inclination angle of 70° forming a butterfly configuration together with LAGEOS-1.Although the primary objective of LARES-2 is verification of the Lense-Thirring effect emerging from general relativity, the satellite also has a substantial impact on the geodetic parameters derived from SLR observations. We process 18 months of LARES-2 data and compare the LAGEOS-1/2 solutions with the combined LAGEOS-1/2+LARES-2 solutions. We show the impact of LARES-2 on the (1) SLR station coordinates, (2) pole coordinates, (3) length-of-day excess, (4) low-degree gravity field parameters focusing on C20 and C30 coefficients, (5) scale of the reference frame, (6) geocenter motion. We show that LARES-2 can especially improve the Z component of the geocenter coordinates and de-correlate C20 from the length-of-day parameter. The secular drifts of the ascending nodes for LARES-1 and LAGEOS-1 caused by C20 are the same in terms of absolute values but with opposite signs. This allows us to successfully separate the measurements of length-of-day excess (or the UT rate) from the C20-induced changes. We also analyze the empirical accelerations acting on LARES-2 which result from unmodeled non-gravitational orbit perturbations, such as thermal effects, and compare them to those observed for LAGEOS satellites. The observation geometry of LARES-2 is especially beneficial for stations located at high and medium latitudes, which allows it to improve the estimation of station coordinates provided by LAGEOS-1/2. Therefore, LARES-2 substantially contributes not only to general relativity and fundamental physics but also to space geodesy improving the future realizations of the international terrestrial reference frames.