© ESO, 2015

1. Introduction

With geodetic Very Long Baseline Interferometry (VLBI) precise group delays – the difference in arrival times of radio waves at two radio telescopes (Schuh & Behrend 2012) – can be measured, and very accurate radio source positions can be obtained (~40 microarcsec, Ma et al. 2009). There are numerous astrometric and geodetic effects that must be accounted for when reducing high-precision VLBI data. One of these effects is gravitational time delay, which is caused by the gravitational fields of massive objects along the observer’s line-of-sight (Shapiro 1964, 1967) and is responsible for the measured group delays. These massive objects also bend the path of light from distant objects, causing an apparent shift of their position in the sky, as predicted by general relativity (Einstein 1916). Within the solar system both the Sun and Jupiter are massive enough to produce these effects. The conventional equation for calculating gravitational delay is formulated in terms of the positions of the radio telescopes within the barycentric reference frame of the solar system (Kopeikin 1990; Eubanks et al. 1991; Soffel et al. 1991; Klioner 1991), rather than the baseline length between the radio telescopes. Moreover, the total VLBI delay formula includes a coordinate term to transform from the geocentric to the barycentric reference frame.

We propose an alternate gravitational delay formula using a Taylor series expansion. In Sect. 2 we show how the conventional formula can be split into a sum of several terms. All the terms are expressed as a ratio of the baseline length to the Earth barycentric distance. One of the terms links the gravitational delay and the well-known formula for the light deflection angle at an arbitrary elongation from the Sun. Another term cancels the coordinate term explicitly presented at the equation of the total VLBI delay. It is more convenient to formulate the effect of general relativity so that it is free from this complicating coordinate term.

In Sect. 3 a numerical comparison between the new and conventional formulae is presented for the cases of the gravitational field of the Sun and Jupiter. Kaplan (1998) compared the astrometric difference in general relativity effects between the group VLBI delay and optical observations, where the observables are measured in angular quantities. It was claimed that the difference is less than 1 microarcsec using some numerical algorithms. Our analytical approach shows that the two approaches are absolutely equivalent for all elongation angles (the difference is equal to zero), except for the small area around the gravitational body, which is a special case.

In Sect. 4 we show that the formulae for the post-post-Newtonian effect can be developed using the same approach as for the standard monopole light deflection.

In Sect. 5 we discuss the application of the new formula to estimate the proper motion of an extragalactic radio source induced by the solar gravitational field.

2. Developing the conventional formula for gravitational delay

In addition to the three classical tests of general relativity (GR), a fourth test – the delay of a signal propagating in the solar gravitational field – was proposed by Shapiro (1964) and is known as the Shapiro delay. The positions of astronomical instruments on Earth are referenced to the solar system barycentre, and the VLBI delay is equal to the terrestrial time (TT) coordinate time interval between two events of the signal arrival at the first and second radio telescopes (Petit et al. 2010). The difference between the two Shapiro delays as measured with two radio telescopes at the barycentric distances r₁ and r₂ from a body of mass M gives a gravitational delay, τ_grav, which must be taken into account during the standard reduction of the high-precision geodetic VLBI data (Petit et al. 2010), (1)where G is the gravitational constant, s is the unit vector in the direction of observed radio source, c is the speed of light, and γ is the parameter of the parametrised post-Newtonian (PPN) formalism (Will 1993), equal to unity in GR. The classical Shapiro delay formula includes the positions of the remote body and the receiver, but the differential delay (1) only depends on the positions of the two receivers in the adopted reference system. Therefore, the formula (1) is valid for all objects regardless of their distance unless the curvature of the observed light front cannot be ignored. All extragalatic radio sources meet this condition.

We introduce the baseline vector, b, as the difference between the two barycentre radius-vectors of two antennas b = r₂ − r₁. Then formula (1) may be re-written as follows (we recall that γ = 1): (2)For typical baselines, | b | ≪ | r₂ | , and we can use the Taylor series expansion for x ≪ 1 to modify (2) as follows: (3)Formula (3) may be re-written using the approximation ln(1 + y) ≈ y − y²/ 2 for y ≪ 1 because the baseline length is much smaller than the barycentric distance to the points on Earth (4)One has to keep the barycentric coordinates r₂ for the second station in the formulae below to provide consistency with the original gravitational delay model (1).

	Fig. 1 Angles ϕ, ψ, θ, and A, originated from the position of gravitational mass (B), quasar (Q), and baseline vector (b). If the Sun plays the role of the gravitational mass, then the point B in Fig. 1 is also the position of the solar system barycentre.
