## Observations of radio sources near the Sun

O. Titov, S. Lambert, B. Soja, F . Shu, A. Melnikov, J. McCallum, L. McCallum, M. Schartner, A. de Witt, D. Ivanov, A. Mikhailov, S. O. Yi, W. Chen, B. Xia, M. Ishigaki, S. Gulyaev, T. Natusch, S. Weston

Abstract Geodetic Very Long Baseline Interferometry (VLBI) data are capable of measuring the light deflection caused by the gravitational field of the Sun and

Oleg Titov

Geoscience Australia, PO Box 378 Canberra 2601 Australia

S´ ebastien Lambert

SYRTE, Observatoire de Paris, Universit´ e PSL, CNRS, Sorbonne Universit´ e, LNE, Paris, France

Benedikt Soja

Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA

Fengchun Shu, Wen Chen abd Bo Xia

Shanghai Astronomical Observatory, 80 Nandan Road, Shanghai, 200030, China large planets with high accuracy. The parameter g 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 g 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.

Alexei Melnikov, Dmitrii Ivanov and Andrei Mikhailov

Institute of Applied Astronomy, Kutuzov Embankment, 10, Saint-Petersburg, 191187, Russia

Jamie McCallum and Lucia McCallum

University of Tasmania, Private Bag 37, Hobart, Tasmania, 7001, Australia

Aletha de Witt

Hartebeesthoek Radio Astronomy Observatory, PO Box 443, Krugersdorp, 1740, South Africa

Matthias Schartner

Department of Geodesy and Geoinformation, Research Group Advanced Geodesy, TU Wien, Gusshausstraße 27-29/E120.4, Wien-1040, Austria

Sang Oh Yi

National Geographic Information Institute, Space Geodetic Observatory, Sejong, PO Box 30060, South Korea

Masafumi Ishigaki

Geospatial Information Authority of Japan 1, Kitasato, Tsukuba 305-0811, Japan

Sergei Gulyaev, Tim Natusch and Stuart Weston

Institute for Radio Astronomy and Space Research, Auckand University of Technology, Auckland, 1010, New Zealand

Keywords VLBI, general relativity, ionosphere

## 1 Introduction

In accordance with General Relativity the radio waves slow down due to the gravitational potential of the Sun (the so-called Shapiro effect; see Shapiro, 1964, 1967), making very long baseline interferometry (VLBI) a useful tool for testing General Relativity by means of the parameterized post-Newtonian (PPN) formalism (Will, 1993). Nevertheless, the accuracy of the PPN parameter g obtained from absolute or differential VLBI observations (Fomalont et al, 2009; Lambert and Le Poncin-Lafitte, 2009, 2011) remains worse than the current best limit of ( 2 1 . ± 2 3 . ) × 10 -5 based on Cassini radio science experiment (Bertotti et al, 2003) by an order of magnitude. The upper limits on the parameter g have been improved substantially in the

past 30 years (Robertson and Carter, 1984; Robertson et al, 1991; Lebach et al, 1995; Fomalont and Kopeikin, 2003), but some authors (Shapiro et al, 2004; Lambert and Le Poncin-Lafitte, 2009) found degradation in the estimates of g with elongation, and suggested that this systematic effect may limit the improvement in the VLBI-derived g upper limits, despite the dramatic growth in the number of observations in recent decades.

The current paper focuses on radio source approaches at angular distances less than three degrees from the centre of the Sun in order to measure the light deflection effect at the highest magnitude and, thus, to avoid a possible bias caused by observations at larger elongations. We report on two special VLBI sessions, AUA020 (May 2017) and AOV022 (May 2018), on the single-session estimates of g .

## 2 Data

A dedicated geodetic VLBI experiment (AUA020, 0102 May, 2017, part of AUSTRAL program) was scheduled to probe the gravitational delay effect using a network of seven radio telescopes (Svetloe, Zelenchukskaya, Badary, HartRAO, Seshan25, Sejong, and Hobart26). Two radio sources 0229+131 and 0235+164 were observed at range of angular distances from 1.15 ◦ to 2.6 ◦ from the Sun. The position of both radio sources with respect to the Sun at the start of the experiment is shown on Fig. 1. A serious issue in such a configuration is the solar thermal noise that penetrates to the signal through the side lobes, and could cause loss of data due to striking the signal-to-noise ratio. To overcome the problem, one has to

1) select strong radio sources with larger correlated flux density in both frequency bands,

2) use large radio telescopes with narrow side lobes and better sensitivity, and

3) use the highest possible data rate recording (e.g., 1 Gbps) to gain a better signal-to-noise ratio during the same integration time.

