Free Access
Volume 529, May 2011
Article Number A70
Number of page(s) 4
Section Celestial mechanics and astrometry
Published online 04 April 2011

© ESO, 2011

1. Introduction

In an earlier paper (Lambert & Le Poncin-Lafitte 2009, hereafter LL09), we estimated the post-Newtonian parameter γ by analyzing group delays recorded by astrometric and geodetic very long baseline interferometry (VLBI) at 8 GHz with an accuracy of 1.5 × 10-4.

In a parallel study, Fomalont et al. (2009) found agreement with General Relativity of 3 × 10-4 using phase-referenced VLBA observations of 3C 279 and 3C 273 at 15, 23, and 43 GHz. Though the accuracy of their measurement does not compete with other studies, the authors pointed out that the precision of the determination of γ could be improved by scheduling specially designed VLBA experiments.

Since then, a number of efforts in the geodesy and astronomy communities led to the inclusion of substantially improved geophysical and astrometric models in the VLBI data reduction. Moreover, we completed the Paris Observatory VLBI analysis center observational database by adding a large number of new observations that had not been processed in LL09, and by upgrading the analysis software to new models. For these reasons, we thought that it was worthwhile to re-launch the LL09 analyses to appreciate the impact of all the improvements listed above.

thumbnail Fig. 1

The main plot displays the observational history of the sources at less than 30° from the Sun (black: observations treated in LL09; red: additional observations of routine experiments not processed in LL09 excluding VLBA+ and VLBA; blue: VLBA+; green: VLBA). The upper plot gives the Sun spot number (SSN, Clette et al. 2007). The right plot displays the deflection angle predicted by General Relativity.

2. Analysis configuration

The session list now includes a number of sessions that were not processed in LL09 (Fig. 1). Most of them are generally designed for geodetic application such as Earth orientation parameter monitoring or station coordinate determination. After 2009, most of the additional sessions are the IVS rapid turn around experiments R1 and R4, scheduled every Monday and Thursday, respectively, that are designed to provide twice weekly Earth orientation parameters and continuity with the NEOS and CORE sessions that were operated before 2002. A relatively small number of additional sessions after 1994 used of the 10-station North American Very Long Baseline Array (VLBA). The observing configuration of the VLBA network allows one to image sources and determine highly accurate station and source positions (see, e.g., Petrov et al. 2009). The VLBA can be used either alone (these sessions will be referred to as VLBA sessions in the following) or together with additional overseas antennas (denoted by VLBA+ sessions in the following) that push baseline lengths to more than 10 000 km. The former category includes the VLBA Calibrator Survey (VCS) programs 1 to 6 (Beasley et al. 2002; Fomalont et al. 2003; Petrov et al. 2005, 2006; Kovalev et al. 2007; Petrov et al. 2008) that were scheduled between 1994 and 2007. They contain observations as close as 1.4° to the Sun, which are indicated as green, vertical bands in Fig. 1. In the latter category, one finds the sessions known as RDV experiments using the VLBA plus up to ten additional geodetic stations located worldwide (in blue in Fig. 1). It appears that VLBA+ sessions stopped observing at less than 15° from the Sun after 2002, like all other routine VLBI experiments. VLBA and VLBA+ sessions usually have a number of observations larger than 10 000 and a postfit rms delay in the range 5–30 ps.

The most complete solution in this study processed 5055 sessions between 3 August 1979 and 30 August 2010, totalling more than 7.3 millions ionosphere-free group delay measurements at 8 GHz. The calculations used the Calc 10.0/Solve 2010.05.21 geodetic VLBI analysis software package1, developed and maintained at NASA Goddard Space Flight Center, and were carried out at the VLBI analysis center of the Paris Observatory (Gontier et al. 2006), which is part of the International VLBI Service for Geodesy and Astrometry (IVS, Schlüter & Behrend 2007). Earth orientation parameters, their rates, and station coordinates were estimated once per session. Loose no-net rotation (NNR) and translation constraints per session were uniformly applied to the positions of all stations excluding Fort Davis (Texas), Pie Town (New Mexico), Fairbanks (Alaska), and the TIGO antenna at Concepción, Chile because of strong non-linear displacements (the latter two sites experienced post-seismic relaxation effects after large earthquakes on the Denali fault in 2003, and between Talca and Concepción in early 2010, respectively). Antenna thermal deformations were mapped using the values of Nothnagel (2009) whereas antenna axis offsets were estimated as global parameters over the full observational time span for a set of 66 stations. The cut-off elevation angle was set to 5°. Zenith wet delays were estimated as a continuous piecewise linear function at 30-min intervals. Troposphere gradients were estimated as 6-hr east and north piecewise functions at all stations except a set of 110 stations with poor observational history.

