Geodetic Very Long Baseline Interferometry (VLBI) data are capable of measuring the light deflection caused by the gravitational field of the Sun and large planets with high accuracy. The parameter $\gamma$ of the parametrized Post-Newtonian (PPN) formalism estimated using observations of reference radio sources near the Sun should be equal to unity in the general relativity. We have run several VLBI experiments tracking reference radio sources from 1 to 3 degrees from the Sun. The best formal accuracy of the parameter $\gamma$ achieved in the single-session mode is less than 0.01 percent, or better than the formal accuracy obtained with a global solution included all available observations at arbitrary elongation from the Sun. We are planning more experiments starting from 2020 using better observing conditions near the minimum of the Solar activity cycle.

Geodetic Very Long Baseline Interferometry (VLBI) measures the group delay in the barycentric reference frame. As the Earth is orbiting around the Solar system barycentre with the velocity $V$ of 30 km/s, VLBI proves to be a handy tool to detect the subtle effects of the special and general relativity theory with a magnitude of $(V/\textrm{c})^2$. The theoretical correction for the second order terms reaches up to 300~ps, and it is implemented in the geodetic VLBI group delay model. The total contribution of the second order terms splits into two effects - the variation of the Earth scale, and the deflection of the apparent position of the radio source. The Robertson-Mansouri-Sexl (RMS) generalization of the Lorenz transformation is used for many modern tests of the special relativity theory. We develop an alteration of the RMS formalism to probe the Lorenz invariance with the geodetic VLBI data. The kinematic approach implies three parameters (as a function of the moving reference frame velocity) and the standard Einstein synchronisation. A generalised relativistic model of geodetic VLBI data includes all three parameters that could be estimated. Though, since the modern laboratory Michelson-Morley and Kennedy-Thorndike experiments are more accurate than VLBI technique, the presented equations may be used to test the VLBI group delay model itself.

Context. Radio-loud active galactic nuclei in the early Universe are rare. The quasars J0906+6930 at redshift z=5.47 and J2102+6015 at z=4.57 stand out from the known sample with their compact emission on milliarcsecond (mas) angular scale with high (0.1-Jy level) flux densities measured at GHz radio frequencies. This makes them ideal targets for very long baseline interferometry (VLBI) observations. Aims. By means of VLBI imaging we can reveal the inner radio structure of quasars and model their brightness distribution to better understand the geometry of the jet and the physics of the sources. Methods. We present sensitive high-resolution VLBI images of J0906+6930 and J2102+6015 at two observing frequencies, 2.3 and 8.6 GHz. The data were taken in an astrometric observing programme involving a global five-element radio telescope array. We combined the data from five different epochs from 2017 February to August. Results. For one of the highest redshift blazars known, J0906+6930, we present the first-ever VLBI image obtained at a frequency below 8 GHz. Based on our images at 2.3 and 8.6 GHz, we confirm that this source has a sharply bent helical inner jet structure within ~3 mas from the core. The quasar J2102+6015 shows an elongated radio structure in the east-west direction within the innermost ~2 mas that can be described with a symmetric three-component brightness distribution model at 8.6 GHz. Because of their non-pointlike mas-scale structure, these sources are not ideal as astrometric reference objects. Our results demonstrate that VLBI observing programmes conducted primarily with astrometric or geodetic goals can be utilized for astrophysical purposes as well.

Context. We highlight the capabilities of the geodetic VLBI technique to test General relativity in the classical astrometric style, i.e., measuring the deflection of light in the vicinity of the Sun. Aims. In previous studies, the parameter was estimated by global analyses of thousands of geodetic VLBI sessions. Here we estimate from a single session where the Sun has approached two strong reference radio sources 0229+131 and 0235+164 at an elongation angle of 1-3 degrees. Methods. The AUA020 VLBI session of 1 May 2017 was designed to obtain more than 1000 group delays from the two radio sources. The Solar corona effect was effectively calibrated with the dual-frequency observations even at small elongation from the Sun. Results. We obtained with a precision better than what is obtained through global analyses of thousands of standard geodetic sessions over decades. Current results demonstrate that the modern VLBI technology is capable of establishing new limits on observational test of General Relativity.

