A&A 423, 717-728 (2004)
DOI: 10.1051/0004-6361:20047169

Kinematics of coronal mass ejections between 2 and 30 solar radii

What can be learned about forces governing the eruption?[*]

B. Vrsnak 1 - D. Ruzdjak 1 - D. Sudar 1 - N. Gopalswamy 2


1 - Hvar Observatory, Faculty of Geodesy, Kaciceva 26, 10000 Zagreb, Croatia
2 - NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA

Received 30 January 2004 / Accepted 21 April 2004

Abstract
Kinematics of more than 5000 coronal mass ejections (CMEs) measured in the distance range 2-30 solar radii is investigated. A distinct anticorrelation between the acceleration, a, and the velocity, v, is found. In the linear form, it can be represented as a=-k1(v-v0), where v0=400 km s-1, i.e., most of CMEs faster than 400 km s-1 decelerate, whereas slower ones generally accelerate.
After grouping CMEs into the width and mean-distance bins, it was found that the slope k1 depends on these two parameters: k1 is smaller for CMEs of larger width and mean-distance. Furthermore, the obtained CME subsets show distinct quadratic-form correlations, of the form $a= -k_2 (v-v_0)\vert v-v_0\vert$. The value of k2 decreases with increasing distance and width, whereas v0 increases with the distance and is systematically larger than the slow solar wind speed by 100-200 km s-1.
The acceleration-velocity relationship is interpreted as a consequence of the aerodynamic drag. The excess of v0 over the solar wind speed is explained assuming that in a certain fraction of events the propelling force is still acting in the considered distance range. In most events the inferred propelling force acceleration at 10 solar radii ranges between $a_{\rm L}=0$ and 10 m s-2, being on average smaller at larger distances. However, there are also events that show $a_{\rm L}>50$ m s-2, as well as events indicating $a_{\rm L}<0$. Implications for the interplanetary motion of CMEs are discussed, emphasizing the prediction of the 1 a.u. arrival time.

Key words: magnetohydrodynamics (MHD) - Sun: corona - Sun: coronal mass ejections (CMEs) - Sun: solar-terrestrial relations - Sun: solar wind - interplanetary medium

1 Introduction

Coronal mass ejection (CME) is a process in which the free energy contained in the coronal non-potential magnetic field is released in the form of kinetic and potential energy of ejected plasma (Forbes 2000). In large eruptions the magnetic flux of some 1023 Wb is expelled into the interplanetary space at velocities of the order of 1000 km s-1, carrying along 1013 kg of coronal plasma (Gosling 1990; Webb et al. 1994).

The non-potential field of the erupting structure carries the electric current, I, providing the Lorentz force that drives the eruption. The free energy of any current-carrying system can most generally be expressed as W=LI2, where L is the inductance of the structure (Jackson 1975; Landau et al. 1984; for a detailed discussion of the MHD aspect see Chen & Krall 2002).

The inductance depends on the size of the current-carrying structure, implying that during the eruption L increases (see, e.g., Anzer 1978; Vrsnak 1990; Chen 1996, and references therein). On the other hand, L is related to the magnetic flux $\Phi$ associated with the current, $\Phi=LI$. Since it can be taken approximately that the erupting magnetic flux $\Phi$ does not change much, this implies that the current has to decrease as the erupting structure enlarges, $I\propto L^{-1}$. Because of the decreasing current, the Lorentz force weakens, eventually ceasing at large heliocentric distances.

Note that the condition $\Phi\approx {\rm const}.$ could be violated in a limited time interval due to reconnection below the erupting flux tube (e.g., Vrsank 1990; Lin & Forbes 2000; see also Wang et al. 2003), or by an injection of the poloidal flux from the subphotospheric layers (e.g., Chen 1996). The supply of the poloidal flux prolongs the acceleration phase of the eruption.

The free energy of the system $W=LI^2\propto L^{-1}$ decreases, and if only the Lorentz force and gravity are considered, it should be directly transformed into kinetic and potential energy. However, CMEs travel through the ambient coronal/interplanetary (C/IP) plasma, causing the emission of MHD waves and appearance of the "aerodynamic'' drag. The drag can have complex characteristics depending on the orientation of the ambient magnetic field (Vandas et al. 1995, 1996; Cargill et al. 1996; Cargill & Schmidt 2002), but still can be expressed as:

 \begin{displaymath}%
a_{\rm d}= -\gamma_2 (v-w)\vert v-w\vert ,
\end{displaymath} (1)

where v is the speed of the plasmoid moving in the medium flowing at velocity w (Cargill et al. In the text 1996). The parameter $\gamma _2$ depends on the cross-section (A) and the mass (m) of the plasmoid, and on the density ( $\rho_{\rm o}$) of the ambient plasma ( $\gamma_2\propto A\rho_{\rm o}/m$).

The interaction of CMEs and the solar wind plasma is under favourable conditions "visualized'' by the dekameter-to-kilometer wavelength type II radio bursts. Such bursts are excited at the C/IP shocks that are driven by CMEs (e.g., Gopalswamy et al. 2001a).

As the Lorentz force and gravity decrease, the drag could become a dominant force. Indeed, Gopalswamy et al. (2000) inferred that in the IP space the CMEs that are faster than the solar wind decelerate, whereas the slower ones are additionally accelerated by the wind. The net acceleration was found to be anticorrelated with the CME speed; the issue was elaborated in a subsequent paper (Goplaswamy et al. 2001b) where the projection effects were partly eliminated and different distances from the Sun were considered. On the other hand, Gopalswamy et al. (2001a) analysed the relationship between acceleration and speed for fast (v > 900 km s-1) and wide (${\phi }$ > 60$^{\circ}$) CMEs associated with type II bursts, and found that in the upper corona faster CMEs show stronger deceleration. A similar result was obtained by Vrsnak (2001a) who in addition showed that the deceleration rate statistically depends not only on the speed of ejecta, but also on the height. An empirical scaling law was established, describing the decrease of the average drag with increasing radial distance, which provides a better prediction of the Sun-Earth transit time (Vrsnak & Gopalswamy 2002).

In the papers by Vrsnak (2001a) and Gopalswamy et al. (2001b) the acceleration/speed relationship (hereinafter a-v relationship) was studied on relatively small samples of events of specific characteristics. The former one embraced only fast events showing an exponential-like decay of speed, whereas the second one included only fast and wide CMEs associated with type II bursts. For this reason, Vrsnak et al. (2003) and Ruzdjak et al. (2004) extended the investigation of the a-v relationship to a large sample of more than 5000 CMEs observed by the Large Angle and Spectrometric Coronagraph (LASCO) aboard the Solar and Heliospheric Observatory (Brueckner et al. 1995) in the period 1996-2001. Again, the data showed a distinct a-v anticorrelation.

In this paper we perform a detailed analysis of the a-v relationship utilizing the same data set as in the preliminary reports by Vrsnak et al. (2003) and Ruzdjak et al. (2004). The aim of the paper is two-fold. Firstly, it is an extension of previous statistical analyses of kinematics of large samples of CMEs (e.g. Hundhausen 1993; Hundhausen et al. 1994; Harrison 1995; St. Cyr et al. 1999, 2000; Moon et al. 2002; Gopalswamy et al. 2003, and references therein). Secondly, the obtained empirical relationships are used as an insight into the physical background of CME dynamics in the near-Sun IP space. Both of these aspects could be essential in advancing the prediction of the CME arrival time at the Earth (Gopalswamy 2002).

The paper is organized as follows: the data set, the data sampling, and the parameters used are described in Sect. 2. Supplementary information concerning the statistical properties of the data set are shown in Appendix A. In Sects. 3 and 4 we present the empirical results, leaving explanations regarding the procedure for Appendix B. The results are discussed and interpreted in Sect. 5 and summarized in Sect. 6.

2 The data

