Baseline-dependent clock offsets in VLBI data analysis

The primary goal of the geodetic Very Long Baseline Interferometry (VLBI) technique is to provide highly accurate terrestrial and celestial reference frames as well as Earth orientation parameters. In compliance with the concept of VLBI, additional parameters reflecting relative offsets and variations of the atomic clocks of the radio telescopes have to be estimated. In addition, reality shows that in many cases significant offsets appear in the observed group delays for individual baselines which have to be compensated for by estimating so-called baseline-dependent clock offsets (BCOs). For the first time, we systematically investigate the impact of BCOs to stress their importance for all kinds of VLBI data analyses. For our investigations, we concentrate on analyzing data from both legacy networks of the CONT17 campaign. Various aspects of BCOs including their impact on the estimates of geodetically important parameters, such as station coordinates and Earth orientation parameters, are investigated. In addition, some of the theory behind the BCO determination, e.g., the impact of changing the reference clock in the observing network on the BCO estimate is introduced together with the relationship between BCOs and triangle delay closures. In conclusion, missing channels, and here in particular at S band, affecting the ionospheric delay calibration, are identified to be the dominant cause for the occurrence of significant BCOs in VLBI data analysis.

The use of astronomy VLBA campaign MOJAVE for geodesy

We investigated the suitability of the astronomical 15 GHz Very Long Baseline Array (VLBA) observing program MOJAVE-5 for estimation of geodetic parameters, such as station coordinates and Earth orientation parameters. We processed a concurrent dedicated VLBA geodesy program observed at 2.3 GHz and 8.6 GHz starting on September 2016 through July 2020 as reference dataset. We showed that the baseline length repeatability from MOJAVE-5 experiments is only a factor of 1.5 greater than from the dedicated geodetic dataset and still below 1 ppb. The wrms of the difference of estimated Earth orientation parameters with respect to the reference IERS C04 time series are a factor of 1.3 to 1.8 worse. We isolated three major differences between the datasets in terms of their possible impact on the geodetic results, i.e. the scheduling approach, treatment of the ionospheric delay, and selection of target radio sources. We showed that the major factor causing discrepancies in the estimated geodetic parameters is the different scheduling approach of the datasets. We conclude that systematic errors in MOJAVE-5 dataset are low enough for these data to be used as an excellent testbed for further investigations on the radio source structure effects in geodesy and astrometry.

VLBI and GPS inter- and intra-technique combinations on the observation level for evaluation of TRF and EOP

We study the effects of combination on the observation level (COL) of different space-geodetic techniques and of networks of the same technique and present the corresponding improvement for the determination of station positions and earth orientation parameters. Data from the continuous geodetic very long baseline interferometry (VLBI) campaign CONT17 are used in a batch least-squares (LSQ) estimator. This campaign includes 15 days of observations with two legacy S/X networks, namely Legacy-1 (L1) and Legacy-2 (L2). For this study the VLBI L1 network is used as the base and reference solution. Data from the L1 network are combined first with data from co-located Global Positioning System (GPS) stations by estimating common tropospheric parameters. The derived station positions repeatabilities of the VLBI and GPS networks are evaluated with respect to single-technique solutions. In terms of precision, we find a 25% improvement for the vertical repeatability of the L1 network, and a 10% improvement for the horizontal one. The GPS network also benefits by 20% and 10% in the horizontal and vertical components, respectively. Furthermore, a combined solution using data of the L1 and L2 network is performed by estimating common earth orientation parameters. The combined L1&GPS and L1&L2 solutions are compared to the reference solution by investigating UT1 and polar motion estimates. UT1 is evaluated in terms of mean bias and formal errors with respect to the International Earth Rotation Service (IERS) C04 products which were used as a priori values. The L1&GPS solution has the lowest formal error and mean bias for UT1 with a 30% improvement. The weighted root mean square (WRMS) and weighted mean offset (WMO) differences between the obtained polar motion estimates and the ones derived by the International GNSS Service (IGS) are also compared. We find that the L1&GPS solution gives the lowest WRMS and WMO, exhibiting an average 40% improvement with respect to the reference solution. The presented results highlight the potential of COL for ongoing transition to multi-space geodetic analysis, e.g., Global Navigation Satellite Systems (GNSS) with the next-generation VLBI system. Graphic Abstract

Effect of non-tidal station loading on Lunar Laser Ranging observatories and on the estimation of Earth orientation parameters

<p>The distance between the observatories on Earth and the retro-reflectors on the Moon has been regularly observed by the Lunar Laser Ranging (LLR) experiment since 1970. In the recent years, observations with bigger telescopes (APOLLO) and at infra-red wavelength (OCA) are carried out, resulting in a better distribution of precise LLR data over the lunar orbit and the observed retro-reflectors on the Moon, and a higher number of LLR observations in total. Providing the longest time series of any space geodetic technique for studying the Earth-Moon dynamics, LLR can also support the estimation of Earth orientation parameters (EOP), like UT1. The increased number of highly accurate LLR observations enables a more accurate estimation of the EOP. In this study, we add the effect of non-tidal station loading (NTSL) in the analysis of the LLR data, and determine post-fit residuals and EOP. The non-tidal loading datasets provided by the German Research Centre for Geosciences (GFZ), the International Mass Loading Service (IMLS), and the EOST loading service of University of Strasbourg in France are included as corrections to the coordinates of the LLR observatories, in addition to the standard corrections suggested by the International Earth Rotation and Reference Systems Service (IERS) 2010 conventions. The Earth surface deforms up to the centimetre level due to the effect of NTSL. By considering this effect in the Institute of Geodesy (IfE) LLR model (called ‘LUNAR’), we obtain a change in the uncertainties of the estimated station coordinates resulting in an up to 1% improvement, an improvement in the post-fit LLR residuals of up to 9%, and a decrease in the power of the annual signal in the LLR post-fit residuals of up to 57%. In a second part of the study, we investigate whether the modelling of NTSL leads to an improvement in the determination of EOP from LLR data. Recent results will be presented.</p>