We propose an alternative method to detect the secular aberration drift induced by the solar system acceleration due to the attraction to the Galaxy centre. This method is free of the individual radio source proper motion caused by intrinsic structure variation. We developed a procedure to estimate the scale factor directly from very long baseline interferometry (VLBI) data analysis in a source-wise mode within a global solution. The scale factor is estimated for each reference radio source individually as a function of astrometric coordinates (right ascension and declination). This approach splits the systematic dipole effect and uncorrelated motions on the level of observational parameters. We processed VLBI observations from 1979.7 to 2016.5 to obtain the scale factor estimates for more than 4,000 reference radio sources. We show that the estimates highlight a dipole systematics aligned with the direction to the centre of the Galaxy. With this method we obtained a Galactocentric acceleration vector with an amplitude of 5.2 $\pm$ 0.2 \mu as/yr and direction $\alpha_G = 281\deg \pm 3\deg$ and $\delta_G = -35\deg \pm 3\deg$.

The geodetic VLBI technique is capable of measuring the Sun's gravity light deflection from distant radio sources around the whole sky. This light deflection is equivalent to the conventional gravitational delay used for the reduction of geodetic VLBI data. While numerous tests based on a global set of VLBI data have shown that the parameter 'gamma' of the post-Newtonian approximation is equal to unity with a precision of about 0.02 percent, more detailed analysis reveals some systematic deviations depending on the angular elongation from the Sun. In this paper a limited set of VLBI observations near the Sun were adjusted to obtain the estimate of the parameter 'gamma' free of the elongation angle impact. The parameter 'gamma' is still found to be close to unity with precision of 0.06 percent, two subsets of VLBI data measured at short and long baselines produce some statistical inconsistency.

In extending our spectroscopic program, which targets sources drawn from the International Celestial Reference Frame (ICRF) Catalog, we have obtained spectra for about 160 compact flat-spectrum radio sources and determined redshifts for 112 quasars and radio galaxies. A further 14 sources with featureless spectra have been classified as BL Lac objects. Spectra were obtained at three telescopes: the 3.58-m ESO New Technology Telescope (NTT), and the two 8.2-m Gemini telescopes in Hawaii and Chile. While most of the sources are powerful quasars, a significant fraction of radio galaxies is also included from the list of non-defining ICRF radio sources.

High-redshift radio-loud quasars are used to, among other things, test the predictions of cosmological models, set constraints on black hole growth in the early universe and understand galaxy evolution. Prior to this paper, 20 extragalactic radio sources at redshifts above 4.5 have been imaged with very long baseline interferometry (VLBI). Here we report on observations of an additional ten z>4.5 sources at 1.7 and 5 GHz with the European VLBI Network (EVN), thereby increasing the number of imaged sources by 50%. Combining our newly observed sources with those from the literature, we create a substantial sample of 30 z>4.5 VLBI sources, allowing us to study the nature of these objects. Using spectral indices, variability and brightness temperatures, we conclude that of the 27 sources with sufficient information to classify, the radio emission from one source is from star formation, 13 are flat-spectrum radio quasars and 13 are steep-spectrum sources. We also argue that the steep-spectrum sources are off-axis (unbeamed) radio sources with rest-frame self-absorption peaks at or below GHz frequencies and that these sources can be classified as gigahertz peaked-spectrum (GPS) and megahertz peaked-spectrum (MPS) sources.

Some models of the expanding Universe predict that the astrometric proper motion of distant radio sources embedded in space-time are non-zero as the radial distance from observer to the source grows. Systematic proper motion effects would produce a predictable quadrupole pattern on the sky that could be detected using Very Long Baseline Interferometry (VLBI) technique. This quadrupole pattern can be interpreted either as an anisotropic Hubble expansion, or as a signature of the primordial gravitational waves in the early Universe. We present our analysis of a large set of geodetic VLBI data spanning 1979--2015 to estimate the dipole and quadrupole harmonics in the expansion of the vector field of the proper motions of quasars in the sky. The analysis is repeated for different redshift zones.