The following analysis is based on the CME distance-time measurements, R(t), where $R=r/r_{\odot}$ is the plane-of-sky heliocentric distance r normalized with respect to the solar radius $r_{\odot}$. The measurements are provided in the LASCO CME catalogue (Yashiro and Michalek, http://cdaw.gsfc.nasa.gov/CME_list/). We utilize the data from 5012 events recorded in the 2<R<30 range from 1996 to 2001. From this sample we extracted two subsamples: one consists of events which had been measured in $n\ge 4$ instants, and the other with $n\ge 6$ (hereinafter n4-set and n6-set, containing 4609 and 3697 events, respectively). Such a sampling was applied to check how the accuracy of the acceleration and velocity estimates affects the results. Demanding a still larger n would certainly improve the accuracy of individual kinematic data, but it would lead to a systematic exclusion of the fastest CMEs due to the limited image acquisition rate of the LASCO instruments.

The basic parameters treated in the following are:

Statistical properties of the data set are given in Appendix A, where we show the distribution of $R_{\rm b}$, $\Delta R$, $R_{\rm m}$, ${\phi }$, $v_{\rm m}$, and $a_{\rm m}$. There is a weak positive correlation between the parameters $R_{\rm m}$, $v_{\rm m}$, and ${\phi }$ (also shown in Appendix A), whereas the acceleration is not correlated with the mean-distance or width.

In the preliminary reports by Vrsnak et al. (2003) and Ruzdjak et al. (2004) it was shown that CME accelerations and velocities are anticorrelated, indicating that the aerodynamic drag plays an important role in the CME dynamics. Since the drag acceleration depends on the ambient solar wind speed and density, as well as on the CME dimensions and density, it would be ideal if all events were recorded in the same range of heights. However, mean distances are found in the range from $R_{\rm m}=3$ to 20 (see Fig. 10 in Appendix A). Therefore, it is instructive to consider the velocity "normalized'' to a particular height so that the ambient wind speed and density are the same for all events. Since a vast majority of events were measured from $R_{\rm b}\approx 3$, we included in the analysis the velocity v3 as a kind of "initial velocity'' at a fixed (plane-of-sky) height. Other reasons for using v3 are explained in more detail and discussed in Sect. 5.1 and Appendix B.

The velocity v5 represents (statistically) the velocity at R somewhere in between the "initial distance'' R=3 and the average value of mean-distances, $\overline R_{\rm m}$, which amounts to $\overline R_{\rm m}\approx 9$ (see Appendix A). We also considered $v(R_{\rm m})$ obtained from the 2nd degree polynomial fit, but since $v(R_{\rm m})\approx v_{\rm m}$ the difference in the outcome for  $v(R_{\rm m})$ and $v_{\rm m}$ turned out to be statistically insignificant. Hereafter, v3, v5, $v_{\rm m}$, and  $v(R_{\rm m})$are denoted as "velocity options''.

An important consequence of employing v3 is that the events with $a_{\rm m}>0$ for which the fitted R(t) parabolas do not have an intercept with the R=3 level dropped out from the analysed sample. Similarly, considering v5 one excludes the events that show a large deceleration and are observed only across a relatively small range of heights, so that the fitted R(t) parabola has no intercept with R=5. Such a "silent elimination'' additionally filters out unreliable/odd data. In this way we got four subsamples which we denote n4v3-subset, n6v3-subset, n4v5-subset, and n6v5-subset, containing 4463, 3581, 4570, and 3681 events, respectively.

3 Raw a-v relationship


  \begin{figure}
\par\includegraphics[width=8.3cm,clip]{0169f1.eps}
\end{figure} Figure 1: The acceleration-speed correlation: a) 2D distribution of events in a-v space for n4vm-subset. N is the number of events in the [$\Delta a$$\times $ [$\Delta v$] = [5 m s-2$\times $ [20 km s-1] bins. b) The n4v3 data. c) n4v3 bin-averaged values $\overline a_{\rm m}$ versus $\overline v_3$ for $\Delta v=100$ km s-1 bins (beyond v>1000 km s-1, five successive bins are merged into $\Delta v=500$ km s-1 bins). Full lines in  b) and  c) show the linear least square fits, whereas 99% confidence lines are drawn dashed. d) Quadratic fit for n6v3 data (bold-dotted line) compared with the linear fit (thin-dashed line).
Open with DEXTER

In Fig. 1 we present the general relationship between CME accelerations and velocities. The four presented graphs are chosen to show basic characteristics of the a-v correlation, and various statistical aspects of the relationship. Different subsets are chosen to illustrate that the general outcome is qualitatively the same, whatever sampling and velocity options are employed.

Since it is impossible to transparently present several thousands of data points in a single graph, we first present in Fig. 1a the dependence  $a_{\rm m}(v_{\rm m})$ by showing a 2-dimensional distribution of data. An example of the "all-data'' $a_{\rm m}(v_{\rm m})$ graph, analogous to the $v_3(a_{\rm m})$ graph that is shown in Fig. 1b, can be found in the paper by Moon et al. (2002), or Ruzdjak et al. (2004).

The complete n4v3-subset data are shown in the $v_3(a_{\rm m})$ graph in Fig. 1b, and in Fig. 1c we show the corresponding bin-averaged values with the standard deviations (for mean errors $M=\sigma/\sqrt{N}$ see Vrsnak et al. 2003). The fit in Fig. 1c is obtained by attributing equal weights to all bin-averaged values in order to eliminate the effect of data grouping. The quadratic fit of the form $a_{\rm m}= -k_2 (v-v_0)\vert v-v_0\vert$ is illustrated in Fig. 1d, where the n6v3-subset is used.

Note that the a(v) graphs are presented in the most convenient form, where the acceleration is expressed in m s-2 and the velocity in km s-1. So, the slope k1 of the linear regression line in such a presentation has the unit [km-1 m s-1] = [10-3 s-1], thus the units-adjusted slope has the value $\gamma_1=10^{-3}k_1$. Analogously, the parameter k2, which appears in the quadratic fit and has the unit [m km-2] = [10-3 km-1], is related to the unit adjusted slope as $\gamma_2=10^{-3}k_2$, where $\gamma _2$ is expressed in km-1.

The same regression analysis as the one shown in Fig. 1 is performed for all other velocity options and data-subsets. The results are summarized in the upper part of Table 1, where in the first two columns the velocity option and the data-subset label are indicated. In the 3rd, and 4th column we present the unit-adjusted slope, $\gamma_1=10^{-3}k_1$ s-1, and the x-axis intercept, $v_0\equiv v_{a=0}$, of the linear least squares fit, $a=-\gamma _1(v-v_0)$. In the last column the correlation coefficients C are presented. According to the F-test the statistical significance $P_{\%}$ is larger than 99% for all correlations listed in Table 1 ( $P^*_{\%}=100-P_{\%}$ gives the probability for C=0).

Table 1: Summary of linear fits $a=-\gamma _1(v-v_0)$.

In order to check how much the accuracy of measurements influences the obtained results, we successively excluded the events with the acceleration-error Ma larger than 30, 20, 10 m s-2. It was found that the results do not depend significantly on such a subsampling.

In the lower part of Table 1 we show the results obtained using the bin-averaged values of  $\overline a$ and  $\overline v$ (Fig. 1c). The linear least squares fits presented in Table 1 exclude the highest velocity bins, 1500<v<2000 km s-1 and 2000<v<2500 km s-1, where only a small number of events is found (e.g., in the n6v3-subset these two bins contain only 21 and 1 event, respectively).

Figure 1 and Table 1 reveal a distinct anticorrelation of accelerations and speeds. The largest correlation coefficients are found for the n6v3-subset and n4v3-subset, which will be addressed in Sect. 5.1 and Appendix B.

Inspecting Table 1, one finds that the slopes are scattered in the range $\gamma_1\approx
0.01$-0.03 $\times $ 10-3 s-1. The acceleration $a_{\rm m}=0$ is found between $v_0\approx$ 350 and 450 km s-1 in most of the data subsets considered. The bin-averaged data give systematically lower values of v0 than the complete data sets ( $\overline v_0=380\pm 40$ and $\overline
v_0=430$ $\pm$ 20 km s-1, respectively).

The results shown in Table 1 can be summarized in general terms as:

\begin{displaymath}%
a_{\rm\left[m~s^{-2}\right]}=-k_1(v-v_0)_{\rm\left[km~s^{-1}\right]}= -0.02_{\pm 0.01} \left(v-400_{\pm 50}\right),
\end{displaymath} (2)

where the acceleration a is expressed in m s-2 and the velocity v in km s-1. Such an anticorrelation of accelerations and speeds is indicative of the aerodynamic drag (Cargill et al. 1996; Vrsnak 2001a; Cargill & Schmidt 2002). Indeed, the velocity $v_0\equiv v_{a=0}\approx 400$ km s-1 is relatively close to the typical solar wind speed in the considered distance range (Sheeley et al. 1997). So, roughly speaking, CMEs faster than the solar wind decelerate, whereas slower ones accelerate. However, note that v0 is in fact somewhat larger than the values measured within the LASCO field-of-view (Sheeley et al. 1997); the issue is addressed in Sect. 5.3.

Finally, we emphasize that Figs. 1a, b, and d expose a void of data-points in the range v<150 km s-1, a<0, revealing that there are practically no slow events that decelerate.

Results presented by Gopalswamy et al. (2001b) and Vrsnak (2001a) indicate that the a-v relationship might be better approximated by the quadratic dependence of the form defined by Eq. (1)[*]. The quadratic fit parameters, given in Table 2, show a much larger scatter than the linear-fit ones (to be discussed in Sect. 5.3). Statistically, the most reliable quadratic fit is obtained for the n6v3-subset which is shown in Fig. 1d. For this sample the least squares fit gives $\gamma_2=0.036$ $\times $ 10-6 km-1 and v0=500 km s-1. The values averaged over all considered data sets are $\overline \gamma_2=0.019$ $\pm$ 0.012 $\times $ 10-6 km-1 and $\overline v_0=270$ $\pm$ 170 km s-1. A more detailed treatment of quadratic fits is presented in Sect. 4.2, where it will be shown that they become statistically relevant after the data are sorted into the width and mean-distance bins.

