Free Access
Issue
A&A
Volume 519, September 2010
Article Number A66
Number of page(s) 8
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201014199
Published online 14 September 2010
A&A 519, A66 (2010)

Reconstructing the solar integrated radial velocity using MDI/SOHO

N. Meunier - A.-M. Lagrange - M. Desort

Laboratoire d'Astrophysique de Grenoble, Observatoire de Grenoble, Université Joseph Fourier, CNRS, UMR 5571, 38041 Grenoble Cedex 09, France

Received 4 February 2010 / Accepted 12 May 2010

Abstract
Context. Searches for exoplanets with radial velocity techniques are increasingly sensitive to stellar activity. It is therefore crucial to characterize how this activity influences radial velocity measurements in their study of the detectability of planets in these conditions.
Aims. In a previous work we simulated the impact of spots and plages on the radial velocity of the Sun. Our objective is to compare this simulation with the observed radial velocity of the Sun for the same period.
Methods. We use Dopplergrams and magnetograms obtained by MDI/SOHO over one solar cycle to reconstruct the solar integrated radial velocity in the Ni line 6768 Å. We also characterize the relation between the velocity and the local magnetic field to interpret our results.
Results. We obtain a stronger redshift in places where the local magnetic field is larger (and as a consequence for larger magnetic structures): hence we find a higher attenuation of the convective blueshift in plages than in the network. Our results are compatible with an attenuation of this blueshift by about 50% when averaged over plages and network. We obtain an integrated radial velocity with an amplitude over the solar cycle of about 8 m s-1, with small-scale variations similar to the results of the simulation, once they are scaled to the Ni line.
Conclusions. The observed solar integrated radial velocity agrees with the result of the simulation made in our previous work within 30%, which validates this simulation. The observed amplitude confirms that the impact of the convective blueshift attenuation in magnetic regions will be critical to detect Earth-mass planets in the habitable zone around solar-like stars.

Key words: techniques: radial velocities - Sun: activity - Sun: surface magnetism - stars: early-type

1 Introduction