The goal of this presentation is to report the latest progress in creation of the next generation of VLBI-based International Celestial Reference Frame, ICRF3. Two main directions of ICRF3 development are improvement of the S/X-band frame and extension of the ICRF to higher frequencies. Another important task of this work is the preparation for comparison of ICRF3 with the new generation optical frame GCRF expected by the end of the decade as a result of the Gaia mission.

The Sun's gravitational field deflects the apparent positions of close objects in accordance with the formulae of general relativity. Optical astrometry is used to test the prediction, but only with the stars close to the Sun and only during total Solar eclipses. Geodetic Very Long Baseline Interferometry (VLBI) is capable of measuring the deflection of the light from distant radio sources anytime and across the whole sky. We show that the effect of light deflection is equivalent to the gravitational delay calculated during the reduction of VLBI data. All reference radio sources display an annual circular motion with the magnitude proportional to their ecliptic latitude. In particular, radio sources near the ecliptic pole draw an annual circle with magnitude of 4 mas. This effect could be easily measured with the current precision of the geodetic VLBI data.

We present an alternate formula for calculating gravitational time delay. We use this formula to reduce geodetic Very Long Baseline Interferometry (VLBI) data, taking into account gravitational effects within the solar system, and to test general relativity. The alternate formula was obtained by expanding the conventional formula in a Taylor series. We show that the gravitational delay can be split into several terms including a term due to the coordinate transformation and terms that are explicitly linked to the light deflection angle. Our formula is compared numerically with the conventional formula, and difference in arrival times within 1 ps are found at 1$^\circ$ from the Sun for a full range of baseline lengths. We conclude that the standard reduction of geodetic VLBI data for the effects of general relativity is equivalent to displacing the reference radio sources from their original catalogue positions in accordance with the classical light deflection formula across the whole sky.

We examine the relationship between source position stability and astrophysical properties of radio-loud quasars making up the International Celestial Reference Frame. Understanding this relationship is important for improving quasar selection and analysis strategies, and therefore reference frame stability. We construct light curves for 95 of the most frequently observed ICRF2 quasars at both the 2.3 and 8.4 GHz geodetic VLBI observing bands. Because the appearance of new quasar components corresponds to an increase in quasar flux density, these light curves alert us to potential changes in source structure before they appear in VLBI images. We test how source position stability depends on three astrophysical parameters: (1) Flux density variability at X-band; (2) Time lag between flares in S and X-bands; (3) Spectral index rms, defined as the variability in the ratio between S and X-band flux densities. We find that small time lags between S and X-band light curves, and low spectral index variability, are good indicators of position stability. On the other hand, there is no strong dependence of source position stability on flux density variability in a single frequency band. These findings can be understood by interpreting the time lag between S and X-band light curves as a measure of the size of the source structure. Monitoring of source flux density at multiple frequencies therefore appears to provide a useful probe of quasar structure on scales important to geodesy. We show how multi-frequency flux density monitoring may allow the dependence on frequency of the relative core positions along the jet to be elucidated. Knowledge of the position-frequency relation has important implications for current and future geodetic VLBI programs, as well as the alignment between the radio and optical celestial reference frames. (Abridged)

Very Long Baseline Interferometry is capable of measuring the gravitational delay caused by the Sun and planet gravitational fields. The post-Newtonian parameter $\gamma$ is now estimated with accuracy of $\sigma_{\gamma}=2\cdot 10^{-4}$ using a global set of VLBI data from 1979 to present (Lambert, Gontier, 2009), and $\sigma_{\gamma}=2\cdot10^{-5}$ by the Cassini spacecraft (Bertotti et. al, 2003). Unfortunately, VLBI observations in S- and X-bands very close to the Solar limb (less than 2-3 degrees) are not possible due to the strong turbulence in the Solar corona. Instead, the close approach of big planets to the line of sight of the reference quasars could be also used for testing of the general relativity theory with VLBI. Jupiter is the most appropriate among the big planets due to its large mass and relatively fast apparent motion across the celestial sphere. Six close approaches of Jupiter with quasars in 2014-2016 were found using the DE405/LE405 ephemerides, including one occultation in 2016. We have formed tables of visibility for all six events for VLBI radio telescopes participating in regular IVS programs. Expected magnitudes of the relativistic effects to be measured during these events are discussed in this paper.