Table 2: Summary of quadratic fits $a= -\gamma_2 (v-v_0)\vert v-v_0\vert$.

4 Dependencies ${\gamma (\phi )}$ and ${\gamma }$($R_{\rm m}$)

The a-v anticorrelation indicates that the motion of CMEs in the R=2-30 range is strongly influenced by the aerodynamic drag. If so, it can be expected that the CME width ${\phi }$ is an important parameter since the drag acceleration depends on the dimensions of the moving body. On the other hand, the drag depends on the external density and flow speed, so one can suspect that another relevant parameter is the height range in which a particular CME is measured. To investigate the role of the parameters ${\phi }$ and $R_{\rm m}$ we rely mainly on the n4-set, which is large enough to provide a further sub-grouping of the data into a number of ${\phi }$ and $R_{\rm m}$ bins. Furthermore, we limit the analysis only to $\phi<200^{\circ}$, to exclude halo CMEs, being predominantly Earth- or anti-Earth- directed, i.e., to exclude events where projection effects are large.

4.1 Simple ${\phi }$-bins


  \begin{figure}
\par\includegraphics[width=8cm,clip]{0169f2.eps}
\end{figure} Figure 2: Slopes $k_1=10^3\gamma _1$ of the a(v) linear fits for the n4v3-subset data sorted into $\Delta \overline \phi =10^{\circ }$ bins: a) k1 presented as a function of bin-averaged widths $\overline \phi $; b) k1 shown versus bin-averaged distances  $\overline R_{\rm m}$. The power-law least squares fits are presented together with the correlation coefficients C. In the insets analogous results are presented for the $\Delta \phi =20^{\circ }$ bins using n4v3-subset (dots), n4v5-subset (crosses), and n4vm-subset (triangles).
Open with DEXTER

In the first step we tried to sort the data into the $\Delta \phi=10^{\circ}$ bins, and for each bin we have drawn the a(v) graphs independently to check if, and how, the slopes $\gamma _1$depend on ${\phi }$. The procedure was applied to all velocity options and it turned out that the $\Delta \phi=10^{\circ}$ bins are generally too narrow, giving a clear $\gamma_1(\phi)$ relationship only in the case of the n4v3-subset. The outcome is shown in Fig. 2a, where the $100^{\circ}<\phi<150^{\circ}$ and $150^{\circ}<\phi<200^{\circ}$ data are merged, since for $\phi>100^{\circ}$ the $\Delta \phi=10^{\circ}$ bins are not abundant enough (see Fig. 10d in Appendix A) to provide statistically significant a-v correlations.

Figure 2a shows that the parameter $\gamma _1$ is smaller for broader CMEs, exposing a distinct power-law-like dependence. On the other hand, as the bin-averaged values  $\overline \phi $increase, the bin-averaged distances  $\overline R_{\rm m}$ also increase because wider CMEs are usually traced to larger heights, i.e., ${\phi }$ and $R_{\rm m}$ are correlated (see Fig. 11c in Appendix A). Consequently, one finds also a distinct  $\gamma_1(R_{\rm m})$ power-law-like dependence which is shown in Fig. 2b.

The decreasing trend of $\gamma_1(\phi)$ and $\gamma_1(R_{\rm m})$ can be recognized also when other velocity options and data-subsets are utilized. However, statistically significant results (P>95% confidence level) are achieved only if the bins are broadened to $\Delta \phi\ge
20^{\circ}$ bins. In the insets in Fig. 2 we show the outcome for n4v3, n4v5, and n4vm subsets, split into the $\Delta \phi =20^{\circ }$ bins. The power-law fit parameters are summarized in Table 3.

Table 3 and the insets in Fig. 2 illustrate that the fit parameters depend considerably on the velocity option - the errors of individual fit parameters are far smaller than the difference between different subsets. The power-law exponents $\beta_R$ received from the $\gamma_1(R)=\alpha_R R^{-\beta_R}$ fit vary from $\beta_R=1.5$ to 3.1, whereas the exponents of the $\gamma_1(\phi)=\alpha_\phi \phi^{-\beta_\phi}$ fit vary from $\beta_\phi=0.36$ to 0.69.

Table 3: Summary of power-law fit parameters ( $k_1=\alpha _\phi \phi ^{-\beta _\phi }$, $k_1=\alpha _R R^{-\beta _R}$) for $\Delta \phi =20^{\circ }$-bin sampling.

4.2 Decoupling the Rm${\phi }$ crosstalk


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{0169f3.eps}
\end{figure} Figure 3: The a-v correlations shown for three $R_{\rm m}$-bins of the n4v3, $\overline \phi =60^{\circ }$, subset: a) $5<R_{\rm m}<7$ ( $\overline R_{\rm m}=6.01$); b)  $7<R_{\rm m}<9$ ( $\overline R_{\rm m}=7.96$); c) $11<R_{\rm m}<15$ ( $\overline R_{\rm m}=13.25$). The number of data points is indicated (N).
Open with DEXTER

Due to the crosstalk of ${\phi }$ and $R_{\rm m}$, the presented simple-bin analysis does not really specify how $\gamma _1$ depends on ${\phi }$ at a given height, and vice versa, how $\gamma _1$ depends on $R_{\rm m}$ for CMEs of a given width class. So, in the next step we try to decouple the $\gamma
(R,\phi)$ dependence. The straightforward procedure would be to select the data within a narrow interval of $R_{\rm m}$, and then to divide the data into a number of ${\phi }$-bins. Analogously, one could try to select CMEs having widths within a certain interval, and then to separate the data into $R_{\rm m}$-bins. However, in the former case the results turned out to be statistically significant only if the $R_{\rm m}$-interval includes at least $6<R_{\rm m}<12$, and analogously in the later case, if the width-interval embraces $30^{\circ}<\phi<90^{\circ}$. Thus, although in this way the $R_{\rm m}$${\phi }$ crosstalk is reduced to a certain degree, the effect of the $R_{\rm m}(\phi )$ dependence is still present because of too wide  $\Delta R_{\rm m}$, or  $\Delta\phi$ intervals used.

So we apply another, less stringent procedure, where we start by creating $\Delta \phi =20^{\circ }$bins, and then proceed by eliminating in each bin CMEs with largest (smallest) values of $R_{\rm m}$, until $\overline R_{\rm m}$ in all ${\phi }$-bins attains a previously chosen value. This provides an insight into the ${\gamma (\phi )}$ dependence at a given height. Analogously, we group the data into $\Delta R_{\rm m}=2$ bins, and then successively eliminate wide (narrow) CMEs until all $R_{\rm m}$-bins show the same mean width  $\overline \phi $. This provides an insight into the $\gamma(R)$ dependence for a chosen mean CME width.

Figure 3 illustrates how the a-v relationship changes with increasing $\overline R_{\rm m}$. We show the n4v3-subset data in three different $\Delta R_{\rm m}=2$ bins, with the average width adjusted to the value $\overline \phi =60^{\circ }$ in each bin. The data are fitted by the linear and quadratic dependencies, and it turned out that the quadratic fits show a much larger statistical significance than in the case of complete data sets (Sect. 3) and the $R_{\rm m}$${\phi }$ coupled data (Sect. 4.1). Note that the linear, as well as the quadratic fits, become succesively less steep, whereas the x-axis intercepts v0 shift to larger velocities. A similar behaviour can be recognized if one focuses on the lower boundary of the data (the "void'' mentioned in Sect. 3)

In Fig. 4 the outcome for n4v3-subset is summarized. The dependencies $\gamma_1(\phi)$ and $\gamma_2(\phi)$ are shown in Fig. 4a, where the ${\phi }$-bin averaged value of $R_{\rm m}$ is adjusted to the value $\overline R_{\rm m}=9$. The top panel depicts the linear-fit results, $\gamma_1(\phi)$, whereas in the bottom panel we present the quadratic-fit outcome, $\gamma_2(\phi)$. Similar results are obtained for $\overline R_{\rm m}=8$. For $\overline R_{\rm m}>9$ and $\overline R_{\rm m}<8$ the statistical significances of a(v) fits in some ${\phi }$-bins become lower than P<95% due to a small number of data, so the ${\gamma (\phi )}$ dependencies become unreliable.

In Fig. 4b we show the $\gamma _{1,2}(R)$ dependence for $\overline \phi =30^{\circ }$ and $\overline \phi =60^{\circ }$. For $\overline \phi<30^{\circ}$ and $\overline \phi>60^{\circ}$ the results become unreliable due to the same reason as in the $\overline R_{\rm m}<8$ and $\overline R_{\rm m}>9$ case.

In Fig. 4c we show the x-axis intercept v0 as a function of  $\overline R_{\rm m}$, for $\overline \phi =30^{\circ }$ and $\overline \phi =60^{\circ }$, as found from the linear-fitting (top) and the quadratic-fitting (bottom). The values of v0 show a distinct increase with the radial distance. The obtained velocities are compared with the empirical solar wind model by Sheeley et al. (1997):

 \begin{displaymath}%
w(R) = w_{\infty} \sqrt{1-{\rm e}^{-(R-R_{\rm b})/R_{\rm a}}},
\end{displaymath} (3)

where $R_{\rm b}=2.8$, $R_{\rm a}=8.1$, and  $w_{\infty}$ is the asymptotic value of the wind speed for $R\rightarrow\infty$, commonly taken to be $w_{\infty }=400$ km s-1. Obviously, the values of v0(R) are systematically larger than the solar wind speed w(R), which will be discussed in Sect. 5.3.

The outcome for the n4vm-subset and the n4v5-subset is statistically less reliable than for the n4v3-subset. These options show a lower significance of a-v correlations, and consequently larger errors of power-law fit coefficients for the $\gamma _{1,2}(\phi )$ and $\gamma _{1,2}(R)$dependencies. Yet, the overall decreasing trend of $\gamma_1(\phi)$ and  $\gamma_1(R)$ is well recognizable, although not as clear as in the $R_{\rm m}$${\phi }$ coupled case (insets in Fig. 2).


  \begin{figure}
\par\includegraphics[width=17.15cm,clip]{0169f4.eps}
\end{figure} Figure 4: Summary of the a-v correlation parameters obtained for the n4v3-subset. a) Slopes $\gamma _1$ ( top) and $\gamma _2$ ( bottom) shown as a function of bin-averaged widths $\overline \phi $ for $\overline R_{\rm m}=9$; b) slopes $\gamma _1$ ( top) and $\gamma _2$ ( bottom) shown as a function of bin-averaged distances $\overline R_{\rm m}$ for $\overline \phi =30^{\circ }$ (circles and fit shown in gray) and $\overline \phi =60^{\circ }$ (crosses and fit shown in black). c) The corresponding x-axis intercepts v0 presented as a function of R and compared with the empirical solar wind model by Sheeley et al. (1997) for $w_{\infty }=400$ km s-1 (bold-gray line). In a) and b) the power-law least squares fits for $\gamma _{1,2}(\phi )$ and $\gamma _{1,2}(R)$ are shown together with the correlation coefficients C.
Open with DEXTER