Searches for exoplanets with radial velocity (RV) techniques will be increasingly sensitive to stellar activity as the sensitivity of instruments is improved (such as Espresso on the VLT or Codex on the E-ELT D'Odorico & the CODEX/ESPRESSO team, 2007). It is therefore crucial to understand how this activity influences the RV measurements. In our previous works, we simulated the influence of spots (Lagrange et al., 2009, hereafter Paper I) and plages (Meunier et al., 2010, hereafter Paper II) on the RV for the Sun seen equator-on. The simulation covered a solar cycle with an excellent temporal cadence. Our study considered the photometric contribution of spots and plages and also the attenuation of the convective blueshift in magnetic regions (Paper II). This last effect appeared to be dominant, with a long-term amplitude over the solar cycle of about 8 m s-1, while the photometric contributions produced a rms of 0.08 m s-1 during cycle minimum and of 0.42 m s-1 during cycle maximum (0.33 m s-1 over the whole cycle), with peak-to-peak amplitudes of a few m s-1. We studied how these RV variations influence the detectability of an Earth-mass planet in the habitable zone and found that it made this detection very difficult.

Most existing observations of the solar RV do not cover our data set. Neither are the investigated lines well characterized in terms of convective blueshift, which makes them difficult to compare with simulations. McMillan et al. (1993) measured the solar RV with the moon light observed in the violet part of the spectrum for the period 1987-1992. They found solar RV variations that are significantly smaller than our simulation. For the period 1983-1992, Deming & Plymate (1994) measured the solar RV in the 2.3 $\mu $m domain while Jimenez et al. (1986) measured it in the potassium line at 769.9 nm (Brookes et al., 1978) during cycle 21 (1976-1984). Both have found an amplitude over the cycle of about 30 m s-1. The potassium line is a strong line however, for which a particularly high convective blueshift is not expected. The GOLF (Gabriel et al., 1995) instrument on SOHO provides integrated measurements of the solar RV in the Na D1 and D2 lines. Unfortunately, the planned calibration procedure could not be applied to the GOLF observations because of an instrumental failure early in the mission. Therefore the long-term variations of the RV signal and its amplitudes are unreliable (Boumier, priv. comm.). We note that after a few years of operation, a comparison has been made between the MDI and GOLF signals (Henney et al., 1998) but focused on the high temporal frequency component and not on the long-term signal. Some attempts to calibrate the long term variations a posteriori, mostly due to the spacecraft orbit, have been made (Ulrich et al., 2000) but the resulting signal still shows some strong instrumental residuals that are not compatible with the amplitude we wish to test.

This is why precise RV observations of the Sun in integrated light are necessary for testing our simulations, in particular for the convective blueshift contribution. It would be best to directly observe the Sun in integrated light over an extended domain of the spectrum with an instrument and data analysis similar to those used to measure stellar RV variations. This is a long-term project though. The aim of this paper is to characterize the solar RV variations with spatially resolved observations by MDI/SOHO in a single spectral line (Ni line at 6768 Å). The use of an imaging instrument also allows us to identify the contribution of solar activity to the solar integrated velocities, in contrast with an instrument such as GOLF. We reconstruct the integrated RV variations (neglecting the photometric component of spots and plages simulated in Papers I and II) and characterize the behavior of the velocity variations measured in the Ni line with respect to the magnetic field. Our observations are compared to the simulations performed in Paper II scaled to the expected behavior of that line.

The input data and data analysis are described in Sect. 2. We describe how we determine the zero velocity for each data point in the time series and estimate the uncertainty of our measurements. The local velocity-magnetic field relation derived from the analysis of the Ni line is presented in Sect. 3. This relation allows the interpretation of the integrated RV variations. We also discuss the origin of the velocity associated to magnetic regions in term of convective blueshift attenuation and downflows. Then we present the integrated RV in Sect. 4 and compare it to the results of the simulation of the convective blueshift attenuation made in Paper II scaled to the Ni line expected behavior. The comparison is made at different time-scales. We conclude in Sect. 5.

2 Observation and data analysis

2.1 Input data

We use a series of MDI/SOHO Dopplergrams and magnetograms (Scherrer et al., 1995) covering a solar cycle between 1996 May 5 (Julian day 2 450 209) and 2007 October 7 (Julian day 2 454 380), for the Ni 6768 line. The intensity in the line core is about 0.36 times the continuum intensity. Dopplergrams are maps of the line-of-sight velocity (positive velocities correspond to a redshift), and are dominated by the large-scale motions (mostly the rotation) and small-scale motions mostly due to granulation and supergranulation. An example is shown in Fig. 1 (uppper left panel). Magnetograms are maps of the line-of-sight magnetic field averaged over each pixel, i.e. flux maps. We selected Dopplergrams obtained on the same days as the magnetograms analyzed in Paper II and separated in time by less than six hours, which leads to 3181 data points (out of which 74% corresponds to a separation less than one hour).

The MDI data observations are made only for the Ni 6768 line. To derive a convective blueshift for our entire spectral range, we used the empirical relation between the convective blueshift and the line depths (Gray, 2009) in Paper II. The average convective blueshift was also weighted by the relative contributions of various lines to the correlation function used to compute radial velocities, to produce a convective blueshift valid for the entire wavelength range. Given the depth of the Ni line, and an attenuation by a factor two-thirds (Brandt & Solanki, 1990) of the convective blueshift in magnetic regions as in Paper II, we expect a vertical velocity $\Delta V$ in the magnetic regions for that line of 227 m s-1 (instead of 190 m s-1 for the entire spectrum). This value is used in Sect. 4 to scale the results of the simulation of Paper II to compare it with our specific conditions.

2.2 Dopplergram analysis

\begin{figure}
\par\includegraphics{14199fig1.ps}
\end{figure} Figure 1:

Upper left panel: Dopplergram for 2000 May 16 (in m s-1). Upper right panel: magnetogram for the same day (in Gauss). Lower left panel: Dopplergram after interpolation (in m s-1). Lower right panel: mask derived from the magnetogram (on the disk, 0 where the Dopplergram is interpolated).

Open with DEXTER

Radial velocities measured by MDI do not present a good long-term stability (e.g. Beck et al., 1998; Evans et al., 1998), mostly due to the orbit of the spacecraft (amplitude of about 500 m s-1) and semi-annual changes in the wave plate angles, producing a typical drift of 1 m s-1 per day and steps of several hundreds of m s-1 typically twice a year. The zero of MDI Dopplergrams is not sufficiently good to directly derive an integrated radial velocity from the observations and we therefore need to define this zero using the data themselves.

The principle of our approach is the following (the details will be given below). We assume that the quiet Sun provides an appropriate zero velocity map. Using the magnetogram observed at a given time, we define a mask by selecting the magnetic regions. This mask is applied to the corresponding (in time) Dopplergram to identify regions where the velocity is affected by magnetic activity. The velocity in those regions will have to be interpolated from the unmasked areas (which represent the quiet Sun). Then we average the velocities over the disk, which produces a reference for the radial velocity variations.

The masks were derived from each magnetogram. First, we corrected the magnetogram for the rotation: if the Dopplergram and magnetogram are not obtained strictly simultaneously, the rotation induces a difference in position between a given region on the two images. A threshold of 50 G was then applied on the slightly smoothed magnetograms to define magnetic regions, producing a binary mask. This mask was then enlarged (10 pixels around each structure). This allowed us to eliminate most flows caused by magnetic regions without eliminating too many pixels. An example of such a mask is shown in Fig. 1 (lower right panel).

The mask was then used to define zones in the Dopplergram where the velocity is interpolated from quiet Sun values around the mask. Because we aim to determine a zero velocity that would correspond to the quiet Sun, we built a Dopplergram that would be similar to a quiet Sun Dopplergram. For each pixel we used a linear interpolation between four pixels after smoothing smoothing the Dopplergrams. These four pixels were the closest pixels outside the mask that are respectively on the right and left (same line) and above and below (along the same column) the considered pixel. Figure 1 (lower left panel) shows a Dopplergram interpolated in this way.

The integrated velocity is derived from the actual Dopplergram (thus providing $V_{\rm obs}$) and from the interpolated reference Dopplergram (thus providing $V_{\rm zero}$) after a weighting by the center-to-limb darkening shown in Fig. 2. If we were actually observing the Sun in integrated light directly, the velocity corresponding to each pixel on the Sun would indeed be weighted by the corresponding intensity. Figure 3 shows the complete time series for  $V_{\rm zero}$ and  $V_{\rm obs}$. Note that in principle it should be possible to use the spacecraft velocity as reference. However, the use of the spacecraft radial velocity shows that some strong residuals, which are related to the instrument, are still present.

2.3 Uncertainties

\begin{figure}
\par\includegraphics{14199fig2.ps}
\end{figure} Figure 2:

Limb-darkening deduced from 78 MDI intensity maps during a very quiet period.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=7.7cm,clip]{14199fig3.ps}
\end{figure} Figure 3:

Upper panel: zero velocity ( $V_{\rm zero}$) determined on interpolated Dopplergrams. Lower panel: velocity ( $V_{\rm obs}$) integrated over Dopplergrams.

Open with DEXTER
\begin{figure}
\par\includegraphics[width=7.7cm,clip]{14199fig4.ps}
\end{figure} Figure 4:

Upper panel: estimated errors of the velocity versus time. Lower panel: distribution of the estimated errors.

Open with DEXTER

The identified sources of uncertainties are the statistical noise in the data caused by the actual motions at the granular and supergranular scales, the Evershed effect around sunspots, the interpolation made to estimate the zero velocity in the regions of magnetic activity and the exact definition of the mask. We study the relative importance of these factors.

To test the zero of the velocity, we considered a very quiet Dopplergram, for which we assumed the observed velocity to be our reference  $V_{\rm ref}$(i.e. we assume that the magnetic regions do not significantly influence the average velocity in that case). Then we applied to that same Dopplergram the mask derived for each data point in our time series, interpolated, and computed the radial velocity V(t) of the corrected Dopplergrams. Ideally the two should be identical. Figure 4 shows $V(t)-V_{\rm ref}$ for the entire time series (3415 data points) and the distribution of the values. The rms is 0.4 m s-1, and is lower during cycle minimum (0.20-0.25 m s-1) than during cycle maximum (0.5 m s-1). Part of this dispersion is due to the interpolation. The dispersion in velocity over the solar disk is also a source of uncertainty (and is included in the previous measurement). After removing large-scale motions (solar rotation and meridional circulation), the rms is 220 m s-1, which leads to an average uncertainty of about 0.3 m s-1 considering the number of pixels. This is not much below 0.4 m s-1, which shows that the contribution of the interpolation to the uncertainty is small (probably not more than 0.2-0.3 m s-1 during high-activity periods). Averaging the Dopplergrams in time (for example one hour) does not help in reducing this noise because supergranulation dominates on longer time scales so that the rms remains large.

We also built similar time series using other quiet Dopplergrams: they lead to a similar dispersion. However, the exact values are different (two of these series are not correlated), so this time series cannot be used to correct the RV signal for instance, only the dispersion is of interest. Note that part of the uncertainty may arise because even during the quiet periods that we consider, there are still a few magnetic structures: depending on the mask some may be eliminated or not, which adds an uncertainty to this error estimate itself.