We propose new estimates of the secular aberration drift, mainly due to the rotation of the Solar System about the Galactic center, based on up-to-date VLBI observations and and improved method of outlier elimination. We fit degree-2 vector spherical harmonics to extragalactic radio source proper motion field derived from geodetic VLBI observations spanning 1979-2013. We pay particular attention to the outlier elimination procedure to remove outliers from (i) radio source coordinate time series and (ii) the proper motion sample. We obtain more accurate values of the Solar system acceleration compared to those in our previous paper. The acceleration vector is oriented towards the Galactic center within 7 deg. The component perpendicular to the Galactic plane is statistically insignificant. We show that an insufficient cleaning of the data set can lead to strong variations in the dipole amplitude and orientation, and statistically biased results.

We examine the relationship between astrometric stability and astrophysical properties in flat-spectrum radio-loud quasars making up the celestial reference frame. We use position determinations from geodetic VLBI measurements, and develop a new metric for source position stability. We then compare this quantity to two measures of source compactness: structure index, which probes structure on milliarcsecond scales; and the presence of interstellar scintillation associated with the quasar, which probes scales of tens of microarcseconds. We find that persistent scintillators have greater position stability than episodic scintillators, which are in turn more stable than non-scintillators. Scintillators are also more likely to be compact on milliarcsecond scales, as measured by the structure index. Persistent scintillators are therefore excellent candidates for inclusion in any future realisation of the celestial reference frame. A list of these sources is presented in the Appendix. We find that slow (characteristic timescale >3 days) scintillators have more stable positions than rapid scintillators, suggesting they are more compact. High-cadence, long term monitoring is therefore required to identify other members of this population of compact, high brightness temperature quasars.

Continuing our program of spectroscopic observations of ICRF sources, we present redshifts for 120 quasars and radio galaxies. Data were obtained with five telescopes: the 3.58m ESO New Technology Telescope (NTT), the two 8.2m Gemini telescopes, the 2.5m Nordic Optical Telescope (NOT), and the 6.0m Big Azimuthal Telescope (BTA) of the Special Astrophysical Observatory in Russia. The targets were selected from the International VLBI Service for Geodesy & Astrometry (IVS) candidate International Celestial Reference Catalogue which forms part of an observational VLBI program to strengthen the celestial reference frame. We obtained spectra of the potential optical counterparts of more than 150 compact flat-spectrum radio sources, and measured redshifts of 120 emission-line objects, together with 19 BL Lac objects. These identifications add significantly to the precise radio--optical frame tie to be undertaken by Gaia, due to be launched in 2013, and to the existing data available for analysing source proper motions over the celestial sphere. We show that the distribution of redshifts for ICRF sources is consistent with the much larger sample drawn from FIRST and SDSS, implying that the ultra-compact VLBI sources are not distinguished from the overall radio-loud quasar population. In addition, we obtained NOT spectra for five radio sources from the FIRST and NVSS catalogs, selected on the basis of their red colors, which yielded three quasars with z>4.

The AuScope geodetic Very Long Baseline Interferometry array consists of three new 12 m radio telescopes and a correlation facility in Australia. The telescopes at Hobart (Tasmania), Katherine (Northern Territory) and Yarragadee (Western Australia) are co-located with other space geodetic techniques including Global Navigation Satellite Systems (GNSS) and gravity infrastructure, and in the case of Yarragadee, Satellite Laser Ranging (SLR) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) facilities. The correlation facility is based in Perth (Western Australia). This new facility will make significant contributions to improving the densification of the International Celestial Reference Frame in the Southern Hemisphere, and subsequently enhance the International Terrestrial Reference Frame through the ability to detect and mitigate systematic error. This, combined with the simultaneous densification of the GNSS network across Australia will enable the improved measurement of intraplate deformation across the Australian tectonic plate. In this paper we present a description of this new infrastructure and present some initial results, including telescope performance measurements and positions of the telescopes in the International Terrestrial Reference Frame. We show that this array is already capable of achieving centimetre precision over typical long-baselines and that network and reference source systematic effects must be further improved to reach the ambitious goals of VLBI2010.