5 Interpretation and discussion

The acceleration-velocity relationship analysed in Sects. 3 and 4 indicates that the aerodynamic drag plays an important role in the dynamics of CMEs in the upper corona and the near-Sun interplanetary space. In this section we interpret the obtained empirical relationships, relying on basic physical characteristics of the drag.

However, first we have to clarify some important issues considering the a-v correlation as a statistical superposition of extracted kinematical parameters of different-width CMEs measured at different mean heights. Then we proceed by discussing the qualitative aspects of results that do not depend on the procedure and the velocity option. Eventually, we compare the results with some previous studies, and discuss the motion beyond R=30.

   
5.1 Meaning of the a-v relationship


  \begin{figure}
\par\includegraphics[width=8cm,clip]{0169f5.eps}
\end{figure} Figure 5: a) A sketch showing three $a= -\gamma_2 (v-v_0)\vert v-w\vert$ curves of different $\gamma _2$ and w, and the statistical outcome for a set of nine fictive-events that were traced during a limited time interval corresponding to bold segments on the a(v)-curves. The mean acceleration and velocity of a given event is marked by circle, and the velocity interval by thin horizontal line. b) Enlarged low-velocity part of the a-v graph showing n4v3-subset data (gray), where the events with $R_{\rm m}<5$ are enhanced by black dots. Two quadratic a(v)-curves are drawn provisionally to depict the border of the "void''.
Open with DEXTER

When the dynamics of an individual CME is considered, Eq. (1) relates the instantaneous acceleration with the instantaneous velocity, whereas actual CME measurements provide only the mean acceleration estimated over a wide distance range. Furthermore, CMEs have different widths ${\phi }$, and are measured at different $R_{\rm m}$, across different $\Delta R$. This means that the interaction of CMEs with the solar wind is (quantitatively) different for various "categories'' of CMEs, as demonstrated in Sect. 4. For this reason it is instructive to pay attention to the meaning of the statistical outcome, taking into account the nature of individual measurements. Furthermore, one should bear in mind that CMEs move in various directions (different projection effects), the solar wind is highly inhomogeneous (different ambient densities and flow velocities for different CMEs), and the relative contribution of other forces might be different.

In Fig. 5a we present schematically the physical background of the statistical relationship between the mean accelerations and velocities for an imaginary set of CMEs whose motion is dominated by the drag. We have drawn three $a=a_{\rm d}= -\gamma_2 (v-w)\vert v-w\vert$ curves, hereinafter abbreviated as a(v)-trajectories. The curves are characterized by a successively larger w and smaller $\gamma _2$, representing the drag acceleration at different heights. Note that at a given ambient flow speed w, CMEs of different widths ${\phi }$ have different $\gamma _2$, which should be represented by a fan of curves crossing the x-axis at the given w. Furthermore, the presence of other forces shifts the curves up or down (see Sect. 5.3).

On the three depicted quadratic curves we have drawn the "real'' a(v)-trajectories (bold segments of the curves) of several CMEs traveling at various velocities in different solar wind environments. CMEs that are faster than the solar wind ( $v>w \Rightarrow a_{\rm d}<0$) move in the left-up direction, whereas those with v<w, i.e., $a_{\rm d}>0$, move right-down. Note that the instantaneous acceleration of CMEs is changing. The measured mean acceleration and velocity are indicated by circles, whereas the change of the velocity is indicated by thin horizontal bars.

Figure 5a explains several aspects of the statistical results obtained in Sects. 3 and 4. First of all, it clarifies why the quadratic fits turned out to be worse than linear ones when the complete data set is considered. Measuring CMEs of different dimensions and at different heights smears out the basic a(v) form into an inclined "cloud'' of data points which is naturally better fitted by the straight line.

On the other hand, the analysed CMEs were measured only above the occulting disc of the LASCO-C2 coronagraph. So, the smearing effect is reduced below some a(v)-curve that is defined by the solar wind speed w and $\gamma _2$ appropriate for $R_{\rm m}$ of CMEs measured at low heights across relatively short distance interval. Indeed, one finds a sharp-edged void at a<0, v<150 km s-1 (see Fig. 1) for any of the subsets and velocity options used. In Fig. 5b we show this part of the n4v3-subset a-v graph enlarged, where two (provisionally chosen) quadratic a(v)-curves are inserted to depict roughly the limits within which lies the border of the void (pluses and crosses represent $\gamma_2=1$ and 3 $\times $ 10-6 km-1, respectively; in both cases w=100 km s-1). Going back to Fig. 3 one finds that as $R_{\rm m}$ increases, the edge of the void successively expands towards higher velocities and becomes less steep.

   
5.2 Qualitative behaviour of the drag

Let us now focus on the results that qualitatively do not depend on the procedure (the velocity option or sampling):

i)
slopes of the a-v correlations, $\gamma _1$ and $\gamma _2$, decrease with height and are smaller for wider CMEs;
ii)
x-axis intercept of the a-v correlation ( v0=va=0) increases with height;
iii)
v0(R) is systematically larger than the solar wind speed w(R) for some 100-200 km s-1.
The parameter $\gamma _2$ in the drag acceleration, Eq. (1), depends on the CME cross-section (A) and mass (m), as well as on the ambient density ( $\rho_{\rm o}$):

 \begin{displaymath}%
\gamma_2\propto \frac{\rho_{\rm o} A}{m}
\end{displaymath} (4)

(cf. Cargill et al. 1996). Employing a characteristic CME dimension, d, one can write for a given R most generally:

\begin{displaymath}%
\gamma_2\propto \frac{\rho_{\rm o}}{\rho_{\rm i}}\frac{d^2}{d^3}
=\frac{\rho_{\rm o}}{\rho_{\rm i}d},
\end{displaymath} (5)

where $\rho_{\rm i}$ stands for the internal CME density averaged over the CME body, and we substituted $m\propto\rho_{\rm i} d^3$ and $A\propto d^2$ (see Cargill 2004 for a more specific approach where the cylindrical shape is presumed). Assuming that d is related to the angular width (e.g., $d\approx r\phi$; $\phi\approx {\rm const}.$; Bothmer & Schwenn 1994), one finds out that for a fixed distance r, $\gamma _2$ should be smaller for larger ${\phi }$. Indeed, the relationships exposed in Figs. 2c and 4a indicate such a dependence.

Here we have implicitly assumed that at a given height there is no statistical correlation between the CME density  $\rho_{\rm i}$ and d. However, wider CMEs seem to be denser: they can be traced to larger distances (see Fig. 11c in Appendix A), implying that they are brighter, meaning also denser. This effect tends to steepen the slope of the $\gamma_2(\phi)$ dependence.