More details about the schedule design are publsihed in (Titov et al, 2018).

The target radio source 0229+131 is a defining source of the ICRF2 whose position is given with an accuracy close to the ICRF2 noise floor of 40 m as. The position of the second target 0235+164 is less accurate

Fig. 1 Geometry of the radio sources close to the sun at the start (Left) and at the end (Right) of VLBI session AUA020 with respect to a LASCO C2 image of the solar corona. The Sun is hidden behind the occultation disc of the coronagraph, with the white inner circle representing the limb of the Sun. The field-ofview is 1.5 degrees elongation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00973v1-figure-24769231-39db-4bdc-b70c-373950dab85c.jpg)

by a factor of five but still at the level of the ICRF2 median error and largely below the millisecond of arc. Both sources are compact and their structure indices measured at the time of the ICRF2 work were of 2.4 and 1.3, respectively, ensuring a structure delay lower than 2 ps (Fey and Charlot, 1997).

## 3 Analysis and Results

For purpose of cross-checking the results and testing their robustness, we processed the VLBI session AUA020 within two independent teams with two independent geodetic VLBI analysis software packages. The duplication of the analyses with two software packages also allows to use some specific options that are available on only one of them. The first analysis package is OCCAM (Titov et al., 2004) that implements the least-squares collocation method (Titov, 2000) for calibrating the wet troposphere fluctuations, and to account for the mutual correlations between observables. The second one is Calc/Solve (Ma et al., 1986), developed and maintained by the geodetic VLBI group at NASA GSFC, that uses classical least-squares. More details about the data analysis design are discussed in (Titov et al, 2018)

Table 1 shows the results of the AUA020 experiment data analysis. Uncertainties on g lie between 0 9 . × 10 -4 and 4 × 10 -4 . Our estimates appear therefore as precise as that obtained from global solutions using thousands of VLBI experiments (Lambert and Le Poncin-Lafitte, 2009, 2011). The formal error is about

two times lower when g is fitted to the observations of the radio source that is two times closer to the centre of the Sun (0235+164) than to that of its counterpart (0229+131). Using all scans returns a result similar to using only scans relevant to 0229+131 and 0235+164, confirm that only sources a low elongation can efficiently constrain the PPN parameter. Solutions from both software packages are consistent within the standard errors. The difference of postfit rms between OCCAM and Calc/Solve might find its origin in the different modeling of the nuisance parameters (stochastic versus CPWL function). No large systematics are detected except a 2 7 . s deviation in the Calc/Solve solution when only 0229+131 is used and whose origin is unclear: as both solutions started from the same a priori, the issue could rather be in the estimation method or in the handling of troposphere/clock parameters.

It appears that during session AUA020, data in three channels at Sejong station were lost due to technical reasons. Therefore, we reprocessed the previous analyses after downweighting (but not suppressing) Sejong data. (We could test this option with OCCAM only since Calc/Solve does not handle downweighting.) The postfit rms of the solution is significantly lowered. The formal error on g is marginally lowered down to 9 × 10 -5 .

For purpose of comparison of the AUA020 session with other standard geodetic VLBI sessions, we estimated g with Calc/Solve using the parameterization described above for each of sessions of the full geodetic VLBI data base made available by the International VLBI Service for geodesy and astrometry (IVS) since 1979 (at the exclusion of intensive sessions). The median postfit rms is 27 ps that is close to the postfit rms of the AUA020 session. The distribution of the obtained values of g -1 is shown in Fig. 2 along with distributions of errors and normalized estimates. The distribution of errors in log-scale is slightly asymmetric, exhibiting a 'tail' on its right side that might traduce results from sessions not designed for precise astrometry. Nevertheless, assuming a Gaussian shape, the logscaled distribution peaks at 10 -2 with a s of ∼ 0 5. . This makes the error estimate from AUA020, that is two orders of magnitude less, somewhat 'outstanding'. The bottom-right panel of Fig. 2 shows that the major part of the sessions does not bring severe systematics, the estimates of g being unity within the error bars; session AUA020 is part of the session group that presents the lowest systematics.

Table 1 Estimates of g -1 for the session AUA020, in unit of 10 -4 , along with the session c 2 and the postfit rms delay r in ps.

