# Radial-velocity variations due to meridional flows in the Sun and solar-type stars: impact on exoplanet detectability

Nadège Meunier, A.-M. Lagrange

► To cite this version:

Nadège Meunier, A.-M. Lagrange. Radial-velocity variations due to meridional flows in the Sun and solar-type stars: impact on exoplanet detectability. Astronomy and Astrophysics - A&A, 2020, 638, pp.A54. 10.1051/0004-6361/201937354. hal-02903835

HAL Id: hal-02903835

https://hal.science/hal-02903835

Submitted on 21 Jul 2020

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
# Radial-velocity variations due to meridional flows in the Sun and solar-type stars: impact on exoplanet detectability

N. Meunier and A.-M. Lagrange

Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
e-mail: nadege.meunier@univ-grenoble-alpes.fr

Received 19 December 2019 / Accepted 7 April 2020

## ABSTRACT

**Context.** Stellar variability due to magnetic activity and flows at different spatial scales strongly impacts radial velocities. This variability is seen as oscillations, granulation, supergranulation, and meridional flows. The effect of this latter process is currently poorly known but could affect exoplanet detectability.

**Aims.** We aim to quantify the amplitude of the meridional flow integrated over the disc and its temporal variability, first for the Sun, as seen with different inclinations, and then for other solar-type stars. We then want to compare these amplitudes with low-mass exoplanetary amplitudes in radial velocity.

**Methods.** We used long time series (covering two 11-yr cycles) of solar latitudinal meridional circulation to reconstruct its integrated contribution and study its properties. We then used scaling laws from hydrodynamical simulations relating the amplitude of the meridional flow variability with stellar mass and rotation rate to estimate the typical amplitude expected for other solar-type stars.

**Results.** We find typical rms of the order of 0.5–0.7 m s⁻¹ (edge-on) and 1.2–1.7 m s⁻¹ (pole-on) for the Sun (peak-to-peak amplitudes are typically 1–1.4 m s⁻¹ and 2.3–3.3 m s⁻¹ resp.), with a minimal jitter for an inclination of 45–55°. This signal is significant compared to other stellar activity contributions and is much larger than the radial-velocity signal of the Earth. The variability is strongly related to the activity cycle, with maximum flows during the descending phase of the cycle, and possible variability on timescales lower than the cycle period. Extension to other solar-type stars shows that the variability due to meridional flows is dominated by the amplitude of the cycle of those stars (compared with mass and rotation rate), and that the peak-to-peak amplitudes can reach 4 m s⁻¹ for the most variable stars when seen pole-on. The meridional flow contribution sometimes represents a high fraction of the convective blueshift inhibition signal, especially for quiet, low-mass stars. For fast-rotating stars, the presence of multi-cellular patterns should significantly decrease the meridional flow contribution to the radial-velocity signal.

**Conclusions.** Our study shows that these meridional flows could be critical for exoplanet detection. Low inclinations are more impacted than edge-on configurations, but these latter still exhibit significant variability. Meridional flows also degrade the correlation between radial velocities due to convective blueshift inhibition and chromospheric activity indicators. This will make the correction from this signal challenging for stars with no multi-cellular patterns, such as the Sun for example, although there may be some configurations for which the line shape variations may be used if the precision is sufficient.

**Key words.** Sun: activity – Sun: photosphere – techniques: radial velocities – stars: kinematics and dynamics – stars: activity – planetary systems

## 1. Introduction

Stellar variability strongly impacts exoplanet radial-velocity (RV) detection and mass characterisation in different ways. Firstly, it is affected by stellar magnetic activity (e.g. Saar & Donahue 1997; Hatzes 2002; Saar et al. 2003; Wright 2005; Desort et al. 2007; Lagrange et al. 2010; Meunier et al. 2010a; Boisse et al. 2012; Dumusque et al. 2014; Borgniet et al. 2015; Dumusque 2016). Small-scale dynamics such as oscillations, granulation, and supergranulation also affect RVs (Dumusque et al. 2011; Cegla et al. 2013, 2018; Meunier et al. 2015; Meunier & Lagrange 2019a; Cegla et al. 2019). All of these effects have been studied. Much less attention has been paid to the effect of large-scale flows such as meridional circulation. Meridional flows, as in differential rotation, arise from the redistribution of angular momentum caused by turbulence. The only estimation of integrated RV that we are aware of was performed by Makarov (2010). This latter author simulated the solar RV over a cycle, including magnetic activity (and in particular the convective blueshift inhibition in plages) and meridional flows.

Makarov (2010) implemented this simulation for the Sun seen edge-on only, and the amplitude of the contribution of the meridional flows is not entirely clear since it was superposed to the other contributions. He derived a total amplitude (including convective blueshift inhibition) of about 7 m s⁻¹. We expect that the contribution of the meridional flows he simulated is lower than 1–2 m s⁻¹ assuming the full amplitude he found (7 m s⁻¹) is compatible with the results of Meunier et al. (2010b), who did not take meridional flows into account.

Here, we aim to quantify the amplitude and timescale of the meridional flow contribution to the RV signal, and to compare the former with the amplitude of an Earth-like planet. We also expect a reversal at a certain inclination because the meridional flows are predominantly poleward, and therefore taking stellar inclination into account is crucial. Long time series of high-quality solar measurements can be used to estimate the effect of meridional flows for the Sun seen with different inclinations. In addition, we use recent results from 3D hydrodynamical simulations to extrapolate the solar values to stars other than the Sun with different spectral types and activity levels.
<table><caption>**Table 1.** Input meridional circulation latitudinal profiles.</caption><thead><tr><th>Reference</th><th>Name</th><th>Temporal coverage</th><th>Nb of points</th><th>Analysis</th></tr></thead><tbody><tr><td>Ulrich (2010)</td><td>U10</td><td>1986-2009</td><td>24</td><td>Dopplergrams</td></tr><tr><td>Basu & Antia (2010)</td><td>BA10</td><td>1996-2009</td><td>15</td><td>Ring diagram</td></tr><tr><td>Hathaway & Rightmire (2010)</td><td>HR10</td><td>1996-2009</td><td>166</td><td>Magnetic feature tracking</td></tr><tr><td>Hathaway & Rightmire (2011)</td><td>HR11</td><td>1996-2010</td><td>180</td><td>Magnetic feature tracking</td></tr></tbody></table>

The outline of this paper is as follows. We first consider the solar case in Sect. 2 for various inclinations and solar latitudinal profiles. The resulting RV contribution is compared to a planetary signal in the habitable zone and other solar activity contributions. These results are then extrapolated to other solar-type stars in Sect. 3, and we conclude in Sect. 4.

## 2. The solar case

We use published solar meridional flow profiles versus latitude and time to build the integrated meridional circulation (MC) at each time step, mc(t), for different solar inclinations. We first describe the latitudinal profiles used as input and then describe the procedure used to compute the disc integration and the results.

### 2.1. Input observed meridional circulation

There is a long history of solar MC measurements using various techniques (e.g. Duvall 1979; Labonte & Howard 1982; Howard & Gilman 1986; Ulrich et al. 1988; Komm et al. 1993; Snodgrass & Dailey 1996; Hathaway 1996; Nesme-Ribes et al. 1997; Meunier 1999), with predominantly poleward flows. Meunier (1999) observed weaker flows during the ascending phase of the cycle compared to solar minimum, already seen by Komm et al. (1993) with larger uncertainties, as well as a converging flow pattern around the activity belt which is not present during the minimum. This pattern was later observed in other studies using different techniques (e.g. Zhao & Kosovichev 2004; Lin & Chou 2018), while some diverging patterns have also been observed (Snodgrass & Dailey 1996; Chou & Dai 2001).

Here, we select monitorings of MC with good temporal coverage of typically about a solar cycle or more so as to have a complete view of its temporal variability along the cycle (see Table 1). We focus on the time series obtained by Ulrich (2010) because MC derived from full-disc Dopplergrams is the most interesting for our purposes, as it corresponds to flows formed at an altitude in the photosphere that is very similar to the altitude we are interested in when observing stars in the optical. The different time series obtained with the other inputs are described and compared in Appendix A to show the robustness of the results. Another relatively long time series by Gizon et al. (2010) is available but we do not include it here because it is slightly shorter than the others and is the combination of two short time series obtained with different techniques.

Ulrich (2010; hereafter U10) measured the meridional flows with a very good temporal coverage (1986–2009, i.e. covering two solar cycles and including three activity minima) from Dopplergrams obtained at the Mount Wilson Observatory. Each yearly latitudinal profile has been digitised (the uncertainty on this numerisation is much smaller than the error bars on the meridional flow estimates) and is shown in Fig. 1. Some profiles

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Latitudinal profiles of meridional flows from U10 for cycles 22 (upper panel) and 23 (middle panel). Positive flows are northwards. The colour code indicates time, from orange (beginning of each cycle) to blue (end of each cycle), equally spread in cycle phase. Lower panel: uncertainties vs. latitude (one curve per year).</figcaption></figure>

show sign reversals at high latitudes, indicating the presence of a second cell in the circulation. Ulrich (2010) argues that they are of solar origin because he took effects such as the limb convective blueshift into account in the analysis, although this was discussed by Hathaway (2012) who casts some doubt on the reality of the flows. Differences between the results of these latter
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Examples of maps (**left panels**) and corresponding velocity distributions (**right panels**) over the solar disc for inclinations of 0°, 30°, 60°, and 90° (**from top to bottom**), for the 1986 latitudinal profile of Ulrich (2010). The levels are in metres per second, and are positive for redshifts.</figcaption></figure>

two authors could be due to the fact that they use different methods probing different layers. Small features across the equator are also present for a few years, but they are not likely to be of solar origin: U10 argues that these could be an artefact coming from the data processing (lack of equatorial symmetry at higher latitude). Uncertainties on the latitudinal profiles are typically around 0.5–1 m s⁻¹ at medium latitudes and can reach several metres per second at high latitudes (and up to 2 m s⁻¹ close to the equator).

## 2.2. Building of the integrated meridional circulation

We subsequently consider the MC latitudinal profile at a given time step. We build a map of the flows over the solar disc according to this profile, taking projection effects into account, for a given solar inclination. Figure 2 shows examples of such maps for different inclinations, together with the velocity distributions, which have a complex shape. We then integrate the flows over the disc, after weighting them with the limb-darkening function of Claret & Hauschildt (2003). The resulting integrated mc(t) is then obtained for inclinations between 0° (pole-on configuration)

and 90° (edge-on configuration), with a step of 5°, for each date. We do not weight the flows with the difference in flux in spots and plages because the impact should be small, especially as we consider long-term variability; in addition, they are located in converging flow regions which should partially compensate each other, hence the lower contribution to the integrated MC variability from these areas.

The estimation of solar meridional circulation is usually not possible at high latitude because of strong projection effects. The maximum latitude depends on the technique but is typically in the 50–80° range. For higher latitudes, we interpolated the published flows between the MC value at the highest latitude and a MC value of zero at the pole (separately in each hemisphere). This is justified by the necessity to have a single turning point for the flows, and by the results of Hathaway (2012) who observed significant meridional flows up to the poles.