Equation (4) can be rewritten also as $\gamma_2\propto \rho_{\rm o} d^2/m$. Let us assume that $m(R)\approx {\rm const}.$, i.e., that the mass loss from the prominence body and the mass pile-up at the leading edge of CME, at least partly compensate each other. Furthermore, utilizing $d\propto r$, and taking into account that in the considered distance range the ambient density decreases with the height as, say, $\rho_{\rm o}\approx r^{-4}$ (see, e.g. Leblanc et al. 1998, or Vrsnak et al. 2003) one finds roughly $\gamma_2\propto r^{-2}$. In other words, $\gamma _2$ has to decrease with the distance, consistent with Fig. 4b. At larger distances, where $\rho_{\rm o}\approx r^{-2}$, it should become approximately constant, except for the mass pile-up effect (for details see Cargill 2004).

   
5.3 Influence of other forces


  \begin{figure}
\par\includegraphics[width=8cm,clip]{0169f6.eps}
\end{figure} Figure 6: a) The parameter $v_0\equiv v_{a=0}$ obtained from n4v3, n6v3, n4vm, and n6vm subsets presented as a function of R. The bin-averaged widths are adjusted to $\overline \phi =30^{\circ }$ and $\overline \phi =60^{\circ }$. The outcome is compared with the empirical solar wind model by Sheeley et al. (1997) for the asymptotic speed of $w_{\infty }=400$ and 600 km s-1 (bold-gray and bold-black curve, respectively). b) The inferred propelling force acceleration shown in the log-log graph for the v3 (squares) and $v_{\rm m}$ (crosses) based subsets that are shown in a). The power-law fits are indicated in black and gray, respectively. The acceleration of gravity is indicated by thin-dashed line.
Open with DEXTER

The slopes of the empirical $\gamma _{1,2}(R)$ and $\gamma _{1,2}(\phi )$ dependencies obtained in Sects. 3 and 4 are quite sensitive on the sampling and velocity option. In such a situation a reliable quantitative comparison with a more detailed theoretical consideration would be exaggerated. On the other hand, the x-axis intercept v0 in the a-v correlation shows a much more coherent behaviour, so in the following we focus on this parameter.

In the absence of other forces the CME acceleration would be governed by Eq. (1), or in linear approximation by $a=\gamma_1(v-w)$. Consequently, CMEs of a given $\overline R_{\rm m}$ would form the a-v correlation that crosses the x-axis around $w(\overline R_{\rm m})$. Here we find a distinct quantitative discrepancy: Fig. 4c shows that the x-axis intercepts v0(R) are systematically higher than the corresponding solar wind speeds w(R).

To check this issue in more detail, in Fig. 6a we compare the v0(R) dependencies obtained from different data sets and velocity options, with the model values evaluated by Eq. (3). One finds out that velocities v0 closely follow the solar wind model for the asymptotic speed of $w_{\infty}\approx 600$ km s-1, which is considerably higher than the widely accepted value $w_{\infty }=400$ km s-1. Going back to Table 1 and comparing the listed values with $w\approx 250$ km s-1 evaluated from Eq. (3) for $\overline R_{\rm m}\approx 9$, one finds again that the values in Table 1 are larger for $\approx$150 km s-1.


  \begin{figure}
\par\includegraphics[width=6.6cm,clip]{0169f7.eps}
\end{figure} Figure 7: Schematic presentation showing how the a-v correlation (bold line) is shifted in the presence of a force which gives rise to the acceleration ax. The resulting dependence (gray-dashed line) intercepts x-axis at v0>w if ax>0, and at v0<w if ax<0.
Open with DEXTER

Such an excess can be interpreted in terms of a non-negligible propelling force[*]. In Fig. 7 we schematically draw how the a-v correlation is modified if the drag is not the only force acting on CMEs in the studied distance range. We consider only the simplest situation where the force (or several forces), giving rise to an acceleration ax, does not statistically depend on v. Such a force does not affect the slope of the a-v correlation - it only changes the x-axis intercept to $v_0\ne
w$. If ax<0 (like gravity), the intercept value decreases (v0<w), and if ax>0the intercept value increases to v0>w. The acceleration ax can be determined from the difference $\Delta v=v_0-w$ as $a_{x}=\gamma_1\Delta v$ (Fig. 7).

In the following we consider the influence of gravity and the propelling (presumably Lorentz) force, giving rise to accelerations g and $a_{\rm L}$. The total acceleration (force per unit mass) can be written as $a=a_{\rm d}-g+a_{\rm L}$, where $a_{\rm d}$ represents the drag acceleration. So, the additional net acceleration $a_{x}=a_{\rm L}-g$ can be estimated as $a_{\rm L}-g=\gamma_1\Delta v$. Since the values $\gamma_1(R)$ are known from observations (e.g., Fig. 4b) and  $\Delta v(R)$ can be estimated using some solar wind model (e.g., Eq. (3)), the Lorentz force acceleration can be evaluated as $a_{\rm L}=\gamma_1\Delta v+g$, where  g(R)=274/R2.

Note that by using g(R)=274/R2 we neglect the buoyancy - assuming $\rho_{\rm o}\approx\rho_{\rm i}$ one gets for the net acceleration of gravity $a_{\rm g}=0$. The value of $a_{\rm L}$ inferred by using g can be considered as an upper limit. Another drawback of the procedure lies in the fact that faster CMEs have on average larger $R_{\rm m}$. So, the slope of the a-v correlation would be altered since g(R) decreases: the influence of g is, statistically, smaller at larger velocities and the a-v correlation becomes less steep ($\gamma _1$ decreases). Analogously if  $a_{\rm L}(R)$ decreases, the a-v correlation becomes steeper. In principle, it would be possible to perform an iterative procedure to estimate more precisely  $a_{\rm L}(R)$, but bearing in mind the accuracy of the observational input and other drawbacks of the applied procedure (see Appendix B), such a treatment would be exaggerated.