|       |                          | g - 1 10 - 4             | s g 10 - 4               | c 2                      | r                        |
|-------|--------------------------|--------------------------|--------------------------|--------------------------|--------------------------|
| OCCAM | All stations             | All stations             | All stations             | All stations             | All stations             |
| OCCAM | All scans                | 0.56                     | 1.15                     | 0.34                     | 28                       |
| OCCAM | 0235+164                 | 1.34                     | 1.58                     | 0.34                     | 28                       |
| OCCAM | 0229+131                 | -1.54                    | 3.41                     | 0.34                     | 28                       |
| OCCAM | Both                     | 0.53                     | 1.14                     | 0.34                     | 28                       |
| OCCAM | With Sejong downweighted | With Sejong downweighted | With Sejong downweighted | With Sejong downweighted | With Sejong downweighted |
| OCCAM | All scans                | 0.91                     | 0.94                     | 0.27                     | 21                       |
| OCCAM | 0235+164                 | 1.64                     | 1.29                     | 0.27                     | 21                       |
| OCCAM | 0229+131                 | 0.32                     | 2.83                     | 0.27                     | 21                       |
| OCCAM | Both                     | 0.89                     | 0.94                     | 0.27                     | 21                       |
|       | All sources              | -0.22                    | 1.10                     | 0.84                     | 26                       |
|       | 0235+164                 | 1.85                     | 1.48                     | 0.84                     | 26                       |
|       | 0229+131                 | -6.84                    | 2.53                     | 0.84                     | 26                       |
|       | Both                     | -0.26                    | 1.09                     | 0.84                     | 26                       |

Fig. 2 Distributions of (Top) estimates of g -1, (Bottom-left) their formal errors, and (Bottom-right) normalized estimates of g -1 for all of the geodetic VLBI sessions. The vertical, red bar stands for the results of the AUA020 session.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00973v1-figure-8548a6b1-8138-4765-a881-4a2b7d0fa370.jpg)

Another experiment (AOV022) was undertaken on 01-02 May, 2018 with ten radio telescopes (Svetloe, Zelenchukskaya, Badary, Hobart26, Seshan25, Kunming, Ishioka, Yarragadee, Katherine, Warkworth). The same radio sources (0229+131 and 0235+164) were scheduled with the same strategy. The statistics of the result was found to be 2-3 times worse than from AUA020, presumably, due to severe source structure delay effect. Fig. 3 shows the post-fit residuals of for radio source 0229+131 for two baselines, IshiokaSeshan25 and Badary-Zelenchk. The variations of the

## 0229+131, Ishioka Seshan25

Fig. 3 Post-fit residuals of 0229+131 with baselines IshiokaSeshan25 and Badary-Zelenchk.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00973v1-figure-ab7017e9-913c-4933-8c65-52b1d661b2bf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00973v1-figure-a9039657-e546-4efd-95d2-60852c8f5a57.jpg)

residuals are consistent to the variations induced by the source structure (Titov and Lopez, 2018) therefore, we believe, that more detailed analysis is required to obtain a better statistic for this experiment.

## 4 Conclusions and outlook

In this paper we discuss our results on testing of general relativity with geodetic VLBI using the close approach of the Sun to the reference radio sources. It was a general misconception in the past that the effect of the plasma of the solar corona completely disturbs the interferometric responce for light rays passing within several degrees from the Sun. We proved that these per- turbations are at an acceptable level unless the signals pass active streamers in the solar corona. Therefore, the standard dual frequency calibration facilitates the stochastic noise induced by charged particles the solar corona by the same way as for the ionosphere around the Earth.

While the systematic effects based on radial or dipole models of the corona appear to be negligible, individual group delay observations are affected by random scatter caused by small-scale coronal structure and temporal variations thereof. Since these perturbations do not systematically affect the observations, we assume that they cancel out over the period of observations (17 hours with observations angularly close to the Sun). Since the ray paths to the radio sources 0235+164 and 0229+131 within small solar elongation happened to be in quiet regions (cf. Fig. 1), the scatter was small enough that precise group delays could be successfully determined at such small elongations.

The major source of stochastic noise in VLBI measurements resides in the unknown wet troposphere delay. The difference between VLBI estimates of the wet troposphere delay and independent radiometer data appears to stay within 3 mm, or 10 ps (Titov and Stanford, 2013) suggesting that the impact of the wet troposphere delay on the astrometric light deflection angle estimation near the Sun is negligible.

Overall, a total improvement of the uncertainty on g by a factor of ten is expected, enabling to challenge the current limit imposed by the Cassini radio science experiment of Bertotti et al (2003), although the Gaia astrometry on Solar system objects is expected to deliver an accuracy of 10 -6 (Mignard and Klioner, 2009).

## 5 Acknowledgments

This paper is published with the permission of the CEO, Geoscience Australia. B. Soja's research was supported by an appointment to the NASA Postdoctoral Program, administered by Universities Space Research Association, at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with National Aeronautics and Space Administration. We are grateful to D. Gordon (GSFC) for post-processing reduction of the AUA020 data.

## References

