The Tropospheric Products of the International VLBI Service for Geodesy and Astrometry

The IVS runs two tropospheric products: The IVS tropospheric parameter rapid combination monitors the zenith wet delay (ZWD) and zenith total delay (ZTD) of the rapid turnaround sessions R1 and R4. Goal of the combination is the identification and the exclusion of outliers by comparison and the assessment of the precision of current VLBI solutions in terms of tropospheric parameters. The rapid combination is done on a weekly basis four weeks after the observation files are released on IVS Data Centers. Since tropospheric and geodetic parameters, such as vertical station components, can significantly correlate, the consistency of the ZTD can be a measure of the consistency of the corresponding TRF as well. The ZWD mainly rely on accurate atmospheric pressure data. Thus, besides estimation techniques, modeling and analyst s noise, ZWD reflects differences in the atmospheric pressure data applied to the VLBI analysis. The second product, called tropospheric parameter long-term combination, aims for an accurate determination of climatological signals, such as trends of the atmospheric water vapor observed by VLBI. Therefore, the long-term homogeneity of atmospheric pressure data plays a crucial role for this product. The paper reviews the methods applied and results achieved so far and describes the new maintenance through DGFI

Estimating integrated water vapor trends from VLBI, GPS,and numerical weather models: sensitivity totropospheric parameterization

©2018. American Geophysical UnionIn this study, we estimate integrated water vapor (IWV) trends from very long baseline interferometry (VLBI) and global navigation satellite systems (GNSS) data analysis, as well as from numerical weather models (NWMs). We study the impact of modeling and parameterization of the tropospheric delay from VLBI on IWV trends. We address the impact of the meteorological data source utilized to model the hydrostatic delay and the thermal deformation of antennas, as well as the mapping functions employed to project zenith delays to arbitrary directions. To do so, we derive a new mapping function, called Potsdam mapping functions based on NWM data and a new empirical model, GFZ‐PT. GFZ‐PT differs from previous realizations as it describes diurnal and subdiurnal in addition to long‐wavelength variations, it provides harmonic functions of ray tracing‐derived gradients, and it features robustly estimated rates. We find that alternating the mapping functions in VLBI data analysis yields no statistically significant differences in the IWV rates, whereas alternating the meteorological data source distorts the trends significantly. Moreover, we explore methods to extract IWV given a NWM. The rigorously estimated IWV rates from the different VLBI setups, GNSS, and ERA‐Interim are intercompared, and a good agreement is found. We find a quite good agreement comparing ERA‐Interim to VLBI and GNSS, separately, at the level of 75%.DFG, 255986470, GGOS-SIM-2: Simulation des "Global Geodetic Observing System

BKG/DGFI Combination Center Annual Report 2012

This report summarizes the activities of the Federal Agency for Cartography and Geodesy (Bundesamt fuer Kartographie und Geodaesie, BKG) and the German Geodetic Research Institute (Deutsches Geodaetisches Forschungsinstitut, DGFI)BKG/DGFI Combination Center in 2011 and outlines the planned activities for the year 2012. The main focus was to stabilize outlier detection and to update the Web presentation of the combined products

The source structure of 0642+449 detected from the CONT14 observations

The CONT14 campaign with state-of-the-art VLBI data has observed the source 0642+449 with about one thousand observables each day during a continuous observing period of fifteen days, providing tens of thousands of closure delays---the sum of the delays around a closed loop of baselines. The closure delay is independent of the instrumental and propagation delays and provides valuable additional information about the source structure. We demonstrate the use of this new "observable" for the determination of the structure in the radio source 0642+449. This source, as one of the defining sources in the second realization of the International Celestial Reference Frame (ICRF2), is found to have two point-like components with a relative position offset of -426 microarcseconds in right ascension and -66 microarcseconds in declination. The two components are almost equally bright with a flux-density ratio of 0.92. The standard deviation of closure delays for source 0642+449 was reduced from 139 ps to 90 ps by using this two-component model. Closure delays larger than one nanosecond are found to be related to the source structure, demonstrating that structure effects for a source with this simple structure could be up to tens of nanoseconds. The method described in this paper does not rely on a priori source structure information, such as knowledge of source structure determined from direct (Fourier) imaging of the same observations or observations at other epochs. We anticipate our study to be a starting point for more effective determination of the structure effect in VLBI observations.Comment: 14 pages, 13 figures, Accepted by Astronomical Journal on 12 Jul, 201