For the matter of illustration, let us first take $\gamma_1\approx 0.02$ $\times $ 10-3 s-1 and v0=350-400 km s-1 ("average'' values from Table 1). Using $w\approx 250$ km s-1 at $\overline R_{\rm m}=9$ as obtained from Eq. (3), one finds $\Delta
v=v_0-w=100$-150 km s-1. Utilizing these values one receives $a_{\rm L}-g=\gamma_1\Delta v=
2$-3 m s-2, i.e., $a_{\rm L}\approx 6$ m s-2, since gR=9=3.4 m s-2.

Let us now evaluate $a_{\rm L}(\overline R_{\rm m})=\gamma_{1}\Delta v + g$ using the differences $\Delta v(\overline R_{\rm m})=v_{0}-w(\overline R_{\rm m})$ and slopes $\gamma_1(\overline R_{\rm m})$ for all individual R-bins considered in Sects. 3. and 4.

The outcome is shown in Fig. 6b for all subsets that are utilized in Fig. 6a. Bearing in mind all drawbacks of the procedure, one can say that the propelling-force acceleration, i.e., the propelling force per unit mass, at $R\approx 10$ amounts to $a_{\rm L}\approx$ 5-10 m s $^{-2}\equiv 5$-10 N kg-1. The data show a decreasing trend which can be roughly expressed as $\approx$R-1. The values based on the velocity v3 are larger than those obtained by using $v_{\rm m}$. This is primarily a result of the steeper slope $\gamma _1$ and somewhat larger v0 of the a-v3 correlation (Sect. 5.1).

Let us note that using the isothermal (T=106 K) solar wind model by Parker (1958; see also Mann et al. 2003) instead of Eq. (3), one gets a 30-50% larger ax(resulting in 10-30% larger $a_{\rm L}$), but the slope of ax(R) remains approximately R-1.

We emphasize the statistical nature of Fig. 6b, i.e., it does not represent an average behaviour of the propelling force  $a_{\rm L}(R)$ in a "typical'' CME. It rather shows that the propelling force is on average weaker at larger distances, being consistent with the anticipated behaviour of the Lorentz force (Sect. 1).

There is also another way to estimate the quantity $a_{\rm L}=a+g-a_{\rm d}$. The drag acceleration for a given CME can be evaluated from $a_{\rm d}\approx -\gamma_1(v-w)$ by substituting the velocity $v=v_{\rm m}$, v3, or v5, and the model-value $w=w(R_{\rm m})$. For estimating $\gamma _1$ one can take the empirical scaling law $\gamma_1=1.16~R^{-1.35}\times10^{-3}$ s-1 found by a different procedure by Vrsnak (2001a). Finally, using the measured values $a_{\rm m}$ and $g(R_{\rm m})=274/R_{\rm m}^2$ one can estimate $a_{\rm L}$ for each event individually.

The calculated values $a_{\rm d}$ are shown in Fig. 8a as a function of v3 (a similar outcome is also found using $v_{\rm m}$ or v5). The overall anticorrelation pattern in Fig. 8a is similar to that exposed by Fig. 1, justifying the proposed interpretation in terms of the drag force, and on the other hand confirming reliability of the empirical scaling law $\gamma(R)$obtained by Vrsnak (2001a). The slopes of the linear fits in Figs. 8a and 1 are similar, whereas there is a difference in the x-axis intercept v0. The linear fit in Fig. 8a crosses the x-axis at $v_0\approx 250$ km s-1, i.e., it is lower for $\Delta v_0\approx 150$ km s-1. Again, this can be attributed to the existence of ax>0. Taking for the slope of the $a_{\rm d}(v)$ correlation $\gamma_1=-0.026$ $\times $ 10-3 s-1 (see the fit in Fig. 8a) one finds $a_{x}=\gamma_1\Delta v_0=+3.9$ m s-2, and assuming $a_{x}=a_{\rm L}-g$ one receives $a_{\rm L}(R_{\rm m})=7.3$ m s-2, consistent with the previous estimates.

Another difference found in Figs. 1 and 8a is a considerably larger scatter of the data in Fig. 1. Firstly, this can be attributed to a variety of the propelling force strengths. Furthermore, the differences in CME widths are not included in the evaluation of $a_{\rm d}$ since the used scaling law by Vrsnak (2001a) does not provide that. Finally, note that in the low-velocity regime, v<150 km s-1, the accelerations appear in Fig. 8a down to $a_{\rm d}\approx -10$ m s-2, whereas in Fig. 1 (see also Fig. 5) the upper edge of the void is at a=0. Such a difference is also consistent with the previous estimates of the propelling force.

In Fig. 8b we compare the distributions $N(a_{\rm d})$ and $N(a_{\rm L})$ for the n6v3-subset (see also the the distribution of measured accelerations, $N(a_{\rm m})$ shown in Fig. 10f in Appendix A). In addition, we present in the inset the distribution N(g), where $g=g(R_{\rm m})$. Mean values are: $\overline a_{\rm d}= -5.4$ $\pm$ 9.7 m s-2, $\overline a_{\rm L}= 10.0$ $\pm$ 13.9 m s-2, and $\overline g= -5.2$ $\pm$ 4.6 m s-2, which explains why CMEs in the upper corona and near-Sun interplanetary space show on average $\overline a_{\rm m}\approx 0$(Appendix A). We emphasize that very similar distributions are obtained utilizing any other subset or velocity option. The drag acceleration is dominantly negative, i.e., most events are faster than the solar wind. On the other hand, in the majority of events, $a_{\rm L}$ lies between 0 and 20 m s-2, consistent with Fig. 6b. The distribution indicates that there is a considerable fraction (18%) of events where $a_{\rm L}$ appears to be larger than 20 m s-2 (see also Vrsnak 2001b). In extreme cases (2%) $a_{\rm L}$ is larger than 50 m s-2 (to be studied in detail in a separate paper). In 17% of events the inferred values of $a_{\rm L}$ are negative, with 5% showing $a_{\rm L}<-10$ m s-2. Note that the $a_{\rm L}<0$ regime is revealed by the "fall-back'' events described by Wang & Sheeley (2002) and by the existence of an upper equilibrium observed in some events (Vrsnak et al. 1990; Vrsnak 2001b). Such a behaviour is anticipated in the flux rope model by Vrsnak (1990).

Finally, let us emphasize one potentially important side-result of the presented analysis. None of the presented graphs (including those in Appendix A) shows any kind of grouping of CMEs that would confirm the existence of two (or more) distinct classes of events regarding the velocity and/or acceleration behaviour, as proposed by some authors who studied small samples of preferentially large CMEs (see, e.g., MacQueen & Fisher 1983, or Andrews & Howard 2001). Rather, the distribution of CMEs in the a-v space indicates that there is a continuum of events over the actual acceleration and velocity range.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{0169f8.eps}
\end{figure} Figure 8: a) Drag accelerations estimated using $v_{\rm m}$ for n6v3-subset, shown in the ad-v3 space. The corresponding linear least squares fit is indicated. b) The n6v3-subset distribution of estimated drag accelerations (black) and the corresponding inferred propelling force accelerations $a_{\rm L}$ (gray). The bin width is 5 m s-2. 23 (0.6%) events in $a_{\rm d}$ distribution and 47 (1.3%) in $a_{\rm L}$ distribution are out of the scale. In the inset the distribution of the acceleration of gravity g is shown for bins of 1 m s-2.
Open with DEXTER

  
5.4 Comparison with previous studies

Values of the exponent $\beta_R$ in the $\gamma_1\propto R^{-\beta_R}$ scaling received in Sect. 4.1, can be directly compared with the result by Vrsnak (2001a) since the $R_{\rm m}$${\phi }$ crosstalk was implicitly present therein too. The obtained values $\beta_R=1.5$-3.1 are larger than $\beta_R=1.35$ found by Vrsnak (2001a), indicating a considerably steeper decrease of $\gamma _1$. However, it should be noted that the results by Vrsnak (2001a) were based on a limited set of only 12 events, all of which showed an exponential-like decay of velocity, indicating that the drag is a dominant force. Furthermore, seven of these events were measured at low heights by tracing H$\alpha$ features, i.e., the deceleration rates could be higher for leading edges of CMEs.

In Fig. 9 we compare the slopes of the a-v correlations obtained herein with the scaling law $\gamma_1(R)$ established by Vrsnak (2001a), and with the results for the interplanetary space obtained by Gopalswamy et al. (2000) and Gopalswamy et al. (2001b). The values for the complete n4v3-subset and n6v3-subset are positioned at $R=\overline R_{\rm m} = 9$ and overlap each other. The slopes from Gopalswamy et al. (2000) and Gopalswamy et al. (2001b) are drawn at R=40 and R=46, respectively, corresponding to the geometrical mean (see the argumentation in Vrsnak 2001a) of the considered distance ranges (0.76 and 1 au, respectively). The fit through the presented data-points matches quite closely to the scaling k1 = 2 R-1.5 that was inferred by Vrsnak & Gopalswamy (2002) from the Sun-Earth transit times.


  \begin{figure}
\par\includegraphics[width=7.6cm,clip]{0169f9.eps}
\end{figure} Figure 9: The slopes $k_1=10^3\gamma _1$ of the linear a-v correlations, presented as a function of radial distance: circle - n6v3-subset; triangle - Goplaswamy et al. (2001b); square - Gopalswamy et al. (2000); crosses connected by the dotted line - Vrsnak (2001a). The power-law fit to these data-points reads k1=1.2 R-1.5 and is drawn by the thin-black line. The gray dashed line shows the scaling k1 = 2 R-1.5 inferred by Vrsnak & Gopalswamy (2002). The bold-black and bold-gray lines show power-law fits from Fig. 2b and from Fig. 4b-top for $\overline \phi =60^{\circ }$, respectively.
Open with DEXTER

5.5 Motion of CMEs beyond R=30

Figure 8b indicates that in the R=2-30 range the majority of events has already more or less adjusted to the solar wind speed since the $N(a_{\rm d})$ distribution (Fig. 8) is sharply peaked at $a_{\rm d}=0$. A similar conclusion can be drawn from Fig. 1a, since most events are grouped around a=0, $v\approx 400$ km s-1.

Let us now inspect the distribution of speeds $v_{\rm e}$ estimated at the time of last measurement (see Sect. 2) for the events that were traced to $R_{\rm e}\ge25$. We find that only 374 out of 4463 n4v3-events (8%) had $v_{\rm e}>500$ km s-1. The number of events with velocities above some referent value, $v_{\rm e}>v^*$, decreases as we increase the limit v*. The events with $v_{\rm e}>v^*$ are on average much broader than the rest of the sample. Furthermore, they show a positive mean acceleration, being statistically larger for faster (and broader) CMEs. The n4v3-subset statistics is summarized in Table 4 for different values of v*. The number (N) and the percentage (%) of the v>v* events is presented in the 2nd and 3rd column. The values in brackets provide the number (and relative percentage) of $a_{\rm m}>0$ events. The average values of widths and accelerations are given in the last two columns. Note that more than half of the $v_{\rm e}>v^*$ events are characterized by $a_{\rm m}>0$. In some cases it is larger than $a_{\rm m}=50$ m s-2.

Table 4: Statistics of CMEs faster than solar wind at $R\ge 25$.

Such a behaviour has important implications for the prediction of the 1 au transit. In most CMEs the propelling force ceases, and the velocity becomes adjusted to the solar wind speed, already below R=30. The arrival of these events to the Earth (if they preserve their integrity) can be estimated to a sufficient accuracy utilizing the $v={\rm const}.=400$ km s-1 approximation. However, for very fast/wide CMEs, which in fact are the most geo-effective events (Zhang et al. 2003), the prediction based on LASCO observations is more difficult since the action of the propelling force is prolonged in such events. To be successful in predicting the 1 au transit, a higher degree polynomial fit to the R(t) LASCO-data should be tried, to estimate if possible, how the acceleration changes. This might give a clue whether the event had entered the drag-dominated regime or not. If yes, the interplanetary motion of the CME should be calculated using the velocity at R=30 as the initial velocity and utilizing the scaling  $\gamma_\phi(R)$ for the corresponding ${\phi }$-bin. Then, the equation of motion given by Eq. (1) should be solved numerically, in the manner demonstrated by Vrsnak & Gopalswamy (2002). Note that events of this kind should have the smallest deceleration rates since they are broad, which implies a relatively small $\gamma _2$.