Station heights were corrected for atmospheric pressure and oceanic tidal loading. The relevant loading quantities were deduced from surface pressure grids from the U. S. NCEP/NCAR reanalysis project atmospheric global circulation model (Kalnay et al. 1996; Petrov & Boy 2004) and from the FES 2004 ocean tide model (Lyard et al. 2004). A critical change with respect to LL09 is the use of atmospheric pressure loading coefficients at diurnal and semi-diurnal frequencies computed with a non-inverted barometer (NIB) hypothesis. In this hypothesis, oceans do not react to any atmospheric pressure variations and any increment in the atmospheric sea level pressure is fully and instantaneously transmitted to the ocean bottom. The NIB ocean is assumed to be static on time scales around and below a few days (e.g., Willebrand et al. 1980; Wunsch & Stammer 1997). In LL09, we assumed that the ocean reacted to balance atmospheric pressure variations at the sea level.

Table 1

Characteristics of the solutions and estimates of γ. Details of the various solutions are given in the text.

We used a priori source coordinates of the second realization of the International Celestial Reference Frame (ICRF2, Fey et al. 2010), adopted by the International Astronomical Union (IAU) in August 2009 as the fundamental realization of the International Celestial Reference System (Feissel & Mignard 1997) and a replacement of the ICRF (Ma et al. 1998). The ICRF2 contains the coordinates of 3414 extragalactic sources determined after VLBI observations at 8 GHz during 1979–2009. With respect to the first ICRF, the noise floor (40 μas) has improved by a factor of five, and the axis stability (10 μas) by a factor of two. Nevertheless, in LL09, we mentioned that the a priori radio source catalogue does not influence the results at a significant level since source coordinates are estimated during the analysis. Therefore, we do not expect the use of the ICRF2 to be the source of the improvement yielded in this paper. As in LL09, we loosely constrained the sources to the ICRF2 catalogue. With respect to LL09, three sources identified as gravitational lenses were removed: 0218+357 and a component 0218+35A (observed in one session each), and PKS 1830–211 (observed in two sessions).

In the new software release, the a priori precession and nutation comply with the IAU 2000/2006 resolutions, which incorporate the nutation model of Mathews et al. (2002) in a way that is consistent with the precession of Capitaine et al. (2003b) and the non-rotating origin-based coordinate transformation between the terrestrial and celestial coordinate systems (Capitaine et al. 2003a).

As a major change with respect to LL09, the Niell (1996) mapping functions (NMF) were replaced by the Vienna mapping functions 1 (VMF1, Böhm et al. 2006), where two of the three coefficients of the continued fraction expressing the hydrostatic mapping functions are computed from the 40-yr reanalysis (ERA–40) of the European Center for Medium-Range Weather Forecast (ECMWF) data. Unlike the previously available mapping functions (Niell 1996; Niell 2000; Böhm & Schuh 2004), they depend on the day of the year, and are no longer symmetric with respect to the equator. Tesmer et al. (2007) showed that VMF1 lead to station height repeatability that is 5.7% better than other mapping functions. They also mentioned that, compared with NMF, the effect on the radio source declinations is  ~10 μas for declinations below  − 20°.

In the next section, we present various solutions designed to test the influence of the changes with respect to LL09 as listed above.

3. Analysis and results

Table 1 summarizes the solution information and displays estimates of γ and the reduced χ2 of the solution. The first line of Table 1 reproduces the most complete solution of LL09. In solution S1, we considered the same session list as in LL09 to check the effects of the improved geophysical and astronomical modeling and the slight changes in the analysis configuration, as listed in Sect. 2. The number of delays in S1 appears to be higher by less than 0.005% than in LL09 because of the use of different versions for a few sessions. It turns out that the formal errors in γ are expected to be comparable, at about 1.5 × 10-4. However, it appears that γ is closer to unity: |γ − 1| was 2 × 10-4 in LL09, and less than 10-4 in S1.