The impacts of source structure on geodetic parameters demonstrated by the radio source 3C371

Closure quantities measured by very long baseline interferometry (VLBI) observations are independent of instrumental and propagation instabilities and antenna gain factors, but are sensitive to source structure. A new method is proposed to calculate a structure index based on the median values of closure quantities rather than the brightness distribution of a source. The results are comparable to structure indices based on imaging observations at other epochs and demonstrate the flexibility of deriving structure indices from exactly the same observations as used for geodetic analysis and without imaging analysis. A three-component model for the structure of source 3C371 is developed by model-fitting closure phases. It provides a real case of tracing how the structure effect identified by closure phases in the same observations as the delay observables affects the geodetic analysis, and investigating which geodetic parameters are corrupted to what extent by the structure effect. Using the resulting structure correction based on the three-component model of source 3C371, two solutions, with and without correcting the structure effect, are made. With corrections, the overall rms of this source is reduced by 1 ps, and the impacts of the structure effect introduced by this single source are up to 1.4 mm on station positions and up to 4.4 microarcseconds on Earth orientation parameters. This study is considered as a starting point for handling the source structure effect on geodetic VLBI from geodetic sessions themselves.Comment: 5 figures, 15 pages, accepted by Journal of Geodesy at 19 Dec., 201

Deriving Ionospheric Total Electron Content by VLBI Global Observing System Data Analysis During the CONT17 Campaign

This article focuses on the new generation of geodetic very long baseline interferometry (VLBI), the VLBI global observing system (VGOS), and measurements carried out during the CONT17 campaign. It uses broadband technology that increases both the number and precision of observations. These characteristics make VGOS a suitable tool for studying the atmosphere. This study focuses on the effects of the ionosphere on VGOS signals using a model that incorporates and extends ideas originally published in Hobiger et al. (2006, ). Our investigation revealed that the differential total electron content (dTEC) data product calculated with the VGOS post-processing software had a sign error that fortunately, does not change the final values of the phase and group delay. Therefore, this study was a way to identify this problem within the dTEC product. After diagnosing and solving this problem, the underlying model was modified such that instead of considering a single unknown for the latitude gradient of the ionosphere, a time series of latitude gradients were considered that enhanced the resulting vertical total electron content (VTEC) estimates. For evaluation purposes, time series of VTEC at each station during the CONT17 campaign were compared with VTEC obtained from the global navigation satellite system (GNSS). The final agreement between VGOS and GNSS was between 1.1 and 5.9 TEC units (TECU).Peer reviewe

Inter-Comparison of UT1-UTC from 24-Hour, Intensives, and VGOS Sessions during CONT17

This work focuses on the assessment of UT1-UTC estimates from various types of sessions during the CONT17 campaign. We chose the CONT17 campaign as it provides 15 days of continuous, high-quality VLBI data from two legacy networks (S/X band), i.e., Legacy-1 (IVS) and Legacy-2 (VLBA) (having different network geometry and are non-overlapping), two types of Intensive sessions, i.e., IVS and Russian Intensives, and five days of new-generation, broadband VGOS sessions. This work also investigates different approaches to optimally compare dUT1 from Intensives with respect to the 24 h sessions given the different parameterization adopted for analyzing Intensives and different session lengths. One approach includes the estimation of dUT1 from pseudo Intensives, which are created from the 24 h sessions having their epochs synchronized with respect to the Intensive sessions. Besides, we assessed the quality of the dUT1 estimated from VGOS sessions at daily and sub-daily resolution. The study suggests that a different approach should be adopted when comparing the dUT1 from the Intensives, i.e., comparison of dUT1 value at the mean epoch of an Intensive session. The initial results regarding the VGOS sessions show that the dUT1 estimated from VGOS shows good agreement with the legacy network despite featuring fewer observations and stations. In the case of sub-daily dUT1 from VGOS sessions, we found that estimating dUT1 with 6 h resolution is superior to other sub-daily resolutions. Moreover, we introduced a new concept of sub-daily dUT1-tie to improve the estimation of dUT1 from the Intensive sessions. We observed an improvement of up to 20% with respect to the dUT1 from the 24 h sessions.The publication costs are supported within the funding programme “Open Access Publikationskosten” Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project Number 491075472. K.B. is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)– Project-ID 434617780—SFB 1464 (TerraQ). S.B. was partially supported by Generalitat Valenciana (SEJIGENT/2021/001), the European Union—NextGenerationEU (ZAMBRANO 21-04) and Ministerio de Ciencia e Innovación (Spanish Project PID2020-119383GB-I00)

