Very Long Baseline Interferometry Analysis Research Articles

Hartebeesthoek Radio Astronomy Observatory (HartRAO) antenna axis offset determined by geodetic VLBI analysis and ground survey

In the Very Long Baseline Interferometry (VLBI) space geodetic technique, various station-specific error sources corrupt the observable VLBI delay. An antenna axis offset (AO) model is applied in the VLBI data analysis for antennas with non-intersecting rotational axes, such as the 26-m and 15-m antennas for the Hartebeesthoek Radio Astronomy Observatory (HartRAO). The a priori AO values recommended by the International VLBI Service for Geodesy and Astrometry (IVS) for use in geodetic VLBI data analysis are taken, where possible, from values measured in ground surveys. The a priori AO values used for the HartRAO antennas in geodetic VLBI analysis have been identified as possible sources of error. The a priori AO value of 6695.3 mm for the 26-m antenna originates from a 2003 co-locational ground survey, conducted before a major bearing repair in 2008, which could have changed the AO. The a priori AO value of 1495.0 mm for the 15-m antenna was determined in 2007 in only a preliminary GPS survey. In this study, the respective AO values of the HartRAO 26-m and 15-m antennas were estimated from a VLBI analysis using the Vienna VLBI and Satellite Software (VieVS) and compared with measurements from co-locational ground surveys. It was found that the VLBI estimated values do not agree within the formal margins of error with the ground survey values, in that they differ by up to eight millimetres (8 mm) for the 26-m antenna and up to five millimetres (5 mm) for the 15-m antenna. As the ground survey values are considered to be more accurate than the VLBI estimated values, a further investigation of the site-specific error sources that may be contaminating the accuracy of VLBI results is required.

Improving VLBI analysis by tropospheric ties in GNSS and VLBI integrated processing

Tropospheric delay modeling is challenging in high-precision Very Long Baseline Interferometry (VLBI) analysis due to the rapid water vapor variation and imperfect observation geometry, where observations from Global Navigation Satellite Systems (GNSS) co-locations can enhance the VLBI analysis. We investigate the impact of tropospheric ties in the VLBI and GNSS integrated processing during the CONT05–CONT17 campaigns, and present a method that automatically handles the systematic tropospheric tie biases. Applying tropospheric ties at VLBI–GNSS co-locations enhances the observation geometry and improves the solution reliability. The VLBI network is stabilized, with station coordinate repeatability improved by 12% horizontally and by 28% vertically, and the network scale improved by 32%. The Earth Orientation Parameters (EOP) improve by up to 20%. Both zenith delay and gradient ties contribute to the improvement of EOP, whereas the gradient ties contribute mainly to the improvement of length of day and celestial pole offsets.

PSR J2222−0137

Context. The PSR J2222−0137 binary system has a set of features that make it a unique laboratory for tests of gravity theories. Aims. To fully exploit the system’s potential for these tests, we aim to improve the measurements of its physical parameters, spin and orbital orientation, and post-Keplerian parameters, which quantify the observed relativistic effects. Methods. We describe an improved analysis of archival very long baseline interferometry (VLBI) data, which uses a coordinate convention in full agreement with that used in timing. We have also obtained much improved polarimetry of the pulsar with the Five hundred meter Aperture Spherical Telescope (FAST). We provide an improved analysis of significantly extended timing datasets taken with the Effelsberg, Nançay, and Lovell radio telescopes; this also includes previous timing data from the Green Bank Telescope. Results. From the VLBI analysis, we have obtained a new estimate of the position angle of the ascending node, Ω = 189−18+19 deg (all uncertainties are 68% confidence limits), and a new reference position for the pulsar with an improved and more conservative uncertainty estimate. The FAST polarimetric results, and in particular the detection of an interpulse, yield much improved estimates for the spin geometry of the pulsar, in particular an inclination of the spin axis of the pulsar of ∼84 deg. From the timing, we obtain a new ∼1% test of general relativity (GR) from the agreement of the Shapiro delay parameters and the rate of advance of periastron. Assuming GR in a self-consistent analysis of all effects, we obtain much improved masses: 1.831(10) M⊙ for the pulsar and 1.319(4) M⊙ for the white dwarf companion; the total mass, 3.150(14) M⊙, confirms this as the most massive double degenerate binary known in the Galaxy. This analysis also yields the orbital orientation; in particular, the orbital inclination is 85.27(4) deg – indicating a close alignment between the spin of the pulsar and the orbital angular momentum – and Ω = 187.7(5.7) deg, which matches our new VLBI estimate. Finally, the timing also yields a precise measurement of the variation in the orbital period, Ṗb = 0.251(8) × 10−12 ss−1; this is consistent with the expected variation in the Doppler factor plus the orbital decay caused by the emission of gravitational waves predicted by GR. This agreement introduces stringent constraints on the emission of dipolar gravitational waves.