Taken individually, none of the analysis configuration changes listed in the previous section are able to fully explain the small shift in the central value of γ toward unity. It results from the combined effects of all the geophysical and astronomical model improvements.

Solution S2 included additional VLBI sessions between 1979 and 2008 that had not been processed in LL09 (excluding the VLBA and VLBA+ sessions). It therefore checked the effect of completing the observational database, but not extending the observing time span. The next three solutions were based on processed sessions until the end of August 2010. Solution S3 included additional VLBI sessions between 2008 and 2010, and checked the ability of the current geodetic VLBI observations (excluding the VLBA and VLBA+ sessions) to constrain γ if the current observing strategy were continued in the future. As for S1, the estimated values of |γ − 1| in S2 and S3 are well below 10-4. Solution S3 processed about 12% more observations than S2. The formal error decreased by 6%. Routine VLBI sessions scheduled every three days on average are only able to decrease the error in γ at a rate of about 10-6 per year.

Solution S4 is similar to S3, except that it included an additional 72 VLBA+ sessions (scheduled bi-montly in average). The error was reduced to 1.3 × 10-4. Finally, solution S5, which is the most complete solution considered in this study, included both VLBA+ and VLBA sessions. We obtained γ = 0.99992 ± 0.00012. Although the lengths of the S3 and S5 session lists differ by less than 3%, the addition of the rather sparse VLBA sessions increases the number of delays by 30% and the number of sources by a factor of four. Compared with LL09 and other solutions that used routine VLBI experiments only (S1 to S3), solution S5 has the formal error that is smaller by  ~15%, where |γ − 1| is below 10-4, and shows that the VLBA experiments, although scheduled sparsely, have a good potential for General Relativity tests.

4. Concluding remarks

Although we do not challenge the results of Bertotti et al. (2003) from Cassini spacecraft measurements (γ = 1.00002 ± 0.00002), the small improvement presented in this study is notable because it illustrates the capability of certain geodetic/astrometric VLBI networks and observing configurations to increase the sensitivity to γ. The VLBI observational data base provides a very flexible way to test General Relativity in the Solar System at the level of 10-4 thanks to the public availability of the data and the low CPU time taken by the solutions.

This data base is still increasing with new observations of very good quality thanks to a great, joint effort of worldwide radio astronomical observatories and space agencies. The upcoming VLBI 2010 will be designed in particular to reduce systematic errors, including possible source structure corrections thanks to faster antennas, larger networks, and higher data rates resulting in a uv coverage that is much better than in the current geodetic experiments (Petrachenko et al. 2008). This new VLBI network will likely lead to improved ground-based tests of General Relativity (Heinkelmann & Schuh 2010). We therefore encourage VLBI observing program committees to schedule observations of sources close to the Sun as in the VLBA calibrator survey sessions.