Investigating the Relationship Between Length of Day and El-Niño Using Wavelet Coherence Method

The relationship between the length of day (LOD) and El-Niño Southern Oscillation (ENSO) has been well studied since the 1980s. LOD is the negative time-derivative of UT1-UTC, which is directly proportional to Earth Rotation Angle (ERA), one of the Earth Orientation Parameters (EOP). The EOP can be determined using Very Long Baseline Interferometry (VLBI), which is a space geodetic technique. In addition, satellite techniques such as the Global Navigation Satellite System (GNSS), Satellite Laser Ranging (SLR), Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) can provide Earth Rotation Parameters, i.e., polar motion and LOD. ENSO is a climate phenomenon occurring over the tropical eastern Pacific Ocean that mainly affects the tropics and the subtropics. Extreme ENSO events can cause extreme weather like flooding and droughts in many parts of the world. In this work, we investigated the effect of ENSO on the LOD from January 1979 to April 2022 using the wavelet coherence method. This method computes the coherence between the two non-stationary time-series in the time-frequency domain using the real-valued Morlet wavelet. We used the Multivariate ENSO index version 2 (MEI v.2) which is the most robust series as the climate index for the ENSO, and LOD time-series from IERS (EOP 14 C04 (IAU2000A)). We also used Oceanic Niño and Southern Oscillation index in this study for comparison. The results show strong coherence of 0.7 to 0.9 at major ENSO events for the periods 2–4 years between LOD and MEI.v2.Kyriakos Balidakis is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project-ID 434617780—SFB 1464 (TerraQ). Santiago Belda was partially supported by Generalitat Valenciana (SEJIGENT/2021/001), the European Union—NextGenerationEU (ZAMBRANO 21-04) and Ministerio de Ciencia e Innovación (Spanish Project PID2020-119383GB-I00). Chaiyaporn Kitpracha acknowledges funding from Deutscher Akademischer Austauschdienst (DAAD) under grant number 91650950

A new method to improve the prediction of the celestial pole offsets

Knowledge of the Earth’s changing rotation is fundamental to positioning objects in space and on the planet. Nowadays, the Earth’s orientation in space is expressed by five Earth Orientation Parameters (EOP). Many applications in astronomy, geosciences, and space missions require accurate EOP predictions. Operational predictions are released daily by the Rapid Service/Prediction Centre of the International Earth Rotation and Reference Systems Service (IERS). The prediction procedures and performances differ for the three EOP classes: polar motion, rotation angle (UT1-UTC), and the two celestial pole offsets (CPO), dX and dY. The IERS Annual Report 2016 shows Rapid Service CPO predictions errors with respect to IERS 08 C04 observations in 2016 ranging from 120 to 140 μas in 40 days for dX, and 100–160 μas for dY. We test a new method for the CPO prediction based on the recent availability of sophisticated empirical models for the Free Core Nutation, a main component of the CPO variations. We found it allows predicting both CPO with error estimates for the period 2000–2016 lower than the 2016 Rapid Service products, reaching about 85 μas after 40 days and near 90 μas after a year. These results would represent a 35–40% improvement.SB and JF work were partially supported by projects AYA2016–79775-P (AEI/FEDER, UE), APOSTD/2026/079 and ERC-2017-STG SENTIFLEX (755617)