# The Sun as a planet-host star: proxies from SDO images for HARPS radial-velocity variations

R. D. Haywood, A. Collier Cameron, Y. C. Unruh, C. Lovis, A. F. Lanza, J.

Llama, M. Deleuil, R. Fares, M. Gillon, C. Moutou, et al.

► To cite this version:

R. D. Haywood, A. Collier Cameron, Y. C. Unruh, C. Lovis, A. F. Lanza, et al. The Sun as a planet-host star: proxies from SDO images for HARPS radial-velocity variations. Monthly Notices of the Royal Astronomical Society, 2016, 457 (4), pp.3637-3651. 10.1093/mnras/stw187. hal-01440671

HAL Id: hal-01440671

https://hal.science/hal-01440671v1

Submitted on 4 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# The Sun as a planet-host star: proxies from SDO images for HARPS radial-velocity variations

R. D. Haywood,¹,²* A. Collier Cameron,¹ Y. C. Unruh,³ C. Lovis,⁴ A.F. Lanza,⁵

J. Llama,¹,⁶ M. Deleuil,⁷ R. Fares,¹,⁵ M. Gillon,⁸ C. Moutou,⁷ F. Pepe,⁴ D. Pollacco,⁹

D. Queloz⁴ and D. Ségransan⁴

¹SUPA, School of Physics and Astronomy, University of St Andrews, St Andrews, Fife KY16 9SS, UK

²Harvard–Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

³Astrophysics Group, Blackett Laboratory, Imperial College London, London SW7 2AZ, UK

⁴Observatoire de Genève, 51 Ch. des Maillettes, CH-1290 Sauverny, Switzerland

⁵INAF-Osservatorio Astrofisico di Catania, via S. Sofia 78, I-95123 Catania, Italy

⁶Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ 86001, USA

⁷Aix Marseille Université, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, F-13388 Marseille, France

⁸Institut d'Astrophysique et de Géophysique, Université de Liège, Allée du 6 août 17, Bat. B5C, B-4000 Liège, Belgium

⁹Department of Physics, University of Warwick, Coventry CV4 7AL, UK

Accepted 2016 January 20. Received 2016 January 15; in original form 2015 August 27

## ABSTRACT

The Sun is the only star whose surface can be directly resolved at high resolution, and therefore constitutes an excellent test case to explore the physical origin of stellar radial-velocity (RV) variability. We present HARPS observations of sunlight scattered off the bright asteroid 4/Vesta, from which we deduced the Sun's activity-driven RV variations. In parallel, the Helioseismic and Magnetic Imager instrument on board the *Solar Dynamics Observatory* provided us with simultaneous high spatial resolution magnetograms, Dopplergrams and continuum images of the Sun in the Fe I 6173 Å line. We determine the RV modulation arising from the suppression of granular blueshift in magnetized regions and the flux imbalance induced by dark spots and bright faculae. The rms velocity amplitudes of these contributions are 2.40 and 0.41 m s⁻¹, respectively, which confirms that the inhibition of convection is the dominant source of activity-induced RV variations at play, in accordance with previous studies. We find the Doppler imbalances of spot and plage regions to be only weakly anticorrelated. Light curves can thus only give incomplete predictions of convective blueshift suppression. We must instead seek proxies that track the plage coverage on the visible stellar hemisphere directly. The chromospheric flux index $R_{HK}$ derived from the HARPS spectra performs poorly in this respect, possibly because of the differences in limb brightening/darkening in the chromosphere and photosphere. We also find that the activity-driven RV variations of the Sun are strongly correlated with its full-disc magnetic flux density, which may become a useful proxy for activity-related RV noise.

**Key words:** techniques: radial velocities – Sun: activity – Sun: faculae, plages – Sun: granulation – sunspots – planets and satellites: detection.

## 1 INTRODUCTION

The surface of a star is constantly bustling with magnetic activity. Starspots and faculae/plage¹ inhibit convective motions taking place

at the stellar surface, thus suppressing part of the blueshift naturally resulting from granulation. In addition, dark starspots (and bright faculae, to a lesser extent) coming in and out of view as the star rotates induce an imbalance between the redshifted and blueshifted halves of the star, which translates into a radial-velocity (RV) variation. Stellar activity, through the perturbation of photospheric convection, induces RV variations on much longer time-scales of the order of years, in tune with their magnetic cycles (e.g. see Santos et al. 2010; Dumusque et al. 2011; Díaz et al. 2016).

Short-term activity-induced RV variations are quasi-periodic: they are modulated by the star's rotation and change as active

* E-mail: rhaywood@cfa.harvard.edu

¹Plages are formed in the upper photosphere, chromosphere and upper layers of the stellar atmosphere (Lean 1997; Murdin 2002). They are not part of the lower photosphere, where the continuum absorption lines originate, from which the RV of a star is measured; however, plage regions do map closely to faculae and sunspots (Schrijver 2002; Hall 2008).
regions (starspots and faculae) emerge, evolve and disappear. The amplitude of these variations is $1-2 \text{ m s}^{-1}$ in 'quiet' stars (Isaacson & Fischer 2010), but they are often larger than this and can mimic or conceal the Doppler signatures of orbiting planets. This has resulted in several false detections (see Queloz et al. 2001; Bonfils et al. 2007; Huélamo et al. 2008; Boisse et al. 2009, 2011; Gregory 2011; Haywood et al. 2014; Robertson et al. 2014; Santos et al. 2014 and many others).

Understanding the RV signatures of stellar activity, especially those at the stellar rotation time-scale, is essential to develop the next generation of more sophisticated activity models and further improve our ability to detect and characterize (super-)Earths and even small Neptunes in orbits of a few days to weeks. In particular, identifying informative and reliable proxies for the activity-driven RV variations is crucial. Desort et al. (2007) found that the traditional spectroscopic indicators [the bisector span (BIS) and full width at half-maximum (FWHM) of the cross-correlation profile] and photometric variations do not give enough information for slowly rotating, Sun-like stars (low $v \sin i$) to disentangle stellar activity signatures from the orbits of super-Earth-mass planets.

The Sun is a unique test case as it is the only star whose surface can be resolved at high resolution, therefore allowing us to investigate directly the impact of magnetic features on RV observations. Early attempts to measure the RV of the integrated solar disc did not provide quantitative results about the individual activity features responsible for RV variability. Jiménez et al. (1986) measured integrated sunlight using a resonant scattering spectrometer and found that the presence of magnetically active regions on the solar disc led to variations of up to $15 \text{ m s}^{-1}$. They also measured the disc-integrated magnetic flux but did not find any significant correlation with RV at the time due to insufficient precision. At about the same time, Deming et al. (1987) obtained spectra of integrated sunlight with an uncertainty level below $5 \text{ m s}^{-1}$, enabling them to see the RV signature of supergranulation. The trend they observed over the 2-yr period of their observations was consistent with suppression of convective blueshift from active regions on the solar surface. A few years later, Deming & Plymate (1994) confirmed the findings of both Jiménez et al. (1986) and Deming et al. (1987), only with a greater statistical significance. Not all studies were in agreement with each other; however, McMillan et al. (1993) recorded spectra of sunlight scattered off the Moon over a 5-yr period and found that any variations due to solar activity were smaller than $4 \text{ m s}^{-1}$.

More recently, Molaro & Centurión (2010) obtained HARPS spectra of the large and bright asteroid Ceres to construct a wavelength atlas for the Sun. They found that these spectra of scattered sunlight provided precise disc-integrated solar RVs, and proposed using asteroid spectra to calibrate high-precision spectrographs used for planet hunting, such as HIRES and HARPS (see a recent paper by Lanza et al. 2016).

In parallel, significant discoveries were made towards a precise quantitative understanding of the RV impact of solar surface features. Lagrange, Desort & Meunier (2010) and Meunier, Desort & Lagrange (2010a) used a catalogue of sunspot numbers and sizes and Michelson Doppler Imager/Solar and Heliospheric Observatory (MDI/SOHO) magnetograms to simulate integrated-Sun spectra over a full solar cycle and deduce the impact of sunspots and faculae/plage on RV variations. The work of Meunier et al. (2010a) also relied on the amplitude of convection inhibition derived by Brandt & Solanki (1990), based on spatially resolved magnetogram observations of plage and quiet regions on the Sun (i.e. independently of full-disc RV measurements).

Sunspot umbrae and penumbrae are cooler and therefore darker than the surrounding quiet photosphere, producing a flux deficit

at the local rotational velocity. Faculae, which tend to be spatially associated with spot groups, are slightly brighter than the quiet photosphere, producing a spectral flux excess at the local rotational velocity. Spots have large contrasts and small area coverage, while faculae have lower contrasts but cover large areas. We thus expect their rotational Doppler signals to be of opposite signs (due to the opposite sign of their flux) and of roughly similar amplitudes. However, facular emission also occurs independently of spot activity in the more extended solar magnetic network (Chapman et al. 2001), so their contributions do not cancel out completely. Meunier et al. (2010a) found that the residual signal resulting from the rotational Doppler imbalance of both sunspots and faculae is comparable to that of the sunspots on their own. Lagrange et al. (2010) estimated the rotational perturbation due to sunspot flux deficit to be of the order of $1 \text{ m s}^{-1}$. In a complementary study, Meunier et al. (2010a) also investigated the effect of sunspots and faculae on the suppression of convective blueshift by magnetic regions. Sunspots occupy a much smaller area than faculae, and as they are dark, they contribute little flux, so their impact on the convective blueshift is negligible. Facular suppression of granular blueshift, however, can lead to variations in RV of up to $8-10 \text{ m s}^{-1}$ (Meunier et al. 2010a). Meunier, Lagrange & Desort (2010b) estimated the disc-integrated solar RV variations expected from the suppression of convective blueshift, directly from MDI/SOHO Dopplergrams and magnetograms (in the 6768 Å Ni I line). Their reconstructed RV variations, over one magnetic cycle, agree with the simulations of Meunier et al. (2010a), establishing the suppression of convective blueshift by magnetic features as the dominant source of activity-induced RV variations. Meunier et al. (2010b) also found that the regions where the convective blueshift is most strongly attenuated correspond to the most magnetically active regions, as was expected.

Following the launch of the *Solar Dynamics Observatory* (SDO; Pesnell, Thompson & Chamberlin 2012) in 2010, continuous observations of the solar surface brightness, velocity and magnetic fields have become available with image resolution finer than the photospheric granulation pattern. This allows us to probe the RV variations of the Sun in unprecedented detail. In the present paper, we deduce the activity-driven RV variations of the Sun based on HARPS observations of the bright asteroid Vesta (Section 2). In parallel, we use high spatial resolution continuum, Dopplergram and magnetogram images of the Fe I 6173 Å line from the Helioseismic and Magnetic Imager (HMI/SDO; Schou et al. 2012) to model the RV contributions from sunspots (and pores), faculae and granulation via inhibition of granular blueshift and flux blocking (Section 3). This allows us to create a model which we test against the HARPS observations (Section 4). Finally, we discuss the implications of our study for the effectiveness of various proxy indicators for activity-driven RV variations in stars. We show that the disc-averaged magnetic flux could become an excellent proxy for activity-driven RV variations on other stars (Section 5).

## 2 HARPS OBSERVATIONS OF SUNLIGHT SCATTERED OFF VESTA

### 2.1 HARPS spectra

The HARPS spectrograph, mounted on the European Southern Observatory 3.6 m telescope at La Silla, was used to observe sunlight scattered from the bright asteroid 4/Vesta (its average magnitude during the run was 7.6). Two to three measurements per night were made with simultaneous Thorium exposures for a total of 98 observations, spread over 37 nights between 2011 September 29 and December 7. The geometric configuration of the Sun and Vesta
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Schematic representation of the Sun, Vesta and Earth configuration during the period of observations (not to scale).</figcaption></figure>

relative to the observer is illustrated in Fig. 1. At the time of the observations, the Sun was just over three years into its 11-yr magnetic cycle; the SDO data confirm that the Sun showed high levels of activity.

The spectra were reprocessed using the HARPS Data Reduction Software (DRS) pipeline (Baranne et al. 1996; Lovis & Pepe 2007). Instead of applying a conventional barycentric correction, the wavelength scale of the calibrated spectra was adjusted to correct for the Doppler shifts due to the relative motion of the Sun and Vesta, and the relative motion of Vesta and the observer (see Section 2.3).

The FWHM and BIS (as defined in Queloz et al. 2001) and log R'_{HK} index were also derived by the pipeline. The median, minimum and maximum signal-to-noise ratio of the reprocessed HARPS spectra at central wavelength 556.50 nm are 161.3, 56.3 and 257.0, respectively. Overall, HARPS achieved a precision of 75 ± 25 cm s⁻¹ (see Table A1). The reprocessed HARPS cross-correlation functions and derived RV measurements are available in the Supplementary Files that are available online.

We account for the RV modulation induced by Vesta's rotation in Section 2.4.1, and investigate sources of intra-night RV variations in Section 2.4.2. We selected the SDO images in such a way as to compensate for the different viewing points of Vesta and the SDO spacecraft: Vesta was trailing SDO, as shown in Fig. 1. This is taken into account in Section 2.5.

## 2.2 Solar rest frame

The data reduction pipeline for HARPS assumes that the target observed is a distant point-like star, and returns its RV relative to the Solar system barycentre. In order to place the observed RVs of Vesta in the solar rest frame, we perform the following operation:

$$RV = RV_{\text{bary,Earth}} + v_{\text{sv}} + v_{\text{ve}}, \qquad (1)$$

where RV_{bary,Earth} is the barycentric RV of the Earth, i.e. the component of the observer's velocity relative to the Solar system barycentre, towards the apparent position of Vesta. It can be found in the fits header for each observation. The two components v_{sv} and v_{ve}, retrieved from the JPL Horizons data base,<fnref n="2" /> correspond to

(1) v_{sv}: the velocity of Vesta relative to the Sun at the instant that light received at Vesta was emitted by the Sun;

(2) v_{ve}: the velocity of Vesta relative to Earth at the instant that light received by HARPS was emitted at Vesta.

This correction accounts for the RV contribution of all bodies in the Solar system and places the Sun in its rest frame. The values of RV_{bary,Earth}, v_{sv} and v_{ve} are given in the Supplementary Files that are available online.

## 2.3 Relativistic Doppler effects

The only relativistic corrections made by JPL Horizons are for gravitational bending of the light and relativistic aberration due to the motion of the observer (Giorgini, private communication). We therefore must correct for the relativistic Doppler shifts incurred. The wavelength correction factor to be applied is given by Lindegren & Dravins (2003) as

$$\frac{\lambda_e}{\lambda_o} = \frac{\sqrt{1 - \frac{v^2}{c^2}}}{1 + \frac{v \cos \theta_o}{c}}, \qquad (2)$$

where λ_e is the wavelength of the light at emission, λ_o is the wavelength that is seen when it reaches the observer, and θ_o is the angle between the direction of the emitter at emission and the observed direction of the light at reception. v is the total magnitude of the velocity vector of the observer relative to the emitter and c is the speed of light.

We apply this correction twice:

(1) The light is emitted by the Sun and received at Vesta. In this case, v is the magnitude of the velocity of Vesta with respect to the Sun (decreasing from approximately 19.9 to 19.4 km s⁻¹ over the duration of the HARPS run), and the radial component v cos θ_o is equal to v_{sv} (defined in Section 2.2, starting at about 1.66 km s⁻¹ and reaching 1.72 km s⁻¹ at opposition near the middle of the run).

(2) Scattered sunlight is emitted from Vesta and received at La Silla. v is the magnitude of the velocity of Vesta with respect to an observer at La Silla (increasing from 20 to 36 km s⁻¹ over the run), and v cos θ_o is v_{ve} (starting at 19.8 km s⁻¹ and reaching 23.3 km s⁻¹ at opposition).

The two wavelength correction factors are then multiplied together in order to compute the total relativistic correction factor to be applied to the pixel wavelengths in the HARPS spectra, from which we derive the correct RVs, shown in Fig. 2(a) (see also column 2 of Table A1). All velocities are measured at the flux-weighted mid-exposure times of observation (MJD_{mid\_UTC}).

The reader may also refer to appendix A of Lanza et al. (2016) for further details on these Doppler corrections. All quantities used to calculate these effects are listed in the Supplementary Files that are available online.

## 2.4 Sources of intra-night RV variations

### 2.4.1 Vesta's axial rotation

Vesta rotates every 5.34 h (Stephenson 1951), so any significant inhomogeneities in its shape or surface albedo will induce an RV modulation. Vesta's shape is close to a spheroid (Thomas et al. 1997), and Lanza & Molaro (2015) found that the RV modulation expected from shape inhomogeneities should not exceed 0.060 m s⁻¹.

<footnote n="2">² Solar System Dynamics Group, Horizons On-Line Ephemeris System, 4800 Oak Grove Drive, Jet Propulsion Laboratory, Pasadena, CA 91109 USA - Information: http://ssd.jpl.nasa.gov/, Jon.Giorgini@jpl.nasa.gov</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) HARPS RV variations in the solar rest frame, corrected for relativistic Doppler effects (but not yet corrected for Vesta's axial rotation). (b) HARPS RV variations of the Sun as a star (after removing the RV contribution of Vesta's axial rotation). (c) Nightly binned HARPS RV variations of the Sun as a star – note the change in scale on the y-axis. The values for each time series are provided in the Supplementary Files that are available online.</figcaption></figure>

Stephenson (1951) presented a photometric study of the asteroid, and suggested that its surface brightness is uneven. He reported brightness variations $\delta m = 0.12$ mag. To make a rough estimate of the amplitude of the RV modulation, we can assume that the brightness variations are due to a single dark equatorial spot on the surface of Vesta, blocking a fraction $\delta f$ of the flux $f$. $\delta m$ and $\delta f$ are related as follows:

$$ \delta m = \frac{2.5 d(\ln f)}{\log(e)} \sim 1.08 \frac{\delta f}{f}, \qquad (3) $$

The fractional flux deficit caused by a dark spot can thus be approximated as

$$ \frac{\delta f}{f} \sim \delta m / 1.08 \sim 0.11. \qquad (4) $$

When the dayside of Vesta is viewed fully illuminated, this spot will give an RV modulation equal to

$$ \Delta RV_{\text{vesta}} = -\frac{\delta f}{f} v_{\text{eq}} \cos \theta \sin \theta, \qquad (5) $$

where $\theta$ is the angle between the spot on the asteroid and our line of sight, and increases from $-\pi/2$ to $+\pi/2$ as it traverses the visible daylight hemisphere. Due to foreshortening, the RV contribution is decreased by a factor $\cos \theta$. The line-of-sight velocity varies with $\sin \theta$. The asteroid's equatorial velocity $v_{\text{eq}}$ is given by

$$ v_{\text{eq}} = 2\pi \frac{R_{\text{vesta}}}{P_{\text{rot}}}. \qquad (6) $$

Using a mean radius $R_{\text{vesta}} = 262.7$ km (Russell et al. 2012) and the rotational period $P_{\text{rot}} = 5.34$ h, we obtain $v_{\text{eq}} = 85.8$ m s$^{-1}$. The maximum RV amplitude of Vesta's rotational modulation, expected at $\theta = \pi/4$, is thus approximately 4.7 m s$^{-1}$. The RV modulation due to surface brightness inhomogeneity should therefore dominate strongly over shape effects.

We find that this RV contribution is well modelled as a sum of Fourier components modulated by Vesta's rotation period:

$$ \Delta RV_{\text{vesta}}(t) = C \cos(2\pi - \lambda(t)) + S \sin(2\pi - \lambda(t)), \qquad (7) $$

where $\lambda(t)$ is the apparent planetographic longitude of Vesta at the flux-weighted mid-times of the HARPS observations and can be retrieved via the JPL Horizons data base (the values of $\lambda$ are listed in Table A1). C and S are scaling parameters, which we determine via a global optimal scaling procedure, performed with the solar activity contributions (see Section 4.1). Since the phase-folded light curve of Vesta shows a double-humped structure (Stephenson 1951), we also tested adding further Fourier terms modulated by the first harmonic of the asteroid's rotation. The improvement to the fit was negligible, so we preferred the simpler model of equation (7).

Fig. 2(b) shows the RV observations obtained after subtracting Vesta's rotational signature, with coefficients C and S of equation (7) derived from the global fit of Section 4.1. The night-to-night scatter has been reduced, even though much of it remains in the first block of observations; this is discussed in the following section.

### 2.4.2 Additional intra-night scatter in first half of HARPS run

The RV variations in the first part of the HARPS run (nights 0 to 11 in Fig. 2) contain some significant scatter, even after accounting for Vesta's rotation. This intra-night scatter does not show in the solar FWHM, BIS or $\log(R'_{HK})$ variations. We investigated the cause of this phenomenon and excluded changes in colour of the asteroid or instrumental effects as a potential source of additional noise. Vesta was very bright (7.6 mag), so we deem the phase and proximity of the Moon unlikely to be responsible for the additional scatter observed.

Solar p-mode oscillations dominate the Sun's power spectrum at a time-scale of about 5 min. Most of the RV oscillations induced by p-mode acoustic waves are therefore averaged out within the 15-min HARPS exposures. Granulation motions result in RV signals of several m s$^{-1}$, over time-scales ranging from about 15 min to several hours. Taking multiple exposures each night and averaging them together (as plotted in panel (c) of Fig. 2) can help to significantly reduce granulation-induced RV variations.

The velocity measurements are sensitive to any displacement of the image of Vesta from the centre of the 1-arcsec fibre. Vesta had a finite angular diameter of 0.49 arcsec at the start of the run, and 0.32 arcsec at the end. Light reflected from Vesta is blueshifted by 172 m s$^{-1}$ at the approaching limb and redshifted by 172 m s$^{-1}$ at the receding limb. We simulated the effects of seeing on an ensemble of photons originating from different points on Vesta's disc, applying random angular deviations with a Gaussian seeing distribution. For a given displacement of Vesta from the centre of the fibre, we computed the mean rotational velocity displacement of those photons falling within the fibre. We found an approximately linear dependence of the mean rotational velocity on mean offset from the centre of the fibre. Our numerical simulations suggest that a mean displacement by 0.1 arcsec of the image from the centre of the fibre in the direction orthogonal to Vesta's rotation axis, averaged over the exposure, gives a velocity error that is closely approximated by the empirical expression

$$ \Delta v = 4.0 \left( \frac{\theta}{0.49 \text{ arcsec}} \right) / \sqrt{\frac{\sigma_{\text{seeing}}}{1.0 \text{ arcsec}}} \text{ m s}^{-1}, \qquad (8) $$

where $\sigma_{\text{seeing}}$ represents the FWHM of the Gaussian seeing distribution. This is of the same order as the excess RV scatter observed during the first part of the run, when successive observations within
a night were widely separated at different airmasses. We conclude that small airmass-dependent guiding errors provide a plausible explanation. During the latter part of the run, the observations were contiguous and were thus conducted at similar airmass, leading to more consistent guiding and smaller scatter. Other factors such as asymmetric image quality, coma of the telescope varying with elevation (and sky rotation on an equatorial telescope), atmospheric dispersion, etc. may also contribute to this effect.

The remaining variations, of the order of 7–10 m s⁻¹, are modulated by the Sun’s rotation and are caused by the presence of magnetic surface markers, such as sunspots and faculae. These variations are the primary focus of this paper, and we model them using SDO/HMI data in Section 3.

## 2.5 Time lag between Vesta and SDO observations

At the time of the observations, the asteroid Vesta was trailing the SDO spacecraft, which orbits the Earth (see Fig. 1). In order to model the solar hemisphere facing Vesta at time t, we used SDO images recorded at t + Δt, where Δt is proportional to the difference in the Carrington longitudes of the Earth/SDO and Vesta at the time of the HARPS observation. These longitudes were retrieved from the JPL Horizons data base. The shortest delay, at the start of the observations, was ~2.8 d, while at the end of the observations it reached just over 6.5 d (see Table A1). We cannot account for the evolution of the Sun’s surface features during this time, and must assume that they remain frozen in this interval. The emergence of sunspots can take place over a few days, but in general large magnetic features (sunspots and networks of faculae) evolve over time-scales of weeks rather than days. Visual inspection of animated sequences of SDO images obtained during the campaign revealed no major flux emergence events on the visible solar hemisphere.

## 3 PIXEL STATISTICS FROM SDO/HMI IMAGES

In the second part of this analysis, we aim to determine the RV contribution from granulation, sunspots and facular regions. We use high-resolution full-disc continuum intensity, line-of-sight Doppler images and line-of-sight magnetograms from the HMI instrument on board SDO.<fnref n="3" /> These were retrieved for the period spanning the HARPS observations of Vesta at times determined by the time lags detailed in Section 2.5 (the exact date and time stamps of the images are listed in the Supplementary Files that are available online). SDO/HMI images the solar disc in the Fe I 6173 Å line at a distance of 45 s, with a spatial resolution of 1 arcsec using a CCD of 4096 × 4096 square pixels. We first converted the SDO/HMI images from pixel coordinates to heliographic coordinates, i.e. to a coordinate system centred on the Sun. This coordinate system is fixed with respect to the Sun’s surface and rotates in the side- real frame once every 25.38 d, which corresponds to a Carrington rotation period (Carrington 1859). A surface element on the Sun, whose image falls on pixel *ij* of the instrument detector, is at position ($w_{ij}$, $n_{ij}$, $r_{ij}$) relative to the centre of the Sun, where *w* is westwards, *n* is northwards and *r* is in the radial direction (see Thompson 2006 for more details on the coordinate system used). The spacecraft is at position (0, 0, $r_{sc}$). The *w*, *n*, *r* components of the spacecraft’s position relative to each element *ij* can thus be

<table><caption>**Table 1.** Solar differential rotation profile parameters from Snodgrass & Ulrich (1990).</caption><thead><tr><td>Parameter</td><td>Value (deg d<sup>-1</sup>)</td></tr></thead><tbody><tr><td>α<sub>1</sub></td><td>14.713</td></tr><tr><td>α<sub>2</sub></td><td>-2.396</td></tr><tr><td>α<sub>3</sub></td><td>-1.787</td></tr></tbody></table>

written as

$$ 
\begin{aligned}
\delta w_{ij} &= w_{ij} - 0 \\
\delta n_{ij} &= n_{ij} - 0 \\
\delta r_{ij} &= r_{ij} - r_{sc}.
\end{aligned} 
\qquad (9) $$