We are grateful to Prof. Harald Schuh for his review that helped in improving the paper. This study could not have been carried out without the work of the International VLBI Service for Geodesy and Astrometry (IVS) community that coordinates observations and correlates and stores geodetic VLBI data.


  1. Beasley, A. J., Gordon, D., Peck, A. B., et al. 2002, ApJSS, 141, 13 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bertotti, B., Iess, L., & Tortora, P. 2003, Nature, 425, 374 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  3. Böhm, J., & Schuh, H. 2004, Geowiss. Mitt. 68, Inst. Geod. Geophys., Vienna [Google Scholar]
  4. Böhm, J., Werl, B., & Schuh, H. 2006, J. Geophys. Res., 111, B02406 [NASA ADS] [CrossRef] [Google Scholar]
  5. Capitaine, N., Chapront, J., Lambert, S. B., & Wallace, P. T. 2003a, A&A, 400, 1145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Capitaine, N., Wallace, P. T., & Chapront, J. 2003b, A&A, 412, 567 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Clette, F., Berghmans, P., Vanlommel, P., et al. 2007, Adv. Sp. Res., 40, 919 [Google Scholar]
  8. Feissel, M., & Mignard, F. 1997, A&A, 331, 33 [Google Scholar]
  9. Fey, A. L., Gordon, D. G., Jacobs, C. S., et al. 2010, International Earth Rotation and Reference Systems Service (IERS) Technical Note 35, Bundesamts für Kartographie und Geodäsie, Frankfurt am Main [Google Scholar]
  10. Fomalont, E. B., Petrov, L., MacMillan, D. S., et al. 2003, AJ, 126, 2562 [NASA ADS] [CrossRef] [Google Scholar]
  11. Fomalont, E., Kopeikin, S., Lanyi, G., & Benson, J. 2009, ApJ, 699, 1395 [Google Scholar]
  12. Gontier, A.-M., Lambert, S. B., & Barache, C. 2006, in Proc. Semaine de l’Astrophysique Française, Journées SF2A, ed. D. Barret et al., 27 [Google Scholar]
  13. Heinkelmann, R., & Schuh, H. 2010, in Proc. IAU 261, ed. S. A. Klioner et al., 286 [Google Scholar]
  14. Kalnay, E., Kanamitsu, M., Kistler, R., et al. 1996, BAMS, 77, 437 [Google Scholar]
  15. Kovalev, Y. Y., Petrov, L., Fomalont, E. B., & Gordon, D. 2007, AJ, 133, 1236 [NASA ADS] [CrossRef] [Google Scholar]
  16. Lambert, S. B., & Le Poncin-Lafitte, C. 2009, A&A, 499, 331 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Lyard, F., Lefèvre, F., Letellier, T., & Francis, O. 2006, Ocean Dyn., 56, 394 [NASA ADS] [CrossRef] [Google Scholar]
  18. Ma, C., Arias, E. F., Eubanks, T. M., et al. 1998, AJ, 116, 516 [NASA ADS] [CrossRef] [Google Scholar]
  19. Mathews, P. M., Herring, T. A., & Buffett, B. A. 2002, J. Geophys. Res., 107 [Google Scholar]
  20. Niell, A. E. 1996, J. Geophys. Res., 101, 3227 [Google Scholar]
  21. Niell, A. E. 2000, Earth Plan. Sp., 52, 699 [NASA ADS] [Google Scholar]
  22. Nothnagel, A. 2009, J. Geod., 83, 787 [NASA ADS] [CrossRef] [Google Scholar]
  23. Petrachenko, W. T., Charlot, P., Collioud, A., et al. 2008, in Proc. Journées 2007 Systèmes de Référence Spatio-Temporels, Observatoire de Paris, France, ed. N. Capitaine [Google Scholar]
  24. Petrov, L., & Boy, J.-P. 2004, J. Geophys. Res., 109, 3405 [Google Scholar]
  25. Petrov, L., Kovalev, Y. Y., Fomalont, E. B., & Gordon, D. 2005, AJ, 129, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  26. Petrov, L., Kovalev, Y. Y., Fomalont, E. B., & Gordon, D. 2006, AJ, 131, 1872 [NASA ADS] [CrossRef] [Google Scholar]
  27. Petrov, L., Kovalev, Y. Y., Fomalont, E. B., & Gordon, D. 2008, AJ, 136, 580 [NASA ADS] [CrossRef] [Google Scholar]
  28. Petrov, L., Gordon, D., Gipson, J., et al. 2009, J. Geod., 83, 859 [NASA ADS] [CrossRef] [Google Scholar]
  29. Schlüter, W., & Behrend, D. 2007, J. Geod., 81, 479 [CrossRef] [Google Scholar]
  30. Shapiro, S. S., Davis, J. L., Lebach, D. E., & Gregory, J. S. 2004, Phys. Rev. Lett., 92, 121101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  31. Tesmer, V., Böhm, J., Heinkelmann, R., & Schuh, H. 2007, J. Geod., 81, 409 [NASA ADS] [CrossRef] [Google Scholar]
  32. Willebrand, J., Philander, S., & Pacanowski, R. 1980, J. Phys. Oceanograp., 10, 411 [NASA ADS] [CrossRef] [Google Scholar]
  33. Wunsch, C., & Stammer, D. 1997, Rev. Geophys., 35, 79 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Characteristics of the solutions and estimates of γ. Details of the various solutions are given in the text.

All Figures

thumbnail Fig. 1

The main plot displays the observational history of the sources at less than 30° from the Sun (black: observations treated in LL09; red: additional observations of routine experiments not processed in LL09 excluding VLBA+ and VLBA; blue: VLBA+; green: VLBA). The upper plot gives the Sun spot number (SSN, Clette et al. 2007). The right plot displays the deflection angle predicted by General Relativity.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.