A priori zenith wet delays in the analysis of VLBI CONT sessions

Tropospheric delay is one of the major sources of error in VLBI (Very Long Baseline Interferometry) analysis. The principal component of this error can be accurately computed through reliable surface pressure data —hydrostatic delay— yet there is also a small but volatile component —wet delay— which is difficult to be modelled a priori. In VLBI analysis, troposphere delay is typically modelled in the theoretical delays using Zenith Hydrostatic Delays (ZHD) and a dry mapping function. Zenith Wet Delay (ZWD) is not modelled but estimated in the analysis process. This work studies inter alia the impact of including external GNSS estimates to model a priori ZWD in VLBI analysis, as well as other models of a priori ZWD.In a first stage, two different sources of GNSS troposphere products are compared to VLBI troposphere estimates in a period of 5 years. The solution with the best agreement to VLBI results is injected in the VLBI analysis as a priori ZWD value and is compared to other options to model a priori ZWD. The dataset used for this empirical analysis consists of the six CONT campaigns.It has been found that modelling a priori ZWD has no significant impact either on baseline length and coordinates repeatabilities. Nevertheless, modelling a priori ZWD can change the magnitude of the estimated coordinates a few millimeters in the up component with respect to the non-modelling approach. In addition, the influence of a priori ZWD on Earth Orientation Parameters (EOP) and troposphere estimates —Zenith Total Delays (ZTD) and gradients—has also been analysed, resulting in a small but significant impact on both geodetic products.

Benefits of non-tidal loading applied at distinct levels in VLBI analysis