Open with DEXTER

Converting the three dot products with angles ϕ,ψ,θ (Fig. 1), where (5)and changing | r₂ | to r₂ and | b | to b for simplicity, one can obtain from Eq. (4) (6)Angle θ is given here with respect to the second station, that is, θ = θ₂ for the sake of simplicity. The difference between θ₁ and θ₂ is discussed in Sect. 3. Similar expressions (without terms that include ) have been obtained for geodetic VLBI (e.g. Finkelstein et al. 1983; Hellings 1986) and for SIM interferometry mission model delays (Turyshev 2009).

For the general case of no coplanar vectors (Fig. 1) (7)and formula (6) for arbitrary angles ϕ and θ is given by (8)Using the small-angle approximation (R ≪ r₂), where R is the impact parameter and is equal to the solar radius for the case of grazing light. Therefore (9)and with a “typical” baseline of ~6000 km the first term ~160 ns, the second term is ~1.2 ns, and the third term is ~3 ps for grazing light. As we show below, the third term is less than 1 ps for observations at θ> 1° from the Sun. This term and other small terms in (8) may be neglected because VLBI observations are not undertaken closer than 4° from the Sun. Then, the gravitational delay is given by (10)or, finally, (11)This equation is sufficient for a picosecond level of accuracy. The first term of Eq. (11) can be excluded because it cancels out a similar term from the total geocentric delay (Petit et al. 2010, Eq. (11.9)). Indeed, this formula (11.9) from the IERS Conventions includes the term (12)Recalling that γ = 1, the resulting effect of general relativity, τ_GR, in the total VLBI delay model is therefore given by the sum of the gravitational delay τ_grav and the coordinate term τ_coord(13)The gravitational delay (11) includes the barycentric distance of the second station r₂ in the denominator, whereas the term (12) includes the barycentric distance of the geocentre r. For the case of the Sun, the coordinate term may be as large as 400 ps (see Sect. 2.1) and should be considered carefully, but for Jupiter it does not exceed 1 ps.

The difference between two terms (12) and (11) is less than 0.1 ps, and it may be ignored for the current level of accuracy. Nonetheless, the angles θ and ψ are to be referred to the vector r₂ for correct calculation.

In the original form the gravitational delay is expressed in the BCRS system, whereas the baseline length in (12) is given in the GCRS system. This means that theoretically, one has to convert the baseline length b and r₂ into Eq. (11) from BCRS to GCRS before comparing the coordinate terms from two different parts. However, in practice, the corresponding correction in time delay does not exceed 1 ps for grazing light and a baseline of 6000 km.

After all reductions, the total contribution of the relativistic effects to the total VLBI delay, τ_GR, can be written as follows: (14)For an arbitrary θ the gravitational delay (1) can be split into several terms in Eq. (13), and the main deflection term is given by (15)where α₂ is the deflection angle at arbitrary θ (Shapiro 1967; Ward 1970), and is given by (16)The minor terms in (14) may be ignored for a sufficiently large θ. Therefore (15) may be used as an alternative to the conventional total delay formula for all sources that are far from to a gravitational body. Thus, formulae (15) and (16) present the deflection angle in a more convenient form for estimation from geodetic VLBI data than previously (Treuhaft & Lowe 1991; Kaplan 1998). These formulae also prove that the two approaches – based on angular observations and on the VLBI group delay – are equivalent for all ranges of elongation angles, therefore the difference between two types of quantities is equal to zero.

2.1. Alternative equation for the gravitational fields of the Sun and Jupiter

We consider two realistic cases and compare the alternative formula τ_GR (14) with the conventional formula τ_grav (1). There are two bodies in the solar system that are able to cause a significant impact on the observations – the Sun and Jupiter. The Sun is the more massive of the two bodies, but observations of radio sources closer than 4° are very rare. Jupiter is three orders of magnitude less massive, but observations at its limb (i.e., at an extremely small impact parameter of 70 000 km) are possible. Other massive planets have a weaker effect (Klioner 2003).

The alternative formula (14) consists of three terms: the first term, t₁, provides the main contribution to the GR effects, while the two other terms, t₂ and t₃, contribute much less than t₁. However, we show that for very small angles (θ less than 30′′) all three terms have a similar contribution. At θ = 180° the relativistic contribution, τ_GR, to the total delay is equal to zero, whereas the gravitational delay, τ_grav, is up to 400 ps at a “standard” baseline of ~6000 km as a result of the coordinate term (12) .