The centrifugial acceleration of the Solar system, resulting from the gravitational attraction of the Galaxy centre, causes a phenomenon known as 'secular aberrration drift'. This acceleration of the Solar system barycentre has been ignored so far in the standard procedures for high-precision astrometry. It turns out that the current definition of the celestial reference frame as epochless and based on the assumption that quasars have no detectable proper motions, needs to be revised. In the future, a realization of the celestial reference system (realized either with VLBI, or GAIA) should correct source coordinates from this effect, possibly by providing source positions together with their proper motions. Alternatively, the galactocentric acceleration may be incorporated into the conventional group delay model applied for VLBI data analysis.

We present the results of spectroscopic observations of the optical counterparts of 47 southern radio sources from the candidate International Celestial Reference Catalogue (ICRC), as part of a Very Long Baseline Interferometry (VLBI) program to strengthen the celestial reference frame, especially in the south. The observations were made with the 3.58-meter European Southern Observatory New Technology Telescope (NTT). We obtained redshifts for 30 quasars and one radio galaxy, with a further 7 objects being probable BL Lac objects with featureless spectra. Of the remainder, four were clear misidentifications with Galactic stars and five had low signal-to-noise spectra and could not be classified. These results, in combination with new VLBI data of the radio sources with redshifts more than 2, add significantly to the existing data needed to refine the distribution of source proper motions over the celestial sphere.

While analyzing decades of very long baseline interferometry (VLBI) data, we detected the secular aberration drift of the extragalatic radio source proper motions caused by the rotation of the Solar System barycenter around the Galactic center. Our results agree with the predicted estimate to be 4-6 micro arcseconds per year ({\mu}as/yr) towards {\alpha} = 266\circ and {\delta} = -29\circ. In addition, we tried to detect the quadrupole systematics of the velocity field. The analysis method consisted of three steps. First, we analyzed geodetic and astrometric VLBI data to produce radio source coordinate time series. Second, we fitted proper motions of 555 sources with long observational histories over the period 1990-2010 to their respective coordinate time series. Finally, we fitted vector spherical harmonic components of degrees 1 and 2 to the proper motion field. Within the error bars, the magnitude and the direction of the dipole component agree with predictions. The dipole vector has an amplitude of 6.4 \pm 1.5 {\mu}as/yr and is directed towards equatorial coordinates {\alpha} = 263\circ and {\delta} = -20\circ. The quadrupole component has not been detected. The primordial gravitational wave density, integrated over a range of frequencies less than 10-9 Hz, has a limit of 0.0042 h-2 where h is the normalized Hubble constant is H0/(100 km s-1).

Many of the compact extragalactic radio sources that are used as fiducial points to define the celestial reference frame are known to have proper motions detectable with long-term geodetic/astrometric Very Long Baseline Interferometry (VLBI) measurements. These changes can be as high as several hundred micro-arcseconds per year for certain objects. When imaged with VLBI at milli-arcsecond (mas) angular resolution, these sources (radio-loud active galactic nuclei) typically show structures dominated by a compact, often unresolved "core" and a one-sided "jet". The positional instability of compact radio sources is believed to be connected with changes in their brightness distribution structure. For the first time, we test this assumption in a statistical sense on a large sample, not only for individual objects. We investigate a sample of 62 radio sources for which reliable long-term time series of astrometric positions as well as detailed 8-GHz VLBI brightness distribution models are available. We compare the characteristic direction of their extended jet structure and the direction of their apparent proper motion. We present our data and analysis method, and conclude that there is indeed a correlation between the two characteristic directions. However, there are cases where the ~1-10-mas scale VLBI jet directions are significantly misaligned with respect to the apparent proper motion direction.

The gravitational attraction of the Galactic centre leads to the centrifugal acceleration of the Solar system barycentre. It results in secular aberration drift which displaces the position of the distant radio sources. The effect should be accounted for in high-precision astrometric reductions as well as by the corresponding update of the ICRS definition.

From observations with the 6-m BTA telescope at SAO RAS, we have determined spectroscopic redshifts of seven optical objects whose coordinates coincide with those of radio sources from the list of IVS (International VLBI Service for Geodesy and Astrometry). When compared to radio data, the obtained spectra and redshifts provide evidence for reliable identification of four observed objects; the other three require further study. The distances to the sources derived from our measurements will make it possible to refine the current estimates for parameters of cosmological models based on proper motions of these objects, which are determined from geodetic VLBI observations.