In addition to the integrated meridional flows mc(t), we study the line shape variability over time. For this purpose, we attribute to each pixel of the disc a spectrum consisting in one spectral line which is shifted according to the projected MC flow at this pixel and projected rotational velocity. The spectral lines from all pixels are then combined (after weighting by the limb-darkening function). We use the synthetic solar optical spectra used in the SAFIR software (Galland et al. 2005) from Kurucz (1993), and use only a single line here for simplification, at 5300 Å, i.e. in approximately the middle of the HARPS spectrograph wavelength coverage. We then compute the bisector of this line after integration over the disc to check whether meridional flows leave a signature in the line shape. To estimate the variation of the bisector shape, we compute the variable BIS at each time step, which is the difference between two quantities, $\lambda_1$ and $\lambda_2$, after conversion to a velocity. These quantities are computed as follows. We use $I$ to denote the intensity at the line core and $I'$ the line depth (1-I). Values of $\lambda_1$ and $\lambda_2$ are then the wavelengths of the middle of the line profile for intensities between $I$ and $I+0.2I'$ and between $I+0.7I'$ and $I+0.9I'$ respectively.

Finally, the time series will be compared to two activity indicators: the spot number from SIDC (Solar Influences Data analysis Center) and the S-index from the Sacramento Peak Observatory. The objective is to compare the phase of the different series during the cycle.

## 2.3. Results

In this section, we first present the time series obtained from this reconstruction. We characterise the amplitude of the variability for different inclinations. Finally, we compare the expected MC signal with an Earth-like signal in the habitable zone and with the contribution due to stellar magnetic activity and other stellar flows.

### 2.3.1. Time series

The time series for input latitudinal profiles from U10 and two configurations (edge-on and pole-on) are shown in Fig. 3. The integrated meridional flow is on average negative when the star is seen pole-on, and positive when it is seen edge-on, because the flows are predominantly poleward, with an average of respectively –8.0 m s⁻¹ and 4.3 m s⁻¹. The MC time series are globally linked to the cycle, but with a complex relationship: from Fig. 3, the comparison of the pole-on time series with activity indicators shows that the meridional circulation is lower at cycle maximum compared to cycle minimum; furthermore, the meridional circulation (in absolute value) is larger during the descending phase
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Time series of integrated meridional flows for pole-on configuration (first panel) and edge-on configuration (second panel) for U10 reconstruction. Two types of uncertainties are shown (see text): a lower limit (red) and an upper limit (black). Third panel: S-index from Sacramento Peak Observatory between 1985 and 2011 (the solid line corresponds to a smoothing over 100 days). Last panel: monthly spot number (stars), superposed on a smoothed time series as a solid line.</figcaption></figure>

of the cycle, as illustrated in Fig. 4. There is also a significant difference between the two cycles: for cycle 22, the maximum is closer to the maximum of the cycle, while for cycle 23 the maximum is closer to the end of the cycle and their amplitudes appear to be different. We attempt to estimate the shift between the activity cycle and MC (derived from U10) by computing the cross-correlation between the S-index and mc(t) shown in the lower panel of Fig. 4: the maximum of the correlation and the minimum are observed for lags different from zero; however the distance between the two peaks is much lower than the cycle period, meaning that a possible time lag is ill-defined.

The uncertainties on each mc(t) value are computed using two assumptions, providing upper and lower limits. The procedure is described in Appendix B. It is likely that the true uncertainties lie between these two estimates, and both are

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Meridional flow vs. S-index (upper panel) for pole-on configuration for cycle 22 (solid line) and cycle 23 (dashed line): ascending phases are shown in red and descending phases in green. The uncertainties in RV are the upper limits. Lower panel: cross-correlation functions between S-index and meridional flows for cycles 22 (solid line) and 23 (dashed line).</figcaption></figure>

shown in Fig. 3. Even when considering the upper uncertainties, we conclude that the long-term variation (cycle) in Fig. 3 is significant, while the variability on smaller timescales may not be. The rms (root-mean-square) of each time series for all input latitudinal profiles is shown in Table 2. The most interesting time series (i.e. from U10) has rms values between 0.46 and 0.71 m s⁻¹ (edge-on). They are higher for the pole-on configuration (1.02–1.59 m s⁻¹). We note that although cycles 22 and 23 are of similar amplitude (from the chromospheric or spot number point of view), the variability of the RV due to meridional flows is larger by about 50% during cycle 23, meaning the relationship between amplitude of the integrated meridional component and the cycle amplitude is complex, because this difference is significant. The amplitude Δ is the amplitude of the sinusoidal fit, and the corresponding peak-to-peak amplitude, 2Δ, is used in most of the paper. We conclude that using other latitudinal profiles for MC gives RV amplitudes of similar magnitude, although they are either slightly larger (BA10) or slightly smaller (HR10, not shown here because very similar to HR11, and HR11).

An important question is whether the variability on timescales smaller than the cycle period is real. From a theoretical point of view, numerical magneto-hydrodynamical simulations
<table><caption>**Table 2.** Solar rms and amplitude of integrated meridional flows.</caption><thead><tr><th>MC time series</th><th>RMS edge-on (m s<sup>-1</sup>)</th><th>RMS pole-on (m s<sup>-1</sup>)</th><th>Δ edge-on (m s<sup>-1</sup>)</th><th>Δ pole-on (m s<sup>-1</sup>)</th></tr></thead><tbody><tr><td><b>U10 cycle 22</b></td><td><b>0.46 (0.01;0.07)</b></td><td><b>1.02 (0.02;0.17)</b></td><td><b>0.50 (0.01;0.21)</b></td><td><b>1.17(0.01;0.39)</b></td></tr><tr><td><b>U10 cycle 23</b></td><td><b>0.71 (0.02;0.08)</b></td><td><b>1.60 (0.02;0.14)</b></td><td><b>0.70 (0.01;0.20)</b></td><td><b>1.67 (0.20;0.38)</b></td></tr><tr><td>BA10 cycle 23</td><td>0.80 (0.05;0.06)</td><td>1.97 (0.10;0.10)</td><td>0.98 (0.07;0.19)</td><td>2.39 (0.14;0.29)</td></tr><tr><td>HR11 cycle 23</td><td>0.59 (0.02;0.02)</td><td>1.11 (0.03;0.04)</td><td>0.64 (0.03;0.07)</td><td>1.15 (0.06;0.11)</td></tr></tbody></table>

**Notes.** The rms is computed on each time series, which have different durations depending on input latitudinal profiles. Δ, also in m s⁻¹, is the amplitude of a sinusoidal fit (and 2Δ the peak-to-peak amplitude). Values between parentheses represent the lower and upper uncertainties (see text). Lines in bold correspond to the main time series studied in this paper.

suggest variability on long timescales only (A. Strugarek, priv. comm.), although short-timescale variability is present due to the activity pattern (see e.g. the two lower panels of Fig. 3). Such variability has been reported in the past (Hathaway 1996) from the latitudinal profiles, and our reconstructions shown in Appendix A also seem to exhibit some year-to-year variability as well as below 1 yr for the reconstructions from HR10 and HR11. However, a first issue is the lack of correspondence between the different reconstructions on those timescales (except between HR10 and HR11 since they correspond to very similar analyses), which casts some doubt on the reality of such variability. A second difficulty arises from the uncertainties on the integrated RV: in Appendix B, we estimate the errors based on two extreme assumptions (correlated and uncorrelated between latitudinal bins), and we estimate that the true uncertainties lie between the two. Unfortunately, the lower bound uncertainties are compatible with a significant short-timescale variability while the upper bound uncertainties are not. In conclusion, we cannot conclude whether the variability at short timescales is significant.

### 2.3.2. Dependence on stellar inclination

The rms over each time series is computed for all inclinations and is shown in Fig. 5 for the two cycles available from U10. The rms is maximum for the pole-on configuration (0°) and is minimal (below 0.3 m s⁻¹) around 45° for cycle 22 and 55° for cycle 23, which also corresponds to the reversal shown in the lower panel. Even with the upper limit for the uncertainties, the differences between the two cycles are significant. Slightly different shapes of the latitudinal profiles are likely to be responsible for the difference in the reversal inclination. The second panel shows the peak-to-peak amplitudes, with similar conclusions about the reversal for inclinations around 50°, the difference between the cycles, and the difference between the pole-on and edge-on configurations. The reversal in sign is illustrated on the lower panel of Fig. 5, which shows the correlation between the meridional time series at different inclinations with the pole-on time series.

### 2.3.3. Effect on the line shape

The bisector span (BIS; see definition and computation method in Sect. 2.2), computed without noise on the spectra in these simulations, varies over time. The results are shown in Fig. 6. The correlation between BIS and RV is close to -1 for pole-on (-0.99) and edge-on (-0.94) configurations but becomes positive for medium inclinations, with a maximum of 0.86 for 50°. The rms of the BIS increases from pole-on to edge-on with