We illustrate this in Fig. 2, where the coordinate term (12) and the main deflection term (first term in (14)) are presented as a function of the impact parameter for A = 180°(cosA = 1) and ϕ = 45° (), b = 6000 km, and r = 1 a.u. The coordinate term is always equal to 280 ps because it does not depend on the variable impact parameter θ. The gravitational delay model (1) does not contain the angle θ explicitly, the solid curve in Fig. 2 is calculated numerically. It is equal to zero at θ = 90° and –280 ps at θ = 180°. The main deflection term, t₁ in (14) , is equal to zero at θ = 180, as was previously noted for the relativistic contribution, τ_GR.

	Fig. 2 Gravitational delay (1), coordinate term, and main deflection term (13) with respect to the angle θ (see more details in the text).
Open with DEXTER

	Fig. 3 Difference between formulae (14) and (17) for a baseline of 10 000 km length.
Open with DEXTER

To investigate the accuracy of the alternative formula (14), one has to sum the original gravitational delay formula (1) and the term (12) (17)Figure 3 shows the difference between the delays τ_GR (14) and τ_GR2 (17) for the Sun and a baseline of 10 000 km: the difference between the two formulae only exceeds 1 ps at . Thus, the difference between the two formulae is negligible for the light propagation in the solar system as measured by geodetic VLBI at the current level of precision.

While the two minor terms, t₂ and t₃, contribute less than t₁, they cannot be ignored completely. The contribution of these terms grows rapidly with baseline length (τ ~ b²) and is 72 ps at θ = 1° for b = 10 000 km (Fig. 4).

	Fig. 4 Contribution of the minor terms to the total GR delay for different baseline lengths.
Open with DEXTER

Herewith, we consider two observational VLBI sessions to separately investigate the effect of t₁ separately during a standard 24 h session for the Sun and Jupiter. The VLBI experiment R&D1208 was carried out on 2–3 October, 2012 to track several reference radio sources in close vicinity to the Sun. While the Sun travelled along the ecliptic during the 24 h session, quasar 1243–072 was tracked in the range of angular distances from of to . The four curves in the left-hand of plot of Fig. 5 reflect the variation in t₁ for four selected VLBI baselines of differing length: Kokee – Tsukub32 (KT; 5755 km), HartRAO – Wettzell (HtW; 7832 km), Onsala60 – Wettzell (OW; 920 km), and HartRAO – Onsala60 (HtO; 8525 km) with respect to the angle θ of the approach of the Sun to the radio source 1243-072. This varies steadily due to the relatively slow apparent motion of the Sun.

	Fig. 5 Main term t1 for 1243-072 (left) and for 1922-224 (right) vs. θ.
Open with DEXTER

The largest solar system planet Jupiter approached the radio source 1922-224 on 18–19 November, 2008. This event was observed during the 24 h experiment OHIG60. Four stations – Hobart26, Kokee, Tsukub32, and Parkes participated in tracking the radio source over only 12 h because of observing constraints. The angle θ between the radio source and Jupiter changed from 1.́4 to about 5′ within 12 h. The right plot of Fig. 5 shows the variations in t₁ as a function of θ during this event. This term reaches its maximum of about 400 ps for the longest available baseline (Hobart26 – Tsukub32 (HT) of b = 8088 km) and becomes negligible after a short period of time. For the shorter baselines (Parkes – Tsukub32 (PT; 7233 km) and Hobart26 – Parkes (HP; 1089 km)), this term is lower in magnitude.

	Fig. 6 Variations of the coordinate term for 1243-072 (left) and for 1922-224 (right).
Open with DEXTER

The coordinate term calculated for both events (Fig. 6) is significant for the Sun (reaching 300 ps) but is negligible for Jupiter (less than 0.1 ps). In accordance with Eq. (17), the coordinate term does not depend on the impact parameter and, therefore, it does not grow at small angles, as seen in Fig. 2. All variations in Fig. 6 are caused by the variations in angle ϕ and baseline length.

	Fig. 7 Comparison of the two models for the approach of the Sun to 1243-072 (left) and of Jupiter to 1922-224 (right). The difference for baseline Hobart26 – Parkes in the right plot is not distinctive for this scale.
Open with DEXTER

Figure 7 illustrates the residuals between Eqs. (14) and (17) for these two events. The difference does not exceed 1 ps for either case. This proves that Eq. (14), proposed in this paper, can be used to model the relativistic effects along with the conventional model (17).

3. Small-angle approximation