We also tested the impact of the mask. As described in the previous section, the main parameters to derive the mask are the threshold used to define magnetic regions, and the dilatation of the resulting mask by 10 pixels around each structure. This last effect dominates the first one because the magnetic regions are defined with a much better precision. It is difficult to use a lower value, because the Evershed effect would then significantly affect the zero, especially for large active regions. We therefore tested slightly different values above 10 pixels, keeping in mind that when considering too high values the resulting pixels will become too few to perform a reliable interpolation. Then we compute the rms of the difference between pairs of series. It appears to be difficult to derive an exact value of the uncertainties, as the rms varies depending on the pair, but values lie typically between 0.2 and 0.5 m s-1 for differences between 1 and 4 pixels for the enlargment, which gives an estimate on the order of magnitude of the uncertainty. This estimation is conservative because part of this rms is caused by the interpolation, which is already taken into account. As the mask is increased, we also expect the uncertainty due to the interpolation to be larger than the 0.4 m s-1 uncertainty previously estimated.

Finally, we note that the radial velocity measurements by MDI in regions where the magnetic field is strong may be biased because of the use of a (quiet Sun) reference line to determine the velocity from the intensity measurements in the line, as lines are broadened by the Zeeman effect in magnetic regions. Wachter et al. (2006) found that the velocity may be underestimated by about 2% for plages, and by up to 5% in spot umbra. The effect on our reconstructed radial velocity is therefore expected to be small.

2.4 Convective blueshift in magnetic structures

The convective blueshift associated to the Ni line observed by MDI has never been studied. Therefore, these observations also constitute an opportunity of characterizing the velocity as a function of the local magnetic field in that particular case. This way, we also check the assumption made in Paper II for $\Delta V$, which was that granulation is abnormal in the network (e.g. Morinaga et al., 2008). For this specific study we considered image pairs separated by less than 12 min (so that the rotation is less than a pixel at disk center), to avoid the additional uncertainty that would be introduced by the rotation correction. This concerns 1901 image pairs. We considered two approaches, a study of the velocity in each pixel of magnetic regions, and a study of averaged velocities over all magnetic structures. In addition, only structures with $\mu $ (cosinus of the angle between the line-of-sight and the vertical to the surface) larger than 0.8 were considered here (except precised otherwise), i.e. very close to disk center, to eliminate as many projection effects as possible.

We first defined magnetic regions using the threshold of Paper II, to be able to make a direct comparison (i.e. 100 G). We considered structures above 7 part-per-million of the solar disk (hereafter ppm). This led to about 18 millions pixels (10 millions for $\mu $ above 0.9). The 100 G threshold also provided 623 808 structures (and 337 808 for $\mu $ above 0.9). We used a lower threshold (40 G) to study the behavior of weaker flux areas, leading to about 37 million pixels (21 millions for $\mu $ above 0.9), while considering structures above 4 ppm. In addition to the magnetic field, velocity field, and position of the pixel, we also retrieved the size of the region (defined according to the same threshold) that the pixel belongs to.

3 Velocity in magnetic regions

3.1 Velocity versus local magnetic field

\begin{figure}
\par\includegraphics{14199fig5.ps}
\end{figure} Figure 5:

Velocity (solid line) and velocity divided by $\mu $ (dashed line) as a function of $\mu $, in m s-1, for magnetic fields larger than 100 G.

Open with DEXTER

\begin{figure}
\par\includegraphics{14199fig6.ps}
\end{figure} Figure 6:

Upper panel: $V/\mu $ versus the magnetic field, for $\mu >0.9$: for all pixels (solid line), pixels belonging to structures larger than 100 ppm (dotted line), and pixels belonging to structures smaller than 100 ppm (dashed line). The horizontal line correspond to 227 m s-1 (Paper II). Middle panel: number of pixels versus the magnetic field for the same pixel selection. Lower panel: relative contribution to the total RV.

Open with DEXTER

3.1.1 Results

Figure 5 shows the velocity V and $V/\mu $ versus $\mu $, averaged over bins of 0.01 in $\mu $. After division by $\mu $ the velocity is expected to be constant if the velocity field is mostly vertical. The observed curve is quite flat, showing that we are mostly sensitive to vertical velocity fields. Strictly speaking, a portion of these pixels should also be associated to horizontal velocity fields, such as the Evershed effect (in active regions) or supergranular velocity field (in the network). Their contribution should be small, however, because we selected pixels very close to the disk center. Below we consider $V/\mu $ only, refered to as V hereafter for simplicity. The velocity is positive, indicative of a redshift, which is consistent with an attenuation of the convective blueshift.