The prolonged action of the propelling force in fast/wide CMEs is consistent with the observations based on the interplanetary scintillation (IPS). For example, Manoharan et al. (2001) have analysed the IP propagation of the CME associated with the large flare of 14 July 2000. Until $R\approx 100$ the speed of ejection was almost constant ( $v\propto R^{-0.08}$), and only beyond this distance it started decelerating more rapidly, as $v\propto R^{-0.72}$. In such cases, the Thompson scattering measurements by the Solar Mass Ejection Imager (SMEI) combined with IPS observations (Jackson et al. 2003; Eyles et al. 2003) could be essential for the understanding of the IP propagation of CMEs and a successful prediction of 1 au transit time.

6 Conclusion

We have shown that there is a well defined anticorrelation between the acceleration and velocity of CMEs measured in the range 2<R<30. The relationship can be comprehended as the statistical outcome for a sample of CMEs whose motion is significantly influenced by the aerodynamic drag: Only 14% of fast (v3>1000 km s-1) CMEs are still speeding up, and 7% of slow CMEs (v3<150 km s-1) show (only a weak) deceleration.

However, the existence of fast events whose velocity is still increasing suggests that in a certain fraction of events the propelling force is not negligible. On the other hand, the deceleration of slow events, and existence of the $a_{\rm L}<0$ events in the $N(a_{\rm L})$ distribution (Fig. 8b) indicates that the deceleration in these events could be caused not only by the gravity, but also by the Lorentz force (see also Wang & Sheeley 2002), which under special conditions might change its sign after passing an upper equilibrium point (Vrsnak 1990).

The overall behaviour of the parameter $\gamma _2$, decreasing with the distance and being of the order of 0.1 $\times $ 10-6 km-1 at $R\approx 10$, is very similar to that obtained by Vrsnak (2001a). In addition, we found herein that the value of $\gamma _2$ is smaller for broader CMEs when measured at the same $R_{\rm m}$. On the other hand, the slope of the linear a-v correlation at $R_{\rm m}\approx 10$ ranges between $\gamma_1=0.01$-0.03 $\times $ 10-3 s-1, which is smaller by roughly a factor of two than found by Vrsnak (2001a).

The results indicate that in the majority of events the $v\approx w$ regime is achieved already in the LASCO field-of-view, i.e., below R<30. Only CMEs with prolonged action of the propelling force come into the $v\approx w$ regime beyond R>30, but the events with, e.g., v>500 km s-1 at R>25 represent less than 10% of the CME population.

Acknowledgements
We are grateful to the LASCO-SOHO team for operating the instruments and performing the basic data reduction, and especially, for the open data policy.

References

 

  
Online Material

Appendix A: Properties of the sample


  \begin{figure}
\par\includegraphics[width=11.25cm,clip]{0169f10.eps}
\end{figure} Figure 10: Distribution of CME parameters in the n4v3-subset: a) initial radial distances $R_{\rm b}$ (bins $\Delta R=1$); b) radial distance range $\Delta R=R_{\rm e}-R_{\rm b}$; c) mean radial distances $R_{\rm m}=(R_{\rm e}+R_{\rm b})/2$; d) angular widths ${\phi }$ (bins $\Delta \phi =20^{\circ }$); e) mean velocities $v_{\rm m}$ (bins $\Delta v_{\rm m}=100$ km s-1); f) mean accelerations $a_{\rm m}$ (bins $\Delta a=5$ m s-2); in the inset the central part of the distribution is presented (bins $\Delta a=2$ m s-2).
Open with DEXTER

The events which were measured at $n\ge 4$ instants, show the distribution centered at $a_{\rm m}\approx 0$, with the standard deviation  $\sigma_{\rm a}=31$ m s-2. The average initial, mean, and final distances amount to $\overline R_{\rm b}=3.4$ $\pm$ 1.0, $\overline R_{\rm m}=9$ $\pm$ 4, and $\overline R_{\rm e}=14$ $\pm$ 7, respectively.

In Fig. 10 we present the distribution of initial radial distances $R_{\rm b}$, the distance range $\Delta R=R_{\rm e}-R_{\rm b}$, mean radial distances $R_{\rm m}$, widths ${\phi }$, mean velocities $v_{\rm m}$, and mean accelerations $a_{\rm m}$, for the n4v3-subset. Note that most often the events start at $3\le R<4$. The n4v3-subset contains 280 CME events (6%) of the angular width larger than  $180^{\circ}$, out of which 184 were classified as $360^{\circ}$-halo CMEs (4%). The basic statistical parameters of the distributions are given in Table 5.

The average error of the accelerations obtained from the 2nd degree polynomial fits is $\overline
M_{\rm a}=7.6$ m s-2 in the n4-set, and it drops to 4.7 m s-2 in the n6-set. The fraction of events with errors $\overline M_{\rm a}>10$ m s-2 is 18% and 10% for the two data-sets, respectively.

The parameters $R_{\rm m}$, ${\phi }$, and $v_{\rm m}$ are mutually correlated (see also, e.g., Hundhausen et al. 1994), which is shown in Fig. 11, where the data are presented together with the linear least squares fits and corresponding correlation coefficients (C). We emphasize that the $\phi(R_{\rm m})$ correlation does not represent the broadening of CMEs with the height, but rather demonstrates that larger CMEs can be traced to larger distances, i.e., that broader CMEs are generally brighter. Similar holds also for the $v_{\rm m}(R_{\rm m})$ relationship - faster CMEs are on average brighter and can be followed to larger heights. No statistically significant correlation between $a_{\rm m}$ and $R_{\rm m}$, or ${\phi }$ is found (C<0.01).

The basic characteristics shown in Figs. 10 and 11 are very similar for all of the considered subsets and velocity options.

Appendix B: CME motion

The main problem included implicitly in the performed analysis is that Eq. (1) relates the instantaneous acceleration, a, with the instantaneous velocity, v, whereas the measurements provide only the mean acceleration $a_{\rm m}$ estimated over a wide distance range, $\Delta R$. Consequently, it would be erroneous to expect ad hoc that the measured mean acceleration is simply related to the mean speed, even if the drag is the dominant force in the considered distance range. Since this issue can cause a systematic error, we have to clarify its influence.

In Fig. 12a the drag-dominated a(v)-trajectory of one a<0 event is sketched in details. The velocity of CME at the beginning of measurement is denoted at the x-axis by "b'' and at the end by "e''. The mean velocity $v_{\rm m}$, defined in the $a={\rm const}.$ approximation as $v_{\rm m}=(v_{\rm b}+v_{\rm e})/2$, is marked by "m''. The mean acceleration, in the $a={\rm const}.$ approximation equivalent to $a_{\rm m}=(v_{\rm e}-v_{\rm b})/(t_{\rm e}-t_{\rm b})$, is depicted by the horizontal dashed line. If the solar wind speed and the parameter $\gamma _2$ are taken in the first approximation to be constant, the CME-trajectory (bold-dotted arrow) would follow the a(v) line defined by $w_{\rm b}$ and  $\gamma_{2{\rm b}}$. However, as the CME advances, the ambient solar wind speed increases whereas $\gamma _2$ decreases (Vrsnak 2001; see also Sects. 4.2 and 5.2). Consequently, the CME "slides'' across a(v)-curves, ending at the one defined by $w_{\rm e}$ and $\gamma_{2{\rm e}}$, as sketched by the bold arrow. Note that the line depicting the mean acceleration intersects the bold arrow at the velocity v different from $v_{\rm m}$. Furthermore, note that the intercept of the $a_{\rm m}$ and $v_{\rm m}$ line is below both depicted a(v)-curves, i.e., it lies on a curve of too large $\gamma _2$, corresponding to some distance smaller than $R_{\rm b}$. On the other hand, the intercept of the $a_{\rm m}$ line with $v_{\rm b}$ is located between the "initial'' and "final'' a(v)-curve, i.e., it is associated with $\gamma _2$corresponding to some distance between $R_{\rm b}$ and $R_{\rm e}$, presumably close to $R_{\rm m}$.

Table 5: Statistical characteristics of the n4v3-subset and (n6v3-subset).

If we want to inspect quantitatively the effect sketched in Fig. 12a, it is instructive to consider firstly what mean acceleration would be measured if an object would move according to Eq. (1) in a homogeneous medium ( $w={\rm const}.$, $\gamma _2={\rm const}.$; depicted by the bold-dashed arrow in Fig. 12a). The mean acceleration of the object, evaluated from the 2nd order polynomial fit to the distance-time data can be expressed, by definition, as:

 \begin{displaymath}%