The monopole effect of the general relativity results in a circular motion of the deflected object (Sazhin et al. 2001; Kopeikin & Makarov 2007). Therefore, the radio sources draw a regular pattern on the sky due to the light deflection. Here we consider an impact of the minor terms in Eq. (14) on the formula for the light deflection, in particular on the intraday time-scale, if a large planet such as Jupiter passes by a reference radio sources very quickly. In the small-angle approximation, Eq. (9), the effect of GR (14) is given by (18)where is the deflection angle for light propagated through a gravitation field (Einstein 1916) for the second station (19)and R₂ = r₂cosθ₂.

Re-arranging formula (18) provides an expression for the deflection angle, , for an arbitrary baseline length in the small-angle approximation as follows: (20)Thus, once in single-telescope mode, only the main angle (19) is presented; in the two-telescope mode, the secondary deflection angle(21)represents the difference between angles and , if and R₁ = r₁cosθ₁. Even though the parallactic shift of the Sun θ₁ − θ₂ is about 8′′ for b = 6000 km, this difference between and , is not negligible in the small-angle approximation. The formula shows that the secondary deflection angle depends on the baseline length, but it reduces to Einstein’s classical expression (19) for a zero-length baseline (i.e. b = 0).

Table 1 shows the meaning of the angle for the Sun and Table 2 for Jupiter with different combinations of the parameters θ₂ and R₂. This effect rapidly grows as b is decreasing and may be detected with a long-baseline radio-interferometer at a very close approach of the Sun (<5°) or large planets (<1′) to a distant radio source.

The angle can change comparatively fast when Jupiter passing by a radio source. An example of the fast variations of is shown in Fig. 8 for the approach of Jupiter and the radio source 1922-224 discussed in Sect. 2.1. The corresponding time delays as shown in the right plot of Fig. 5 are different, but the light deflection angles in Fig. 8. for all baselines are the same.

	Fig. 8 Arc for baselines Hobart26-Tsukub32 (upper left), Kokee-Tsukub32 (upper right), Hobart26-Parkes (bottom left), Kokee-Parkes (bottom right).
Open with DEXTER

4. Post-post-Newtonian effect

As shown by (Klioner 1991), the main term of the post-post-Newtonian effect is given by (22)Using the mathematical transformation of Eqs. (2) and (5), (23)If the baseline length is negligible with respect to the barycentric distance (b ≪ r₂), one could recall and apply the approximation using Eq. (7) to solve the spherical triangle (Fig. 1) and to develop (23) as follows: (24)In the small-angle approximation (R ≪ r₂) (24) is given by (25)For the grazing light the first term in Eq. (25) is about 1 ps, whereas the second one can reach 300 ps at a “standard” baseline of 6000 km (Klioner 1991).

When the first term in Eq. (24) is negligible, we can introduce a new variable (26)and re-write the second term of Eq. (24) as follows: (27)It is obvious that the variable α_PPN plays the role of the post-post-Newtonian deflection angle at arbitrary θ, and converting between the deflection angle and the corresponding part of gravitational delay uses the same factor for the post-Newtonian effect as in Eqs. (15) and (16). It enables us to present the sum of the two delays τ_GR (15) and τ_PPN (27) in the form of two angles, Eqs. (19) and (26), multiplied by , that is, (28)Finally, we recall that the second term of Eq. (24) corresponds to formula (11.14) from the IERS Conventions (Petit et al. 2010).

5. Proper motion of extragalactic radio sources induced by the Sun

The solar gravitational field will result in a proper motion of extragalactic radio sources if the effects of general relativity are not accounted for in VLBI data reduction. The rapid change of the observed positions during a typical 24 h VLBI session may be detected for observations within 10° from the Sun.

The partial derivative is given by (29)For an arbitrary angle θ, the proper motion (where α is from (15)) may be presented by re-arranging Eq. (29) as follows: (30)Figure 9 illustrates Eq. (30) for a range of angles θ from 0° to 10°. It is obvious that all distant radio sources across the sky will display a circular motion if general relativity is not accounted for in VLBI data reduction. The magnitude of the effect is proportional to . In particular, a source near the ecliptic poles θ = 90° will draw a circle of 4 mas in size. The daily displacement between the consecutive days will be equal to . For the Sun, , and the daily displacement is ≈0.07 mas. For all other sources, the daily displacement will be proportional to the factor in Eq. (30).

	Fig. 9 Daily proper motion calculated from Eq. (30) near the Sun.
Open with DEXTER