In the analysis of very long baseline interferometry (VLBI) observations, many geophysical models are used for correcting the theoretical signal delay. In addition to the conventional models described by Petit and Luzum (eds) (IERS Conventions, 2010), we are applying different parts of non-tidal site loading, namely the atmospheric, oceanic, and hydrological ones. To investigate their individual contributions, these parts are considered both separately and combined to a total loading. The application of the corresponding site displacements is performed at two distinct levels of the geodetic parameter estimation process (observation and normal equation level), which turn out to give very similar results in many cases. To validate our findings internally, the site displacements are provided by two different data centres: the Earth-System-Modelling group at the Deutsches GeoForschungsZentrum in Potsdam (ESMGFZ, see Dill and Dobslaw, J Geophys Res Solid Earth, 2013. https://doi.org/10.1002/jgrb.50353)] and the International Mass Loading Service [IMLS, see Petrov (The international mass loading service, 2015)]. We show that considering non-tidal loading is actually useful for mitigating systematic effects in the VLBI results, like annual signals in the station height time series. If the sum of all non-tidal loading parts is considered, the WRMS of the station heights and baseline lengths is reduced in 80–90% of all cases, and the relative improvement is about -,3.5% on average. The main differences between our chosen providers originate from hydrological loading.

Geodetic VLBI for precise orbit determination of Earth satellites: a simulation study

Recent efforts of tracking low Earth orbit and medium Earth orbit (MEO) satellites using geodetic very long baseline interferometry (VLBI) raise questions on the potential of this novel observation concept for space geodesy. Therefore, we carry out extensive Monte Carlo simulations in order to investigate the feasibility of geodetic VLBI for precise orbit determination (POD) of MEO satellites and assess the impact of quality and quantity of satellite observations on the derived geodetic parameters. The MEO satellites are represented in our study by LAGEOS-1/-2 and a set of Galileo satellites. The concept is studied on the basis of 3-day solutions in which satellite observations are included into real schedules of the continuous geodetic VLBI campaign 2017 (CONT17) as well as simulated schedules concerning the next-generation VLBI system, known as the VLBI Global Observing System (VGOS). Our results indicate that geodetic VLBI can perform on a comparable level as other space-geodetic techniques concerning POD of MEO satellites. For an assumed satellite observation precision better than 14.1 mm (47 ps), an average 3D orbit precision of 2.0 cm and 6.3 cm is found for schedules including LAGEOS-1/-2 and Galileo satellites, respectively. Moreover, geocenter offsets, which were so far out of scope for the geodetic VLBI analysis, are close to the detection limit for the simulations concerning VGOS observations of Galileo satellites, with the potential to further enhance the results. Concerning the estimated satellite orbits, VGOS leads to an average precision improvement of 80% with respect to legacy VLBI. In absolute terms and for satellite observation precision of 14.1 mm (47 ps), this corresponds to an average value of 17 mm and 7 mm concerning the 3D orbit scatter and precision of geocenter components, respectively. As shown in this study, a poor satellite geometry can degrade the derived Earth rotation parameters and VLBI station positions, compared to the quasar-only reference schedules. Therefore, careful scheduling of both quasar and satellite observations should be performed in order to fully benefit from this novel observation concept.

Gravitational deformation of ring-focus antennas for VGOS: first investigations at the Onsala twin telescopes project

The receiving properties of radio telescopes used in geodetic and astrometric very long baseline interferometry (VLBI) depend on the surface quality and stability of the main reflector. Deformations of the main reflector as well as changes in the sub-reflector position affect the geometrical ray path length significantly. The deformation pattern and its impact on the VLBI results of conventional radio telescopes have been studied by several research groups using holography, laser tracker, close-range photogrammetry and laser scanner methods. Signal path variations (SPV) of up to 1 cm were reported, which cause, when unaccounted for, systematic biases of the estimated vertical positions of the radio telescopes in the geodetic VLBI analysis and potentially even affect the estimated scale of derived global geodetic reference frames. As a result of the realization of the VLBI 2010 agenda, the geodetic VLBI network is currently extended by several new radio telescopes, which are of a more compact and stiffer design and are able to move faster than conventional radio telescopes. These new telescopes will form the backbone of the next generation geodetic VLBI system, often referred to as VGOS (VLBI Global Observing System). In this investigation, for the first time the deformation pattern of this new generation of radio telescopes for VGOS is studied. ONSA13NE, one of the Onsala twin telescopes at the Onsala Space Observatory, was observed in several elevation angles using close-range photogrammetry. In general, these methods require a crane for preparing the reflector as well as for the data collection. To reduce the observation time and the technical effort during the measurement process, an unmanned aircraft system (UAS) was used for the first time. Using this system, the measurement campaign per elevation angle took less than 30 min. The collected data were used to model the geometrical ray path and its variations. Depending on the distance from the optical axis, the ray path length varies in a range of about pm ,1,hbox {mm}. To combine the ray path variations, an illumination function was introduced as weighting function. The resulting total SPV is about - ,0.5 mm. A simple elevation-dependent SPV model is presented that can easily be used and implemented in VLBI data analysis software packages to correct for gravitational deformation in VGOS radio telescopes. The uncertainty is almost 200,upmu hbox {m} (2sigma ) and is derived by Monte Carlo simulations applied to the entire analysis process.

Comparing atmospheric data and models at station Wettzell during CONT17

Abstract. During the Continuous Very Long Baseline Interferometry (VLBI) Campaign 2017 (CONT17), carried out from 28 November through 12 December 2017, an extensive data set of atmospheric observations was acquired at the Geodetic Observatory Wettzell. In addition to in situ measurements of temperature, humidity, pressure or wind speed at the surface, radiosonde ascents yielded meteorological parameters continually up to 25 km height, and integrated water vapor (IWV) was obtained at several elevations and azimuths from a water vapor radiometer. Troposphere delays estimated from Global Navigation Satellite Systems (GNSS) observations plus comparative values from two different Numerical Weather Models (NWMs) complete the abundance of data. In this presentation, we compare these data sets to parameters of the Vienna Mapping Functions 1 and 3 (VMF1 & VMF3), which are based on NWM data by the ECMWF, and to estimates of VLBI analysis using the Vienna VLBI and Satellite Software (VieVS). On the one hand, we contrast the variety of troposphere delays in zenith direction with each other, while on the other hand we utilize radiosonde data and meteorological observations at the site to create local mapping functions which can then be compared to VMF3 and VMF1 at Wettzell. In general, we thus received very good accordance between the different solutions. Also in terms of the mapping functions, the local radiosonde mapping function is in consistence with VMF1 and VMF3 with differences less than 5 mm at 5∘ elevation.

Source Structure and Measurement Noise Are as Important as All Other Residual Sources in Geodetic VLBI Combined

AbstractThe celestial sources observed by geodetic very‐long‐baseline interferometry (VLBI), most of which are quasars, are resolved—they have a finite angular size as measured by the interferometer—yet modern geodetic VLBI analysis techniques still treat these sources as if they are point‐like. The errors introduced by source structure, the angular distribution of emission, have traditionally been believed to be a small fraction of the total measurement error budget. We analyzed the VLBI data of the CONT14 continuous observing campaign using standard geodetic and imaging techniques, along with an analysis of closures, combinations of interferometric measurements where the error contributions from station‐based effects such as tropospheric and ionospheric errors, clock and cable length errors, station position errors, and so on, cancel exactly, leaving only nonclosing errors, principally source structure and measurement noise. Our statistical analysis of the closure information, based on actual VLBI observations revealed that nearly half (40%) of the total geodetic VLBI error budget (in terms of variance) in CONT14 came from source‐structure effects, nearly the same variance as all station‐based error sources (43%). The remaining error due to other nonsource‐structure, nonclosing errors is also significant (17%). Previous studies, which were primarily simulation‐based, did not account for the full range of effects of source structure on real measurements. Source structure is a major contributor to errors in geodetic VLBI, and source structure must be taken into account in the entire VLBI operational chain, from scheduling to analysis, in order to mitigate its effects on astrometric and geodetic measurements.

Refined discrete and empirical horizontal gradients in VLBI analysis

Missing or incorrect consideration of azimuthal asymmetry of troposphere delays is a considerable error source in space geodetic techniques such as Global Navigation Satellite Systems (GNSS) or Very Long Baseline Interferometry (VLBI). So-called horizontal troposphere gradients are generally utilized for modeling such azimuthal variations and are particularly required for observations at low elevation angles. Apart from estimating the gradients within the data analysis, which has become common practice in space geodetic techniques, there is also the possibility to determine the gradients beforehand from different data sources than the actual observations. Using ray-tracing through Numerical Weather Models (NWMs), we determined discrete gradient values referred to as GRAD for VLBI observations, based on the standard gradient model by Chen and Herring (J Geophys Res 102(B9):20489–20502, 1997. https://doi.org/10.1029/97JB01739) and also for new, higher-order gradient models. These gradients are produced on the same data basis as the Vienna Mapping Functions 3 (VMF3) (Landskron and Böhm in J Geod, 2017. https://doi.org/10.1007/s00190-017-1066-2), so they can also be regarded as the VMF3 gradients as they are fully consistent with each other. From VLBI analyses of the Vienna VLBI and Satellite Software (VieVS), it becomes evident that baseline length repeatabilities (BLRs) are improved on average by 5% when using a priori gradients GRAD instead of estimating the gradients. The reason for this improvement is that the gradient estimation yields poor results for VLBI sessions with a small number of observations, while the GRAD a priori gradients are unaffected from this. We also developed a new empirical gradient model applicable for any time and location on Earth, which is included in the Global Pressure and Temperature 3 (GPT3) model. Although being able to describe only the systematic component of azimuthal asymmetry and no short-term variations at all, even these empirical a priori gradients slightly reduce (improve) the BLRs with respect to the estimation of gradients. In general, this paper addresses that a priori horizontal gradients are actually more important for VLBI analysis than previously assumed, as particularly the discrete model GRAD as well as the empirical model GPT3 are indeed able to refine and improve the results.

Application of ray-traced tropospheric slant delays to geodetic VLBI analysis

The correction of tropospheric influences via so-called path delays is critical for the analysis of observations from space geodetic techniques like the very long baseline interferometry (VLBI). In standard VLBI analysis, the a priori slant path delays are determined using the concept of zenith delays, mapping functions and gradients. The a priori use of ray-traced delays, i.e., tropospheric slant path delays determined with the technique of ray-tracing through the meteorological data of numerical weather models (NWM), serves as an alternative way of correcting the influences of the troposphere on the VLBI observations within the analysis. In the presented research, the application of ray-traced delays to the VLBI analysis of sessions in a time span of 16.5 years is investigated. Ray-traced delays have been determined with program RADIATE (see Hofmeister in Ph.D. thesis, Department of Geodesy and Geophysics, Faculty of Mathematics and Geoinformation, Technische Universität Wien. http://resolver.obvsg.at/urn:nbn:at:at-ubtuw:1-3444, 2016) utilizing meteorological data provided by NWM of the European Centre for Medium-Range Weather Forecasts (ECMWF). In comparison with a standard VLBI analysis, which includes the tropospheric gradient estimation, the application of the ray-traced delays to an analysis, which uses the same parameterization except for the a priori slant path delay handling and the used wet mapping factors for the zenith wet delay (ZWD) estimation, improves the baseline length repeatability (BLR) at 55.9% of the baselines at sub-mm level. If no tropospheric gradients are estimated within the compared analyses, 90.6% of all baselines benefit from the application of the ray-traced delays, which leads to an average improvement of the BLR of 1 mm. The effects of the ray-traced delays on the terrestrial reference frame are also investigated. A separate assessment of the RADIATE ray-traced delays is carried out by comparison to the ray-traced delays from the National Aeronautics and Space Administration Goddard Space Flight Center (NASA GSFC) (Eriksson and MacMillan in http://lacerta.gsfc.nasa.gov/tropodelays, 2016) with respect to the analysis performances in terms of BLR results. If tropospheric gradient estimation is included in the analysis, 51.3% of the baselines benefit from the RADIATE ray-traced delays at sub-mm difference level. If no tropospheric gradients are estimated within the analysis, the RADIATE ray-traced delays deliver a better BLR at 63% of the baselines compared to the NASA GSFC ray-traced delays.

Tropospheric delay determination by Kalman filtering VLBI data

The troposphere is one of the most important error sources for space geodetic techniques relying on radio signals. Since it is not possible to model the wet part of the tropospheric delay with sufficient accuracy, it needs to be estimated from the observational data. In the analysis of very long baseline interferometry (VLBI) data, the parameter estimation is routinely performed using a least squares adjustment. In this paper, we investigate the application of a Kalman filter for parameter estimation, specifically focusing on the tropospheric delays. The main advantages of a Kalman filter are its real-time capability and stochastic approach. We focused on the latter and derived stochastic models for VLBI zenith wet delays, taking into account temporal and location-based differences. Compared to a static noise model, the quality of station coordinates, also estimated in the Kalman filter, increased as a result. In terms of baseline length and station coordinate repeatabilities, this improvement amounted to 2.3 %. Additionally, we compared the Kalman filter and least squares results for VLBI with zenith wet delays derived from GPS (Global Positioning System), water vapor radiometers, and ray tracing in numerical weather models. The agreement of the Kalman filter VLBI solution with respect to water vapor radiometer data was larger than that of the least squares solution by 6–15 %. Our investigations are based on selected VLBI data (CONT campaigns) that are closest to how future VLBI infrastructure is designed to operate. With the aim for continuous and near real-time parameter estimation and the promising results which we have achieved in this study, we expect Kalman filtering to grow in importance in VLBI analysis.

Tropospheric delay ray tracing applied in VLBI analysis

AbstractTropospheric delay modeling error continues to be one of the largest sources of error in VLBI (very long baseline interferometry) analysis. For standard operational solutions, we use the VMF1 elevation‐dependent mapping functions derived from European Centre for Medium‐Range Weather Forecasts data. These mapping functions assume that tropospheric delay at a site is azimuthally symmetric. As this assumption is not true, we have instead determined the ray trace delay along the signal path through the troposphere for each VLBI quasar observation. We determined the troposphere refractivity fields from the pressure, temperature, specific humidity, and geopotential height fields of the NASA Goddard Space Flight Center Goddard Earth Observing System version 5 numerical weather model. When applied in VLBI analysis, baseline length repeatabilities were improved compared with using the VMF1 mapping function model for 72% of the baselines and site vertical repeatabilities were better for 11 of 13 sites during the 2 week CONT11 observing period in September 2011. When applied to a larger data set (2011–2013), we see a similar improvement in baseline length and also in site position repeatabilities for about two thirds of the stations in each of the site topocentric components.

A complete VLBI delay model for deforming radio telescopes: the Effelsberg case

Deformations of radio telescopes used in geodetic and astrometric very long baseline interferometry (VLBI) observations belong to the class of systematic error sources which require correction in data analysis. In this paper we present a model for all path length variations in the geometrical optics of radio telescopes which are due to gravitational deformation. The Effelsberg 100 m radio telescope of the Max Planck Institute for Radio Astronomy, Bonn, Germany, has been surveyed by various terrestrial methods. Thus, all necessary information that is needed to model the path length variations is available. Additionally, a ray tracing program has been developed which uses as input the parameters of the measured deformations to produce an independent check of the theoretical model. In this program as well as in the theoretical model, the illumination function plays an important role because it serves as the weighting function for the individual path lengths depending on the distance from the optical axis. For the Effelsberg telescope, the biggest contribution to the total path length variations is the bending of the main beam located along the elevation axis which partly carries the weight of the paraboloid at its vertex. The difference in total path length is almost $$-$$ 100 mm when comparing observations at 90 $$^\circ $$ and at 0 $$^\circ $$ elevation angle. The impact of the path length corrections is validated in a global VLBI analysis. The application of the correction model leads to a change in the vertical position of $$+120$$ mm. This is more than the maximum path length, but the effect can be explained by the shape of the correction function.

Ray‐traced tropospheric delays in VLBI analysis

We develop a ray‐tracing package for the calculation of path delays of microwave signals in the troposphere based on numerical weather models which we use for the determination of the delays of geodetic Very Long Baseline Interferometry (VLBI) observations. We show results for a two‐week campaign of continuous VLBI sessions in 2008 (CONT08), where we apply those ray‐traced delays and analyze the repeatability of baseline lengths in comparison to a standard approach with zenith delays and mapping functions. We find improvement in baseline length repeatabilities when no tropospheric gradients are estimated in the analysis. Furthermore, ray‐traced delays are applied for Intensive sessions containing the stations Tsukuba (Japan) and Wettzell (Germany) for the determination of Universal Time (UT1). We perform an external validation using GPS‐derived length‐of‐day values and find an improvement for UT1 with ray‐traced delays by up to 4.5%.

Improved geodetic European very-long-baseline interferometry solution using models of antenna gravitational deformation

Very-long-baseline interferometry (VLBI) is used for establishing global geodetic networks where the coordinates attain a 1-mm level of precision. Technique-dependent bias can degrade the VLBI positioning accuracy if it is present and unaccounted for. Among the potential bias, gravitational flexure of VLBI telescopes can vary the path traveled by the incoming radio signal and induce a bias in the height component of the station position. We process here more than 100 European VLBI sessions spanning 1990-2009 with VLBI time delay/Solve software, as the only VLBI analysis package that can be used to correct signal-path variation (SPV) due to gravitational flexure of VLBI telescopes. Currently, SPV models are neglected in VLBI data analysis. To determine the kinematics of the European area over the last 20 years and to assess the effects of telescope gravitational deformation on geodetic VLBI estimates, we perform two VLBI solutions with and without SPV models for telescopes in Medicina (northern Italy) and Noto (southern Italy). The two solutions differ by 8.8 mm and 7.2 mm in their height components, with this bias being one order of magnitude larger than the formal errors of the estimated heights. SPV models impact uniquely on the height component of stations where SPVs are modeled. Velocities are not affected by the use of the Medicina and Noto SPV models, and we show that the crustal kinematics derived from VLBI does not suffer from a lack of information with regard to the flexure of other telescopes.

Atmospheric loading corrections at the observation level in VLBI analysis

This paper investigates whether in very long baseline interferometry (VLBI) analysis atmospheric loading corrections should be applied a priori at the observation level or whether it is sufficient to correct for atmospheric loading effects a posteriori by adding constant values per session to the estimated station coordinates. Simulated observations at single stations corresponding to the precise point positioning approach of global navigation satellite systems show that the atmospheric loading effect can be fully recovered by a posteriori corrections, i.e., the height differences between both approaches stay well below 1 mm. However, real global VLBI network solutions with sessions from 1984 to 2008 reveal that the effect of neglected atmospheric loading corrections at the stations is distributed to the other stations in the network, thus resulting in station height differences between solutions with observation level and with a posteriori corrections which can be as large as 10 mm and a ‘damping’ effect of the corrections. As soon as the terrestrial reference frame and the corresponding coordinate time series are determined, it would be conceptually wrong to apply atmospheric loading corrections at the VLBI stations. We recommend the rigorous application of atmospheric loading corrections at the observation level to all stations of a VLBI network because the seven parameters for translation, rotation, and in particular the network-scale of VLBI networks are significantly affected.

Troposphere gradients from the ECMWF in VLBI analysis

Modeling path delays in the neutral atmosphere for the analysis of Very Long Baseline Interferometry (VLBI) observations has been improved significantly in recent years by the use of elevation-dependent mapping functions based on data from numerical weather models. In this paper, we present a fast way of extracting both, hydrostatic and wet, linear horizontal gradients for the troposphere from data of the European Centre for Medium-range Weather Forecasts (ECMWF) model, as it is realized at the Vienna University of Technology on a routine basis for all stations of the International GNSS (Global Navigation Satellite Systems) Service (IGS) and International VLBI Service for Geodesy and Astrometry (IVS) stations. This approach only uses information about the refractivity gradients at the site vertical, but no information from the line-of-sight. VLBI analysis of the CONT02 and CONT05 campaigns, as well as all IVS-R1 and IVS-R4 sessions in the first half of 2006, shows that fixing these a priori gradients improves the repeatability for 74% (40 out of 54) of the VLBI baseline lengths compared to fixing zero or constant a priori gradients, and improves the repeatability for the majority of baselines compared to estimating 24-h offsets for the gradients. Only if 6-h offsets are estimated, the baseline length repeatabilities significantly improve, no matter which a priori gradients are used.

Radio source instability in VLBI analysis

The source position time-series for many of the frequently observed radio sources in the NASA geodetic very long baseline interferometry (VLBI) program show systematic linear and non-linear variations of as much as 0.5 mas (milli-arc-seconds) to 1.0 mas, due mainly to source structure changes. In standard terrestrial reference frame (TRF) geodetic solutions, it is a common practice to only estimate a global source position for each source over the entire history of VLBI observing sessions. If apparent source position variations are not modeled, they produce corresponding systematic variations in estimated Earth orientation parameters (EOPs) at the level of 0.02–0.04 mas in nutation and 0.01–0.02 mas in polar motion. We examine the stability of position time-series of the 107 radio sources in the current NASA geodetic source catalog since these sources have relatively dense observing histories from which it is possible to detect systematic variations. We consider different strategies for handling source instabilities where we (1) estimate the positions of unstable sources for each session they are observed, or (2) estimate spline parameters or rate parameters for sources chosen to fit the specific variation seen in the position-time series. We found that some strategies improve VLBI EOP accuracy by reducing the biases and weighted root mean square differences between measurements from independent VLBI networks operating simultaneously. We discuss the problem of identifying frequently observed unstable sources and how to identify new sources to replace these unstable sources in the NASA VLBI geodetic source catalog.

Modeling thermal deformation of VLBI antennas with a new temperature model

Temperature variations at very long baseline interferometry (VLBI) sites cause thermal deformations of the VLBI antennas and corresponding displacements of the VLBI reference points. The thermal deformation effects typically contain seasonal and daily signatures. The amplitudes of the annual vertical motion of the antenna reference point can reach several millimeters, depending on the design of the antenna structure, on the material, and on the environmental effects such as global station position, station height and climatology effects. Simple methods to correct this effect use the difference of the environmental temperature with respect to a defined reference temperature, the antenna dimensions, the elevation of the antenna, the material of antenna structure. Applying these simple models for thermal deformation in the VLBI data analysis improves the baseline length repeatability by 3.5%. A comparison of these simple models with local thermal deformation measurements at the antennas in Onsala and Wettzell show that the local measurements and the modeled corrections agree well when the temperature of the antenna structure is used, but agree less good when the surrounding air temperatures are used. To overcome this problem we present a method to model temperature penetration into the antenna structures, that allows to model thermal deformation effects that agree with the observed vertical deformation of the Onsala and Wettzell radio telescopes with a root mean square deviation of 0.07 and 0.13 mm, respectively. Possible implementations in the VLBI analysis are presented, and the definition of an adequate reference temperature is discussed.