| Bertotti B, Iess L, Tortora P (2003) Atest of general relativity us- ing radio links with the Cassini spacecraft. Nature 425:374- 376                                                                                             |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Charlot P, Jacobs CS, Gordon D, et al. (2019) The third realiza- tion of the International Celestial Reference Frame by very long baseline interferometry. Astronomy &Astrophysics (in preparation)                               |
| Fey AL, Charlot P (1997) VLBA Observations of Radio Ref- erence Frame Sources. II. Astrometric Suitability Based on Observed Structure. Astrophysical Journal 111:95-142                                                          |
| Fey AL, Gordon D, Jacobs CS, Ma C, et al. (2015) The Second Realization of the International Celestial Reference Frame by Very Long Baseline Interferometry. Astronomical Jour- nal 150:58                                        |
| Measurements of the Gravitational Bending of Radio Waves Using the VLBA. Astrophysical Journal 699:1395-1402 Fomalont EB, Kopeikin SM (2003) The measurement of the light deflection from jupiter: Experimental results. The As-  |
| trophysical Journal 598(1):704 Lambert SB, Le Poncin-Lafitte C (2009) Determining the rela- tivistic parameter g using very long baseline                                                                                         |
| interferometry. Astronomy &Astrophysics 499:331-335 Lambert SB, Le Poncin-Lafitte C (2011) Improved determina-                                                                                                                    |
| tion of g by VLBI. Astronomy &Astrophysics 529:A70 Lebach DE, Corey BE, Shapiro II, et al. (1995) Measurement of the solar gravitational deflection of radio waves using long-baseline interferometry. Phys Rev Lett 75:1439-1442 |
| very- Ma, C., Clark, T. A., Ryan, J. W., et al. 1986, Astronomical Jour- nal, 92, 1020                                                                                                                                            |
| Mignard F, Klioner S (2009) Gaia: Relativistic modelling and testing. In: Klioner S, Seidelman PK (eds) Relativity in Fun- damental Astronomy - Proceedings IAU Symposium No. 261, Cambridge University Press, vol 5, pp 306-314  |
| Robertson DS, Carter WE(1984) Relativistic deflection of radio signals in the solar gravitational field measured with vlbi. Nature 310:572                                                                                        |
| Robertson DS, Carter WE, Dillinger WH (1991) New ment of solar gravitational deflection of radio signals                                                                                                                          |
| vlbi. Nature 349:768 Shapiro II (1964) Fourth Test of General Relativity. Phys                                                                                                                                                    |
| Shapiro II (1967) New method for the detection of light tion by solar gravity. Science 157(3790):806 Shapiro SS, Davis JL, Lebach DE, Gregory JS (2004)                                                                           |
| surement of the Solar Gravitational Deflection of Waves using Geodetic Very-Long-Baseline                                                                                                                                         |
| Titov, O. 2000, IERS Technical Note, 28, 11 O., Tesmer, V., & Boehm, J. 2004, in International VLBI Service for Geodesy and Astrometry 2004 General Meeting                                                                       |
| Titov, Proceedings, ed. N. R. Vandenberg &K. D. Baver, 267                                                                                                                                                                        |
| Titov O, Stanford L (2013) Comparison of wet troposphere vari- ations estimated from VLBI and WVR. In: Zubko N, Pouta- nen M(eds) 21st Meeting of the European VLBI Group for                                                     |
| Data, 1979 1999. Physical Review Letters 92(12):121101                                                                                                                                                                            |
| Interferometry                                                                                                                                                                                                                    |
| Radio                                                                                                                                                                                                                             |
| Rev Lett 13(26):789-791 deflec- Mea-                                                                                                                                                                                              |
| measure- using                                                                                                                                                                                                                    |

Geodesy and Astronomy, Reports of the Finnish Geodetic Institute, p. 151-154

- Titov O, Girdiuk A, Lambert SB, et al. (2018) Testing general relativity with geodetic VLBI. What a single, specially designed experiment can teach us. Astronomy &amp; Astrophysics 618:A8
- Titov O, Lopez Yu (2018) Two-component structure of the radio source 0014+813 from VLBI observations within the CONT14 program. Astronomy Letters, 44, 139

Will CM (1993) Theory and Experiment in Gravitational Physics
Oleg Titov Sebastien Lambert Benedikt Soja Fengchun Shu Alexei Melnikov Jamie McCallum Lucia McCallum Matthias Schartner Aletha de Witt Dmitrii Ivanov Andrei Mikhailov Sang Oh Yi Wen Chen Bo Xia Masafumi Ishigaki Sergei Gulyaev Tim Natusch Stuart Weston
Observations of radio sources near the Sun
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.