The spacecraft’s motion and the rotation of the Sun introduce velocity perturbations, which we determine in Sections 3.1 and 3.2, respectively. These two contributions are then subtracted from each Doppler image, thus revealing the Sun’s magnetic activity velocity signatures. We compute the RV variations due to the suppression of convective blueshift and the flux blocked by sunspots on the rotating Sun in Sections 3.7.3 and 3.7.2. We show that the Sun’s activity-driven RV variations are well reproduced by a scaled sum of these two contributions in Section 4. Finally, we compute the disc-averaged magnetic flux and compare it as an RV proxy against the traditional spectroscopic activity indicators in Section 5.

### 3.1 Spacecraft motion

The *w*, *n*, *r* components of the velocity incurred by the motion of the spacecraft relative to the Sun, $v_{sc}$, are given in the fits header of each SDO/HMI observation to a precision of 10⁻⁶ m s⁻¹. The magnitude of the spacecraft’s velocity away from pixel *ij* can therefore be expressed as

$$ v_{sc,ij} = - \frac{\delta w_{ij} v_{sc,w_{ij}} + \delta n_{ij} v_{sc,n_{ij}} + \delta r_{ij} v_{sc,r_{ij}}}{d_{ij}}, \qquad (10) $$

where

$$ d_{ij} = \sqrt{\delta w_{ij}^2 + \delta n_{ij}^2 + \delta r_{ij}^2} \qquad (11) $$

is the distance between pixel *ij* and the spacecraft. We note that all relative velocities in this paper follow the natural sign convention that velocity is rate of change of distance.

### 3.2 Solar rotation

The solar rotation as a function of latitude was measured by Snodgrass & Ulrich (1990) in low-resolution full-disc Dopplergrams and magnetograms obtained at the Mount Wilson 150-foot tower telescope between 1967 and 1987. By cross-correlating time series of Dopplergrams and magnetograms, they were able to determine the rate of motion of surface features (primarily supergranulation cells and sunspots) and deduce the rate of rotation of the Sun’s surface as a function of latitude. The solar differential rotation profile $\omega(\phi)$ at each latitude $\phi$ is commonly described by a least-squares polynomial of the form

$$ \omega(\phi) = \alpha_1 + \alpha_2 \sin^2 \phi + \alpha_3 \sin^4 \phi. \qquad (12) $$

The best-fitting parameters found by Snodgrass & Ulrich (1990), used in this analysis, are given in Table 1. We apply this rotation profile in the heliographic frame to determine the *w*, *n*, *r* components induced by the solar rotation velocity along the line of sight to a

<footnote n="3">³ HMI data products can be downloaded online via the Joint Science Operations Center website: http://jsoc.stanford.edu.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** First three panels: SDO/HMI flattened intensity $I_{\text{flat}}$, line-of sight velocity $v$ (km s$^{-1}$) for the non-rotating Sun and unsigned radial magnetic flux $|B_r|$ (G) of the Sun, observed on 2011 November 10 at 00:01:30 UTC. Last panel: our thresholded image, highlighting faculae (red/lighter shade pixels) and sunspots (blue/darker shade pixels). For this set of observations (representative of the whole run), faculae account for 9 per cent of the total pixel count, while sunspots account for less than 0.4 per cent. The remaining ~90 per cent of the pixels on the solar disc are magnetically quiet.</figcaption></figure>

given image pixel, $v_{\text{rot},w}$, $v_{\text{rot},n}$ and $v_{\text{rot},r}$. Normalizing again by $d$, we can write

$$v_{\text{rot}} = - \frac{\delta w v_{\text{rot},w} + \delta n v_{\text{rot},n} + \delta r v_{\text{rot},r}}{d} \quad (13)$$

### 3.3 Flattened continuum intensity

We flatten the continuum intensity images using a fifth-order polynomial function $L_{ij}$ with the limb-darkening constants given in Astrophysical Quantities (Allen 1973), through the IDL subroutine DARKLIMB_CORRECT.PRO.<fnref n="4" /> The flattened and non-flattened continuum intensities are related via the limb-darkening function $L$ as follows:

$$I_{\text{flat},ij} = \frac{I_{ij}}{L_{ij}} \quad (14)$$

We also define a scaling factor $\hat{K}$ which will be applied to $L_{ij}$ later on:

$$\hat{K} = \frac{\sum_{ij} I_{ij} L_{ij} W_{ij}}{\sum_{ij} L_{ij}^2 W_{ij}} \quad (15)$$

where the weighting factor $W_{ij}$ is set to unity for quiet-Sun pixels, and zero elsewhere.

### 3.4 Unsigned magnetic field strength

The SDO/HMI instrument measures the line-of-sight magnetic field strength $B_{\text{obs}}$. The magnetic field of the Sun stands radially out of the photosphere with a strength $B_r$. Due to foreshortening, the observed field $B_{\text{obs}}$ is less than the true (radial) field by a factor

$$\mu_{ij} = \cos \theta_{ij} \quad (16)$$

where $\theta_{ij}$ is the angle between the outward normal to the feature on the solar surface and the direction of the line of sight of the SDO spacecraft.

We can thus recover the full magnetic field strength by dividing by $\mu_{ij}$:

$$B_{r,ij} = B_{\text{obs},ij} / \mu_{ij} \quad (17)$$

The noise level in HMI magnetograms is a function of $\mu$ (Yeo, Solanki & Krivova 2013). It is lowest for pixels in the centre of the CCD, where it is close to 5 G, and increases towards the edges and reaches 8 G at the solar limb. For our analysis, we assume that the noise level is constant throughout the image with a conservative

<footnote n="4">⁴ Source code available at http://hesperia.gsfc.nasa.gov/ssw/gen/idl/solar/</footnote>
<figure><img src="image_6.png" /><figcaption>Figure 4. Observed pixel line-of-sight (unsigned) magnetic field strength, $|B_{\text{obs,ij}}|$ (G), as a function of flattened intensity $I_{\text{flat,ij}}$, for the Sun on 2011 November 10 at 00:01:30 UTC. The top and right histograms show the distributions of $|B_{\text{obs,ij}}|$ and $I_{\text{flat,ij}}$, respectively. The dashed lines represent the cutoff criteria selected to define the quiet photosphere, faculae and sunspots.</figcaption></figure>

value $\sigma_{B_{\text{obs,ij}}} = 8$ G, in agreement with the results of Yeo et al. (2013). We therefore set $B_{\text{obs,ij}}$ and $B_{\text{r,ij}}$ to 0 for all pixels with a line-of-sight magnetic field measurement ($B_{\text{obs,ij}}$) below this value.

## 3.5 Identifying quiet-Sun regions, faculae and sunspots

As is routinely done in solar work, we do not consider pixels that are very close to the limb ($\mu_{ij} < 0.1$), as the limb-darkening model becomes unreliable on the very edge of the Sun. This affects about 1 per cent of the pixels of the solar disc.

The first three panels of Fig. 3 show an SDO/HMI flattened intensitygram, line-of-sight Dopplergram and unsigned radial magnetogram for a set of images taken on 2011 November 10, after removing the contributions from spacecraft motion and solar rotation. We identify quiet-Sun regions, faculae and sunspots by applying magnetic and intensity thresholds.

**Magnetic threshold:** the distribution of pixel unsigned observed magnetic field strength as a function of pixel flattened intensity is shown in Fig. 4. In the top histogram and main panel, we see that the distribution of magnetic field strength falls off sharply with increasing field strength. The vast majority of pixels are clustered close to the noise level: these pixels are part of the quiet-Sun surface. Note that the fragmented distribution of pixels with magnetic field strength less than a few gauss arises from the numerical precision in the SDO/HMI images. This is not an issue in our analysis, however, as these pixels are well below the magnetic noise threshold and we set their field value to zero. We separate active regions from quiet-Sun regions by applying a threshold in unsigned radial magnetic field strength for each pixel. Yeo et al. (2013) investigated the intensity contrast between the active and quiet photosphere using SDO/HMI data, and found an appropriate cutoff at

$$|B_{\text{r,ij}}| > 3 \sigma_{B_{\text{obs,ij}}} / \mu_{ij}, \qquad (18)$$

<figure><img src="image_7.png" /><figcaption>Figure 5. Pixel line-of-sight velocity, $v_{ij}$ (in m s$^{-1}$), as a function of flattened intensity $I_{\text{flat,ij}}$, for the Sun on 2011 November 10 at 00:01:30 UTC. The top and right histograms show the distributions of $v_{ij}$ and $I_{\text{flat,ij}}$, respectively, in bins of 1000. In the case of active pixels (yellow dots), the line-of-sight velocity is invariant with pixel brightness. For quiet-Sun pixels (black dots), however, brighter pixels are blueshifted while fainter pixels are redshifted: this effect arises from granular motions.</figcaption></figure>

where $\sigma_{B_{\text{obs,ij}}}$ represents the magnetic noise level in each pixel (see the last paragraph of Section 3.4). As in Yeo et al. (2013), we exclude isolated pixels that are above this threshold as they are likely to be false positives. We can thus write

$$|B_{\text{r,thresh,ij}}| = 24 \text{ G} / \mu_{ij}. \qquad (19)$$

The combined effects of imperfect limb-darkening correction and the $\mu$-correction for the radial magnetic field strength at the very edge of the solar disc ($\mu_{ij} < 0.3$) result in a rim of dark pixels being identified as sunspot pixels. In order to avoid this, we set all such pixels as quiet-Sun elements. In any case, sunspots become invisible near the edge of the solar disc because of the Wilson depression (Loughhead & Bray 1958), so our cut should not affect the identification of real sunspot pixels.

**Intensity threshold:** the distribution of line-of-sight velocity as a function of pixel flattened intensity is shown in Figs 4 and 5. The main panel allows us to further categorize active-region pixels into faculae and sunspots (umbra and penumbra). We apply the intensity threshold between faculae and spots of Yeo et al. (2013):

$$I_{\text{thresh}} = 0.89 \hat{I}_{\text{quiet}}, \qquad (20)$$

where $\hat{I}_{\text{quiet}}$ is the mean pixel flattened intensity over quiet-Sun regions. It can be calculated by summing the flattened intensity of each pixel:

$$\hat{I}_{\text{quiet}} = \frac{\sum_{ij} I_{\text{flat,ij}} W_{ij}}{\sum_{ij} W_{ij}}, \qquad (21)$$

where the weighting factor $W_{ij}$ is set to 1 if $|B_{\text{r,ij}}| < |B_{\text{r,thresh,ij}}|$, 0 otherwise.

In the main panel of Fig. 5, quiet-Sun pixels are plotted in black, while active-region pixels are overplotted in yellow. The last panel of Fig. 3, which shows the thresholded image according to these
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Zoom-in on the largest sunspot group on the approaching hemisphere of the Sun observed on 2011 November 10 at 00:01:30 UTC. Red pixels represent faculae, while the sunspot pixels are colour-coded according to the direction of their velocity flows, thereby revealing the presence of Evershed flows: dark blue pixels are directed towards the observer ($v < 0$), while light cyan pixels are directed away from the observer ($v > 0$).</figcaption></figure>

$I_{\text{flat,ij}}$ and $|B_{r,ij}|$ criteria, confirms that these thresholding criteria are effective at identifying sunspot and faculae pixels correctly.

## 3.6 Surface velocity flows in sunspot penumbrae and umbrae

The range of velocities in the penumbral pixels (which form the horizontal oval of yellow points just below $I_{\text{thresh}}$ in Fig. 5) is large owing to the Evershed effect (Evershed 1909). The darker umbral pixels have a narrower velocity distribution, allowing the umbral flows that feed the Evershed effect to be resolved into a blueshifted and redshifted velocity component. The much broader range of velocities seen in Fig. 5 for penumbral versus umbral pixels confirms that these velocity flows accelerate with distance from the middle of the spot (Evershed 1909). This effect is highlighted in Fig. 6, which shows a zoom-in on the largest sunspot group in the images of Fig. 3. The Evershed flows are tangential to the surface. They will be most visible for sunspots located away from disc centre, where a larger proportion of them will be directed along our line of sight.

The sunspot group illustrated in Fig. 6 is located near the approaching side of the Sun. The flows in its left half are directed away from the observer, while in its right half they are directed towards the observer.

## 3.7 Decomposing the Sun's RV into individual feature contributions

The total disc-integrated RV of the Sun is the sum of all contributions from the quiet-Sun, sunspot and faculae/plage regions. The 'quiet' parts of the Sun's surface are in constant motion due to granulation, while sunspots and faculae induce RV variations via two processes (cf. introduction).

Photometric effect: as the Sun rotates, the presence of dark spots or bright faculae on the solar surface breaks the Doppler balance between the approaching (blueshifted) and receding (redshifted) hemispheres. We calculated the contribution of dark spots and bright faculae separately and our findings confirm those of Meunier et al. (2010a), who found that this residual signal is approximately equal to the photometric contribution from sunspots. We note, however,

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Correlation diagram showing the relationship between the rotational Doppler imbalances resulting from sunspots (horizontal axis) and faculae (vertical axis), if computed separately. The Spearman correlation coefficient is -0.69. The variations are only partially anticorrelated, reflecting both the tendency of faculae to be spatially associated with sunspot groups and the existence of faculae in bright networks not associated with sunspots. The values of these two basis functions are given in the Supplementary Files that are available online.</figcaption></figure>