a_{\rm m}=\frac{(v_{\rm e}-v_{\rm b})}{(t_{\rm e}-t_{\rm b}...
...rac{(\hat v_{\rm e} - \hat v_{\rm b})}{(t_{\rm e}-t_{\rm b})},
\end{displaymath} (6)

where $v_{\rm b}$ and $v_{\rm e}$ denote the rest-frame velocity at the beginning ($t_{\rm b}$) and at the end ($t_{\rm e}$) of the measurements, whereas $\hat v_{\rm b}$ and  $\hat v_{\rm e}$ represent the corresponding velocities relative to the ambient gas ( $\hat v=v-w=v-{\rm const}.$).


  \begin{figure}
\par\includegraphics[width=7cm,clip]{0169f11.eps}
\end{figure} Figure 11: Relation between different CME parameters in the n4v3-subset: a) $v_{\rm m}(R_{\rm m})$ - black dots and the dashed linear least squares fit represent $\phi =360^{\circ }$ events. The $\phi <360^{\circ }$ events are depicted by gray circles and full-line least squares fit. b) $v_{\rm m}(\phi )$; c) $R_{\rm m}(\phi )$. In b) and c) the $\phi =360^{\circ }$ events (crosses) are excluded from the fit.
Open with DEXTER

The solution of Eq. (1) in the case $v_{\rm b}>w$ reads:

 \begin{displaymath}%
\hat v_{\rm e} = \hat v_{\rm b} {\rm e}^{-\gamma_2 \Delta \hat r},
\end{displaymath} (7)

where $\Delta \hat r=\hat r_{\rm e} -\hat r_{\rm b}$ is the distance traveled relative to the ambient gas in the time $\Delta t=t_{\rm e}-t_{\rm b}$. Integrating d$\hat r$ = $\hat v$ dt one finds another form of the solution:

 \begin{displaymath}%
\hat v_{\rm e} = \frac{\hat v_{\rm b}}{1+\hat v_{\rm b}\gamma_2\Delta t},
\end{displaymath} (8)

(see also Cargill et al. 1996), providing:

 \begin{displaymath}%
\Delta t =\frac{{\rm e}^{\gamma_2 \Delta \hat r}-1}{\gamma_2 V_{\rm b}},
\end{displaymath} (9)

Substituting Eqs. (9) and (7) into Eq. (6) one finds:

 \begin{displaymath}%
a_{\rm m} = -\gamma_2 \hat v_{\rm b} \hat v_{\rm e} = -\gamma_2 (v_{\rm b}-w)(v_{\rm e}-w),
\end{displaymath} (10)

which is different from the ad hoc expected $\gamma_2 (v_{\rm m}-w)(v_{\rm m}-w)$. Using Eq. (7) one can rewrite Eq. (10) as:

 \begin{displaymath}%
a_{\rm m}=\gamma_2 (v_{\rm b}-w)^2 {\rm e}^{-\gamma_2 \Delta \hat r} =\gamma_2^*(v_{\rm b}-w)^2,
\end{displaymath} (11)

where we abbreviated $\gamma_2^*=\gamma_2{\rm e}^{-\gamma_2 \Delta \hat r}$. So, if one would take the measured values $a_{\rm m}$ and $v_{\rm b}$ (or equivalently $v_{\rm b}\approx v_3$ as done in Sect. 4) to estimate $\gamma_2^*$ utilizing Eq. (11), its value would be smaller than the real value $\gamma _2$ by the factor $X={\rm e}^{-\gamma_2 \Delta\hat r}$. Using Eq. (8), the deviation X can be expressed as:

 \begin{displaymath}%
X =\left(\gamma_2 \hat v_{\rm b} \Delta t +1\right)^{-1}.
\end{displaymath} (12)

For a given CME the deviation X depends on the initial velocity and the time interval during which the kinematics was measured, and this varies from one event to another.

Actually, if $\Delta t$ increases, $R_{\rm m}=R_{\rm b}+\Delta R/2$ increases too, which means that $\gamma _2$ decreases (Sect. 4). Thus, the two effects partly compensate each other. In Fig. 12b we show the deviation X as a function of $R_{\rm m}$. In evaluating $X(R_{\rm m})$ we have used $w=w(R_{\rm m})$ as given by Eq. (3), and $\gamma_2=\gamma_2(R_{\rm m})$ estimated utilizing the two scalings displayed in the bottom panel of Fig. 4b ( $\gamma _2=5.9$ $\times $ 10-6R-2 and $\gamma _2=0.76$ $\times $ 10-6R-0.98).

Note that in the considered approximation, in most of events ( $4<R_{\rm m}<15$; see Appendix A) the value of  $\gamma_2^*$ is underestimated for not more than 50% with respect to $\gamma _2$.


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{0169f12.eps}
\end{figure} Figure 12: a) Comparison of the a(v)-trajectory of an event moving in the $w={\rm const}.$ environment (bold-dotted arrow) with the one moving in the solar wind where the density decreases and flow speed w increases (bold arrow). The initial, mean and final velocities are denoted by b, m, and e, respectively. b) The deviation X for the approximation $\gamma _2={\rm const}.$, $w={\rm const}.$ c) The deviation X for the solar wind speed given by Eq. (3) combined with $\gamma _2=0.76$ $\times $ 10-6R-0.98 (black) and $\gamma _2=5.9$ $\times $ 10-6R-2 (gray). The successively thicker lines in  b) and  c) represent v3=1000, v3=1500, and v3=2000 km s-1.
Open with DEXTER

Let us now consider the more realistic situation, where the wind speed and density, as well as the dimensions and density of CME, depend on the radial distance, so w and $\gamma _2$ are functions of R. We presume that w(R) increases, whereas  $\gamma_2(R)$ decreases.

We performed a series of numerical simulations from which we estimated $a_{\rm m}$ to determine how much the approximate value $\gamma_2^*(R_{\rm m})=a_{\rm m}/(v_{\rm b}-w(R_{\rm m}))^2$ differs from the true $\gamma_2(R_{\rm m})$. Again we used the two $\gamma_2(R)$ scalings from the bottom panel of Fig. 4b, combined with the solar wind model given by Eq. (3). The results are illustrated by Fig. 12c, revealing a complex situation. The deviation X is generally X<1, but it can also be X>1, i.e., $\gamma_2^*$ can be overestimated (see also the explanation of Fig. 12a). The later case is more pronounced for larger values of the exponent $\beta_R$ in the $\gamma(R)\propto R^{-\beta_R}$ scaling. When $v_{\rm m}$ is considered instead of $v_{\rm b}$, the value of  $\gamma_2^*$ is always overestimated. For this reason in the previous analysis (Sects. 3 and 4) we relied more on the velocity v3 than on $v_{\rm m}$. Note that weak-acceleration events are not affected much by the choice of v3 or $v_{\rm m}$.

Finally, we stress another important aspect of employing the velocity v3. Considering what happens in the a-v graph when the $a_{\rm m}(v_{\rm m})$ data are replaced by  $a_{\rm m}(v_3)$, one finds that the usage of v3 shifts the $a_{\rm m}>0$ events to lower velocities (to the left in the graph), and  $a_{\rm m}<0$ events to higher velocities (to the right). This would tend to decrease the slope of the a-v correlation if only the drag force dominated events would be present (see the explanation of Fig. 12a). However, a compensating effect appears at large velocities where we find also the events with large $a_{\rm m}>0$ (presumably those whose motion is still governed by the propelling force). These events, which tend to smear the correlation and decrease $\gamma_{1,2}$, are shifted to the left towards v0, which reduces their influence in the a(v) fitting. This is a reason reason why the $a_{\rm m}(v_3)$ relationship shows higher correlation coefficients than  $a_{\rm m}(v_{\rm m})$ relationship and a steeper slope $\gamma _1$.


  \begin{figure}
\par\includegraphics[width=6.9cm,clip]{0169f13.eps}
\end{figure} Figure 13: The fitted $\gamma _2(R,\phi )$ surface obtained by using Eq. (13). In the inset log-scale presentation is shown.
Open with DEXTER

In conclusion, let us note that Eq. (10) provides another possibility to estimate $\gamma_2(\phi,R)$ scaling, different from that performed in Sect. 4. Combining the observed values $a_{\rm m}$, $v_{\rm b}$, and $v_{\rm e}$, with the model values for the solar wind speed  $w(R_{\rm b})$ and  $w(R_{\rm e})$ as given by Eq. (3), one can estimate:

 \begin{displaymath}%
\tilde\gamma_2 = - \frac{a_{\rm m}}{(v_{\rm b}-w_{\rm b})\vert(v_{\rm e}-w_{\rm e})\vert}
\end{displaymath} (13)

for each event. Here we meet a problem if the $v_{\rm b}\approx w_{\rm b}$ events are included - because of the small value in the denominator the outcome becomes highly unreliable. So we excluded the $\vert v_{\rm b}-w_{\rm b}\vert<50$ km s-1 events. The remaining data are then fitted in the $\tilde\gamma_2(R,\phi)$ space using the function $\tilde\gamma_2(R,\phi)=\alpha~\phi^{-\beta_\phi}R^{-\beta_R}$ (the plane in log-log-log presentation). The outcome is shown in Fig. 13: the parameters obtained by the linear least squares fit read $\alpha=0.015$ $\times $ 10-3 km-1, $\beta_\phi=0.42$, and $\beta_R=0.95$. So, the value of  $\tilde\gamma_2$ decreases with R and ${\phi }$. Note that the slopes  $\beta_\phi$ and $\beta_R$ are consistent with those shown in Figs. 4a and 4b.



Copyright ESO 2004