ESA's Earth Orientation Parameter product

<p>The availability of highly accurate Earth Orientation Parameters (EOPs) in near real time is of major importance for any type of positioning and navigation applications on Earth, Sea, Air and also in Space. This is equally true for all ESA missions and the EU space programs Galileo, EGNOS and Copernicus.</p><p>To ensure operational capability, ESA’s Navigation Support Office developed independent EOP products and services.</p><p>The EOPs are estimated based on a rigorous combination of the ESA’s contributions to the International Association of Geodesy (IAG) that are used as an input for the generation of the International Earth Rotation Service (IERS) products. For the ESA/ESOC EOP products, the individual parameters are combined on normal equation level and propagated with the contribution of model-based predicted Effective Angular Momentum (EAM) functions.</p><p>The ESA/ESOC’s EOP product generation is currently running in pre-operational mode.</p><p>This presentation will provide a high-level overview of the methodology and the status of ESA’s EOP products and services. In this context, the accuracy achieved in the test operations and the roadmap for the publication of ESA’s EOP products and services will be outlined.</p>

Earth Orientation Parameter Estimation by Integrating VLBI and GNSS on the Observation Level

<p>The terrestrial and celestial reference frames are linked by the Earth Orientation Parameters (EOP), which describe the irregularities of the Earth's rotation and are determined by the space geodetic techniques, namely, Very Long Baseline Interferometry (VLBI), Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS). The satellite geodetic techniques (SLR, GNSS, and DORIS) cannot determine the UT1-UTC or celestial pole offsets (CPO), rendering VLBI the only technique capable of determining full EOP set. On the other hand, the GNSS technique provides precise polar motion estimates due to the continuous observations from a globally distributed network. Integrating VLBI and GNSS provides the full set of EOP and guarantees a superior accuracy than any single-technique solution.</p><p>In this study we focus on the integrated estimation of the full EOP set from GNSS and VLBI. Using five VLBI continuous observing campaigns (CONT05–CONT17), the GNSS and VLBI observations are processed concurrently in a common least-squares estimator. The impact of applying global ties (EOP), local ties, and tropospheric ties, and combinations thereof is investigated. The polar motion estimates in integrated solution are dominated by the huge GNSS observations, and the accuracy in terms of weighted root mean squares (WRMS) is ~40 μas compared to the IERS 14 C04 product, which is much better than that of the VLBI-only solution. The UT1-UTC and CPO in the integrated solution also show slight improvement compared to the VLBI-only solution. Moreover, the CPO agreement between the two networks in CONT17, i.e., the VLBA and IVS networks, shows an improvement of 20% to 40% in the integrated solution with different types of ties applied.</p>

The IDS Contribution to the ITRF2020

<p>In the context of the realization of the next International Terrestrial Reference Frame (ITRF2020), the International DORIS Service (IDS) is involved in the estimation of DORIS station positions/velocities as well as Earth orientation parameters from DORIS data. Thus, the 4 IDS Analysis Centers  have re-analyzed all the DORIS observations from the fifteen DORIS satellites from January 1993 to December 2020.0.</p><p>The primary objective of this study is to analyze the DORIS contribution to ITRF2020 in terms of (1) geocenter and scale solutions; (2) station positions and week-to-week repeatability; (3) Earth orientation parameters; (4) a cumulative position and velocity solution.</p><p>Comparisons with the IDS contribution to ITRF2014 will address the benefits of the new antenna models, new models, including improved methods to handle non-conservative force model error on the Jason satellites, as well as the addition of data (compared to ITRF2014) from the latest DORIS missions (e.g. Jason-3, Sentinel-3A/B) in the IDS combination.</p>

Earth Orientation Parameters determination by GNSS & VLBI Combination at Normal Equation Level

<p>The Earth orientation parameters (EOP), the regular products of IERS Earth Orientation Centre, are computed at daily bases by combination of EOP solutions using different astro-geodetic techniques. At SYRTE we have developed a strategy of combination of the <strong>Global Navigation Satellite Systems</strong> (GNSS) and <strong>Very Long Baseline Interferometry</strong> (VLBI) techniques at normal equation level using Dynamo software maintained by CNES (France). This approach allows to produce the EOP at the daily bases, which contains polar coordinates (x,y) and their rates (x<sub>r</sub>,y<sub>r</sub>), universal time UT1 and its rate LOD, and corrections from IAU2000A/2006 precession-nutation model (dX,dY), and in the same run station coordinates constituting the terrestrial frame (TRF). The recorded EOP solutions obtained from GNSS and VLBI combination at weekly bases is recently maintained by SYRTE.</p><p>The strategy applied to consistently combine the IGS and IVS solutions provided in Sinex format over the time period 2000-2021 are presented and the resulting EOP, station positions (TRF) are analysed and evaluated, differences w.r.t. the individual solutions and the IERS time-series investigated.</p>