that this is merely a coincidence arising from the specific geometrical configuration and ratio of sunspots to faculae/plage on the Sun. In other words, this assumption may not be valid for other stars with different spot to faculae configurations and/or filling ratios. We therefore decide to include the effect of faculae in this study (even though it is effectively negligible in the case of the Sun). Because the two contributions are correlated (see Fig. 7), we must account for them in a single term $\Delta\hat{v}_{\text{phot}}$, which we describe in Section 3.7.2.

Convective effect: sunspots and faculae are strongly magnetized features that inhibit convective motions. Sunspots, which cover a small area of the solar surface and contribute little flux, have a very small contribution. In the case of faculae, however, this contribution is large and is expected to be the dominant contribution to the total solar RV variations (Meunier et al. 2010a,b). We compute this contribution $\Delta\hat{v}_{\text{conv}}$ in Section 3.7.3.

We can thus write the total disc-integrated RV of the Sun as

$$ \hat{v} = \hat{v}_{\text{quiet}} + \Delta\hat{v}_{\text{phot}} + \Delta\hat{v}_{\text{conv}}, \quad (22) $$

where $\hat{v}_{\text{quiet}}$ represents the velocity contribution of convective motions in quiet-Sun pixels, which we calculate in Section 3.7.1.

### 3.7.1 Velocity contribution of convective motions in quiet-Sun regions

We estimate the average RV of the quiet Sun by summing the intensity-weighted velocity of non-magnetized pixels, after removing the spacecraft motion and the Sun's rotation:

$$ \hat{v}_{\text{quiet}} = \frac{\sum_{ij} (v_{ij} - \delta v_{\text{sc,ij}} - \delta v_{\text{rot,ij}}) I_{ij} W_{ij}}{\sum_{ij} I_{ij} W_{ij}}. \quad (23) $$

For this calculation, the weights are defined as in

$$ \begin{aligned} W_{ij} &= 1 & \text{if } |B_{r,ij}| < |B_{r,\text{thresh},ij}|, \\ W_{ij} &= 0 & \text{if } |B_{r,ij}| > |B_{r,\text{thresh},ij}|. \end{aligned} \quad (24) $$
This velocity field is thus averaged over the vertical motions of convection granules on the solar surface. Hot and bright granules rise up to the surface, while cooler and darker fluid sinks back towards the Sun's interior. This process is visible in the main panel of Fig. 5: quiet-Sun pixels (black dots) are clustered in a tilted ellipse. The area of the upflowing granules is larger than that enclosed in the intergranular lanes, and the granules are carrying hotter and thus brighter fluid. This results in a net blueshift, as seen in Fig. 5.

### 3.7.2 Rotational Doppler imbalance due to dark sunspots and bright faculae

This velocity perturbation can be obtained by summing the line-of-sight velocity of sunspot pixels corrected for the spacecraft's motion, and weighted by the flux deficit or excess produced by a dark spot or a bright plage, respectively:

$$ \Delta \hat{v}_{\text{phot}} = \frac{\sum_{ij} (v_{ij} - \delta v_{\text{sc},ij}) (I_{ij} - \hat{K} L_{ij}) W_{ij}}{\sum_{ij} I_{ij}}. \quad (25) $$

The weights $W_{ij}$ are set to 1 for pixels with $|B_{r,ij}| > |B_{r,thresh,ij}|$ (otherwise they are set to 0).

As illustrated in Fig. 7, the photometric effects of sunspots and faculae are quite similar in amplitude and roughly anticorrelated, due to their opposite flux signs. However, because they are not strictly anticorrelated, they sum into a net signal that has an amplitude similar to the photometric effect of sunspots (Meunier et al. 2010a). The value of $\Delta\hat{v}_{phot}$ at each time of the HARPS observations is listed in Table A1. From the SDO images, we derive an rms amplitude of 0.17 m s⁻¹ for this signal, which is smaller than the observational uncertainties of the HARPS velocities. This value is slightly lower than that found by Meunier et al. (2010a) during the peak of the Sun's activity cycle, of 0.42 m s⁻¹, but remains broadly consistent with their results given the small amplitude of this signal.

### 3.7.3 Suppression of convective blueshift from active regions

The suppression of granular blueshift induced by magnetically active regions ($|B_{r,ij}| > |B_{r,thresh,ij}|$, predominantly faculae) is

$$ \Delta \hat{v}_{\text{conv}} = \hat{v} - \hat{v}_{\text{quiet}} - \hat{v}_{\text{phot}}. \quad (26) $$

We measure the total disc-averaged velocity of the Sun $\hat{v}$ by summing the velocity contribution of each pixel $ij$, weighted by their intensity $I_{ij}$, after subtracting the spacecraft motion and solar rotation:

$$ \hat{v} = \frac{\sum_{ij} (v_{ij} - \delta v_{\text{sc},ij} - \delta v_{\text{rot},ij}) I_{ij}}{\sum_{ij} I_{ij}}. \quad (27) $$

The value of $\Delta\hat{v}_{conv}$ at each time of the HARPS observations is listed in Table A1. The rms amplitude of this basis function (unscaled) is 1.30 m s⁻¹, which is consistent with the value of 1.39 m s⁻¹ computed by Meunier et al. (2010a) during the peak of the solar activity cycle.

# 4 REPRODUCING THE RV VARIATIONS OF THE SUN

## 4.1 Total RV model

We combine our model of Vesta's rotational RV signal (presented in Section 2.4.1) with the two magnetic activity basis functions determined in Sections 3.7.3 and 3.7.2, in order to reproduce the

RV variations seen in the HARPS observations. The final model has the form

$$ \Delta RV_{\text{model}}(t) = A \Delta \hat{v}_{\text{phot}}(t) + B \Delta \hat{v}_{\text{conv}}(t) + \Delta RV_{\text{vesta}}(t) + RV_0. \quad (28) $$

We carry out an optimal scaling procedure in order to determine the scaling factors (A, B, C and S) of each of the contributions, as well as the constant offset $RV_0$ and a constant variance term $s^2$ added in quadrature to the observational errors. This variance term will account for any remaining uncorrelated noise arising from granulation motions. The SDO/HMI images were exposed for 45 s each, while the HARPS observations were exposed for 15 min and binned nightly. We expect these differences to result in a small night-to-night uncorrelated noise contribution. In addition, this variance term will absorb any residual signal of Vesta's axial rotation, which will naturally not be completely sinusoidal.

For each value of s, the basis functions are orthogonalized by subtracting their inverse-variance weighted means prior to performing the scaling. We determine the maximum likelihood via a procedure similar to the one described in Collier Cameron et al. (2006), whereby we maximize the likelihood $\mathcal{L}$ of the solution

$$ \ln \mathcal{L} = -\frac{n}{2} \ln(2\pi) - \frac{1}{2} \chi^2 - \frac{1}{2} \sum_{i=1}^{i=n} \ln(\sigma_i^2 + s^2) \quad (29) $$

with respect to s, where $\chi^2$ is the chi-squared value of the n HARPS data points, with uncertainties $\sigma_i$. This procedure is applied to the unbinned (not nightly averaged) HARPS data set, in order to determine the appropriate scaling coefficients (C and S) for Vesta's axial rotation (see equation 7). The total amplitude of the modulation induced by Vesta's rotation is equal to 2.39 m s⁻¹, which is of the same order as the amplitude we estimated in Section 2.4.1.

## 4.2 Agreement between HARPS observations and SDO-derived model

After all the scaling coefficients were determined, we grouped the observations in each night by computing the inverse-variance weighted average for each night. The final model is shown in Fig. 8.

We list the best-fitting values of the scaling parameters for each of the basis functions in Table 2. We note that the values of A and B, which represent the difference in response of HARPS and SDO to the Doppler imbalance and the suppression of convective blueshift, differ from unity. The HMI/SDO images are based on measurements of a single spectral line, namely the Fe I 6173 Å line, which may not be representative of the several thousand lines from which the HARPS RVs are derived. The value of $\Delta\hat{v}_{phot}$ is determined by the contrast between the magnetic elements (spots, faculae) and the quiet-Sun photosphere. Most of the lines in the HARPS mask lie bluewards of Fe I 6173 Å where the contrast is greater for both types of feature. We therefore expect to have to scale the SDO-derived $\Delta\hat{v}_{phot}$ by a factor of order, but greater than, unity. We measure it to be 2.45 ± 2.02; the amplitude of $\Delta\hat{v}_{phot}$ is so small that A cannot be distinguished from unity when the additional variance s is taken into account. The strength of the convection inhibition changes with line depth (cf. Gray 2009 and Meunier et al. 2010b). We find that the HARPS response to suppression of granular blueshift is 1.85 ± 0.27 times greater than we predict from SDO.

We show a plot of the observed HARPS RVs ($\Delta RV_{Sun}$) versus our SDO-modelled RVs ($\Delta RV_{model}$) in Fig. 9, in a similar fashion to Fig. 9 of Meunier et al. (2010b). We see good agreement between the data and the model. The Spearman correlation coefficient is 0.64 for
<figure><img src="image_10.png" /><figcaption>**Figure 8.** (a) HARPS RV variations of the Sun as a star, ΔRV<sub>Sun</sub>. (b) Scaled basis function for the suppression of convective blueshift, Δv̂<sub>conv</sub>, derived from SDO/HMI images. (c) Scaled basis function for the rotational Doppler imbalance due to spots and faculae, Δv̂<sub>phot</sub>. (d) Total RV model, ΔRV<sub>model</sub> (red/lighter shade, with errors including additional variance s), overlaid on top of the HARPS RV variations (blue/darker shade points). (e) Residuals obtained after subtracting the model from the observations. All RVs are in m s<sup>-1</sup>. Note that the scale of the y-axis is different from that used in Fig. 2. The values of these nightly binned time series are provided in the Supplementary Files available online.</figcaption></figure>

<table><caption>**Table 2.** Best-fitting parameters and rms amplitudes resulting from the optimal scaling procedure.</caption><thead><tr><th>Parameter</th><th>Value</th><th>Basis function</th><th>rms amplitude (unscaled)</th></tr></thead><tbody><tr><td>A</td><td>2.45 ± 2.02</td><td>Δv̂<sub>phot</sub></td><td>0.17 m s<sup>-1</sup></td></tr><tr><td>B</td><td>1.85 ± 0.27</td><td>Δv̂<sub>conv</sub></td><td>1.30 m s<sup>-1</sup></td></tr><tr><td>C</td><td>2.19 ± 0.42</td><td>cos(2π − λ)</td><td>0.74 m s<sup>-1</sup></td></tr><tr><td>S</td><td>0.55 ± 0.47</td><td>sin(2π − λ)</td><td>0.66 m s<sup>-1</sup></td></tr><tr><td>RV<sub>0</sub></td><td>99.80 ± 0.28 m s<sup>-1</sup></td><td></td><td></td></tr><tr><td>s</td><td>2.70 m s<sup>-1</sup></td><td></td><td></td></tr></tbody></table>

the full data set and 0.87 when considering only the second part of the run. This is close to, although not as good as, the correlation coefficient of 0.94 found by Meunier et al. (2010b) between simulated velocities of Meunier et al. (2010a) and MDI/SOHO Dopplergrams.

Panel (e) of Fig. 8 shows the residuals remaining after subtracting the total model ΔRV<sub>model</sub> from the HARPS observations of

<figure><img src="image_11.png" /><figcaption>**Figure 9.** HARPS RV variations of the Sun as a star versus our model derived from SDO/HMI images. Observations from the first part of the run are highlighted in a lighter shade.</figcaption></figure>

the Sun as a star ΔRV<sub>Sun</sub>. We see that the residuals are within the level of the error bars, when considering an extra variance term s of 2.7 m s<sup>-1</sup>. Maximum-likelihood analysis of the additional uncorrelated noise during the first part of the run (nights 0–11) yields an additional variance with rms amplitude 4.0 m s<sup>-1</sup>, while the second part (nights 36–68) has an additional noise rms of 1.5 m s<sup>-1</sup>. As discussed in Section 2.4.2, we attribute the excess scatter in the first block of nights to airmass-dependent guiding errors, which would have been more important in the first part of the run, where the observations within each night were taken at very different airmasses. The additional noise rms of 1.5 m s<sup>-1</sup> in the second part of the run is consistent with the rms due to solar granulation and supergranulation, ranging between 0.28 and 1.12 m s<sup>-1</sup>, as recently found by Meunier et al. (2015).

## 4.3 Relative importance of suppression of convective blueshift and rotational velocity imbalance

We see that the activity-induced RV variations of the Sun are well reproduced by a scaled sum of the two basis functions, v̂<sub>conv</sub> and v̂<sub>phot</sub> (shown in panels (b) and (c) of Fig. 8, respectively). As previously predicted by Meunier et al. (2010a), we find that the suppression of convective blueshift plays a dominant role (rms of 2.40 m s<sup>-1</sup>). This was also found to be the case for CoRoT-7, a main-sequence G9 star with a rotation period comparable to that of the Sun (Haywood et al. 2014). The relatively low amplitude of the modulation induced by sunspot flux blocking (rms of 0.17 m s<sup>-1</sup>) is expected in slowly rotating stars with a low v sin i (Desert et al. 2007).

## 4.4 Zero-point of HARPS