Figure 6 (solid line) shows the dependence of the velocity on the magnetic field for $\mu >0.9$, averaged over all pixels with the magnetic field in a certain range. We observe a stronger velocity field for larger magnetic fields, indicative of a stronger attenuation of the convective blueshift as expected. The average over the entire disk (for the 100 G threshold) is 123 m s-1, which is lower than the expected 227 m s-1 derived from the work of Brandt & Solanki (1990) and Gray (2009) in Paper II. This value is used in the next section to compare our observations with the simulation. The maximum value for the largest magnetic fields is close to the 340 m s-1 (convective blueshift for the Ni line derived from the assumptions made in Paper II) which would correspond to the total suppression of the convective blueshift for that line, which is indeed expected for umbral regions. The percentage of pixels with these strong fields is lower than the proportion of pixels expected in umbral area (a few percent of all pixels, excluding the network), so we expect a significant proportion of the pixels with a strong RV to belong to umbral regions. The amplitude of the curve presents a factor $\sim$4 only between the small and large magnetic fields (while these vary by a factor 25).

When considering pixels down to 40 G, the average velocity over the entire disk is 83 m s-1. Note that an average velocity of 227 m s-1 corresponds to pixels with a magnetic field larger than 500 G. The lower panel of Fig. 6 shows that weak fields contribute very significantly to the total RV.

3.1.2 Convective blueshift attenuation versus downflows

It is well known that the network magnetic structures can be associated to downflows, both at the granular scale (in intergranules) and at the supergranular scale (in regions of strong converging flows, outlining supergranular cells). Part of what we measure in the network could therefore be due to an actual downflow (which should be taken into account when computing the zero) and not to the attenuation of the convective blueshift. This downflow would also produce a redshift. Because of the poor resolution of full-disk images obtained by MDI (2 arcsec, i.e. similar to a typical granule), the small-scale contribution should be negligible because we average the velocity over granules. We now investigate the possibility that part of the measured velocity associated to magnetic regions could be caused by an actual downflow in supergranules.

It is difficult to observationally make the difference between an actual downflow and an attenuation of the convective blueshift. Miller et al. (1984) for example observed velocities associated to network structure. They found that the three Fe lines they studied were redshifted in the network by 75-200 m s-1 (in excellent agreement with our results), but only when measured from the line wings. The signal was not present in the line core, which they interpreted as a sign that this redshift was not caused by an actual downflow, but by an attenuation of the convective blueshift (which leads to line distorsions). Several results have been obtained by various groups (see the review by Solanki, 1993), and they are generally consistent with a modification of granulation due to magnetic fields.

To investigate this question we here selected pixels according to the size of the region they belong to, because we do not expect any significant downflow bias in active regions due to the altered supergranulation there. This is shown in Fig. 6. We expect pixels in the network (i.e. small structures) to have a higher velocity (because they are predominantly in a downward flow which is superimposed on the redshift produced by the attenuation of the convective blueshift) compared to pixels with the same magnetic field in larger structures. We find that this is true only for pixels with a relatively large magnetic field of at least a few 100 G. These pixels represent only a few % of all network pixels, so we consider that this effect can be neglected in our analysis. Also, the velocity-magnetic field relationship starts to change for structures below $\sim$200 ppm: these are already large structures, for which it would be difficult to have all pixels associated to an actual downflow. We also point out that due to the large mask, the local downflows should not significantly impact the estimation of the zero velocity at each time.

We conclude that most of the velocity amplitude we measure is probably due to the attenuation of the convective blueshift, in agreement with Miller et al. (1984).

\begin{figure}
\par\includegraphics[width=7.2cm,clip]{14199fig7.ps}
\end{figure} Figure 7:

Upper panel: velocity averaged over the structure versus the structure size for the 100 G threshold structures. Middle panel: number of structures versus the size. Lower panel: relative contribution of each category of structures to the total RV versus the size.

Open with DEXTER
\begin{figure}
\par\includegraphics[width=8cm,clip]{14199fig8.ps}
\end{figure} Figure 8:

Upper panel: residual velocity ( $V_{\rm obs}-V_{\rm zero}$). Lower panel: simulated RV variations due to the convective blueshift attenuation (from Paper II) scaled to a $\Delta V$ of 227 m s-1.

Open with DEXTER

3.2 Structure analysis

Figure 7 (upper panel) shows the velocity field averaged over each magnetic structure as defined in Paper II. The velocity scales with the size of the structures, as expected from the larger contribution of strong magnetic field pixels. Note the change in slope at 300 ppm. This dependence of the RV on the size has not been taken into account in our simulation so far. This will be done in a forthcoming paper.