values up to 0.14 m s⁻¹. The peak-to-peak amplitude is small in general, but can reach up to 0.4 m s⁻¹ for edge-on configurations, which may be detectable in some cases. This means that MC will not significantly affect this observable in many cases and the BIS cannot be used to correct the RV signal for the meridional flows for low and intermediate inclinations. This is probably due to the fact that, unlike spot or plage crossing across the disc, meridional flow maps are symmetric with respect to the central meridian (i.e. the same at all longitudes). However, there may be a possibility to use it for close-to edge-on configurations if the precision is sufficient (amplitude greater than 0.3 m s⁻¹ for inclinations higher than 70°, representing about one-third of the RV amplitude, although one should keep in mind that it will be superposed to other contributions to BIS variations when the star is active.

### 2.3.4. Comparison with Earth-like planetary signals

We now compare the RV signal due to meridional flows mc(t) with the signal due to an Earth-mass planet in the habitable zone. We consider orbital periods of 274 days (inner side), 505 days (middle of the habitable zone), and 777 days (external side) as in Meunier & Lagrange (2019b). For that purpose, we first interpolate the mc(t) time series over a new sampling (the interpolation is performed using a spline function) consisting of 1000 points and covering 12.5 yr between 1996 and 2009 (i.e. cycle 23), with a gap of 4 months every year to simulate the fact that a star can in general not be observed at all times.

Figure 7 (upper panel) shows the time series derived from U10 (edge-on configuration) and for planets with masses of 1<sub>M</sub>arth and 2<sub>M</sub>arth in the middle of the habitable zone: the planetary signals are much smaller than the meridional flow signal. The periodograms of the different types of signal are then compared for the two masses and the three orbital periods in the habitable zone, which also shows that such planetary signals are much smaller than the MC signal in this range of periods.

A correction of the stellar signal using a direct linear relationship between RV and the S-index is not useful here because of the strong time lag between the two, as shown in Sect. 2.3.1. A correction using a smoothed version of the S-index for cycle 23 (the smoothing is performed using a simple moving average) and taking the lag into account for the 14 points of cycle 23 lead to a rms of 0.58 m s⁻¹ for the residuals instead of 0.71 m s⁻¹, which is a small reduction. Here we use a simple sinusoidal fit, which decreases the rms for the same time series to 0.45 m s⁻¹. On the interpolated series shown in Fig. 7, the rms decreases from 0.57 to 0.26 m s⁻¹. This sinusoidal fit of mc(t) significantly decreases the power due to meridional flows, as shown in the
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Root mean square of meridional flow time series (*upper panel*), peak-to-peak amplitude (*middle panel*), and correlation of meridional time series with pole-on time series (*lower panel*) vs. stellar inclination, for U10 reconstruction: full time series (black), cycle 22 (orange), and cycle 23 (green). The uncertainties with the same colour correspond to the upper limits, and those in red to the lower limit. The vertical dotted lines indicate the position of the reversals for the two cycles.</figcaption></figure>

lower panel, but the residuals still exhibit significant power for periods in the habitable zone. The exact amplitude of the residuals should however depend on the exact variability on timescales shorter than the cycle period (10s and 100s of days), as discussed in Sect. 2.3.1. The other time series (Appendices A.2 and A.3), although providing slightly different periodograms, lead to similar conclusions. Furthermore, when superposed on other strong signals, a simple sinusoidal fit will not be the best model. More sophisticated techniques will have to be implemented to correct for such flows.

This example shows the impact of the edge-on configuration. We conclude that in this case, even after the correction

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Peak-to-peak amplitude (solid line) and rms of BIS (dashed line) vs. inclination (*upper panel*), and correlation between BIS and mc(t) (*lower panel*), to be compared with rms variation in Fig. 5.</figcaption></figure>

described above, the detectability of Earth-like planets in the habitable zone will be challenging, although the performance may be better for such planets on the inner side of the habitable zone. However, specific configurations at inclinations close to the reversal will be more suitable. The impact may also depend on the inclination of the orbital plane of the planet with respect to the equatorial plane, as a 2 $M_{\text{Earth}}$ with an inclination of 30° for example corresponds to a 1 $M_{\text{Earth}}$ planet seen edge-on. More results taking into account this impact are shown in Sect. 3.

### 2.3.5. Superposition to the magnetic activity signal

Finally, we combine mc(t) to the contribution due to solar magnetic activity, which is dominated by the inhibition of the convective blueshift, but also includes the signal due to spots and plages crossing the disc (through their intensity contrast). We first superpose mc(t) to the magnetic activity RV derived from the solar model of Borgniet et al. (2015), which statistically follows cycle 23 general behavior<fnref n="1" />: this allows us to consider configurations different from edge-on. Figure 8 shows the activity signal alone (in black) and after superposition (in red) with

<footnote n="1">¹ The exact magnetic features are different because these RVs result from a model where spots and plages are generated, but they follow cycle 23 amplitude and cycle shape.</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Radial velocity due to edge-on meridional flows vs. time (stars and solid line) and due to 1 (red) and 2 (green) Earth-mass planets in the middle of the habitable zone (first panel). The two following panels show the periodograms for the meridional flow (black; full range, then focusing on small powers). The planet periodograms are shown for different orbital periods: inner side of the habitable zone (solid), middle (dashed), and external side (dotted), and the planet peaks are indicated by arrows. Last panel: similar but for the residuals after the sinusoidal correction (see text).</figcaption></figure>

the meridional flow signal for inclinations between 10 and 90°. Except for inclinations corresponding to the reversal, the shape is significantly altered, especially for low inclinations, because the meridional flow component, although related to the cycle, is shifted in time and is of strong amplitude. The corresponding periodograms for periods longer than 100 days are also shown. The last plot (lower right panel) of Fig. 8 shows the correlation between the two signals, that is, before and after superposition

of the meridional flow signal, versus inclination. The departure from a correlation of one is low when seen edge-on, while very large departures are present for inclinations below 40°. There is no departure from a correlation of one around 55°, where the reversal is occurring.

We now look at the edge-on configuration in more detail. The departure from the original correlation between RV and the S-index is not very large, but it is present. We consider the solar magnetic activity signal reconstructed for cycle 23 in Meunier et al. (2010a) and compare it to the Sacramento Peak Observatory S-index (edge-on case only). The correlation between RV and chromospheric emission is computed over 3586 days covering cycle 23, and is equal to 0.91 without the meridional flow (i.e. for magnetic activity alone), and 0.88 after addition of the meridional flow, hence a small decrease. When binning the time series over 6 months, the correlation (sensitive mostly to long timescales) decreases from 0.99 to 0.94. The impact of such departure on the RV signal is illustrated in Fig. 9, showing the binned RV versus the binned chromospheric emission. As shown in Meunier et al. (2019a), even with no MC, this curve presents an hysteresis along the cycle due to projection effects and the butterfly diagram. When taking meridional flows into account (dashed curve), the effect is amplified in this particular configuration (edge-on). The correction method proposed in Meunier et al. (2019a), which was already taking a departure from a linear correlation between RV and chromospheric emission into account, will therefore have to be adapted and include a time lag introduced by the presence of meridional flows: an additional component of the form $\alpha \times \log R'_{HKsm}(t+\text{lag})$, where $\log R'_{HKsm}$ is a smooth version of $\log R'_{HK}$ (to remove the rotational contribution), could be added to their model.

We have so far compared the MC contribution with the convective blueshift inhibition, which is dominating the long-term variability in the solar case, with a typical amplitude over the cycle of the order of 8 m s⁻¹. We now compare the MC amplitude with other sources of stellar signal which have been identified:

- Spot and plage contrasts: their contribution mostly affects periods close to the rotation period or its harmonics (e.g. Saar & Donahue 1997), with an rms of the order of 0.33 m s⁻¹ for the Sun seen edge-on on average over the solar cycle (Lagrange et al. 2010; Meunier et al. 2010a). Lower values are obtained at lower inclinations due to projection effects (Borgniet et al. 2015). Their contribution at long periods is therefore much smaller than the MC contribution.

- Oscillations: the five-minute oscillation signal due to p-modes has a negligible contribution compared to the other stellar signals at long periods and can be easily averaged out (e.g. Chaplin et al. 2019).

- Granulation: the rms of this contribution is estimated to be between 0.4 and 0.8 m s⁻¹ (Elsworth et al. 1994; Pallé et al. 1999; Meunier et al. 2015). Figure 10 (upper panel) shows a comparison of the periodogram due to granulation after a one-hour binning (the original rms is 0.8 m s⁻¹, and 0.39 m s⁻¹ after binning) for time series similar to those shown in Fig. 7 (MC and planets); the MC contribution is higher than the granulation contribution in the range of periods we are interested in.

- Supergranulation: the rms of this contribution is less well determined, but could be between ~0.3 and 1.2 m s⁻¹ (Meunier et al. 2015); Pallé et al. (1999) also found a rms of 0.78 m s⁻¹. Figure 10 (lower panel) shows a comparison of the periodogram due to supergranulation (rms of 0.7 m s⁻¹, also binned over one hour, which in this case does not significantly change the rms) for time series similar to those shown
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Radial velocity due to activity (black dots) and activity superposed on meridional flows (red dots) vs. time for various inclinations, and corresponding periodograms. Meridional flows alone (reconstruction from U10) are shown as a solid green line after interpolation as in Fig. 7. Offsets are arbitrary (the lowest value for each series has been arbitrarily set to zero). The *last plot* shows the correlation between the two time series (without and with meridional flows) versus inclination.</figcaption></figure>

in Fig. 7 (MC and planets). In this case, although the supergranulation contribution is slightly lower than that from MC, it can lead to peaks of similar amplitudes in certain cases in the range of periods corresponding to the habitable zone.

### 3. From F to K stars

In this section, we estimate the amplitude of the integrated meridional flows for a range of stars similar to the magnetic activity simulations in Meunier et al. (2019b), i.e. covering the spectral range F6-K4 and activity levels corresponding to old

main sequence stars. The range of parameters is shown in Fig. 11. Stellar RVs are expected to depend on spectral type and rotation rate (i.e. in our case on activity level). Since the variability of the solar integrated MC obtained in Sect. 2 is related to the cycle, stellar MC amplitudes should also depend on the cycle amplitude. To our knowledge, there has so far not been any observation of these flows for stars other than the Sun because they are very weak. However, a number of 3D hydrodynamical simulations have been performed in several studies, which allows us to scale the stellar RV to the solar one. We then study the effect of other properties, such as for example the likely presence of
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Radial velocity vs. Sacramento Peak Observatory chromospheric emission index, both binned over 6 months, for activity alone (solid lines) and for activity superposed on meridional flows (dashed lines). Red curves correspond to the ascending phase of cycle 23, and green curves to the descending phase. Curves are arbitrarily shifted in RV.</figcaption></figure>

multi-cellular patterns for certain stars, and the impact of the limb darkening function. The MC amplitude is compared to the amplitude of Earth-like planets in the habitable zone and to the RV amplitude due to the inhibition of the convective blueshift (the dominating part due to magnetic activity over the cycle) to quantify the importance of the MC contribution.

## 3.1. Scaling of the solar amplitudes with spectral type and Prot

In this section, we first consider the effect of rotation rate and mass on the global meridional flow amplitude. Ballot et al. (2007) simulations (made for $1 M_\odot$ only) showed that the amplitude of the meridional flows decreases for smaller rotation periods (i.e. larger values of $\Omega$). For a given mass, stars are strongly affected by stellar activity, which is proportional to the rotation rate: the quietest stars have the longest rotation periods and therefore we expect them to exhibit the highest meridional circulation. This means that in our grid of parameters, for a given mass, the quietest stars should show the greatest meridional circulation. Ballot et al. (2007) modelled its amplitude with a scaling law in $\Omega^\beta$, where the exponent $\beta$ took values between -0.53 and -1.02 depending on the set of simulations. In addition, Matt et al. (2011) obtained larger meridional flows for high-mass stars for a solar rotation rate. Both studies correspond to a limited range of parameters, and these works were subsequently extended to a larger range of parameters in Brun et al. (2017), who obtained similar trends. The range of parameters covered in Brun et al. (2017) is shown in Fig. 11 (middle panel, orange and green dots). Most of the parameters we consider correspond to solar differential rotation, as discussed in Meunier et al. (2019b), which is also apparent in this figure (which indicates which hydrodynamical simulations correspond to solar or anti-solar differential rotation). We model the meridional flow amplitude obtained by Brun et al. (2017) as proportional to $\Omega^\beta M^\alpha$ from the green dots (i.e. solar differential rotation only). Adjusting the power law in $\Omega$ for each of the four ranges in mass leads to $\beta$ varying between -0.19 and -0.36, with an average of -0.26. The dependence on the mass can then be modelled with $\alpha = 3.44$. The results are not very different when considering all simulations (both solar and anti-solar differential rotation).

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

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Periodograms of time series shown in the third panel of Fig. 7 (original time series and planetary RV for 1 and 2 $M_\odot$), added to periodogram due to granulation (upper panel) and supergranulation (lower panel). The colour code for MC and planets is similar to that of Fig. 7, and the granulation and supergranulation periodograms are shown in brown.</figcaption></figure>

We apply these scaling laws (with $\alpha = 3.44$, $\beta = -0.26$) to the solar peak-to-peak amplitudes (the values corresponding to rms RV are smaller by about a factor two; see Table 2) of mc(t) derived in Sect. 2 for the U10 solar profiles (cycle 23) and four inclinations: 0° (pole-on), 30°, 55° (where the variability is the lowest), and 90° (edge-on). The results corresponding to a solar cycle amplitude are shown in the upper panels of Fig. 12. The dependence is dominated by the mass, while the effect of the rotation is weaker, within the range of parameters we consider.

We then scale the MC amplitude with the amplitude of the cycle in our grid of simulations as described in Meunier et al. (2019b), following the prescription shown in Fig. 11 (lower panel) and assuming the amplitude of the meridional flow variability is linearly correlated with the cycle amplitude. These results are presented in the lower panels of Fig. 12, again for four inclinations. In this case, the meridional flow amplitudes are dominated by the cycle amplitude, that is, by the magnetic activity variability. Fine details in the scaling laws describing how the meridional flows vary with mass and rotation rate are therefore not as critical here. More work on how the meridional flows
<figure><img src="image_13.png" /><figcaption>**Fig. 11.** Average log $R'_{HK}$ vs. B-V (upper panel), rotation rate Ω vs. stellar mass (middle panel), and cycle amplitude vs. B-V (lower panel) covered in this work, all extracted from the grid of parameters in Meunier et al. (2019b). The cycle amplitude is the difference between the prescribed maximum and minimum log $R'_{HK}$ over the cycle. The green and orange dots in the middle panels represent the position the MHD simulation parameters in Brun et al. (2017) for solar differential rotation (green) and anti-solar differential rotation (orange).</figcaption></figure>

depend on the magnetic activity variability will be needed in the future however, since we have considered a linear relationship here.

As a consequence, there is no strong trend of the expected MC amplitude versus mass once it is scaled by the activity amplitude, which shows a very large dispersion. This is because each given spectral type has a large range of possible activity levels. For the edge-on configuration, the expected peak-to-peak amplitude takes values up to 1.8 m s⁻¹ (edge-on), with many values in the 0.2–1 m s⁻¹ range. For the pole-on configuration, the amplitude can be as high as 4 m s⁻¹. Considering stellar inclination

where the integrated RV is reversing would lead to much lower amplitudes (lower than 0.5 m s⁻¹). The distributions of peak-to-peak amplitudes for different inclinations is shown in Fig. 13. The maximum of the global distribution is around 0.3 m s⁻¹. For edge-on configurations, the maximum of the distribution is in the range 0.3–0.8 m s⁻¹.

In this section, we directly scale the solar flow (obtained with a limb darkening function corresponding to the solar $T_{\text{eff}}$) using the values obtained in Sect. 2. We check what would be the impact of a limb-darkening function computed for the proper temperature over the range of parameters on the integrated meridional flow. Trends are observed with stellar type as expected, but we find that the effect of the limb-darkening function can be neglected, with a maximum departure of the order of 0.8% (for the edge-on configuration, and B-V around 1.05). We therefore neglect this effect here.

For fast rotating stars (typically faster than the Sun), hydro-dynamical simulations show that the latitudinal profiles could be more complex than the solar one, and exhibit a multi-cellular pattern (Matt et al. 2011; Guerrero et al. 2013, 2016). We investigated the impact of such a pattern and found that this tends to reduce the RV effects induced by MC due to cancellation between different latitudes. However, such fast rotating stars are not very good targets to search for low-mass planets due to their faster rotation than the Sun, because of the lower achievable RV precision.

## 3.2. Comparison with Earth-like planetary signals

Figure 12 also shows the peak-to-peak amplitude of planetary signals for 1 $M_{\text{Earth}}$ and 2 $M_{\text{Earth}}$ planets orbiting in the habitable zone (three values of the orbital periods). This amplitude can be compared with the MC contribution without noise. The amplitude shown here is projected according to the assumption that the orbital plane is similar to the stellar equatorial plane. As a consequence, the amplitude is zero for pole-on configurations (which are not very suitable for detecting planets). However, there could be departures from these curves if there is an inclination between the two planes: in this case, the maximum planetary amplitude is the one shown in the panel corresponding to the edge-on configuration. In most cases, the MC amplitude is larger than an Earth-mass planet in the habitable zone, except for a small number of simulations corresponding to very quiet stars seen with an intermediate inclination. The percentage of simulations above those planetary levels is shown in Table 3 for three values of habitable zone periods (PHZ) and two planet masses, again assuming that the orbital plane and the equatorial plane are similar. These percentages are lowest only for 2 $M_{\text{Earth}}$ and intermediate inclinations. All percentages are shown in Fig. 14; for 1 $M_{\text{Earth}}$, they are never lower than 70% and are their lowest for inclinations around 50° (red curves) and 60° (brown curves). Percentages are lower for 2 $M_{\text{Earth}}$ and can be close to zero for these intermediate inclinations, although the percentage is always above 60% for other configurations.

## 3.3. Comparison to the RV amplitude due to magnetic activity

A detailed analysis of the effect of MC on exoplanet detectability is beyond the scope of this paper because it requires the inclusion of magnetic activity in order to be realistic given the distortion of the variability introduced by this latter (see Sect. 2.3.4). However, it is interesting to estimate how the amplitudes we have found compare with the other processes leading to significant long-term variability, such as the inhibition of the convective
<figure><img src="image_14.png" /><figcaption>**Fig. 12.** Stellar RV peak-to-peak amplitude (scaled to solar meridional flows from U10, cycle 23) vs. *B*–*V* for solar cycle 23 amplitude (*upper panels*) and scaled to cycle amplitude from our grid of parameters (*lower panels*) for four inclinations. Orange symbols represent rotation rates higher than the solar one. The black lines on the left of each panel indicate where the mass is outside the range of hydrodynamical simulation parameters. The peak-to-peak amplitude of the planetary signal, assuming the planet orbital plane is similar to the stellar equatorial plane, is shown for 1 *M*<sub>Earth</sub> (solid lines) and 2 *M*<sub>Earth</sub> (dashed lines) and for different orbital periods in the habitable zone: PHZ<sub>in</sub> (black), PHZ<sub>med</sub> (red), and PHZ<sub>out</sub> (green). The planetary curves for a pole-on configuration therefore follow the *x*-axis.</figcaption></figure>

blueshift obtained in Meunier et al. (2019b), which we know
is very challenging for exoplanet detections. The comparison
between peak-to-peak MC amplitudes and the long-term ampli-
tude of the convective blueshift inhibition obtained in Meunier
et al. (2019b) is shown in Fig. 15 for all inclinations (upper
panel). We find that in general the convective blueshift contri-
bution is larger, except for a few configurations for inclinations
lower than 30°. Only 1.3% of the simulations correspond to MC
amplitude larger than the convective blueshift inhibition ampli-
tude. We compute the proportion of simulations for which the
MC amplitude represents at least a certain fraction (5, 10, and
50%) of the convective blueshift amplitude (the straight lines in
the first panel indicate those proportions). These proportions are
8.9, 53.4, and 80.8% for the 50, 10, and 5% levels. The median
ratio over all realisations is 9.2, with ratios for individual reali-
sations between 0.5 and 208. We also computed the median ratio
separately for each inclination, as shown in the second panel of
Fig. 15: the median ratio at a given inclination is always higher
than 1 (although not very far from it for low inclinations), and is
significantly higher for inclinations around 40–50° (median ratio

higher than 30): those inclinations are the most suitable config-
urations for exoplanet detectability, from the point of view of
meridional flows.

The percentage of simulations where MC represents more
than a certain fraction of the convective blueshift inhibition is
shown in the following panels versus different variables: incli-
nation, *B*−*V*, log *R*'<sub>HK</sub>, and *A*<sub>cyc</sub>. The 50% level is usually low,
with a predominance of low inclinations, and to a lesser extent of
quiet, low-mass stars: the percentage of simulations where MC
represents more than half of the convective blueshift inhibition
can be as high as 40% for an inclination of 0°, showing again
the possibility for strong distortions of the RV-chromospheric
emission relationship. A large proportion of simulations have
MC representing more than 10% of the convective blueshift inhi-
bition amplitude (except around 50–60°), which should lead to
significant distortions of the relationship between RV and activ-
ity indicators: this will therefore affect corrections and exoplanet
detectability. The most suitable configurations are again medium
inclinations, followed by high mass, more active stars (with
log *R*'<sub>HK</sub> above approximately −4.7), at least from the point of
<figure><img src="image_15.png" /><figcaption>**Fig. 13.** Distributions of stellar MC values for all simulations (black) and vs. inclinations (from pole-on in yellow to edge-on in blue, as shown by the colour bar).</figcaption></figure>

<table><caption>**Table 3.** Percentage of simulations with MC larger than planetary signal.</caption><thead><tr><th rowspan="2">PHZ</th><th colspan="3">1 M<sub>Earth</sub></th><th colspan="3">2 M<sub>Earth</sub></th></tr><tr><th>0°</th><th>55°</th><th>90°</th><th>0°</th><th>55°</th><th>90°</th></tr></thead><tbody><tr><td>PHZ<sub>in</sub></td><td>100</td><td>45.9</td><td>87.7</td><td>100</td><td>3.2</td><td>69.5</td></tr><tr><td>PHZ<sub>med</sub></td><td>100</td><td>61.2</td><td>90.9</td><td>100</td><td>10.2</td><td>76.1</td></tr><tr><td>PHZ<sub>out</sub></td><td>100</td><td>67.0</td><td>92.4</td><td>100</td><td>19.1</td><td>80.1</td></tr></tbody></table>

view of the relative contribution of meridional flows compared to the inhibition of the convective blueshift.

Finally, we estimate the distortion of the RV-log $R'_{HK}$ relationship brought about by the presence of MC. For each simulation, we estimate the meridional circulation by smoothing the log $R'_{HK}$ time series, shifting it with a time lag similar to the solar one (0.23 times the cycle period), and scaling its amplitude as in Sect. 3.1. The MC signal is then added to the magnetic activity contribution of Meunier et al. (2019b) and the correlation between RV (before adding MC and after) and log $R'_{HK}$ is computed. The results are shown in Fig. 16. As for the previous comparisons, the main effects are seen for low inclinations (the median ratio, where the ratio is between the correlation with activity+MC and the correlation with activity alone, is down to 0.9 for pole-on stars), especially for quiet (on average), low-mass stars, mostly when exhibiting some variability.

## 4. Conclusion

We reconstructed the integrated solar RV due to meridional flows with very good temporal coverage (two cycles, all inclinations). The MC variability is clearly related to the cycle, although cycles of similar amplitude (cycles 22 and 23) lead to different amplitudes: the difference is of the order of 50%, which is most likely the main source of difficulty that we face when attempting to predict expected amplitudes. The most interesting time series for exoplanet detectability is the reconstruction made from latitudinal profiles obtained from Dopplergram analysis (Ulrich 2010) because they correspond to flows similar to

<figure><img src="image_16.png" /><figcaption>**Fig. 14.** Percentage of simulations for which MC is higher than planetary signal vs. B-V. Each panel corresponds to a different mass and orbital period. The colour code is similar to Fig. 13: inclination from pole-on in yellow to edge-on in blue (see colour bar on Fig. 13). Red corresponds to 50°.</figcaption></figure>

those seen when measuring RV in the visible. The integrated RV exhibits a time lag with respect to the cycle. The other reconstructions (using results obtained with ring diagram helioseismology, and magnetic feature tracking) are more closely anti-correlated with the cycle (maximum at activity minimum and vice versa), and amplitudes are robust. The integrated RV also varies on shorter timescales, but the exact behaviour of the integrated RV due to meridional flows at these scales is still unclear.

Typical values for the rms of the signal (over cycle 23) are of the order of 0.7 m s⁻¹ for the edge-on configuration and 1.6 m s⁻¹ for the pole-on configuration (cycle 23), with a minimum around 0.3 m s⁻¹ for an inclination of 55° (where there is a reversal). Peak-to-peak amplitudes are respectively of 1.4 and 3.3 m s⁻¹ for the edge-on and pole-on configurations. The exact position of the reversal may vary depending on the exact shape of the latitudinal profile over time (45° for cycle 22).
<figure><img src="image_17.png" /></figure>

**Fig. 15.** Meridional circulation amplitude vs. long-term amplitude of convective blueshift inhibition (upper panel) on our grid of parameters (from Meunier et al. 2019b) for different inclinations (from pole-on in yellow to edge-on in blue). Only one point in five is shown for clarity. The straight lines correspond to MC amplitude representing 5% (dotted-dashed line), 10% (solid line), 50% (dashed line), and 100% (dotted line) of the convective blueshift inhibition. Second panel: median ratio between the convective blueshift inhibition component and MC vs. inclination. The four following panels show the percentage of simulations in which MC represents more than 5% (resp. 10, 50%) (in dotted-dashed lines, resp. solid lines, dashed lines) of the convective blueshift amplitude vs. inclination, *B*–*V*, log *R*'HK, and *A*<sub>cyc</sub>.

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

**Fig. 16.** Correlation between activity+MC signal and log *R*'HK vs. correlation with activity alone. The colour code represents inclination as in Fig. 15. The four following panels show the median ratio between correlation for activity+MC divided by correlation for activity alone vs. inclination, *B*–*V*, log *R*'HK, and *A*<sub>cyc</sub> (red lines). The green lines represent the 5th and 95th percentiles, and the dotted lines the 10th and 90th percentiles.

We compared these solar time series with the signal that would be due to an Earth-like planet in the habitable zone around the Sun: the power due to meridional flows is much larger than the planetary signal. Even after a sinusoidal fit to correct for the meridional flow contribution, the power at periods corresponding to the habitable zone is still large and above the level corresponding to Earth-like planets. Furthermore, when added to the magnetic activity signal, the correlation existing between
RV due to magnetic structures (in particular the inhibition of the convective blueshift) and chromospheric emission is degraded. The degradation is moderate for edge-on configurations (but must be taken into account to reach very low detection limits). We find that the hysteresis pattern observed between these two observables (Meunier et al. 2019a) is amplified by the presence of meridional circulation due to the time lag (around 2.5 yr) between RV due to meridional flows and the cycle. The effect is naturally very small when reaching the reversal around 55°, but becomes major at lower inclinations, leading to RV variations that show departure from a good correlation with the chromospheric emission variations. This alteration of the RV-chromospheric emission relationship should be a major difficulty to correct for this signal as well as to correct for magnetic activity in such conditions, especially for configurations close to pole-on. Nevertheless, we note that the pole-on configuration is not the most suitable for exoplanet detection in RV, because the signal of the planet falls to zero if the orbital plane is similar to the stellar equatorial plane. The bisector variations may be significant (up to 0.4–0.45 m s⁻¹) for configurations close to edge-on, which could be used to remove part of the MC contribution in those cases.

Finally, using scaling laws from the numerical simulations of Brun et al. (2017), we extrapolated the amplitude of the integrated meridional flows that we find for the Sun to other solar-type stars. For an amplitude of the activity cycle similar to the solar one (i.e. considering only a scaling in mass and rotation), the variability is dominated by the mass, with a much larger RV variability for F stars compared to K stars. However, if we also scale this variability with the amplitude of the cycle expected for different types of stars while assuming a linear relationship valid for all stars then the variability is in fact dominated by the cycle amplitude. In view of the differences of up to 50% that we see for the Sun, it will be very important to better understand this relationship in the future. For stars seen edge-on, the peak-to-peak amplitude can reach 1.8 m s⁻¹, and up to 4 m s⁻¹ for pole-on stars. At stellar inclinations close to the reversal, the amplitudes are much smaller. Furthermore, some stars may exhibit multi-cellular patterns, as shown from numerical simulations: this should significantly decrease the RV amplitude, and so these configurations may be significantly more suitable. We expect this effect to play a role for F stars and early G stars, but not for late G and K stars, in the range of rotation periods we consider here.

The present study focuses on characterising the amplitude of the MC contribution. We find that its variability is of the same order of magnitude as or slightly larger than other processes (granulation, supergranulation, spots, plages), i.e. typically in the range 0.5–1 m s⁻¹ for a star with a solar-like cycle, and can be much higher for stars with low inclination. It also represents a significant proportion of the amplitude of the signal due to the inhibition of the convective blueshift, which will have a significant impact on exoplanet detectability, as shown in Sects. 2.3.5 and 3.3. However, we have not estimated the performance in terms of detection rates and false-positive levels because this must be done in conjunction with the presence of magnetic activity in order to obtain meaningful results, which is beyond the scope of this paper: this will be done in a future study, following a systematic approach combining MC and magnetic activity contributions using blind tests.

*Acknowledgements.* We thank S. Brun and A. Strugarek for usefull discussions about this work. This work has been funded by the ANR GIPSE ANR-14-CE33-0018. This work was supported by the "Programme National de Physique

Stellaire" (PNPS) of CNRS/INSU co-funded by CEA and CNES. This work was supported by the Programme National de Planétologie (PNP) of CNRS/INSU, co-funded by CNES. The monthly sunspot number has been provided by the SIDC-team, World Data Center for the Sunspot Index, Royal Observatory of Belgium (http://www sidc.be/sunspot-data/). This work made use of the Ca K index provided by the Sacramento Peak Observatory of the US Air Force Phillips Laboratory.

## References

Ballot, J., Brun, A. S., & Turck-Chièze, S. 2007, ApJ, 669, 1190

Basu, S., & Antia, H. M. 2010, ApJ, 717, 488

Boisse, I., Bonfils, X., & Santos, N. C. 2012, A&A, 545, A109

Borgniet, S., Meunier, N., & Lagrange, A.-M. 2015, A&A, 581, A133

Brun, A. S., Strugarek, A., Varela, J., et al. 2017, ApJ, 836, 192

Cegla, H. M., Shelyag, S., Watson, C. A., & Mathioudakis, M. 2013, ApJ, 763, 95

Cegla, H. M., Watson, C. A., Shelyag, S., et al. 2018, ApJ, 866, 55

Cegla, H. M., Watson, C. A., Shelyag, S., Mathioudakis, M., & Moutari, S. 2019, ApJ, 879, 55

Chaplin, W. J., Cegla, H. M., Watson, C. A., Davies, G. R., & Ball, W. H. 2019, AJ, 157, 163

Chou, D.-Y., & Dai, D.-C. 2001, ApJ, 559, L175

Claret, A., & Hauschildt, P. H. 2003, A&A, 412, 241

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

Dumusque, X. 2016, A&A, 593, A5

Dumusque, X., Lovis, C., Ségransan, D., et al. 2011, A&A, 535, A55

Dumusque, X., Boisse, I., & Santos, N. C. 2014, ApJ, 796, 132

Duvall, T. L., J. 1979, Sol. Phys., 63, 3

Elsworth, Y., Howe, R., Isaac, G. R., et al. 1994, MNRAS, 269, 529

Galland, F., Lagrange, A.-M., Udry, S., et al. 2005, A&A, 443, 337

Gizon, L., Birch, A. C., & Spruit, H. C. 2010, ARA&A, 48, 289

Guerrero, G., Smolarkiewicz, P. K., Kosovichev, A. G., & Mansour, N. N. 2013, ApJ, 779, 176

Guerrero, G., Smolarkiewicz, P. K., de Gouveia Dal Pino, E. M., Kosovichev, A. G., & Mansour, N. N. 2016, ApJ, 819, 104

Hathaway, D. H. 1996, ApJ, 460, 1027

Hathaway, D. H. 2012, ApJ, 760, 84

Hathaway, D. H., & Rightmire, L. 2010, Science, 327, 1350

Hathaway, D. H., & Rightmire, L. 2011, ApJ, 729, 80

Hatzes, A. P. 2002, Astron. Nachr., 323, 392

Howard, R., & Gilman, P. A. 1986, ApJ, 307, 389

Komm, R. W., Howard, R. F., & Harvey, J. W. 1993, Sol. Phys., 147, 207

Kurucz, R. L. 1993, CD-ROM, 13, 18 http://kurucz.harvard.edu

Labonte, B. J., & Howard, R. 1982, Sol. Phys., 80, 361

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

Lin, C.-H., & Chou, D.-Y. 2018, ApJ, 860, 48

Makarov, V. V. 2010, ApJ, 715, 500

Matt, S. P., Do Cao, O., Brown, B. P., & Brun, A. S. 2011, Astron. Nachr., 332, 897

Meunier, N. 1999, ApJ, 527, 967

Meunier, N., & Lagrange, A. M. 2019a, A&A, 625, L6

Meunier, N., & Lagrange, A. M. 2019b, A&A, 628, A125

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

Meunier, N., Lagrange, A.-M., & Cuzacq, S. 2019a, A&A, 632, A81

Meunier, N., Lagrange, A. M., Boulet, T., & Borgniet, S. 2019b, A&A, 627, A56

Nesme-Ribes, E., Meunier, N., & Vince, I. 1997, A&A, 321, 323

Pallé, P. L., Roca Cortés, T., Jiménez, A., GOLF Team, & Virgo Team. 1999, ASP Conf. Ser. 173, 297

Saar, S. H., & Donahue, R. A. 1997, ApJ, 485, 319

Saar, S. H., Hatzes, A., Cochran, W., & Paulson, D. 2003, Cool Stars, Stellar Systems, and the Sun, Vol. 12, The Future of Cool-Star Astrophysics: 12th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, eds. A. Brown, G. M. Harper, & T. R. Ayres (Berlin: Springer), 694

Snodgrass, H. B., & Dailey, S. B. 1996, Sol. Phys., 163, 21

Ulrich, R. K. 2010, ApJ, 725, 658

Ulrich, R. K., Boyden, J. E., Webster, L., et al. 1988, Sol. Phys., 117, 291

Wright, J. T. 2005, PASP, 117, 657

Zhao, J., & Kosovichev, A. G. 2004, ApJ, 603, 776
# Appendix A: Integrated meridional-flow time series

<table><caption>**Table A.1.** Integrated meridional flow from Ulrich (2010).</caption><thead><tr><th>Time (yr)</th><th>RV at 0° (m s<sup>-1</sup>)</th><th>σ<sub>high</sub> at 0° (m s<sup>-1</sup>)</th><th>σ<sub>low</sub> at 0° (m s<sup>-1</sup>)</th><th>RV at 55° (m s<sup>-1</sup>)</th><th>σ<sub>high</sub> at 55° (m s<sup>-1</sup>)</th><th>σ<sub>low</sub> at 55° (m s<sup>-1</sup>)</th><th>RV at 90° (m s<sup>-1</sup>)</th><th>σ<sub>high</sub> at 90° (m s<sup>-1</sup>)</th><th>σ<sub>low</sub> at 90° (m s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>1986</td><td>-6.45</td><td>0.51</td><td>0.05</td><td>0.12</td><td>0.02</td><td>0.05</td><td>3.88</td><td>0.25</td><td>0.02</td></tr><tr><td>1987</td><td>-6.05</td><td>0.49</td><td>0.06</td><td>0.07</td><td>0.02</td><td>0.05</td><td>3.51</td><td>0.26</td><td>0.03</td></tr><tr><td>1988</td><td>-5.89</td><td>0.48</td><td>0.05</td><td>-0.52</td><td>0.01</td><td>0.04</td><td>3.38</td><td>0.25</td><td>0.02</td></tr><tr><td>1989</td><td>-6.38</td><td>0.56</td><td>0.07</td><td>-0.20</td><td>0.02</td><td>0.06</td><td>3.72</td><td>0.30</td><td>0.02</td></tr><tr><td>1990</td><td>-5.75</td><td>0.69</td><td>0.08</td><td>-0.09</td><td>0.04</td><td>0.06</td><td>3.39</td><td>0.35</td><td>0.03</td></tr><tr><td>1991</td><td>-8.09</td><td>0.58</td><td>0.07</td><td>0.32</td><td>0.04</td><td>0.05</td><td>4.50</td><td>0.31</td><td>0.02</td></tr><tr><td>1992</td><td>-8.51</td><td>0.53</td><td>0.05</td><td>0.50</td><td>0.04</td><td>0.05</td><td>4.70</td><td>0.27</td><td>0.02</td></tr><tr><td>1993</td><td>-7.66</td><td>0.52</td><td>0.07</td><td>0.32</td><td>0.01</td><td>0.04</td><td>4.22</td><td>0.28</td><td>0.02</td></tr><tr><td>1994</td><td>-7.94</td><td>0.49</td><td>0.05</td><td>-0.06</td><td>0.01</td><td>0.03</td><td>4.29</td><td>0.25</td><td>0.01</td></tr><tr><td>1995</td><td>-8.09</td><td>0.50</td><td>0.04</td><td>-0.33</td><td>0.01</td><td>0.03</td><td>4.34</td><td>0.26</td><td>0.02</td></tr><tr><td>1996</td><td>-7.24</td><td>0.69</td><td>0.09</td><td>-0.06</td><td>0.01</td><td>0.05</td><td>3.96</td><td>0.32</td><td>0.02</td></tr><tr><td>1997</td><td>-6.67</td><td>0.51</td><td>0.06</td><td>-0.19</td><td>0.01</td><td>0.04</td><td>3.78</td><td>0.25</td><td>0.02</td></tr><tr><td>1998</td><td>-6.78</td><td>0.59</td><td>0.06</td><td>-0.10</td><td>0.02</td><td>0.04</td><td>3.69</td><td>0.31</td><td>0.03</td></tr><tr><td>1999</td><td>-8.24</td><td>0.54</td><td>0.06</td><td>0.08</td><td>0.02</td><td>0.04</td><td>4.41</td><td>0.27</td><td>0.03</td></tr><tr><td>2000</td><td>-9.47</td><td>0.52</td><td>0.07</td><td>-0.09</td><td>0.01</td><td>0.05</td><td>5.06</td><td>0.27</td><td>0.02</td></tr><tr><td>2001</td><td>-8.84</td><td>0.53</td><td>0.08</td><td>0.15</td><td>0.02</td><td>0.05</td><td>4.59</td><td>0.28</td><td>0.02</td></tr><tr><td>2002</td><td>-9.64</td><td>0.52</td><td>0.04</td><td>0.15</td><td>0.01</td><td>0.04</td><td>4.99</td><td>0.27</td><td>0.02</td></tr><tr><td>2003</td><td>-8.87</td><td>0.51</td><td>0.06</td><td>0.10</td><td>0.01</td><td>0.04</td><td>4.55</td><td>0.26</td><td>0.02</td></tr><tr><td>2004</td><td>-10.87</td><td>0.49</td><td>0.05</td><td>-0.15</td><td>0.01</td><td>0.04</td><td>5.49</td><td>0.25</td><td>0.02</td></tr><tr><td>2005</td><td>-10.31</td><td>0.54</td><td>0.06</td><td>-0.15</td><td>0.05</td><td>0.05</td><td>5.31</td><td>0.27</td><td>0.02</td></tr><tr><td>2006</td><td>-9.81</td><td>0.45</td><td>0.04</td><td>0.08</td><td>0.03</td><td>0.03</td><td>5.10</td><td>0.23</td><td>0.02</td></tr><tr><td>2007</td><td>-10.45</td><td>0.44</td><td>0.05</td><td>-0.52</td><td>0.07</td><td>0.03</td><td>5.29</td><td>0.14</td><td>0.02</td></tr><tr><td>2008</td><td>-8.84</td><td>0.50</td><td>0.06</td><td>-0.17</td><td>0.02</td><td>0.05</td><td>4.86</td><td>0.25</td><td>0.02</td></tr><tr><td>2009</td><td>-5.57</td><td>0.51</td><td>0.04</td><td>-0.51</td><td>0.00</td><td>0.04</td><td>3.13</td><td>0.26</td><td>0.02</td></tr></tbody></table>

**Notes.** Low and high estimates of uncertainties correspond to two assumptions, see text in Sect. 2.3.1.

## A.1. Table with the reconstructed RV from Ulrich (2010)

Table A.1 provides the integrated RV due to MC that we reconstructed from the latitudinal profiles of Ulrich (2010). Lower and upper estimates of uncertainties are also provided. The results are discussed in Sect. 2.3.1.

## A.2. MC from Basu & Antia (2010)

Basu & Antia (2010; hereafter BA10) performed a ring diagram analysis of Dopplergrams from MDI/SOHO data from between 1996 and 2009. This technique probes flows in different layers below the photosphere down to about 10 Mm. Here we use the flows corresponding to the closest layer to the surface, i.e. 2.8 Mm deep and the OLA (Optimally Localized Average) inversion (two inversions are provided in BA10, but they give very similar results). BA10 decomposed each latitudinal profile into a sum of associated Legendre polynomials: we digitised these coefficients versus time to reconstruct the profiles. Four even coefficients (corresponding to the antisymmetric component of the flows, c2 to c8) were used for a good reconstruction, and are on average equal to 9.5, -3.6, -0.8, and -2.2 m s⁻¹ respectively, with uncertainties on average of the order of 0.5, 0.7, 0.7, and 0.5 m s⁻¹ respectively. The resulting time series mc(t) are given in Table A.2 for 0°, 55°, and 90°.

## A.3. MC from Hathaway & Rightmire (2010, 2011)

Hathaway & Rightmire (2010; hereafter HR10) used a very different method: they estimated the meridional flows by tracking

<table><caption>**Table A.2.** Integrated meridional flow from Basu & Antia (2010).</caption><thead><tr><th>Time (yr)</th><th>RV at 0° (m s<sup>-1</sup>)</th><th>RV at 55° (m s<sup>-1</sup>)</th><th>RV at 90° (m s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>1996</td><td>-9.87</td><td>0.01</td><td>5.10</td></tr><tr><td>1997</td><td>-7.08</td><td>-0.10</td><td>3.91</td></tr><tr><td>1998</td><td>-8.15</td><td>-0.06</td><td>4.45</td></tr><tr><td>1999</td><td>-5.22</td><td>-0.26</td><td>3.31</td></tr><tr><td>2000</td><td>-4.11</td><td>-0.30</td><td>2.86</td></tr><tr><td>2001</td><td>-5.82</td><td>-0.08</td><td>3.26</td></tr><tr><td>2002</td><td>-4.99</td><td>-0.17</td><td>2.94</td></tr><tr><td>2002</td><td>-3.34</td><td>-0.35</td><td>2.51</td></tr><tr><td>2003</td><td>-6.61</td><td>-0.19</td><td>3.84</td></tr><tr><td>2003</td><td>-6.49</td><td>-0.09</td><td>3.56</td></tr><tr><td>2004</td><td>-7.61</td><td>-0.17</td><td>4.27</td></tr><tr><td>2005</td><td>-7.74</td><td>-0.19</td><td>4.39</td></tr><tr><td>2006</td><td>-8.47</td><td>-0.03</td><td>4.44</td></tr><tr><td>2008</td><td>-9.25</td><td>0.03</td><td>4.72</td></tr><tr><td>2009</td><td>-9.41</td><td>-0.03</td><td>4.97</td></tr></tbody></table>

magnetic features in magnetograms from MDI/SOHO from between 1996 and 2009 because they were interested to see the flows actually seen by these magnetic features. This was crucial to make better estimates of the contribution of these flows in flux transport dynamo models, in which they play a very significant role. We use their time average latitudinal profile, which is then scaled by the amplitude versus time, both digitised from
<figure><img src="image_19.png" /><figcaption>Fig. A.1. Time series of integrated meridional flows for pole-on (upper panel) and edge-on configuration (lower panel), for the U10 reconstruction (black), BA10 (brown), HR10 (res), and HR11 (orange).</figcaption></figure>

HR10. Unlike the results from U10 and BA10, which were averaged over one year, the temporal cadence is much better here: the MC latitudinal profiles are averaged over the rotation period instead of one year. The use of the averaged profile means that any converging flow pattern superposed on the relatively smooth profile is not taken into account here. These flows, in addition to being weaker than those derived by U10, also present less polar reversals. This may be due to the fact that they may probe different layers (below the photosphere) and/or may represent how the magnetic flux tubes move across the flows, as discussed by Ulrich (2010).

Hathaway & Rightmire (2011; hereafter HR11) use the same technique, and achieve results that are only slightly different and with a slightly longer temporal coverage (1996–2010). These

authors modelled the latitudinal profiles with two functions in latitude: we digitised the two corresponding coefficients $S_1(t)$ and $S_3(t)$ and reconstructed the profiles using:

$$mc(t, \theta) = S_1(t)(2 \sin(\theta) \cos(\theta)) + S_3(t) \times (7 \sin^2(\theta) - 3 \sin(\theta)) \cos(\theta). \quad (A.1)$$

As in HR10, the description of the flows does not include complex patterns such as the converging flows in the activity belt. The coefficients $S_1(t)$ and $S_3(t)$ are on average equal to 11 and $-2.2 \text{ m s}^{-1}$ respectively (with dispersion over time of the order of $2 \text{ m s}^{-1}$), and the typical uncertainties are 0.9 and $0.4 \text{ m s}^{-1}$ respectively. We only show the table (Table B.1) for Hathaway & Rightmire (2011) since the results are very similar to those from Hathaway & Rightmire (2010), but slightly longer.

## A.4. Comparison between time series

Figure A.1 compares the time series of integrated meridional flows for U10, B10, HR10, and HR11. The amplitudes are of the same order of magnitude, but B10, HR10, and HR11 are more anti-correlated with the cycle phase than U10. This may be due to the fact that they do not trace the same layers.

# Appendix B: Uncertainties on integrated RV

We computed uncertainties on each value of mc(t) based on two different assumptions. First, we consider that for a given latitudinal profile (one of those shown in Fig. 1), the uncertainties between latitudinal bins are uncorrelated: we compute synthetic latitudinal profiles which are the sum of the original one added to a random deviation derived from a Gaussian distribution corresponding to the uncertainty in this latitude bin. We perform 25 such realisations of the time series: the rms of the RV values after integration at each time step is then the uncertainty for a given time step. The resulting uncertainties are very small, of the order of $0.06 \text{ m s}^{-1}$ for pole-on time series and $0.02 \text{ m s}^{-1}$ for the edge-on one.

However, we cannot exclude the fact that part of the uncertainties on the latitudinal profiles are correlated, and correspond to an uncertainty on the global amplitude of the MC. Therefore, we also consider the possibility that all these uncertainties are correlated, i.e. that the uncertainty is mostly due to a bias on the global amplitude of the latitudinal profile, to provide an upper limit. Here we compute two integrated RV values, one for the upper envelope of the latitudinal profile, and the other one for the lower envelope, which provides an upper limit of the uncertainty on the integrated flow at each time step, shown in Fig. 3: the typical uncertainties on integrated RV are of the order of $1 \text{ m s}^{-1}$ for the pole-on time series and $0.5 \text{ m s}^{-1}$ for the edge-on one. Synthetic time series of mc(t) assuming such uncertainties are then generated in order to compute uncertainties on other variables such as the rms, amplitude, or correlation (see Sect. 2.3.2). It is likely that the true values lie between these two estimates.
<table><caption>**Table B.1.** Integrated meridional flow from Hathaway & Rightmire (2011).</caption><thead><tr><th>Time (yr)</th><th>RV at 0° (m s<sup>-1</sup>)</th><th>RV at 55° (m s<sup>-1</sup>)</th><th>RV at 90° (m s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>1996.52</td><td>-5.98</td><td>-0.03</td><td>3.16</td></tr><tr><td>1996.68</td><td>-5.25</td><td>-0.07</td><td>2.87</td></tr><tr><td>1996.75</td><td>-5.51</td><td>-0.03</td><td>2.91</td></tr><tr><td>1996.82</td><td>-6.01</td><td>-0.08</td><td>3.28</td></tr><tr><td>1996.90</td><td>-5.22</td><td>-0.02</td><td>2.73</td></tr><tr><td>1996.97</td><td>-5.42</td><td>-0.05</td><td>2.90</td></tr><tr><td>1997.06</td><td>-4.80</td><td>-0.05</td><td>2.58</td></tr><tr><td>1997.13</td><td>-5.93</td><td>-0.06</td><td>3.20</td></tr><tr><td>1997.20</td><td>-5.53</td><td>-0.03</td><td>2.92</td></tr><tr><td>1997.28</td><td>-5.64</td><td>-0.05</td><td>3.01</td></tr><tr><td>1997.38</td><td>-5.72</td><td>-0.06</td><td>3.09</td></tr><tr><td>1997.45</td><td>-5.12</td><td>-0.01</td><td>2.66</td></tr><tr><td>1997.55</td><td>-5.37</td><td>-0.12</td><td>3.04</td></tr><tr><td>1997.61</td><td>-4.84</td><td>-0.05</td><td>2.62</td></tr><tr><td>1997.68</td><td>-4.51</td><td>-0.04</td><td>2.41</td></tr><tr><td>1997.76</td><td>-4.58</td><td>-0.04</td><td>2.45</td></tr><tr><td>1997.84</td><td>-4.76</td><td>-0.07</td><td>2.62</td></tr><tr><td>1997.90</td><td>-4.62</td><td>-0.02</td><td>2.43</td></tr><tr><td>1997.98</td><td>-5.39</td><td>0.01</td><td>2.74</td></tr><tr><td>1998.06</td><td>-3.46</td><td>-0.04</td><td>1.88</td></tr><tr><td>1998.13</td><td>-3.45</td><td>-0.08</td><td>1.96</td></tr><tr><td>1998.20</td><td>-4.93</td><td>0.03</td><td>2.47</td></tr><tr><td>1998.28</td><td>-4.76</td><td>0.02</td><td>2.40</td></tr><tr><td>1998.36</td><td>-5.33</td><td>0.02</td><td>2.69</td></tr><tr><td>1998.43</td><td>-4.10</td><td>0.01</td><td>2.08</td></tr><tr><td>1998.82</td><td>-4.15</td><td>-0.00</td><td>2.14</td></tr><tr><td>1999.18</td><td>-3.39</td><td>-0.02</td><td>1.79</td></tr><tr><td>1999.26</td><td>-4.66</td><td>0.03</td><td>2.34</td></tr><tr><td>1999.34</td><td>-3.35</td><td>0.07</td><td>1.57</td></tr><tr><td>1999.39</td><td>-4.46</td><td>0.00</td><td>2.29</td></tr><tr><td>1999.48</td><td>-4.92</td><td>0.05</td><td>2.42</td></tr><tr><td>1999.55</td><td>-4.67</td><td>-0.01</td><td>2.43</td></tr><tr><td>1999.63</td><td>-4.54</td><td>-0.02</td><td>2.38</td></tr><tr><td>1999.70</td><td>-5.15</td><td>0.01</td><td>2.64</td></tr><tr><td>1999.77</td><td>-4.66</td><td>-0.04</td><td>2.48</td></tr><tr><td>1999.85</td><td>-4.80</td><td>0.03</td><td>2.40</td></tr><tr><td>1999.93</td><td>-5.51</td><td>-0.00</td><td>2.84</td></tr><tr><td>1999.99</td><td>-4.58</td><td>0.11</td><td>2.09</td></tr><tr><td>2000.07</td><td>-3.95</td><td>0.04</td><td>1.95</td></tr><tr><td>2000.15</td><td>-3.25</td><td>-0.00</td><td>1.68</td></tr><tr><td>2000.22</td><td>-3.90</td><td>-0.05</td><td>2.12</td></tr><tr><td>2000.29</td><td>-4.85</td><td>0.01</td><td>2.47</td></tr><tr><td>2000.38</td><td>-3.78</td><td>-0.01</td><td>1.97</td></tr><tr><td>2000.45</td><td>-4.48</td><td>-0.05</td><td>2.41</td></tr><tr><td>2000.52</td><td>-4.62</td><td>-0.02</td><td>2.42</td></tr><tr><td>2000.60</td><td>-3.96</td><td>-0.03</td><td>2.12</td></tr><tr><td>2000.68</td><td>-4.72</td><td>-0.03</td><td>2.50</td></tr><tr><td>2000.74</td><td>-4.67</td><td>0.02</td><td>2.36</td></tr><tr><td>2000.84</td><td>-4.48</td><td>0.01</td><td>2.28</td></tr><tr><td>2000.89</td><td>-3.61</td><td>-0.06</td><td>1.99</td></tr><tr><td>2000.98</td><td>-4.12</td><td>0.06</td><td>1.99</td></tr><tr><td>2001.05</td><td>-4.50</td><td>0.10</td><td>2.09</td></tr><tr><td>2001.11</td><td>-5.67</td><td>0.02</td><td>2.88</td></tr><tr><td>2001.18</td><td>-3.47</td><td>0.04</td><td>1.70</td></tr><tr><td>2001.28</td><td>-3.92</td><td>0.00</td><td>2.02</td></tr><tr><td>2001.34</td><td>-3.89</td><td>0.03</td><td>1.92</td></tr><tr><td>2001.41</td><td>-5.11</td><td>-0.03</td><td>2.71</td></tr><tr><td>2001.48</td><td>-4.76</td><td>0.01</td><td>2.43</td></tr><tr><td>2001.58</td><td>-4.50</td><td>0.01</td><td>2.30</td></tr></tbody></table>

**Table B.1.** continued.

<table><thead><tr><th>Time (yr)</th><th>RV at 0° (m s<sup>-1</sup>)</th><th>RV at 55° (m s<sup>-1</sup>)</th><th>RV at 90° (m s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>2001.62</td><td>-5.10</td><td>0.07</td><td>2.46</td></tr><tr><td>2001.72</td><td>-4.79</td><td>0.03</td><td>2.39</td></tr><tr><td>2001.80</td><td>-5.65</td><td>0.07</td><td>2.74</td></tr><tr><td>2001.87</td><td>-5.14</td><td>0.00</td><td>2.64</td></tr><tr><td>2001.96</td><td>-5.09</td><td>0.01</td><td>2.61</td></tr><tr><td>2002.01</td><td>-3.27</td><td>-0.01</td><td>1.70</td></tr><tr><td>2002.09</td><td>-4.34</td><td>-0.05</td><td>2.35</td></tr><tr><td>2002.16</td><td>-4.62</td><td>0.10</td><td>2.14</td></tr><tr><td>2002.24</td><td>-5.17</td><td>-0.00</td><td>2.66</td></tr><tr><td>2002.30</td><td>-5.81</td><td>0.01</td><td>2.97</td></tr><tr><td>2002.38</td><td>-5.28</td><td>0.08</td><td>2.54</td></tr><tr><td>2002.46</td><td>-4.56</td><td>0.01</td><td>2.33</td></tr><tr><td>2002.55</td><td>-5.73</td><td>0.03</td><td>2.87</td></tr><tr><td>2002.61</td><td>-4.84</td><td>-0.00</td><td>2.50</td></tr><tr><td>2002.69</td><td>-4.48</td><td>0.01</td><td>2.29</td></tr><tr><td>2002.77</td><td>-4.59</td><td>-0.01</td><td>2.39</td></tr><tr><td>2002.85</td><td>-4.51</td><td>0.01</td><td>2.31</td></tr><tr><td>2002.92</td><td>-5.31</td><td>0.02</td><td>2.68</td></tr><tr><td>2002.99</td><td>-5.14</td><td>0.04</td><td>2.55</td></tr><tr><td>2003.06</td><td>-5.08</td><td>0.01</td><td>2.58</td></tr><tr><td>2003.14</td><td>-5.11</td><td>0.05</td><td>2.52</td></tr><tr><td>2003.22</td><td>-3.95</td><td>0.01</td><td>2.02</td></tr><tr><td>2003.29</td><td>-5.08</td><td>0.07</td><td>2.45</td></tr><tr><td>2003.36</td><td>-6.23</td><td>0.01</td><td>3.18</td></tr><tr><td>2003.43</td><td>-5.15</td><td>-0.03</td><td>2.73</td></tr><tr><td>2003.52</td><td>-5.76</td><td>0.01</td><td>2.95</td></tr><tr><td>2003.57</td><td>-6.10</td><td>0.06</td><td>3.00</td></tr><tr><td>2003.67</td><td>-6.06</td><td>0.02</td><td>3.09</td></tr><tr><td>2003.74</td><td>-6.27</td><td>0.04</td><td>3.13</td></tr><tr><td>2003.80</td><td>-6.35</td><td>0.01</td><td>3.24</td></tr><tr><td>2003.88</td><td>-5.90</td><td>0.01</td><td>3.03</td></tr><tr><td>2003.96</td><td>-5.54</td><td>0.03</td><td>2.78</td></tr><tr><td>2004.04</td><td>-5.33</td><td>-0.08</td><td>2.92</td></tr><tr><td>2004.11</td><td>-4.87</td><td>-0.06</td><td>2.64</td></tr><tr><td>2004.18</td><td>-5.68</td><td>-0.06</td><td>3.06</td></tr><tr><td>2004.25</td><td>-5.89</td><td>-0.03</td><td>3.10</td></tr><tr><td>2004.32</td><td>-5.98</td><td>0.03</td><td>3.01</td></tr><tr><td>2004.41</td><td>-5.78</td><td>-0.02</td><td>3.02</td></tr><tr><td>2004.49</td><td>-7.30</td><td>-0.03</td><td>3.82</td></tr><tr><td>2004.56</td><td>-5.03</td><td>0.02</td><td>2.54</td></tr><tr><td>2004.64</td><td>-8.27</td><td>-0.10</td><td>4.48</td></tr><tr><td>2004.71</td><td>-6.13</td><td>0.02</td><td>3.12</td></tr><tr><td>2004.78</td><td>-5.23</td><td>-0.01</td><td>2.71</td></tr><tr><td>2004.85</td><td>-6.61</td><td>-0.03</td><td>3.47</td></tr><tr><td>2004.94</td><td>-5.74</td><td>-0.02</td><td>3.01</td></tr><tr><td>2005.00</td><td>-7.06</td><td>0.01</td><td>3.60</td></tr><tr><td>2005.07</td><td>-5.93</td><td>-0.01</td><td>3.07</td></tr><tr><td>2005.14</td><td>-4.78</td><td>-0.02</td><td>2.50</td></tr><tr><td>2005.24</td><td>-5.95</td><td>0.00</td><td>3.06</td></tr><tr><td>2005.30</td><td>-6.31</td><td>0.08</td><td>3.06</td></tr><tr><td>2005.38</td><td>-5.41</td><td>-0.07</td><td>2.95</td></tr><tr><td>2005.44</td><td>-7.80</td><td>-0.01</td><td>4.04</td></tr><tr><td>2005.53</td><td>-7.95</td><td>0.01</td><td>4.08</td></tr><tr><td>2005.61</td><td>-7.66</td><td>0.01</td><td>3.92</td></tr><tr><td>2005.67</td><td>-6.76</td><td>-0.03</td><td>3.56</td></tr><tr><td>2005.75</td><td>-7.30</td><td>0.00</td><td>3.75</td></tr><tr><td>2005.82</td><td>-6.43</td><td>0.02</td><td>3.27</td></tr><tr><td>2005.90</td><td>-6.22</td><td>-0.03</td><td>3.27</td></tr><tr><td>2005.97</td><td>-5.30</td><td>-0.00</td><td>2.74</td></tr><tr><td>2006.04</td><td>-5.86</td><td>0.02</td><td>2.98</td></tr></tbody></table>
<table><caption>Table B.1. continued.</caption><thead><tr><th>Time (yr)</th><th>RV at 0° (m s<sup>-1</sup>)</th><th>RV at 55° (m s<sup>-1</sup>)</th><th>RV at 90° (m s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>2006.12</td><td>-6.37</td><td>-0.07</td><td>3.45</td></tr><tr><td>2006.20</td><td>-7.72</td><td>-0.08</td><td>4.17</td></tr><tr><td>2006.26</td><td>-7.07</td><td>0.04</td><td>3.55</td></tr><tr><td>2006.35</td><td>-6.87</td><td>0.02</td><td>3.49</td></tr><tr><td>2006.42</td><td>-7.00</td><td>-0.03</td><td>3.68</td></tr><tr><td>2006.50</td><td>-7.15</td><td>-0.04</td><td>3.77</td></tr><tr><td>2006.57</td><td>-7.42</td><td>-0.03</td><td>3.88</td></tr><tr><td>2006.64</td><td>-7.56</td><td>-0.03</td><td>3.96</td></tr><tr><td>2006.72</td><td>-6.71</td><td>0.07</td><td>3.30</td></tr><tr><td>2006.78</td><td>-6.49</td><td>-0.03</td><td>3.42</td></tr><tr><td>2006.84</td><td>-6.68</td><td>-0.03</td><td>3.50</td></tr><tr><td>2006.92</td><td>-7.67</td><td>-0.02</td><td>3.99</td></tr><tr><td>2006.99</td><td>-7.11</td><td>-0.04</td><td>3.76</td></tr><tr><td>2007.06</td><td>-5.52</td><td>-0.04</td><td>2.94</td></tr><tr><td>2007.13</td><td>-6.37</td><td>-0.08</td><td>3.45</td></tr><tr><td>2007.20</td><td>-6.59</td><td>-0.11</td><td>3.65</td></tr><tr><td>2007.28</td><td>-5.97</td><td>-0.03</td><td>3.14</td></tr><tr><td>2007.34</td><td>-6.72</td><td>0.02</td><td>3.40</td></tr><tr><td>2007.43</td><td>-6.94</td><td>0.01</td><td>3.56</td></tr><tr><td>2007.51</td><td>-6.50</td><td>-0.02</td><td>3.40</td></tr><tr><td>2007.59</td><td>-7.51</td><td>0.02</td><td>3.82</td></tr><tr><td>2007.66</td><td>-6.02</td><td>-0.01</td><td>3.14</td></tr><tr><td>2007.72</td><td>-7.23</td><td>0.02</td><td>3.68</td></tr><tr><td>2007.80</td><td>-6.84</td><td>-0.04</td><td>3.61</td></tr><tr><td>2007.88</td><td>-7.46</td><td>-0.03</td><td>3.92</td></tr><tr><td>2007.96</td><td>-5.97</td><td>-0.03</td><td>3.15</td></tr><tr><td>2008.03</td><td>-4.08</td><td>-0.02</td><td>2.16</td></tr><tr><td>2008.11</td><td>-5.65</td><td>-0.05</td><td>3.03</td></tr><tr><td>2008.18</td><td>-6.01</td><td>-0.01</td><td>3.12</td></tr><tr><td>2008.26</td><td>-6.78</td><td>0.02</td><td>3.44</td></tr><tr><td>2008.33</td><td>-6.64</td><td>0.02</td><td>3.37</td></tr><tr><td>2008.41</td><td>-7.00</td><td>-0.03</td><td>3.66</td></tr><tr><td>2008.47</td><td>-5.83</td><td>-0.02</td><td>3.05</td></tr><tr><td>2008.55</td><td>-6.89</td><td>0.00</td><td>3.54</td></tr><tr><td>2008.64</td><td>-7.34</td><td>-0.03</td><td>3.85</td></tr><tr><td>2008.68</td><td>-6.22</td><td>0.03</td><td>3.13</td></tr><tr><td>2008.78</td><td>-5.92</td><td>0.01</td><td>3.03</td></tr><tr><td>2008.84</td><td>-5.99</td><td>-0.03</td><td>3.15</td></tr><tr><td>2008.93</td><td>-6.23</td><td>-0.04</td><td>3.31</td></tr><tr><td>2009.00</td><td>-6.53</td><td>0.02</td><td>3.32</td></tr><tr><td>2009.07</td><td>-6.15</td><td>-0.02</td><td>3.21</td></tr><tr><td>2009.15</td><td>-6.38</td><td>0.01</td><td>3.27</td></tr><tr><td>2009.22</td><td>-5.71</td><td>-0.06</td><td>3.07</td></tr><tr><td>2009.30</td><td>-6.88</td><td>-0.01</td><td>3.57</td></tr><tr><td>2009.36</td><td>-6.05</td><td>-0.07</td><td>3.28</td></tr><tr><td>2009.44</td><td>-7.52</td><td>0.05</td><td>3.77</td></tr><tr><td>2009.52</td><td>-6.77</td><td>-0.03</td><td>3.55</td></tr><tr><td>2009.59</td><td>-8.21</td><td>0.11</td><td>3.97</td></tr><tr><td>2009.66</td><td>-5.92</td><td>0.02</td><td>3.00</td></tr><tr><td>2009.74</td><td>-7.22</td><td>-0.01</td><td>3.75</td></tr><tr><td>2009.82</td><td>-5.81</td><td>0.00</td><td>2.99</td></tr><tr><td>2009.89</td><td>-6.64</td><td>0.02</td><td>3.38</td></tr><tr><td>2009.97</td><td>-6.85</td><td>0.05</td><td>3.42</td></tr><tr><td>2010.05</td><td>-4.85</td><td>-0.03</td><td>2.57</td></tr><tr><td>2010.11</td><td>-5.97</td><td>0.02</td><td>3.04</td></tr><tr><td>2010.18</td><td>-5.50</td><td>0.00</td><td>2.83</td></tr><tr><td>2010.28</td><td>-6.79</td><td>-0.01</td><td>3.53</td></tr><tr><td>2010.34</td><td>-6.12</td><td>0.02</td><td>3.11</td></tr><tr><td>2010.42</td><td>-6.54</td><td>0.02</td><td>3.33</td></tr><tr><td>2010.48</td><td>-5.75</td><td>-0.02</td><td>3.00</td></tr><tr><td>2010.57</td><td>-5.79</td><td>-0.04</td><td>3.08</td></tr></tbody></table>