The wavelength adjustments that were applied to the HARPS RVs were based on precise prior dynamical knowledge of the rate of change of distance between the Earth and Vesta, and between Vesta and the Sun. The offset RV<sub>0</sub> = 99.80 ± 0.28 m s<sup>-1</sup> thus represents the zero-point of the HARPS instrument, including the mean granulation blueshift for the Sun. A previous study of integrated sunlight reflected by the Moon, by Molaro et al. (2013), determined a value of 102.2 ± 0.86 m s<sup>-1</sup>, which is significantly larger than our
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Correlation plots of the nightly averaged HARPS RV variations of the Sun as a star ΔRV<sub>Sun</sub> and suppression of convective blueshift Δv̂<sub>conv</sub> against (from left to right): the disc-averaged observed magnetic flux |B̂<sub>obs</sub>| (G), filling factor *f* (per cent), FWHM (km s<sup>-1</sup>), BIS (m s<sup>-1</sup>) and log(R'<sub>HK</sub>). Observations from the first part of the run are highlighted in a lighter shade. Spearman correlation coefficients are displayed in the bottom-right corner of each panel: for the full observing run (in bold and black), and for the second part of the run only (in blue).</figcaption></figure>

value. They did not account for the effect of sunspots and faculae,
however, so their result could be affected by activity-induced solar
variations at that time.

5 TOWARDS BETTER PROXIES FOR RV
OBSERVATIONS

**5.1 Spatial distributions of sunspots and faculae**

Aigrain, Pont & Zucker (2012) have shown that it is possible to pre-
dict the rotational Doppler imbalance due to photospheric surface
brightness inhomogeneities from a simultaneous high-precision op-
tical light curve. If one further assumes that faculae/plage regions
are co-spatial with spot groups, then they can also predict the form
of the RV variation caused by suppression of granular blueshift. A
recent analysis of the active host star CoRoT-7 by Haywood et al.
(2014) modelled activity-induced RV variations via the FF' method
of Aigrain et al. (2012). The predicted Doppler imbalance was much
smaller than the observed activity-driven RV variations. The asso-
ciated suppression of convective blueshift was of larger amplitude
than, and partially correlated with, the observed RVs. The resid-
uals, however, had a similar amplitude and shared the covariance
properties of the star's (simultaneous) light curve.

The present study provides a natural explanation of this mis-
match: on the Sun, the faculae are not perfectly co-spatial with
sunspot groups. Indeed, Fig. 7 shows us that the location of sunspot
groups gives an incomplete prediction of the facular coverage. Since
the suppression of granular blueshift is the dominant process at play
in slowly rotating stars such as CoRoT-7 and the Sun, it is there-
fore important to develop proxies that are directly sensitive to the
distribution of faculae on the stellar surface.

**5.2 Correlations between RV and traditional activity indicators**

Fig. 10 presents the correlations between ΔRV<sub>Sun</sub>, Δv̂<sub>conv</sub> and the
following activity indicators: the full-disc magnetic flux |B̂<sub>obs</sub>| and
filling factor *f* computed from the SDO/HMI images, and the ob-
served FWHM, BIS and log(R'<sub>HK</sub>) derived from the HARPS DRS

reduction pipeline. We computed the Spearman correlation coeffi-
cient to get a measure of the degree of monotone correlation between
each variable (the correlation between two variables is not neces-
sarily linear, for example between RV and BIS). The coefficients
are displayed in each panel of Fig. 10, both including and excluding
the observations made in the first part of the run, which show a lot
of intra-night scatter.

We do not show similar correlation plots for Δv̂<sub>phot</sub> because we
do not find any significant correlations with any of the activity
indicators; this is expected since Δv̂<sub>phot</sub> is such that it crosses zero
when the surface covered by spots and/or faculae is at a maximum,
i.e. when they are in the middle of the stellar disc (Δv̂<sub>phot</sub> is of
course still related to |B̂<sub>obs</sub>|, the FWHM, BIS and log(R'<sub>HK</sub>), but the
Spearman coefficient is close to zero).

We note that the relatively weak correlation between the observed
RVs and the BIS is not completely unexpected in the case of the Sun,
which has a vsin *i* of about 2 km s<sup>-1</sup>. The line profile distortions
induced by solar activity, including those measured by the BIS, will
therefore be smaller than the resolution of HARPS that is close to
2.5 km s<sup>-1</sup>. In other words, the resolution of HARPS is not adequate
to fully resolve the BIS variations in a star rotating as slowly as the
Sun, which could reduce the correlation coefficients computed with
the BIS.

We find a relatively weak correlation between the observed RVs
and the chromospheric flux index log(R'_{HK}), with a Spearman cor-
relation coefficient of 0.26 for the second half of the run (0.18 for
the full run).

Although naively one might expect log(R'_{HK}) to be a good pre-
dictor of plage filling factor, and hence of the convective RV com-
ponent, there are several physical factors that might reasonably be
expected to degrade the correlation over short time-scales (our data
set spans two to three solar rotations). Foreshortening and limb
darkening affect the Ca II H&K emission cores and the brightness-
weighted line-of-sight granular velocities in different ways. Near
the limb, the Ca II emission in plages originates in higher, hotter
regions of the chromosphere and remains bright. The limb darken-
ing of facular pixels is less than that of quiet-Sun pixels, but near
the limb the line-of-sight component of the radial motion of bright
granule cores is reduced by foreshortening. The disc-averaged line-
of-sight magnetic field strength is attenuated by both foreshortening
and limb darkening in approximately the same way as the radial motion of the granular flow. This may explain why, even though the RVs and log(R'_{HK}) values were measured simultaneously from the same spectra, the correlation between the variations in log(R'_{HK}) and suppression of granular blueshift appears weak. In their study of long-term solar RV variations spanning over 8 years, Lanza et al. (2016) find a stronger correlation between log(R'_{HK}) and RV variations, with a Spearman coefficient of 0.357. This positive correlation is in agreement with previous studies of quiet late-type stars (Lovis et al. 2011; Gomes da Silva et al. 2012), and shows that the log(R'_{HK}) may be a more useful proxy for long-term RV variations induced by the stellar magnetic cycle.

We note that in Fig. 10, the variations in log(R'_{HK}) look similar to the variations in RV, except that they are shifted by a few days (this is especially noticeable towards the end of the run). The reality and origin of this shift will be the subject of future studies, thanks to the wealth of Sun as a star RV observations that are currently being gathered by the solar telescope at HARPS-N.

## 5.3 Disc-averaged observed magnetic flux |B̂_{obs}|

We compute the full-disc line-of-sight magnetic flux of the Sun, by summing the intensity-weighted line-of-sight unsigned magnetic flux in each pixel:

$$ |\hat{B}_{\text{obs}}| = \frac{\sum_{ij} |B_{\text{obs},ij}| I_{ij}}{\sum_{ij} I_{ij}}. \qquad (30) $$

The values at the time of each HARPS observation are listed in Table A1. The variations in |B̂_{obs}| are shown in panel (b) of Fig. 11, together with the nightly averaged HARPS RV variations of the Sun as a star, in panel (a). We see that the variations in the disc-averaged magnetic flux are in phase with the RV variations, despite the scatter in RV in the first part of the run (discussed in Section 2.4.2). If we only consider the observations in the second part of the run, the Spearman correlation coefficient between ΔRV_{Sun} and |B̂_{obs}| is 0.80 (see Fig. 10). The correlation is stronger between v̂_{conv} and |B̂_{obs}|, with a correlation coefficient of 0.87. This is expected since magnetized areas are known to suppress convective blueshift (see Meunier et al. 2010a,b).

We note that these observations were taken close to the solar cycle maximum, during which the solar photosphere was mostly dominated by a few large sunspot groups (surrounded by facular networks). Meunier et al. (2010b) found that the convective shift attenuation is greater in larger structures, since they contain a stronger magnetic field. The relationship between convective shift attenuation and magnetic field, however, is not linear. We thus expect larger RV variations for a few large active regions than when the Sun is in a phase of lower magnetic activity, when the photosphere might be dominated by several smaller structures, even though they would still give the same total flux. In fact, Lanza et al. (2016) find a much weaker correlation between the mean total magnetic flux measured by the Synoptic Optical Long-term Investigations of the Sun (SOLIS) Vector Spectromagnetograph (VSM) instrument and RV (Spearman correlation coefficient of 0.131) over their 8-yr data set, which spans both active and less active phases of the solar cycle.

When compared against correlations with the traditional spectroscopic activity indicators (the FWHM, BIS and log(R'_{HK})), we see that the disc-averaged magnetic flux |B̂_{obs}| is a much more effective predictor of activity-induced RV variations, over the time-scale of a few rotation periods. The averaged magnetic flux may therefore be a useful proxy for activity-driven RV variations as it should map on to areas of strong magnetic fields, which suppress the Sun's con-

<figure><img src="image_13.png" /><figcaption>Figure 11. (a) HARPS RV variations of the Sun as a star. (b and c) Time series of the disc-averaged line-of-sight magnetic flux |B̂_{obs}| and filling factor f, respectively, determined from the SDO/HMI magnetograms. (d–f) Time series of the FWHM, BIS and log(R'_{HK}), respectively, determined from the HARPS DRS reduction pipeline. The values of these nightly binned time series are provided in the Supplementary Files available online.</figcaption></figure>

vective blueshift. The line-of-sight magnetic flux density and filling factor on the visible hemisphere of a star can be measured from the Zeeman broadening of magnetically sensitive lines (Robinson 1980; Reiners et al. 2013). Their product gives the disc-averaged flux density that we are deriving from the solar images. We note that such measurements are still very difficult to make for other stars than the Sun, because the Zeeman splitting of magnetically sensitive lines is so small that the technique can only be applied to bright, slowly rotating stars. Fortunately, such stars are also the best targets for planet searches.

## 5.4 Magnetic filling factor f

In addition to the disc-averaged observed magnetic flux, we also computed the filling factor of magnetic regions on the solar disc. It
is weighted by the foreshortening at the location of each pixel, and counted as a fraction of the total pixel count:

$$f = \frac{1}{N_{\text{pix}}} \sum_{ij} W_{ij}, \qquad (31)$$

where $N_{pix}$ is the total number of pixels in the solar disc and the weight $W_{ij}$ is set to 1 in magnetically active regions, and 0 in the quiet Sun. The variations of the filling factor are shown in panel (c) of Fig. 11. As expected, they follow the disc-averaged magnetic flux closely. The correlation with the predicted $\Delta RV_{\text{conv}}$ and the observed $\Delta RV_{\text{Sun}}$ is none the less weaker than that found for the brightness-weighted line-of-sight magnetic field $|\hat{B}_{\text{obs}}|$, since no correction is made for limb darkening.

## 6 CONCLUSION

In the present analysis, we decomposed activity-induced RV variations into identifiable contributions from sunspots, faculae and granulation, based on Sun as a star RV variations deduced from HARPS spectra of the bright asteroid Vesta and high spatial resolution images in the Fe I 6173 Å line taken with the HMI instrument aboard the SDO. We find that the RV variations induced by solar activity are mainly caused by the suppression of convective blueshift from magnetically active regions, while the flux deficit incurred by the presence of sunspots on the rotating solar disc only plays a minor role. We further compute the disc-averaged line-of-sight magnetic flux and show that although we cannot yet measure it with precision on other stars at present, it is a very good proxy for activity-driven RV variations, much more so than the FWHM and BIS of the cross-correlation profile, and the Ca II H&K activity index. These findings are in agreement with the previous works of Meunier et al. (2010a,b).

In addition to the existing 2011 HARPS observations of sunlight scattered off Vesta, there will soon be a wealth of direct solar RV measurements taken with HARPS-N, which will be regularly fed sunlight through a small 2-inch telescope developed specifically for this purpose. A prototype for this is currently being commissioned at HARPS-N (see Dumusque et al. 2015; Glenday et al., in preparation). Gaining a deeper understanding of the physics at the heart of activity-driven RV variability will ultimately enable us to better model and remove this contribution from RV observations, thus revealing the planetary signals.

## ACKNOWLEDGEMENTS

We wish to thank the referee for their thoughtful comments, which have greatly helped improve the analysis presented in this paper. RDH gratefully acknowledges STFC studentship grant number ST/J500744/1, and a grant from the John Templeton Foundation. The opinions expressed in this publication are those of the authors and do not necessarily reflect the views of the John Templeton Foundation. ACC and RF acknowledge support from STFC consolidated grants numbers ST/J001651/1 and ST/M001296/1. JL acknowledges support from NASA Origins of the Solar System grant No. NNX13AH79G and from STFC grant ST/M001296/1. The Solar Dynamics Observatory was launched by NASA on 2011 February 11, as part of the Living With A Star programme. This research has made use of NASA's Astrophysics Data System Bibliographic Services. We also used data from Haywood et al. (2016).

## REFERENCES

* Aigrain S., Pont F., Zucker S., 2012, MNRAS, 419, 3147

* Allen C., 1973, Allen: Astrophysical Quantities, 3rd edn. Athlone Press, London

* Baranne A. et al., 1996, A&AS, 119, 373

* Boisse I. et al., 2009, A&A, 495, 959

* Boisse I., Bouchy F., Hébrard G., Bonfils X., Santos N., Vauclair S., 2011, A&A, 528, A4

* Bonfils X. et al., 2007, A&A, 474, 293

* Brandt P. N., Solanki S. K., 1990, A&A, 231, 221

* Carrington R. C., 1859, MNRAS, 19, 81

* Chapman G. A., Cookson A. M., Dobias J. J., Walton S. R., 2001, ApJ, 555, 462

* Collier Cameron A. et al., 2006, MNRAS, 373, 799

* Deming D., Plymate C., 1994, ApJ, 426, 382

* Deming D., Espenak F., Jennings D. E., Brault J. W., Wagner J., 1987, ApJ, 316, 771

* Desort M., Lagrange A. M., Galland F., Udry S., Mayor M., 2007, A&A, 473, 983

* Díaz R. F. et al., 2016, A&A, 585, 134

* Dumusque X. et al., 2011, A&A, 535, A55

* Dumusque X. et al., 2015, ApJ, 814, L21

* Evershed J., 1909, MNRAS, 69, 454

* Gomes da Silva J., Santos N. C., Bonfils X., Delfosse X., Forveille T., Udry S., Dumusque X., Lovis C., 2012, A&A, 541, A9

* Gray D. F., 2009, ApJ, 697, 1032

* Gregory P. C., 2011, MNRAS, 415, 2523

* Hall J. C., 2008, Living Rev. Sol. Phys., 5, 2

* Haywood R. D. et al., 2014, MNRAS, 443, 2517

* Haywood R. D. et al., 2016, Data underpinning - The Sun as a planet-host star: Proxies from SDO images for HARPS radial-velocity variations, St Andrews. Available at: http://dx.doi.org/10.17630/bb43e6a3-72e0-464c-9fdd-fbe5d3e56a09

* Huélamo N. et al., 2008, A&A, 489, L9

* Isaacson H., Fischer D., 2010, ApJ, 725, 875

* Jiménez A., Pallé P. L., Régulo C., Roca Cortes T., Isaak G. R., 1986, COSPAR and IAU, 6, 89

* Lagrange A. M., Desort M., Meunier N., 2010, A&A, 512, A38

* Lanza A. F., Molaro P., 2015, Exp. Astron., 39, 461

* Lanza A. F., Molaro P., Monaco L., Haywood R. D., 2016, A&A, preprint (arXiv:1601.05646)

* Lean J., 1997, ARA&A, 35, 33

* Lindegren L., Dravins D., 2003, A&A, 401, 1185

* Loughhead R. E., Bray R. J., 1958, Aust. J. Phys., 11, 177

* Lovis C., Pepe F., 2007, A&A, 468, 1115

* Lovis C., Dumusque X., Santos N. C., Udry S., Mayor M., 2011, AAS/Division for Extreme Solar Systems Abstracts, 2, 0202

* McMillan R. S., Moore T. L., Perry M. L., Smith P. H., 1993, ApJ, 403, 801

* Meunier N., Desort M., Lagrange A. M., 2010a, A&A, 512, A39

* Meunier N., Lagrange A.-M., Desort M., 2010b, A&A, 519, A66

* Meunier N., Lagrange A.-M., Borgniet S., Rieutord M., 2015, A&A, 583, A118

* Molaro P., Centurión M., 2010, A&A, 525, A74

* Molaro P., Monaco L., Barbieri M., Zaggia S., 2013, MNRAS, 429, L79

* Murdin P., 2002, Murdin P., ed., Encyclopedia of Astronomy and Astrophysics. Inst. Phys. Pub., Bristol

* Pesnell W. D., Thompson B. J., Chamberlin P. C., 2012, Sol. Phys., 275, 3

* Queloz D. et al., 2001, A&A, 379, 279

* Reiners A., Shulyak D., Anglada-Escudé G., Jeffers S. V., Morin J., Zechmeister M., Kochukhov O., Piskunov N., 2013, A&A, 552, 103

* Robertson P., Mahadevan S., Endl M., Roy A., 2014, Science, 345, 440

* Robinson R. D., Jr, 1980, ApJ, 239, 961

* Russell C. T. et al., 2012, Science, 336, 684

* Santos N. C., Gomes da Silva J., Lovis C., Melo C., 2010, A&A, 511, A54

* Santos N. C. et al., 2014, A&A, 566, 35

* Schou J. et al., 2012, Sol. Phys., 275, 229
Schrijver C. J., 2002, Astron. Nachr., 323, 157
Snodgrass H. B., Ulrich R. K., 1990, ApJ, 351, 309
Stephenson C. B., 1951, ApJ, 114, 500
Thomas P. C., Binzel R. P., Gaffey M. J., Zellner B. H., Storrs A. D., Wells
E., 1997, Icarus, 128, 88
Thompson W. T., 2006, A&A, 449, 791
Yeo K. L., Solanki S. K., Krivova N. A., 2013, A&A, 550, A95

Solar Dynamics Observatory. The HARPS observations, together with tables for the results presented in this paper, are available in electronic format at: The SDO/HMI images can be downloaded from http://jsoc.stanford.edu

http://dx.doi.org/10.17630/
bb43e6a3-72e0-464c-9fdd-fbe5d3e56a09

Please note: Oxford University Press are not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.

SUPPORTING INFORMATION

Additional Supporting Information may be found in the online version of this article:

Based on observations made with the HARPS instrument on the 3.6 m telescope under the programme ID 088.C-0323 at Cerro La Silla (Chile) and the Helioseismic and Magnetic Imager on board the

APPENDIX A:

<table><caption>**Table A1.** HARPS 2011–2012 data for the asteroid Vesta, processed by the HARPS pipeline with the correct barycentric RV and accounting for the relativistic correction, and SDO-derived quantities resulting from our analysis. From left to right are given: Julian date (flux-weighted mid-exposure times of observation), RV, the estimated error $\sigma_{RV}$ on the RV measurement, the full width at half-maximum (FWHM) and the line bisector of (BIS) of the cross-correlation function (as defined in Queloz et al. 2001), the CaII H&K activity indicator $\log(R'_{HK})$ and its error $\sigma_{\log(R'_{HK})}$, the time lag $\Delta t$ between the HARPS observations and the time at which SDO observed the same hemisphere of the Sun, the apparent planetographic longitude of Vesta $\lambda$, the values of the basis functions for the sunspot velocity signal, $\hat{v}_{phot}$ and the suppression of granular blueshift, $\hat{v}_{conv}$, the full-disc magnetic flux $|\hat{B}_{obs}|$ and the filling factor of magnetic regions $f$ both computed from the SDO images. An expanded version of this table, together with all Supplementary Files, is available electronically at http://dx.doi.org/10.17630/bb43e6a3-72e0-464c-9fdd-fbe5d3e56a09.</caption><thead><tr><th>Julian date<br>(MJD<sub>mid-UTC</sub>)</th><th>RV<br>(km s<sup>-1</sup>)</th><th>σ<sub>RV</sub><br>(km s<sup>-1</sup>)</th><th>FWHM<br>(km s<sup>-1</sup>)</th><th>BIS<br>(km s<sup>-1</sup>)</th><th>log(R'<sub>HK</sub>)</th><th>σ<sub>log(R'<sub>HK</sub>)</sub></th><th>Δt<br>(d)</th><th>λ<br>(deg)</th><th>ν<sub>phot</sub><br>(m s<sup>-1</sup>)</th><th>ν<sub>conv</sub><br>(m s<sup>-1</sup>)</th><th>|B̂<sub>obs</sub>|<br>(G)</th><th>f<br>(per cent)</th></tr></thead><tbody><tr><td>2455834.48296192</td><td>0.085 70</td><td>0.000 43</td><td>7.089 02</td><td>-0.024 17</td><td>-4.9411</td><td>0.0024</td><td>2.8365</td><td>143.99</td><td>-0.122 949</td><td>6.612 382</td><td>13.390 075</td><td>6.442 092</td></tr><tr><td>2455834.61321887</td><td>0.099 36</td><td>0.000 39</td><td>7.0883</td><td>-0.025 59</td><td>-4.9435</td><td>0.0019</td><td>2.8451</td><td>351.75</td><td>-0.122 949</td><td>6.612 382</td><td>13.390 075</td><td>6.442 092</td></tr><tr><td>2455835.53355275</td><td>0.094 77</td><td>0.000 49</td><td>7.089 03</td><td>-0.024 78</td><td>-4.9513</td><td>0.0028</td><td>2.8901</td><td>39.73</td><td>-0.372 071</td><td>6.728 378</td><td>13.461 877</td><td>6.553 660</td></tr><tr><td>2455835.66225392</td><td>0.099 36</td><td>0.000 60</td><td>7.090 69</td><td>-0.024 77</td><td>-4.9498</td><td>0.0044</td><td>2.9002</td><td>247.49</td><td>-0.369 973</td><td>6.800 018</td><td>13.456 769</td><td>6.559 029</td></tr><tr><td>2455836.50488464</td><td>0.084 31</td><td>0.000 46</td><td>7.098 78</td><td>-0.024 77</td><td>-4.9450</td><td>0.0027</td><td>2.9396</td><td>171.93</td><td>-0.148 130</td><td>6.458 323</td><td>13.693 823</td><td>6.750 327</td></tr><tr><td>2455836.63463460</td><td>0.102 54</td><td>0.000 47</td><td>7.087 21</td><td>-0.025 89</td><td>-4.9506</td><td>0.0029</td><td>2.9487</td><td>19.69</td><td>-0.149 324</td><td>6.451 065</td><td>13.612 058</td><td>6.682 767</td></tr><tr><td>2455837.58225198</td><td>0.099 16</td><td>0.000 43</td><td>7.086 32</td><td>-0.024 37</td><td>-4.9575</td><td>0.0026</td><td>2.9941</td><td>112.57</td><td>-0.062 027</td><td>7.066 728</td><td>13.760 690</td><td>6.785 311</td></tr><tr><td>2455837.66801992</td><td>0.101 48</td><td>0.000 54</td><td>7.089 11</td><td>-0.026 32</td><td>-4.9597</td><td>0.0041</td><td>3.0056</td><td>252.95</td><td>-0.062 027</td><td>7.066 728</td><td>13.760 690</td><td>6.785 311</td></tr><tr><td>2455838.56557661</td><td>0.095 52</td><td>0.000 65</td><td>7.093 28</td><td>-0.026 98</td><td>-4.9373</td><td>0.0043</td><td>3.0525</td><td>261.60</td><td>-0.104 478</td><td>7.448 517</td><td>13.797 603</td><td>6.798 914</td></tr><tr><td>2455838.67384272</td><td>0.099 25</td><td>0.000 77</td><td>7.0939</td><td>-0.027 17</td><td>-4.9475</td><td>0.0061</td><td>3.0623</td><td>81.28</td><td>-0.104 478</td><td>7.448 517</td><td>13.797 603</td><td>6.798 914</td></tr><tr><td>2455839.56369186</td><td>0.106 00</td><td>0.000 49</td><td>7.090 53</td><td>-0.023 31</td><td>-4.9480</td><td>0.0031</td><td>3.1099</td><td>78.70</td><td>0.282 247</td><td>6.290 580</td><td>13.538 572</td><td>6.711 449</td></tr><tr><td>2455839.66500316</td><td>0.101 44</td><td>0.000 54</td><td>7.089 58</td><td>-0.025 47</td><td>-4.9404</td><td>0.0038</td><td>3.1058</td><td>241.53</td><td>0.285 453</td><td>6.290 580</td><td>13.538 572</td><td>6.711 449</td></tr><tr><td>2455840.54096382</td><td>0.101 81</td><td>0.000 42</td><td>7.090 45</td><td>-0.025 73</td><td>-4.9519</td><td>0.0024</td><td>3.1535</td><td>216.48</td><td>0.492 760</td><td>6.355 797</td><td>13.337 171</td><td>6.545 630</td></tr><tr><td>2455841.51821591</td><td>0.097 43</td><td>0.000 44</td><td>7.0908</td><td>-0.024 34</td><td>-4.9443</td><td>0.0023</td><td>3.2110</td><td>359.87</td><td>0.548 444</td><td>6.054 414</td><td>13.267 717</td><td>6.419 010</td></tr><tr><td>2455842.54067654</td><td>0.101 56</td><td>0.000 58</td><td>7.093 98</td><td>-0.023 31</td><td>-4.9419</td><td>0.0041</td><td>3.2579</td><td>210.63</td><td>0.285 453</td><td>5.964 030</td><td>13.453 311</td><td>6.502 014</td></tr><tr><td>2455842.64940682</td><td>0.102 21</td><td>0.000 65</td><td>7.096 29</td><td>-0.022 21</td><td>-4.9521</td><td>0.0051</td><td>3.2673</td><td>30.31</td><td>0.285 453</td><td>5.964 030</td><td>13.453 311</td><td>6.502 014</td></tr><tr><td>2455843.53446953</td><td>0.101 35</td><td>0.000 42</td><td>7.097 55</td><td>-0.024 79</td><td>-4.9529</td><td>0.0023</td><td>3.3128</td><td>16.47</td><td>0.289 295</td><td>6.042 467</td><td>14.255 929</td><td>7.096 244</td></tr><tr><td>2455843.66473254</td><td>0.097 22</td><td>0.000 57</td><td>7.097 89</td><td>-0.0264</td><td>-4.9598</td><td>0.0043</td><td>3.3214</td><td>229.83</td><td>0.289 295</td><td>6.042 467</td><td>14.255 929</td><td>7.096 244</td></tr><tr><td>2455844.53056656</td><td>0.096 56</td><td>0.000 58</td><td>7.090 06</td><td>-0.024 79</td><td>-4.9538</td><td>0.0039</td><td>3.3722</td><td>187.91</td><td>0.338 963</td><td>7.427 578</td><td>15.036 871</td><td>7.600 587</td></tr><tr><td>2455844.66179197</td><td>0.102 92</td><td>0.000 50</td><td>7.095 19</td><td>-0.0258</td><td>-4.9513</td><td>0.0033</td><td>3.3799</td><td>41.27</td><td>0.412 640</td><td>7.942 586</td><td>16.157 705</td><td>7.742 510</td></tr><tr><td>2455845.54569082</td><td>0.101 43</td><td>0.000 53</td><td>7.093 54</td><td>-0.0214</td><td>-4.9472</td><td>0.0032</td><td>3.4265</td><td>33.04</td><td>0.214 037</td><td>9.072 355</td><td>15.656 077</td><td>8.091 010</td></tr><tr><td>2455845.65206487</td><td>0.098 30</td><td>0.000 54</td><td>7.097 19</td><td>-0.023 43</td><td>-4.9456</td><td>0.0036</td><td>3.4313</td><td>201.48</td><td>0.147 598</td><td>9.104 808</td><td>15.530 288</td><td>7.938 438</td></tr><tr><td>2455870.51420322</td><td>0.108 43</td><td>0.000 59</td><td>7.098 49</td><td>-0.019 02</td><td>-4.9355</td><td>0.0032</td><td>4.7844</td><td>86.70</td><td>-0.230 017</td><td>9.870 595</td><td>18.187 723</td><td>9.696 133</td></tr><tr><td>2455870.52215421</td><td>0.107 01</td><td>0.000 53</td><td>7.101 07</td><td>-0.021 57</td><td>-4.9399</td><td>0.0035</td><td>4.7834</td><td>97.93</td><td>-0.230 017</td><td>9.870 595</td><td>18.187 723</td><td>9.696 133</td></tr><tr><td>2455870.53060607</td><td>0.103 79</td><td>0.000 53</td><td>7.098 26</td><td>-0.022 32</td><td>-4.9419</td><td>0.0036</td><td>4.7819</td><td>109.16</td><td>-0.230 017</td><td>9.870 595</td><td>18.187 723</td><td>9.696 133</td></tr><tr><td>2455871.50280525</td><td>0.102 50</td><td>0.000 61</td><td>7.093 61</td><td>-0.018 83</td><td>-4.9341</td><td>0.0033</td><td>4.8375</td><td>241.14</td><td>-0.247 643</td><td>9.628 828</td><td>17.731 211</td><td>9.520 507</td></tr><tr><td>2455871.51026399</td><td>0.101 90</td><td>0.000 53</td><td>7.097 64</td><td>-0.020 95</td><td>-4.9334</td><td>0.0033</td><td>4.8370</td><td>252.37</td><td>-0.247 643</td><td>9.628 828</td><td>17.731 211</td><td>9.520 507</td></tr><tr><td>2455871.51752074</td><td>0.102 79</td><td>0.000 53</td><td>7.098 29</td><td>-0.020 04</td><td>-4.9307</td><td>0.0033</td><td>4.8366</td><td>269.21</td><td>-0.247 643</td><td>9.628 828</td><td>17.731 211</td><td>9.520 507</td></tr><tr><td>2455872.49503353</td><td>0.106 38</td><td>0.001 77</td><td>7.103 66</td><td>-0.023 63</td><td>-4.8856</td><td>0.0193</td><td>4.8869</td><td>46.81</td><td>-0.217 090</td><td>8.539 591</td><td>16.708 896</td><td>8.931 849</td></tr><tr><td>2455872.50361680</td><td>0.107 06</td><td>0.000 63</td><td>7.098 57</td><td>-0.0206</td><td>-4.9459</td><td>0.0044</td><td>4.8853</td><td>63.65</td><td>-0.212 479</td><td>8.339 659</td><td>16.718 279</td><td>8.928 224</td></tr><tr><td>2455872.51065080</td><td>0.104 33</td><td>0.000 66</td><td>7.093 74</td><td>-0.021 82</td><td>-4.9417</td><td>0.0046</td><td>4.8852</td><td>74.88</td><td>-0.232 778</td><td>8.259 415</td><td>16.704 289</td><td>8.916 693</td></tr><tr><td>2455872.51832551</td><td>0.104 31</td><td>0.000 66</td><td>7.094 74</td><td>-0.020 18</td><td>-4.9297</td><td>0.0044</td><td>4.8914</td><td>86.10</td><td>-0.218 351</td><td>8.044 784</td><td>16.715 321</td><td>8.926 018</td></tr><tr><td>2455873.50406872</td><td>0.098 70</td><td>0.000 50</td><td>7.091 24</td><td>-0.022 93</td><td>-4.9400</td><td>0.0030</td><td>4.9404</td><td>240.54</td><td>-0.159 657</td><td>7.999 471</td><td>15.740 116</td><td>8.388 542</td></tr><tr><td>2455873.51131777</td><td>0.099 37</td><td>0.000 50</td><td>7.092 52</td><td>-0.021 97</td><td>-4.9425</td><td>0.0031</td><td>4.9401</td><td>251.76</td><td>-0.149 055</td><td>8.115 433</td><td>15.749 984</td><td>8.389 838</td></tr><tr><td>2455873.51863681</td><td>0.099 33</td><td>0.000 51</td><td>7.096 35</td><td>-0.022 45</td><td>-4.9380</td><td>0.0030</td><td>4.9397</td><td>262.99</td><td>-0.137 447</td><td>8.065 520</td><td>15.735 783</td><td>8.370 678</td></tr><tr><td>2455874.50355099</td><td>0.101 33</td><td>0.000 77</td><td>7.096 52</td><td>-0.023 86</td><td>-4.9467</td><td>0.0051</td><td>4.9964</td><td>57.42</td><td>-0.289 008</td><td>8.449 905</td><td>15.011 817</td><td>7.870 854</td></tr><tr><td>2455874.51074342</td><td>0.102 94</td><td>0.000 74</td><td>7.092 61</td><td>-0.023 14</td><td>-4.9402</td><td>0.0057</td><td>4.9962</td><td>68.65</td><td>-0.305 588</td><td>8.072 379</td><td>14.993 994</td><td>7.858 010</td></tr><tr><td>2455874.51805876</td><td>0.101 32</td><td>0.000 89</td><td>7.097 96</td><td>-0.025 44</td><td>-4.9379</td><td>0.0069</td><td>4.9958</td><td>79.87</td><td>-0.305 588</td><td>8.072 379</td><td>14.993 994</td><td>7.858 010</td></tr><tr><td>2455875.50408302</td><td>0.098 39</td><td>0.000 58</td><td>7.089 79</td><td>-0.024 18</td><td>-4.9544</td><td>0.0040</td><td>5.0515</td><td>234.30</td><td>-0.160 578</td><td>7.051 506</td><td>14.539 260</td><td>7.545 515</td></tr><tr><td>2455875.51155021</td><td>0.096 71</td><td>0.000 57</td><td>7.091 51</td><td>-0.023 84</td><td>-4.9472</td><td>0.0039</td><td>5.0509</td><td>245.52</td><td>-0.160 578</td><td>7.051 506</td><td>14.539 260</td><td>7.545 515</td></tr><tr><td>2455875.51886866</td><td>0.096 45</td><td>0.000 58</td><td>7.092 36</td><td>-0.019 25</td><td>-4.9500</td><td>0.0040</td><td>5.0506</td><td>256.75</td><td>-0.160 578</td><td>7.051 506</td><td>14.539 260</td><td>7.545 515</td></tr><tr><td>2455878.50606288</td><td>0.104 81</td><td>0.001 23</td><td>7.093 06</td><td>-0.023 84</td><td>-4.9368</td><td>0.0129</td><td>5.2162</td><td>44.91</td><td>0.239 367</td><td>7.862 843</td><td>14.083 289</td><td>7.241 949</td></tr><tr><td>2455878.51386630</td><td>0.103 96</td><td>0.001 25</td><td>7.087 15</td><td>-0.020 17</td><td>-4.9521</td><td>0.0137</td><td>5.2153</td><td>56.14</td><td>0.219 175</td><td>7.727 767</td><td>14.087 698</td><td>7.262 066</td></tr><tr><td>2455878.52076438</td><td>0.103 24</td><td>0.001 20</td><td>7.092 72</td><td>-0.026 94</td><td>-4.9529</td><td>0.0128</td><td>5.2153</td><td>67.36</td><td>0.215 755</td><td>7.792 477</td><td>14.102 572</td><td>7.257 334</td></tr></tbody></table>
<table><caption>Table A1 – continued.</caption><thead><tr><th>Julian date<br>(MJD<sub>mid</sub>_UTC)</th><th>RV<br>(km s<sup>-1</sup>)</th><th>σ<sub>RV</sub><br>(km s<sup>-1</sup>)</th><th>FWHM<br>(km s<sup>-1</sup>)</th><th>BIS<br>(km s<sup>-1</sup>)</th><th>log(R'<sub>HK</sub>)</th><th>σ<sub>log(R'<sub>HK</sub>)</sub></th><th>Δt<br>(d)</th><th>λ<br>(deg)</th><th>ν<sub>phot</sub><br>(m s<sup>-1</sup>)</th><th>ν<sub>conv</sub><br>(m s<sup>-1</sup>)</th><th>|B̂<sub>obs</sub>|<br>(G)</th><th>f<br>(per cent)</th></tr></thead><tbody><tr><td>2455879.50702795</td><td>0.097 25</td><td>0.001 12</td><td>7.094 57</td><td>-0.025 22</td><td>-4.9197</td><td>0.0103</td><td>5.2777</td><td>227.38</td><td>0.232 115</td><td>7.128 572</td><td>14.466 632</td><td>7.447 019</td></tr><tr><td>2455879.51427708</td><td>0.097 05</td><td>0.001 07</td><td>7.089 51</td><td>-0.028 34</td><td>-4.9360</td><td>0.0100</td><td>5.2774</td><td>238.61</td><td>0.230 304</td><td>6.983 342</td><td>14.483 479</td><td>7.457 781</td></tr><tr><td>2455879.52152229</td><td>0.097 81</td><td>0.001 05</td><td>7.094 96</td><td>-0.025 36</td><td>-4.9146</td><td>0.0093</td><td>5.2771</td><td>249.84</td><td>0.232 562</td><td>7.438 171</td><td>14.526 323</td><td>7.492 319</td></tr><tr><td>2455886.55971131</td><td>0.096 68</td><td>0.000 85</td><td>7.099 98</td><td>-0.021 19</td><td>-4.9537</td><td>0.0073</td><td>5.6625</td><td>109.52</td><td>0.232 784</td><td>6.416 208</td><td>14.173 657</td><td>6.983 742</td></tr><tr><td>2455886.56696080</td><td>0.094 96</td><td>0.000 86</td><td>7.094 83</td><td>-0.023 48</td><td>-4.9481</td><td>0.0075</td><td>5.6622</td><td>120.75</td><td>0.232 784</td><td>6.416 208</td><td>14.173 657</td><td>6.983 742</td></tr><tr><td>2455886.57392844</td><td>0.092 57</td><td>0.000 93</td><td>7.100 74</td><td>-0.0237</td><td>-4.9477</td><td>0.0086</td><td>5.6622</td><td>131.98</td><td>0.232 784</td><td>6.416 208</td><td>14.173 657</td><td>6.983 742</td></tr><tr><td>2455887.52342646</td><td>0.096 87</td><td>0.000 60</td><td>7.090 76</td><td>-0.023 92</td><td>-4.9636</td><td>0.0042</td><td>5.7196</td><td>224.60</td><td>0.254 368</td><td>6.578 471</td><td>13.681 434</td><td>6.611 017</td></tr><tr><td>2455887.53095415</td><td>0.097 45</td><td>0.000 67</td><td>7.091 91</td><td>-0.026 85</td><td>-4.9726</td><td>0.0052</td><td>5.7190</td><td>235.83</td><td>0.254 368</td><td>6.578 471</td><td>13.681 434</td><td>6.611 017</td></tr><tr><td>2455887.53806579</td><td>0.096 16</td><td>0.000 72</td><td>7.093 53</td><td>-0.026 17</td><td>-4.9578</td><td>0.0056</td><td>5.7189</td><td>247.06</td><td>0.254 368</td><td>6.578 471</td><td>13.681 434</td><td>6.611 017</td></tr><tr><td>2455888.52593242</td><td>0.101 10</td><td>0.000 60</td><td>7.090 39</td><td>-0.026 13</td><td>-4.9661</td><td>0.0043</td><td>5.7727</td><td>47.04</td><td>0.157 666</td><td>6.316 140</td><td>13.525 647</td><td>6.498 637</td></tr><tr><td>2455888.54049630</td><td>0.101 75</td><td>0.000 62</td><td>7.090 83</td><td>-0.025 48</td><td>-4.9577</td><td>0.0047</td><td>5.7720</td><td>69.50</td><td>0.157 666</td><td>6.316 140</td><td>13.525 647</td><td>6.498 637</td></tr><tr><td>2455889.52578074</td><td>0.097 55</td><td>0.000 60</td><td>7.091 31</td><td>-0.021 11</td><td>-4.9642</td><td>0.0045</td><td>5.8284</td><td>223.87</td><td>0.183 892</td><td>6.607 664</td><td>13.455 927</td><td>6.430 930</td></tr><tr><td>2455889.53310116</td><td>0.096 01</td><td>0.000 60</td><td>7.093 05</td><td>-0.024 11</td><td>-4.9633</td><td>0.0046</td><td>5.8280</td><td>235.10</td><td>0.183 892</td><td>6.607 664</td><td>13.455 927</td><td>6.430 930</td></tr><tr><td>2455889.54034195</td><td>0.094 61</td><td>0.000 65</td><td>7.093 77</td><td>-0.022 42</td><td>-4.9796</td><td>0.0054</td><td>5.8277</td><td>246.32</td><td>0.183 892</td><td>6.607 664</td><td>13.455 927</td><td>6.430 930</td></tr><tr><td>2455890.50917526</td><td>0.100 97</td><td>0.000 59</td><td>7.088 86</td><td>-0.028 67</td><td>-4.9505</td><td>0.0040</td><td>5.8867</td><td>12.62</td><td>0.079 914</td><td>6.820 067</td><td>13.946 109</td><td>6.891 653</td></tr><tr><td>2455890.51649149</td><td>0.100 74</td><td>0.000 58</td><td>7.087 63</td><td>-0.026 59</td><td>-4.9564</td><td>0.0040</td><td>5.8863</td><td>23.85</td><td>0.079 914</td><td>6.820 067</td><td>13.946 109</td><td>6.891 653</td></tr><tr><td>2455890.52382222</td><td>0.099 96</td><td>0.000 59</td><td>7.087 79</td><td>-0.024 85</td><td>-4.9611</td><td>0.0042</td><td>5.8859</td><td>35.08</td><td>0.079 914</td><td>6.820 067</td><td>13.946 109</td><td>6.891 653</td></tr><tr><td>2455891.52246765</td><td>0.099 41</td><td>0.000 60</td><td>7.091 51</td><td>-0.025 35</td><td>-4.9591</td><td>0.0044</td><td>5.9428</td><td>211.90</td><td>0.079 914</td><td>6.820 067</td><td>13.946 109</td><td>6.891 653</td></tr><tr><td>2455891.52978471</td><td>0.099 21</td><td>0.000 63</td><td>7.0924</td><td>-0.025 49</td><td>-4.9649</td><td>0.0049</td><td>5.9424</td><td>223.12</td><td>-0.128 128</td><td>7.724 229</td><td>15.031 973</td><td>7.761 893</td></tr><tr><td>2455891.53703495</td><td>0.098 40</td><td>0.000 63</td><td>7.093 53</td><td>-0.026 16</td><td>-4.9482</td><td>0.0049</td><td>5.9421</td><td>234.35</td><td>-0.128 128</td><td>7.724 229</td><td>15.031 973</td><td>7.761 893</td></tr><tr><td>2455892.52735539</td><td>0.103 98</td><td>0.000 61</td><td>7.086 42</td><td>-0.022 54</td><td>-4.9442</td><td>0.0045</td><td>5.9935</td><td>34.33</td><td>0.250 746</td><td>8.037 871</td><td>15.638 486</td><td>7.970 436</td></tr><tr><td>2455892.53453612</td><td>0.103 58</td><td>0.000 79</td><td>7.091 45</td><td>-0.024 41</td><td>-4.9370</td><td>0.0068</td><td>5.9932</td><td>45.55</td><td>0.239 745</td><td>8.154 407</td><td>15.602 640</td><td>7.956 911</td></tr><tr><td>2455892.54214391</td><td>0.102 47</td><td>0.000 69</td><td>7.0878</td><td>-0.023 79</td><td>-4.9422</td><td>0.0054</td><td>5.9995</td><td>62.40</td><td>0.239 745</td><td>8.154 407</td><td>15.602 640</td><td>7.956 911</td></tr><tr><td>2455893.51437423</td><td>0.100 22</td><td>0.000 59</td><td>7.095 16</td><td>-0.024 93</td><td>-4.9440</td><td>0.0040</td><td>6.0551</td><td>194.30</td><td>0.178 788</td><td>9.501 673</td><td>16.492 246</td><td>8.578 127</td></tr><tr><td>2455893.52626540</td><td>0.102 39</td><td>0.000 57</td><td>7.091 83</td><td>-0.025 14</td><td>-4.9578</td><td>0.0040</td><td>6.0501</td><td>211.14</td><td>0.178 788</td><td>9.501 673</td><td>16.492 246</td><td>8.578 127</td></tr><tr><td>2455893.53359078</td><td>0.102 28</td><td>0.000 60</td><td>7.096 83</td><td>-0.024 59</td><td>-4.9607</td><td>0.0044</td><td>6.0497</td><td>222.37</td><td>0.178 788</td><td>9.501 673</td><td>16.492 246</td><td>8.578 127</td></tr><tr><td>2455894.50541095</td><td>0.102 65</td><td>0.000 59</td><td>7.094 96</td><td>-0.023 71</td><td>-4.9499</td><td>0.0041</td><td>6.0988</td><td>354.27</td><td>-0.143 119</td><td>10.188 700</td><td>16.843 670</td><td>8.978 122</td></tr><tr><td>2455894.51266019</td><td>0.104 00</td><td>0.000 62</td><td>7.093 07</td><td>-0.023 55</td><td>-4.9416</td><td>0.0044</td><td>6.1054</td><td>5.49</td><td>-0.143 119</td><td>10.188 700</td><td>16.843 670</td><td>8.978 122</td></tr><tr><td>2455894.52011844</td><td>0.105 48</td><td>0.000 68</td><td>7.094 66</td><td>-0.023 04</td><td>-4.9486</td><td>0.0052</td><td>6.1049</td><td>16.72</td><td>-0.143 119</td><td>10.188 700</td><td>16.843 670</td><td>8.978 122</td></tr><tr><td>2455895.51190055</td><td>0.098 55</td><td>0.000 58</td><td>7.0944</td><td>-0.0217</td><td>-4.9341</td><td>0.0038</td><td>6.1617</td><td>182.30</td><td>-0.219 694</td><td>9.459 038</td><td>16.678 771</td><td>8.985 445</td></tr><tr><td>2455895.51928533</td><td>0.100 14</td><td>0.000 61</td><td>7.093 07</td><td>-0.022 48</td><td>-4.9339</td><td>0.0041</td><td>6.1613</td><td>193.53</td><td>-0.219 694</td><td>9.459 038</td><td>16.678 771</td><td>8.985 445</td></tr><tr><td>2455895.52668202</td><td>0.103 85</td><td>0.000 60</td><td>7.094 09</td><td>-0.020 13</td><td>-4.9456</td><td>0.0041</td><td>6.1608</td><td>204.76</td><td>-0.219 694</td><td>9.459 038</td><td>16.678 771</td><td>8.985 445</td></tr><tr><td>2455896.51221054</td><td>0.103 11</td><td>0.000 62</td><td>7.095 55</td><td>-0.0246</td><td>-4.9082</td><td>0.0042</td><td>6.2170</td><td>359.11</td><td>-0.118 240</td><td>9.149 870</td><td>16.324 338</td><td>8.893 371</td></tr><tr><td>2455896.51945075</td><td>0.104 02</td><td>0.000 64</td><td>7.0943</td><td>-0.021 61</td><td>-4.9219</td><td>0.0044</td><td>6.2167</td><td>10.34</td><td>-0.118 240</td><td>9.149 870</td><td>16.324 338</td><td>8.893 371</td></tr><tr><td>2455896.52676887</td><td>0.104 50</td><td>0.000 65</td><td>7.0945</td><td>-0.0215</td><td>-4.9086</td><td>0.0045</td><td>6.2163</td><td>21.56</td><td>-0.118 240</td><td>9.149 870</td><td>16.324 338</td><td>8.893 371</td></tr><tr><td>2455897.51313557</td><td>0.098 05</td><td>0.000 57</td><td>7.0929</td><td>-0.021 41</td><td>-4.9279</td><td>0.0041</td><td>6.2716</td><td>175.91</td><td>0.310 313</td><td>8.521 387</td><td>15.731 010</td><td>8.518 003</td></tr><tr><td>2455897.52017185</td><td>0.101 99</td><td>0.000 61</td><td>7.090 74</td><td>-0.020 29</td><td>-4.9207</td><td>0.0045</td><td>6.2715</td><td>187.14</td><td>0.310 313</td><td>8.521 387</td><td>15.731 010</td><td>8.518 003</td></tr><tr><td>2455897.52797234</td><td>0.101 03</td><td>0.000 66</td><td>7.094 31</td><td>-0.022 93</td><td>-4.9178</td><td>0.0049</td><td>6.2776</td><td>203.98</td><td>0.310 313</td><td>8.521 387</td><td>15.731 010</td><td>8.518 003</td></tr><tr><td>2455898.51251923</td><td>0.102 34</td><td>0.000 76</td><td>7.094 42</td><td>-0.025 29</td><td>-4.9139</td><td>0.0059</td><td>6.3278</td><td>352.71</td><td>0.208 771</td><td>7.346 064</td><td>15.023 078</td><td>7.994 100</td></tr><tr><td>2455898.51955262</td><td>0.102 92</td><td>0.000 74</td><td>7.095 72</td><td>-0.020 88</td><td>-4.9290</td><td>0.0060</td><td>6.3277</td><td>3.94</td><td>0.193 525</td><td>7.375 161</td><td>15.016 566</td><td>7.987 264</td></tr><tr><td>2455898.52692661</td><td>0.101 53</td><td>0.000 78</td><td>7.091 94</td><td>-0.018 18</td><td>-4.9312</td><td>0.0067</td><td>6.3272</td><td>15.17</td><td>0.193 006</td><td>7.321 089</td><td>15.020 403</td><td>7.988 140</td></tr><tr><td>2455899.51280365</td><td>0.094 28</td><td>0.000 77</td><td>7.092 21</td><td>-0.021 72</td><td>-4.9295</td><td>0.0065</td><td>6.3830</td><td>169.51</td><td>0.355 292</td><td>6.343 507</td><td>14.470 741</td><td>7.472 408</td></tr><tr><td>2455899.51984723</td><td>0.095 35</td><td>0.000 75</td><td>7.088 03</td><td>-0.025 29</td><td>-4.9478</td><td>0.0066</td><td>6.3829</td><td>180.74</td><td>0.360 379</td><td>6.548 910</td><td>14.477 587</td><td>7.468 980</td></tr><tr><td>2455899.52750987</td><td>0.096 53</td><td>0.000 76</td><td>7.093 12</td><td>-0.021 94</td><td>-4.9296</td><td>0.0063</td><td>6.3822</td><td>191.96</td><td>0.340 881</td><td>6.579 109</td><td>14.478 097</td><td>7.470 433</td></tr><tr><td>2455900.51316244</td><td>0.099 20</td><td>0.000 71</td><td>7.089 69</td><td>-0.024 18</td><td>-4.9457</td><td>0.0060</td><td>6.4382</td><td>346.30</td><td>0.270 016</td><td>6.616 161</td><td>13.979 841</td><td>7.170 614</td></tr><tr><td>2455900.52048381</td><td>0.097 37</td><td>0.000 71</td><td>7.091 26</td><td>-0.026 17</td><td>-4.9389</td><td>0.0059</td><td>6.4378</td><td>357.53</td><td>0.290 577</td><td>6.656 814</td><td>13.981 226</td><td>7.170 553</td></tr><tr><td>2455900.52779857</td><td>0.099 82</td><td>0.000 73</td><td>7.0925</td><td>-0.023 89</td><td>-4.9388</td><td>0.0062</td><td>6.4444</td><td>14.37</td><td>0.287 024</td><td>6.723 330</td><td>13.985 253</td><td>7.170 031</td></tr><tr><td>2455901.51376528</td><td>0.094 12</td><td>0.000 84</td><td>7.087 55</td><td>-0.026 95</td><td>-4.9244</td><td>0.0074</td><td>6.4932</td><td>163.10</td><td>0.527 470</td><td>5.870 234</td><td>13.607 637</td><td>6.908 622</td></tr><tr><td>2455901.52080259</td><td>0.094 32</td><td>0.000 79</td><td>7.094 11</td><td>-0.023 79</td><td>-4.9370</td><td>0.0068</td><td>6.4931</td><td>174.32</td><td>0.526 638</td><td>5.968 377</td><td>13.626 008</td><td>6.919 914</td></tr><tr><td>2455901.52818692</td><td>0.097 43</td><td>0.000 84</td><td>7.089 43</td><td>-0.019 41</td><td>-4.9387</td><td>0.0078</td><td>6.4996</td><td>191.16</td><td>0.528 856</td><td>5.908 731</td><td>13.682 738</td><td>6.957 918</td></tr><tr><td>2455902.51342158</td><td>0.098 49</td><td>0.001 00</td><td>7.095 81</td><td>-0.022 93</td><td>-4.9453</td><td>0.0099</td><td>6.5491</td><td>339.89</td><td>0.080 486</td><td>5.991 561</td><td>13.461 406</td><td>6.723 588</td></tr><tr><td>2455902.52087765</td><td>0.098 11</td><td>0.000 94</td><td>7.087 12</td><td>-0.023 49</td><td>-4.9608</td><td>0.0097</td><td>6.5555</td><td>356.73</td><td>0.072 490</td><td>6.135 274</td><td>13.455 276</td><td>6.727 012</td></tr><tr><td>2455902.52825791</td><td>0.099 28</td><td>0.000 96</td><td>7.086 18</td><td>-0.024 03</td><td>-4.9391</td><td>0.0096</td><td>6.5551</td><td>7.96</td><td>0.084 431</td><td>6.272 675</td><td>13.472 634</td><td>6.750 239</td></tr></tbody></table>

This paper has been typeset from a *TEX/LATEX* file prepared by the author.