The bottom panel of Fig. 7 also shows the relative contribution to the total velocity field for each category of structures. This relative contribution is the product defined by $N \times S \times V$ (normalized so that the sum of all contributions is equal to 1), where N is the number of structures of size S (middle panel) and V the velocity averaged over the structures of size S (upper panel). The largest contribution comes from large structures (67% for structures larger than 100 ppm of the solar disk, i.e. $\sim$150 Mm2). Small structures still contribute significantly however: for example, structures smaller than 26 ppm (about 80 Mm2, typical of the quiet network, Wang, 1988) represent 14% of the total velocity. The shape of the relative contribution into two components (for sizes above and below 200 ppm respectively) is related to the size distribution (middle panel of Fig. 7). Although the effect is small, there is a change in slope around 200 ppm of this size distribution, which creates the gap for this size: the product $S \times N$ also shows that shape.

4 Integrated radial velocity: results and analysis

4.1 General results

Figure 8 shows the complete time series of the deduced RV signal $V_{\rm obs}-V_{\rm zero}$ together with the simulated RV variations scaled to 227 m s-1 (we recall that this value is derived from the assumptions made in Paper II). The periodogram of the whole time series is very similar to the periodogram for the simulation computed in Paper II (see the bottom of their Fig. 10) and therefore variations on various time-scales agree well. Figure 9 shows these simulated RV as a function of the observed ones, whenever Dopplergrams and simulations are separated in time by less than one hour (corresponding to 2328 points). They are correlated well (correlation of 0.94). The observed long-term variations are significantly above the uncertainties we estimated in Sect. 2.

Note that the observations show very few values below zero (less than 4%), out of which 80% are above -0.2 m s-1. This is compatible with our uncertainty estimate made in Sect. 2.3.

\begin{figure}
\par\includegraphics[width=8cm,clip]{14199fig9.ps}
\end{figure} Figure 9:

Simulated velocity (from Paper II, scaled to 227 m s-1) versus the observed velocity, whenever Dopplergrams and simulations are separated in time by less than one hour.

Open with DEXTER
\begin{figure}
\par\includegraphics[width=7.5cm,clip]{14199fig10.ps}
\end{figure} Figure 10:

Upper panel: simulated RV smoothed over 30 days (from Paper II) scaled to 227 m s-1 (dashed line) and observation (solid), whenever Dopplergrams and simulations are separated in time by less than one hour. Lower panel: ratio between the two (simulation divided by observation). The horizontal dashed line represents the median value (1.42).

Open with DEXTER

4.2 Long-time scale variations

We now investigate the long-term variations of the RV. Below, we smooth the time series over 30 days to facilitate the comparison of the various levels, and focus on the long-term variations (i.e. above the solar rotation period). The smoothing is performed using a running mean (on these plots a boxcar is used, but the use of triangular mean gives similar results). Figure 10 compares the observed RV with the results of the simulation obtained in Paper II for the convection and scaled to 227 m s-1, once smoothed over 30 days. The two curves are quite similar within 30%. The observed long-term amplitude is about 70% of the amplitude obtained in the simulation, which is a good consensus. It shows that our simulation produces the correct order of magnitude. The variations at various scales are also well reproduced.

We now compare the values deduced from our simulation scaled to $\Delta V$ derived from our present MDI analysis (Sect. 3.1). Figure 11 shows the observations compared to the simulation scaled to 123 m s-1 (this corresponds to the average velocity for structures defined as in Paper II) instead of 227 m s-1 as in Fig. 10. Our simulation then predicts a smaller radial velocity signal than what is observed. The difference is because we considered in Paper II only pixels above 100 G. However, pixels below 100 G also contribute to the signal. This is shown by the dotted line, representing the RV simulated in Paper II scaled to 83 m s-1 and to the correct filling factor for this threshold of 40 G (there is a factor 2 between the number of pixels above 100 G and the number of pixels above 40 G): it is very close to the observed velocity.

In conclusion, the simulation in Paper II overestimated the $\Delta V$ that should be taken into account for the considered structures by about 30%; but it underestimated the surface covered by magnetic regions.

Finally, the $\Delta V$ of 227 m s-1 represented 2/3 of the expected convective blueshift for that line (340 m s-1, derived from the Ni line depth). The simulation would reproduce the RV for the structures defined by the 100 G threshold for an average $\Delta V$ of 70% of the 227 m s-1. This corresponds to an average attenuation of the convective blueshift for these structures by a factor 0.47 (instead of 0.67). We insist that the actual value of the $\Delta V$ to be used on a given data set of magnetic structures depends on their definition, because there is a crosstalk between $\Delta V$ and the size of the structures (only the product is important).

\begin{figure}
\par\includegraphics[width=8cm,clip]{14199fig11.ps}
\end{figure} Figure 11:

Smoothed simulation (from Paper II) scaled to 123 m s-1 (dashed line) and to 166 m s-1 (dotted line), and observation (solid line), whenever Dopplergrams and simulations are separated in time by less than one hour.

Open with DEXTER

4.3 Small-time scale variations

In addition to this general agreement for the long-term variations, it is interesting to compare the behavior of the observed and simulated short-term variations. We consider now the residuals after substraction of the smoothed time series defined in the previous section: the rms is 1.14 m s-1 for the observation, compared to 1.26 m s-1 for the simulation. If we multiply this last rms by 0.7 to take into account the difference in scale between observations and the simulation as shown in the previous section, the rms associated to the simulation becomes 0.88 m s-1. Taking into account the uncertainties due to the interpolation (0.4 m s-1), the quadratic difference leads to an extra rms for the observations of 0.6 m s-1. It is confirmed by the rms of the observation minus the simulation scaled by 70%, which leads to an rms of 0.75 m s-1 (and the quadratic difference with 0.4 m s-1 also gives 0.6 m s-1). The additionnal uncertainty due to the mask definition is less precise, but a conservative estimate shows that it lies between 0.2 and 0.5 m s-1. It could therefore be compatible with the remaining 0.6 m s-1, although a small additionnal dispersion may be present. This means that oursimulations reproduce the observations on small-time scales quite well, but also that other effects not taken into account in our simulation such as supergranulation motions and the Evershed effect contribute very little. Their contribution is therefore much lower than that of the convective blueshift attenuation.

5 Conclusion

We reconstructed the integrated RVs of the Sun from MDI Dopplergrams using the Ni line at 6768 Å$\:$ over a solar cycle. We obtained a long-term amplitude of the variations of about 8 m s-1. The resulting RV variations agree very well with the simulated RV (Paper II) once scaled to the Ni line: the small-scale and long term variations are well reproduced within 30%. We therefore validate the simulations made in Paper II.

We also analyzed in detail the variations of the relative velocity (with respect to the quiet Sun) in magnetic regions, which we interpreted as an attenuation of the convective blueshift in most cases. We characterized this velocity as a function of the local magnetic field and as a function of the structure size. Our results are compatible with an average attenuation of the convective blueshift in plages and network regions by about 50%, compared to a factor 2/3 in Brandt & Solanki (1990) for plages. The difference probably arises because we included weaker field regions in our analysis.

The variations of the RV due to the convective blueshift attenuation versus the structure size could be implemented in the simulations like that made in Paper II, because the structure size is one of the input parameters defining magnetic regions. The dependence on the local magnetic field is important to understand and test our RV variations. However, if we want to extrapolate this work to other stars, this relation may not be easy to take into account, because the magnetic field distribution on other stars is not well characterized.

Acknowledgements
SOHO is a mission of international cooperation between the European Space Agency (ESA) and NASA.

References

  1. Beck, J. G., Scherrer, P. H., Bush, R. I., & Tarbell, T. D. 1998, in Structure and Dynamics of the Interior of the Sun and Sun-like Stars, ed. S. Korzennik, ESA SP, 418, 105 [Google Scholar]
  2. Brandt, P. N., & Solanki, S. K. 1990, A&A, 231, 221 [NASA ADS] [Google Scholar]
  3. Brookes, J. R., Isaak, G. R., & van der Raay, H. B. 1978, MNRAS, 185, 1 [NASA ADS] [Google Scholar]
  4. Deming, D., & Plymate, C. 1994, ApJ, 426, 382 [Google Scholar]
  5. D'Odorico, V., & the CODEX/ESPRESSO team 2007, Mem. Soc. Astron. Ital., 78, 712 [Google Scholar]
  6. Evans, S., Ulrich, R. K., Scherrer, P. H., Bush, R. I., & Tarbell, T. D. 1998, in Structure and Dynamics of the Interior of the Sun and Sun-like Stars, ed. S. Korzennik, ESA SP, 418, 157 [Google Scholar]
  7. Gabriel, A. H., Grec, G., Charra, J., et al. 1995, Sol. Phys., 162, 61 [NASA ADS] [CrossRef] [Google Scholar]
  8. Gray, D. F. 2009, ApJ, 697, 1032 [NASA ADS] [CrossRef] [Google Scholar]
  9. Henney, C. J., Ulrich, R. K., Bertello, L., et al. 1998, in Structure and Dynamics of the Interior of the Sun and Sun-like Stars, ed. S. Korzennik, ESA SP, 418, 213 [Google Scholar]
  10. Jimenez, A., Palle, P. L., Regulo, C., Roca Cortes, T., & Isaak, G. R. 1986, Adv Space Res., 6, 89 [Google Scholar]
  11. Lagrange, A. M., Desort, M., & Meunier, N. 2009, A&A, 512, A38 [Google Scholar]
  12. McMillan, R. S., Moore, T. L., Perry, M. L., & Smith, P. H. 1993, ApJ, 403, 801 [NASA ADS] [CrossRef] [Google Scholar]
  13. Meunier, N., Desort, M., & Lagrange, A. M. 2010, A&A, 512, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Miller, P., Foukal, P., & Keil, S. 1984, Sol. Phys., 92, 33 [NASA ADS] [CrossRef] [Google Scholar]
  15. Morinaga, S., Sakurai, T., Ichimoto, K., et al. 2008, A&A, 481, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129 [Google Scholar]
  17. Solanki, S. K. 1993, Space Sci. Rev., 63, 1 [NASA ADS] [CrossRef] [Google Scholar]
  18. Ulrich, R. K., Boumier, P., Robillot, J., et al. 2000, A&A, 364, 816 [NASA ADS] [Google Scholar]
  19. Wachter, R., Schou, J., & Sankarasubramanian, K. 2006, ApJ, 648, 1256 [NASA ADS] [CrossRef] [Google Scholar]
  20. Wang, H. 1988, Sol. Phys., 117, 343 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

  \begin{figure}