In the small-angle approxiation, Eq. (30) may be written as follows: (31)For background radio sources that are co-planar with the plane of the body, the apparent motion, μ, rapidly grows from 0.5 μas/year at θ = 10°, to 1.6 μas/year at θ = 5° and 40 μas/year at θ = 1°. Therefore, a systematic pattern of proper motion in a particular area of the sky may serve as an indicator of a hidden mass in this direction.

6. Conclusion

The effects of general relativity are explicitly contained in both components of the total VLBI delay model – the gravitational delay and the coordinate term with a factor of . While the former component uses the individual barycentre positions of the radio telescopes to calculate the effect, the latter component is expressed in terms of the baseline between the radio telescopes. Coupling between both parts has not been investigated until now.

The Shapiro effect as measured by radars and the deflection of light measured with traditional astronomical instruments are considered two independent tests of general relativity. We showed that the total group delay model joins the two tests within one observational technique – geodetic VLBI. The gravitational delay that traditionally originated from the Shapiro effect is linked to the light deflection angle. Therefore, the two approaches, VLBI delay and angular, are absolutely equivalent.

For almost all realistic situations this angle does not depend on the baseline length, thus, a standard geodetic VLBI interferometer acts as a traditional astronomical instrument. In addition, the coordinate term explicitly presented in the conventional total delay model ceases to exist because it is compensated by the same effect in the gravitational delay with opposite sign. Thus, the proposed alternative version of the general relativity contribution to the total VLBI group delay model is free of the coordinate effects. Therefore, the two approaches, time delay and angular, are absolutely equivalent.

The final equation of the general relativity contribution expressed in angular terms also comprises two smaller terms that are significant at very small angular separation between

the deflecting body and distant radio source. These terms may be considered as an increment in the light deflection angle due to the additional time delay for the propagation of the light from station 1 to station 2. This effect becomes significant at .

It is unlikely that the effect will be detected with ground-based interferometers, but it may be detected with space-baseline interferometers, for example, RadioAstron with b up to 300 000 km (Kardashev et al. 2013).

Acknowledgments

We thank Craig Harrison for proofreading the manuscript and valuable comments. We are also grateful to Slava Turyshev, Sergei Kopeikin, Laura Stanford, Sergei Klioner, Benedikt Soja, Natalia Pavlovskaya, Nicole Capitaine, Sebastien Lambert, Adrien Bourgoine, Pierre Teyssandier, Aurelien Hees, Christophe Le Ponchin-Lafitte, Vladimir Zharov and Mihail Sazhin for useful discussion on the mathematical details of general relativity. Oleg Titov was supported by the Russian fund “Dynasty” for two short-term working visits to Saint-Petersburg University in 2011 and 2013. A special “hand made” observing schedule was prepared by Dirk Behrend (GSFC) to track the radio source 1922-224 on 18–19 November, 2008. This paper is published with the permission of the CEO, Geoscience Australia.

References

All Tables

All Figures

	Fig. 1 Angles ϕ, ψ, θ, and A, originated from the position of gravitational mass (B), quasar (Q), and baseline vector (b). If the Sun plays the role of the gravitational mass, then the point B in Fig. 1 is also the position of the solar system barycentre.
Open with DEXTER
In the text

	Fig. 2 Gravitational delay (1), coordinate term, and main deflection term (13) with respect to the angle θ (see more details in the text).
Open with DEXTER
In the text

	Fig. 3 Difference between formulae (14) and (17) for a baseline of 10 000 km length.
Open with DEXTER
In the text

	Fig. 4 Contribution of the minor terms to the total GR delay for different baseline lengths.
Open with DEXTER
In the text

	Fig. 5 Main term t1 for 1243-072 (left) and for 1922-224 (right) vs. θ.
Open with DEXTER
In the text

	Fig. 6 Variations of the coordinate term for 1243-072 (left) and for 1922-224 (right).
Open with DEXTER
In the text

	Fig. 7 Comparison of the two models for the approach of the Sun to 1243-072 (left) and of Jupiter to 1922-224 (right). The difference for baseline Hobart26 – Parkes in the right plot is not distinctive for this scale.
Open with DEXTER
In the text

	Fig. 8 Arc for baselines Hobart26-Tsukub32 (upper left), Kokee-Tsukub32 (upper right), Hobart26-Parkes (bottom left), Kokee-Parkes (bottom right).
Open with DEXTER
In the text

	Fig. 9 Daily proper motion calculated from Eq. (30) near the Sun.
Open with DEXTER
In the text