We report on recent VLBA/VERA/IVS observational tests of General Relativity. First, we will summarize the results from the 2005 VLBA experiment that determined gamma with an accuracy of 0.0003 by measuring the deflection of four compact radio sources by the solar gravitational field. We discuss the limits of precision that can be obtained with VLBA experiments in the future. We describe recent experiments using the three global arrays to measure the aberration of gravity when Jupiter and Saturn passed within a few arcmin of bright radio sources. These reductions are still in progress, but the anticipated positional accuracy of the VLBA experiment may be about 0.01 mas.

A new list of physical characteristics of 4261 astrometric radio sources, including all 717 ICRF-Ext.2 sources has been compiled. Comparison of our data of optical characteristics with the official International Earth Rotation and Reference Systems Service (IERS) list showed significant discrepancies for about half of 667 common sources. We also found that asymmetry in the radio sources distribution between hemispheres could cause significant correlation between the vector spherical harmonics, especially if the case of sparse distribution of the sources with high redshift. We identified radio sources having many-year observation history and lack redshift. This sources should be urgently observed at large optical telescopes. The list of optical characteristics created in this paper is recommended for use as a supplement material for the next International Celestial Reference Frame (ICRF) realization. It can be also effectively used for cosmological studies and planning of observing programs both in radio and optics.

A new list of physical characteristics of 3914 astrometric radio sources, including all 717 ICRF-Ext.2 sources, observed during IVS and NRAO VCS sessions have been compiled. The list includes source type, redshift and visual magnitude (if available). In case of doubt detailed comment is provided. The list of sources with their positions was taken from the Goddard VLBI astrometric catalog with addition of two ICRF-Ext.2 sources. At this stage the source characteristics were mainly taken from the NASA/IPAC Extragalactic Database (NED). 667 sources from our list are included into the IERS list. Comparison has shown a significant difference in characteristics for about half of these 667 common sources. We compiled a list of frequently observed sources without known physical characteristics for urgent optical identification and spectrophotometric observations with large optical telescopes. This presented list of physical characteristics can be used as a supplement material for the ICRF-2, as well as a database for kinematic studies of the Universe and other related works, including scheduling of dedicated IVS programs.

Some models of the expanding Universe predict that the astrometric proper motion of distant radio sources embedded in space-time are non-zero as radial distance from observer to the source grows. Systematic effects due to this proper motion can even increase with distance making possible to measure them with high precision astrometric techniques like VLBI. We analyzed a large set of geodetic VLBI data spanning from 1979 till 2008 to estimate the dipole and the quadrupole harmonics in the expansion of the vector field of the proper motions of quasars in the sky. We estimated the vector spherical harmonics (three parameters for the dipole and ten - for the quadrupole systematic) directly from the VLBI group delays without intermediate calculation of the individual proper motion. The estimates have been obtained separately for different red shift zones. It was shown that the dipole harmonic does not vary significantly, whereas the amplitude of the quadrupole gradually increases with the red shift. This quadrupole pattern can be interpreted either as an anisotropic Hubble expansion, or as an indication of the primordial gravitational waves in the early Universe. However, more prosaic explanations may also be possible.

The motion of relativistic jets from the active extragalactic nuclei can reach several hundred microseconds per year and mimic proper motion of the distant radio sources observed by VLBI. Such motion of individual quasar is not correlated and its magnitude exceeds the small systematic effects induced by the rotation of the Solar system around the centre of the Galaxy. In this paper we search for the cause of the systematic effect and discuss the results.

The galactocentric rotation of the Solar system generates the systematic effect in proper motions known as 'secular aberration drift'. This tiny effect (about five microseconds per year) in the quasar proper motion can be measured by VLBI. However, the motions of relativistic jets from the active extragalactic nuclei can reach several hundred microseconds per year and mimic the proper motion of the observed radio sources. These apparent motions exceed the secular aberration drift by factor of 10-100. In this paper we search the ways to overcome the difficulties and discuss our estimates of the secular aberration drift using OCCAM software.