\par\includegraphics{14199fig1.ps}
\end{figure} Figure 1:

Upper left panel: Dopplergram for 2000 May 16 (in m s-1). Upper right panel: magnetogram for the same day (in Gauss). Lower left panel: Dopplergram after interpolation (in m s-1). Lower right panel: mask derived from the magnetogram (on the disk, 0 where the Dopplergram is interpolated).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics{14199fig2.ps}
\end{figure} Figure 2:

Limb-darkening deduced from 78 MDI intensity maps during a very quiet period.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{14199fig3.ps}
\end{figure} Figure 3:

Upper panel: zero velocity ( $V_{\rm zero}$) determined on interpolated Dopplergrams. Lower panel: velocity ( $V_{\rm obs}$) integrated over Dopplergrams.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{14199fig4.ps}
\end{figure} Figure 4:

Upper panel: estimated errors of the velocity versus time. Lower panel: distribution of the estimated errors.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics{14199fig5.ps}
\end{figure} Figure 5:

Velocity (solid line) and velocity divided by $\mu $ (dashed line) as a function of $\mu $, in m s-1, for magnetic fields larger than 100 G.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics{14199fig6.ps}
\end{figure} Figure 6:

Upper panel: $V/\mu $ versus the magnetic field, for $\mu >0.9$: for all pixels (solid line), pixels belonging to structures larger than 100 ppm (dotted line), and pixels belonging to structures smaller than 100 ppm (dashed line). The horizontal line correspond to 227 m s-1 (Paper II). Middle panel: number of pixels versus the magnetic field for the same pixel selection. Lower panel: relative contribution to the total RV.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.2cm,clip]{14199fig7.ps}
\end{figure} Figure 7:

Upper panel: velocity averaged over the structure versus the structure size for the 100 G threshold structures. Middle panel: number of structures versus the size. Lower panel: relative contribution of each category of structures to the total RV versus the size.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{14199fig8.ps}
\end{figure} Figure 8:

Upper panel: residual velocity ( $V_{\rm obs}-V_{\rm zero}$). Lower panel: simulated RV variations due to the convective blueshift attenuation (from Paper II) scaled to a $\Delta V$ of 227 m s-1.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{14199fig9.ps}
\end{figure} Figure 9:

Simulated velocity (from Paper II, scaled to 227 m s-1) versus the observed velocity, whenever Dopplergrams and simulations are separated in time by less than one hour.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{14199fig10.ps}
\end{figure} Figure 10:

Upper panel: simulated RV smoothed over 30 days (from Paper II) scaled to 227 m s-1 (dashed line) and observation (solid), whenever Dopplergrams and simulations are separated in time by less than one hour. Lower panel: ratio between the two (simulation divided by observation). The horizontal dashed line represents the median value (1.42).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{14199fig11.ps}
\end{figure} Figure 11:

Smoothed simulation (from Paper II) scaled to 123 m s-1 (dashed line) and to 166 m s-1 (dotted line), and observation (solid line), whenever Dopplergrams and simulations are separated in time by less than one hour.

Open with DEXTER
In the text


Copyright ESO 2010

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

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

Initial download of the metrics may take a while.