Free Access
Issue
A&A
Volume 635, March 2020
Article Number A112
Number of page(s) 10
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/201936016
Published online 17 March 2020

© ESO 2020

1. Introduction

Coronal mass ejections (CMEs) were discovered in the 1970s by the observations of the white-light coronograph experiment on board the OSO-7 mission (Tousey et al. 1973). Since then, a large number of data have been accumulated that help us understand the basic characteristics of CMEs, but many questions still remain unanswered (see e.g., Gopalswamy et al. 2004; Webb & Howard 2012; Gopalswamy 2016, and references therein). In particular, the Large Angle and Spectrometric Coronograph (LASCO, Brueckner et al. 1995) on board the Solar and Heliospheric Observatory (SOHO) spacecraft has been very successful in contributing considerably to this understanding of CME characteristics. In this work, we use the high number of LASCO observations and the CME characteristics, measured and made available online in the CDAW database (Gopalswamy et al. 2009), to explore the statistical properties of CMEs. In particular, we are interested in consecutive CMEs, which may help in the understanding of the energy storage and its release via a triggering mechanism, and the early stages of the CME dynamics. Coronal mass ejections and flares consume only part of the free magnetic energy available in active regions (Feng et al. 2013). The released energy (considering the flare and CME energy) may be, on average, a third of the free magnetic energy (Emslie et al. 2012). Furthermore, CMEs carry on average 7% of the total dissipated energy during the eruptive event (Aschwanden et al. 2017). Therefore, in terms of the stored energy, more than one CME can be ejected from the same active region (Gopalswamy et al. 2005; Michalek & Yashiro 2013; Gopalswamy 2017), and more importantly, within short waiting times (WT), which is the time elapsed between the launch of two successive CMEs (see, Sect. 3). The triggering mechanisms remain unsolved (Forbes 2000). Recently, it has been proposed that CMEs may destabilize nearby active regions and trigger “sympathetic” CMEs (see, Lugaz et al. 2017 and references therein), in such cases short WT are expected, although, to support this causality, the WT distribution should be different from one of a Poisson random process. Therefore, the WT distribution gives valuable information about the so-called “memory” of the system. For memory-less systems (stochastic process) the distribution is exponential or gamma. On the other hand, the so called “fat-tailed” distributions (e.g., Power Law, Pareto, Weibull, etc.) are often associated with memory systems (see, Barabási 2005; Samorodnitsky 2007; Sheng et al. 2011, for a general discussion of this subject). As evidenced by coronal EUV dimmings, the footprints of CMEs (Thompson et al. 2000), the corona behind CMEs is “evacuated” (Krista & Reinard 2017) during a few hours (López et al. 2017). Therefore, a second CME launched in the same region within a short WT encounters a different coronal medium than the first CME, and its dynamics are different according to the drag force models of CME propagation (Cargill 2004; Borgazzi et al. 2009; Vršnak et al. 2010). If the second CME is faster than the first one, and the source region is similar, there is a high probability of interaction at any heliospheric distance. This interaction causes CME cannibalism (Gopalswamy et al. 2001), and the merged regions have a different dynamics than a single CME (Niembro et al. 2015), making the prediction of the travel time and arrival speed at 1 AU difficult (Gopalswamy et al. 2013). The paper is organized as follows: The SOHO/LASCO data used in this work is presented in Sect. 2. We determine the waiting time (WT) between consecutive CMEs, and present statistical characteristics in Sect. 3. The WT has a strong dependence on the solar cycle phase, which we describe in Sect. 4. We study the source region of consecutive CMEs, analyzing the statistical behavior of the difference of position angle (PA) (Sect. 5) and speed (Sect. 6) between consecutive CME pairs. Given that both WT and speed have heavy-tailed distributions, which are associated with long-range dependence, we perform a detrended fluctuation analysis to corroborate the long-range dependence in Sect. 7. This dependence is clearly seen in the difference between the speed of trailing and leading CMEs (Sect. 8). Finally, our discussion and conclusions are presented in Sects. 9 and 10, respectively.

2. CME database

The CME database is maintained by the CDAW team at the CDAW Data Center, NASA Goddard Space Flight Center1. The database includes the basic sky-plane measurements of CMEs such as speed, angular width, central position angle, acceleration, mass, and kinetic energy. The first-appearance time of each CME in the LASCO field of view (FOV) is also given in the catalog, which we use for determining the WT. The observations have been continuously available since 1996, with a major three-month gap in 1998 when the SOHO spacecraft was temporarily disabled (Yashiro et al. 2004, 2008; Gopalswamy et al. 2009). In order to avoid possible instrumental errors during the first stage of the LASCO mission, from the 29 577 CMEs recorded from January 1996 until the end of December 2018, in this work we consider only 28 270 CMEs observed after the second major data gap, starting with the CME observed on February 2, 1999 at 16:35, and finishing with the CME observed on December 31, 2018 at 09:48 UT. Taking into account that our analysis is performed over pairs of consecutive CMEs, using the second CME of each pair as a reference, the number of pairs is equal to the number of CMEs minus one. Furthermore, we skip the first CME after a data gap (which corresponds to the second CME of the registered pair during the data gap), this way discarding the data gaps from the analysis (it is important to note that we do not apply any restriction to select the leading-trailing CME pairs, so that a leading one can be a trailing one of a previous event, and these are treated as an independent pair). During the entire period (1996–2018), 1003 data gaps longer than two hours were reported, from these data gaps, 673 occurred during our period of study. Although, in some cases (specifically, 164), more than one data gap took place in between two observed CMEs. Therefore, our final set contains 27 761 CMEs.

3. Waiting time

We define the CME WT as the time elapsed between the first observation of two consecutive CMEs in the LASCO FOV. It is important to note that the CME initial time, as reported when using coronographic observations, is not exactly the launching time, but the time when a CME first appears in the coronograph FOV above the occulting disk. Therefore, there is a small non-uniform error in assigning the initial time. This error is smaller for limb CMEs and increases when the source region of the CME approaches the center of the disk. This effect can be seen in Fig. 1, where we show the delay time of a CME observed by a coronograph with a two solar radii (R) occulting disk. This is equivalent to the time that it takes for a CME with given constant speed (from 200 to 2000 km s−1, marked by colors in the figure) to travel a distance of 2 R projected in the plane of the sky, as a function of the helio-longitude of the CME source region. The delay time (t) is given as t = [sin(θ) + tan (w/2) cos (θ)]1/2 R/v, where θ is the helio-longitude, and w and v are the CME width and speed, respectively. This analysis shows that for very slow CMEs (V ≤ 200 km s−1), the delay time is < 1 h when the source region is located beyond 45° of helio-longitude, whereas for very fast CMEs (V >  1000 km s−1), the delay time may be considered constant for helio-longitudes higher than 30°. We note that this basic analysis assumes similar visibility for front- and back-sided CMEs. The delay time is also shorter for wider CMEs for a given speed.

thumbnail Fig. 1.

Delayed time of first observation of a CME with constant speed, occulted by a 2 R disk as a function of the helio-longitude of the source region. Continuous and dashed lines correspond to CME width of 30° and 60°, respectively.

As shown in Fig. 1, narrow and slow CMEs launched close to the center of the disk () have a delay of the observed initial time longer than two hours. Although, CMEs may measure any width and may be launched in any longitude (the center/limb launching position is an observational effect, and thus, the probability that a CME of a given speed is launched inside this longitude range is ∼1/3). Therefore, as we have a large number of events (launched with different velocities, widths and longitudes), by selecting a minimum of two hours for the WT, we are able to neglect the error introduced by the slow-narrow CMEs launched close to the center of the disk.

Another source of systematic errors is the cadence of the instrument. In the case of LASCO, this is between 12 and 40 min (depending on the observation mode). Therefore, WTs < 1 h are difficult to quantify properly. On the other hand, only slow (V <  600 km s−1) and central source region (|θ| < 15°) CMEs have delays > 2 h (Fig. 1). Therefore, for the majority of CMEs the delay is under two hours, and this can be taken as the WT resolution. As a consequence, the bin size of the distributions in the next section has been taken as two hours. We note that the cadence of LASCO was doubled after August 2010, this caused a steep elevation of the CME rate reported in automated (CACTus and SEEDS), but not in CDAW catalogs (Wang & Colaninno 2014; Hess & Colaninno 2017). This is aside from the fact that the real number of CMEs was higher for cycle 24 than for cycle 23 (Petrie 2015). The cadence change affected only very narrow CMEs, while wide CMEs such as halos were not affected (Gopalswamy et al. 2013).

The statistical distributions are useful for understanding the physical nature of the underlying process. In particular, the CME WT may help to determine whether the CME process is purely stochastic or if there is a dependence or physical connection between consecutive CMEs. This is an important question in terms of CME triggering, and therefore a major issue in terms of prediction of space weather.

Figure 2 shows the distribution of the WT (black circles) of 27 761 CMEs observed by LASCO from February 1999 up to December 2018.

thumbnail Fig. 2.

Relative frequency (divided by the bin size of 2 h) of observed WT distribution (black circles), along the fit exponential (blue) and Pareto (red) distributions. The inner plot is the so-called “Q–Q plot”, which graphically illustrates that the Pareto distribution (red line) is closer to the “true” distribution (represented by the gray line) than the exponential distribution.

The main characteristics of the WT distribution are: (i) Almost all CMEs (98%) occurred inside a time interval of 25 h. Only 567 CMEs have a time difference longer than 25 h. (ii) ∼85% of the events have a time difference shorter than 10 h. (iii) The WT is under five hours for ∼61% of the events.

Using a maximum likelihood method, we fit an exponential and a Pareto Type IV distribution to the observed WT. These distributions are plotted with blue and red continuous lines, respectively, in Fig. 2. Their probability density functions (PDF) and best fit parameters found through a maximum likelihood method are:

  • Exponential distribution (blue curve):

    f ( x ) = λ exp x λ , $$ \begin{aligned} f(x)=\lambda \exp ^{-x \lambda } , \end{aligned} $$(1)

    with λ = 0.17 and a mean value of 5.75.

  • Pareto Type-IV distribution (red curve):

    f ( x ) = κ 1 / γ α γ [ 1 + ( κ x μ ) 1 / γ ] 1 α ( x μ ) 1 + 1 / γ , $$ \begin{aligned} f(x)= \frac{ \kappa ^{-1/\gamma } \alpha }{\gamma } \left[ 1 + \left( \frac{\kappa }{ x - \mu } \right)^{-1/\gamma }\right]^{-1 - \alpha } (x - \mu )^{-1 + 1/\gamma } , \end{aligned} $$(2)

    with κ = 10.49, α = 2.84, γ = 0.81, and μ = 0.07, and a mean value of 5.79.

The Exponential distribution is memory-less and arises from a “simple” stochastic process like the waiting times of a Poisson process (Marshall 2007); the memory-less property implies that there is no dependence or precondition between consecutive events. On the other hand, the Pareto distribution has a heavy-tail2 (Marshall 2007) and has been associated with long-range dependence or long-memory processes in a large variety of systems (Samorodnitsky 2007). It is a very interesting property often associated with nonstationary processes, and the scaling and fractal behavior of the system and phase transitions (see, Samorodnitsky 2007, for a survey of different points of view of the long-range dependence). As shown in Fig. 2, the Pareto distribution (red line) better follows the observed WT (black circles) in the entire range. This fact is better demonstrated by the quantile-quantile or “Q–Q plot” on the inner frame of Fig. 2.

Q–Q plots are used to graphically assess the quality of the fit between the model and the empirical distribution, which is represented as a straight gray line in this case. To quantify the differences between the proposed and the observed distributions, we use the Bayesian information criterion (BIC, Schwarz 1978), which uses the maximum likelihood to determine the deviation between the empirical and the proposed distributions, with a penalty term for the number of parameters of the distributions. The distribution favored by BIC ideally corresponds to the candidate model which is a posteriori. Therefore, the relevant indicator of the BIC is the relative change. The fitting results are shown in the last row of Table 1, where the Pareto – exponential ΔBIC is shown in the last column. The fit parameters of the Pareto Type IV and exponential are shown in Cols. 2–5 and 7, respectively.

Table 1.

WT distribution parameters.

As expected and clearly seen in Fig. 3, the WT changes appreciably during the different phases of the solar cycle. This causes the nonstationary aspect of the process. Therefore, an analysis taking into account the changes on each phase of the solar cycle is necessary.

thumbnail Fig. 3.

WT between every two consecutive CMEs (gray plus symbols) as a function of time during the analyzed period. The black dots correspond to a smoothed running average of the WT, whereas the cyan line corresponds to the Sunspot number. The different phases of the solar cycle (upper right) are marked with colored plus symbols, which correspond to the 30-day WT average.

4. Waiting time during the solar cycle

The WT varies with time and follows the solar cycle, as seen in Fig. 3, where we have plotted the CME WT as a function of time (gray plus symbols); its smoothed version (running average of 60 points, black dots) and a 30-day WT average (colored plus symbols), which shows a clear differentiation of the WT during the different solar cycle phases. In order to compare these changes with other solar cycle parameters, we plotted the sunspot number with a cyan line (taken from WDC-SILSO, Royal Observatory of Belgium, Brussels). In general, the shortest WTs correspond to periods of maximum solar activity and viceversa. In fact, during the maximum of activity of solar cycle 24, the mean WT was ≤4 h, whereas the mean WT reaches ∼11 h during the descending phase of cycle 24.

There are clear differences of the CME WT rate evolution during the solar cycle, the 30-day average of the WT shows major changes along the time during low activity phases, whereas for high activity phases, it remains relatively constant. Following these changes, we divided the observed data in seven periods as marked in the upper right of Fig. 3 with different colors. The ascending phase of solar cycle 23 corresponds to the first stage of LASCO, is somehow less stable than the other phases (as seen by the magenta symbols in Fig. 3) and also contains the two major LASCO data gaps. Therefore, we did not take into account this phase for the rest of the analysis. The observed WT distributions for each selected phase of the cycle are plotted in Fig. 4 with colored circles. Similarly to the entire time range WT distribution, we fit Pareto Type IV (Eq. (2)) and an exponential (Eq. (1)) distributions.

thumbnail Fig. 4.

Similar to Fig. 2 but for different phases of solar cycles 23 and 24, marked by colors. The continuous and dash lines correspond to Pareto Type IV and exponential distributions, respectively.

There are important differences between the Pareto Type IV distribution parameters at each phase of the solar cycle. The wider WT distribution (mean = 11.23 h) is associated with the descending phase of cycle 24. The minimum phase of the cycle also has a wide WT distribution (mean = 7.88 h). On the other hand, the maximum of 24 has the narrowest distribution (with a mean of 3.90 h). This is in accordance with the increase of the number of CMEs observed in cycle 24 (Petrie 2015). Table 1 shows the parameters of the Pareto Type IV (Cols. 2–5) and exponential (Col. 7) distributions that best fit the WT during each phase of the solar cycles.

The BIC shows that the observed WT distributions follow better the Pareto Type IV distribution. The differences between these distributions are smaller during the phases of low activity (minimum and descending phase of cycle 24), pointing towards a simple stochastic WT process. On the other hand, during high activity periods, the BIC differences clearly show that the WT follows the heavy-tailed Pareto Type IV distribution pointing towards a long-range dependence or memory process. This long-range dependence of the WT time maybe reflected in other characteristics of the CME phenomena. Therefore, we explore the source region and speed of consecutive CMEs in the following sections.

5. Angular (PA) difference of consecutive CMEs

The CME source region may be used to explore the possible relationship between consecutive CMEs. In this analysis, we use the CME position angle (PA) as a proxy of the CME source region (Lara 2008). In particular, we use the central position angle (CPA), which is the mid-angle between the CME edges3. Out of the total number of observed CMEs, we excluded halo CMEs (which, by definition, do not have CPAs) and performed the present analysis over 26 086 CMEs. The rationale is that if there is long-range dependence in the CME WT, it should be statistically reflected in the distribution of the angular difference between the PA of consecutive CMEs (PAdiff). If the source regions of both consecutive events are relatively close, then the PAdiff must approach zero.

The normalized distribution of the observed PAdiff between consecutive CMEs is marked with black circles in Fig. 5. Over-plotted are the Gaussian PDFs,

f ( x ) = 1 2 π σ exp ( ( x μ ) 2 2 σ 2 ) , $$ \begin{aligned} f(x)=\frac{1}{\sqrt{2\pi }\sigma } \exp \left(- \frac{(x-\mu )^2}{2 \sigma ^2}\right), \end{aligned} $$(3)

thumbnail Fig. 5.

Distribution of the difference of the position angle between any two consecutive CMEs with a bin size of 10°. The green, pink, red, and blue continuous lines represent the Gaussian distributions added up (black curve) to fit the observations (black circles). The orange diamonds represent the observed distribution but for narrow (W <  50°) CMEs.

which fit the data. We note that the three wide Gaussian distributions are centered at μ1, 3 = ∼ ± 180° (green and blue curves), and μ2 ∼ 0° (pink curve) representing the PAdiff of randomly distributed events. What is unexpected for a random distribution of PAs is the narrow peak at the center of the distribution (around 0°), which can also be fit by a Gaussian (but narrower) distribution shown by the red line in Fig. 5. This large peak suggests that an important number of consecutive CMEs have less than 30° of PA difference, meaning that these CMEs are produced in a very close source region (even assuming that half of the events in the central peak of Fig. 5 were produced on opposite sides of the Sun, this peak is ∼3 times larger than the others). It is worth noting that we have conserved the PA reference point (north pole) and the sign of the PAdiff to facilitate the fitting process and to make the symmetries of the PA difference distribution clearer. Of course, the real angular distance ranges from 0° to 180°. We are aware of the difficulty of using coronographs to give information about the real position, in the low corona, of the source region of CMEs. As the observations are projected in the plane of the sky, it is difficult to characterize the actual position of the CME source region. Nevertheless, the PA is a good indicator of the radial direction of the CME, which in turn, extrapolated backwards to the solar surface, indicates the CME source region (see also, Gopalswamy et al. 2003; Lara 2008; Murray et al. 2018).

To diminish the effect of the projection, we performed the same analysis using narrow (width < 50°) CMEs only (we note that wherever we restrict the number of events, the imposed condition has to be fulfilled by the two events of the leading-trailer pair of CMEs, if one of them does not meet the condition, the pair is discarded for that particular part of the analysis). Obviously, the number of events is lower (12 477), but the statistical characteristics do not change, as shown by the observed PAdiff distribution of narrow CMEs plotted with diamonds and the associated Gaussian distributions (orange color) in Fig. 5. In fact, the inference of a close source region of consecutive CMEs is more accentuated when we take into account only the narrow CMEs, as suggested by the larger height of the central peak. These results suggest that the subset of CMEs with small PAdiff may be spatially-related, in concordance with the scenario where new emerging flux frequently occurs in an active region which has already emerged (Liggett & Zirin 1985) or in its vicinity forming nests or clusters of flux emergence known as active longitudes (see, van Driel-Gesztelyi & Green 2015 and references therein).

6. Speed difference of consecutive CMEs

Taking into account the short time difference (WT ≤ 10 h) between a vast majority of consecutive CMEs (note that after 10 h, a CME with speed of 500 km s−1 has traveled 25 R, and therefore remains inside the LASCO C3 field of view) and the close spatial relation between the source region (PAdiff ≤ 30°) of a group of them, a natural question arises: Is there any relationship between the speed of the leading (Vlead) and trailing (Vtrail) CMEs? To explore this possibility, we tested several distributions and found that the Vdiff = Vtrail − Vlead best follows a generalized student T-distribution of the form:

f ( x ) = 1 ν σ Beta ( ν 2 , 1 2 ) [ ν ν + ( x μ σ ) 2 ] 1 + ν 2 · $$ \begin{aligned} f(x) = \frac{1}{\sqrt{\nu } \sigma \text{ Beta} \left(\frac{\nu }{2},\frac{1}{2}\right)} \left[ \frac{\nu }{\nu + \left( \frac{x-\mu }{\sigma }\right)^2 }\right]^{\frac{1+\nu }{2}}\cdot \end{aligned} $$(4)

Figure 6 shows the observed (open circles) and the fitted T-distributions (continuous lines) of the Vdiff during the whole period (teal color) and during the individual periods of the solar cycle, plotted with different colors as marked in the Figure. The corresponding parameters are shown in Table 2. It is important to note that the student T-distribution is also a heavy-tailed distribution. In this case, the proper equivalent for an exponential distribution is given by the Laplace distribution (often known as the double-exponential distribution, Forbes 2010) of the form:

f ( x ) = { 1 2 β exp ( ( x μ ) β ) if x 0 1 2 β exp ( ( μ x ) β ) if x < 0 . $$ \begin{aligned} f(x)= {\left\{ \begin{array}{ll} \frac{1}{2\beta }\exp \left(\frac{-(x-\mu )}{\beta }\right)&\mathrm{if } \ x \ge 0 \\ \frac{1}{2\beta }\exp \left(\frac{-(\mu -x)}{\beta }\right)&\mathrm{if } \ x < 0. \end{array}\right.} \end{aligned} $$(5)

thumbnail Fig. 6.

Observed (circles) and fit generalized student T- (continuous lines) and Laplace (dotted lines) distributions of the Vdiff during different phases of the solar cycle and the entire period of study (teal color).

Table 2.

Parameters of the Vdiff distributions.

The BIC shows that the student T-distribution better fits the observed Vdiff than the double exponential distribution (the parameters of both distributions as well as the Δ BIC = BICStudent − BICLaplace are shown in Table 2) and allows us to conclude that Vdiff distribution is heavy-tailed, and therefore, the associated CMEs may have long-range dependence.

7. Long-range dependence

The heavy-tailed distributions, such as the WT and the speed difference of consecutive CMEs, suggest that the underlying process may not be associated with a simple Poisson process, instead, some long-range dependence between the events may be present. In this section, we apply a test, often used to determine the presence of long-range dependence, on the CME variables. Historically, the long-range dependence has been associated with time series where correlation decays slowly (decay exponentially in terms of the lag) or with a spectral density pole at zero frequency (i.e., at the origin). Although, to apply these criteria, the time series under consideration must be stationary, and unfortunately in our case, this condition is not fulfilled. To overcome this restriction, the detrended fluctuation analysis (DFA) was proposed by Peng et al. (1994) to study the long-range dependence of nucleoids of DNA chains. Basically, the goal is the determination of the scaling exponent αs of the (log–log) dependence between the detrended fluctuation function F(s) and the size s of the window where F(s) was computed. The computational details can be found in Kantelhardt et al. (2001).

We applied the DFA to the WT, PAdiff, and the speed difference between the trailing and leading CME (Vdiff) over the entire time range of study and the resulting DFA are shown in Fig. 7, in black, blue, and green, respectively. Valuable information about the range dependence of the data is given by the scaling exponent. For a random process such as Brownian motion, αs = 0.5. This is the case of the PAdiff and the randomly re-arranged WT and Vdiff. The former is in accordance with its (sum of) Gaussian distribution.

thumbnail Fig. 7.

Detrended fluctuation as a function of the window size or scale (s) for the entire time period (upper left) of the WT (black), Vdiff (blue) and PAdiff (green). The filled circles and continuous lines correspond to the actual data, whereas the empty circles and dashed lines correspond to the same, but randomly re-arranged (shuffled), data. The same analysis but for each phase of the cycle applied to the PAdiff (upper right), Vdiff (lower left) and WT (lower right). The αs value during each phase of the cycle is marked with the corresponding color. Open circles and dashed lines correspond to the DFA applied to a random vector created for comparison purposes.

On the other hand, αs >  0.5 implies long-range dependence. In this case, the αs is 0.67 for the speed difference and 0.77 for the WT. To ensure these results, we randomly shifted the values of those series and applied the DFA analysis, the results are plotted with open circles and dashed lines in Fig. 7 and for all cases αs ∼ 0.5.

Similar to the probability distributions of the WT, PAdiff, and Vdiff, the scaling exponent changes with the phase of the solar cycle. These changes are relatively small, for PAdiff, αs ≈ 0.5 during all phases (upper-right panel of Fig. 7). Similarly to during the minimum activity period, αs ≈ 0.5 for both Vdiff and WT (bottom-left and right panels, respectively). Conversely, the αs reaches larger values during the maximum and descending phases of the solar cycle.

In summary, the DFA shows that there is long-range dependence on both the WT and Vdiff differences between successive CMEs. This is evident during the maximum and declining phase of the solar cycle.

8. Trailing and leading CME speed relationship

As suggested by the precedent analysis, consecutive CMEs have some dependence. In this section, we explore the possibility of finding a relationship between the speed of these events. In particular, the possibility that the speed of a CME may be affected by a previous event. For instance, it may be surmised that the leading CME changes the ambient medium where the trailing CME is going to travel, in such a way that the leading CME has an influence on the dynamics of the trailing CME as predicted by drag force models (Cargill 2004; Borgazzi et al. 2009; Vršnak et al. 2010). At this point, it is important to recall that we do not apply any restriction to select the leading-trailing CME pairs, so that a leading can be a trailing of a previous event, and these are treated as an independent pair. In order to investigate the possible statistical differences between the speed of the leading and trailing CMEs, we constructed subsets of spatial and temporally related pairs of CMEs with the following restrictions: (i) A group with no restriction, meaning all events, hereinafter called AllCMEs. (ii) A group of spatially-related consecutive CMEs where the PA difference (PAdiff) between the leading and trailing CMEs must be less than 30°. (iii) For both groups, the temporal association was established creating subsets of time windows where the WT is less than or equal to 2, 3, 4, ..., 25 h.

Then, we computed the mean speed value ( V ¯ $ \overline{V} $) of each subset of leading and trailing CMEs and found a small but persistent difference between the V ¯ $ \overline{V} $ of the trailing and leading CMEs. This is shown in the left panel of Fig. 8 where we plotted the (blue) histogram of the difference ( V diff ¯ $ \overline{V_{\mathrm{diff}}} $) between trailing minus the leading CME means (i.e., V diff ¯ = V trail ¯ V lead ¯ $ \overline{V_{\mathrm{diff}}} = \overline{V_{\mathrm{trail}}} - \overline{V_{\mathrm{lead}}} $) of the AllCMEs subsets. The main part of distribution is inside the −5 to 10 km s−1 range, but with an asymmetrical component that extends towards positive differences up to ∼17 km s−1. The asymmetry is evident for the spatially-related CMEs groups (red histogram), where the difference attains higher values. Gopalswamy et al. (2004) found larger differences for the case of CMEs associated with solar energetic particles (SEP), although these were CMEs from the same active region.

thumbnail Fig. 8.

Differences of trailing minus leading CME mean (left) and median (right) speed of AllCMEs (blue) and spatially-related (red) subsets for all the considered WTs.

The CME speed is well-fit by a lognormal distribution (Lara et al. 2006), which is also a heavy-tailed, and therefore the standard statistical tools, such as sample mean and sample standard deviation, are highly unstable (Pisarenko 2010). In this case, the standard deviations of the speed distributions are ∼200 km s−1, and the standard errors of the mean are of the order of the differences. Then, a value of a few km s−1 used to characterize the differences of the V ¯ $ \overline{V} $ has low statistical significance if we use the mean to characterize the speed distributions. Therefore, we have to use other parameters to quantify the observed differences. The quartiles (particularly the second quartile or median) are better suited to characterize skewed distributions. So, we computed the quartiles of all subsets: for example, Fig. 9 shows the quartiles4 of the speed distributions as a function of the phase of the cycle for four different WT subsets. Blue and red lines correspond to the trailing and leading spatially-related CMEs; whereas pink and light-blue correspond to the trailing and leading AllCMEs groups, respectively. As expected, the solar cycle variations are clear in the quartile behavior (Fig. 9), showing higher values during the maxima, and the lowest values are reached during the minimum activity period. We note that the quartiles clearly show the differences between the distributions of the leading and trailing CMEs.

thumbnail Fig. 9.

Quartiles of speed distributions of the leading (light-blue) and trailing (pink) AllCMEs consecutive pairs, as well as spatially-related leading (blue) and trailing (red) consecutive CME pairs, as a function of the phase of the solar cycle for four different WTs (as marked in each panel).

For comparison (with the V diff ¯ $ \overline{V_{\mathrm{diff}}} $), the right panel of Fig. 8 shows the differences between the medians of the trailing minus the leading speed distributions (i.e., V diff = V trail V lead $ \widetilde{V_{\mathrm{diff}}} = \widetilde{V_{\mathrm{trail}}}- \widetilde{V_{\mathrm{lead}}} $). In this case, the maximum of the V diff $ \widetilde{V_{\mathrm{diff}}} $ for the spatially-related CMEs is shifted towards negative values, but the high asymmetry is conserved and attains higher values than the V diff ¯ $ \overline{V_{\mathrm{diff}}} $. The most important fact shown by both V diff ¯ $ \overline{V_{\mathrm{diff}}} $ and V diff $ \widetilde{V_{\mathrm{diff}}} $ is the asymmetry towards positive values, implying that the Vtrail tends to be higher than the Vlead for consecutive CMEs. In order to quantify these differences and their statistical significance, we performed a sign test with the null hypothesis H0: V trail V lead 0 $ \widetilde{V_{\mathrm{trail}}} - \widetilde{V_{\mathrm{lead}}} \leq 0 $. Therefore, if it is unlikely that H0 is true, the alternative hypothesis, in this case Ha: V trail V lead > 0 $ \widetilde{V_{\mathrm{trail}}} - \widetilde{V_{\mathrm{lead}}} > 0 $, meaning V trail > V lead $ \widetilde{V_{\mathrm{trail}}} > \widetilde{V_{\mathrm{lead}}} $ is likely to be true.

The V diff $ \widetilde{V_{\mathrm{diff}}} $ and V diff ¯ $ \overline{V_{\mathrm{diff}}} $ (upper- and bottom-right panels) as well as the differences of the lower and upper quartiles (upper and bottom-left panels) are shown in Figs. 10 and 11 for cycles 23 and 24, respectively. For the sake of clarity, we selected the subsets whose median sign-test probability is lower than 0.15, 0.1, and 0.05, marked with small, medium, and large symbols, respectively. Moreover, the number of events in each subset with probabilities lower than 0.05 are shown in the lower-right panel. The error bars in these Figures correspond to the median absolute deviation divided by the square root of the number of elements of each subset. Here, it is clear that the positive Vdiff is statistically significative, and that there are excesses of up to ∼100 km s−1 for WT ≤ 2 h during the ascending and maximum phases of cycle 23, and a more modest excess of ∼20 km s−1 during the decreasing phase of cycle 24. These Vtrail excesses decrease rapidly when the WT increases reaching ∼20 km s−1 for cycle 23, and ∼0 km s−1 for cycle 24 when WT <  10 h.

thumbnail Fig. 10.

Difference between location parameters of the Vtrail − Vlead distributions as a function of the WT. The three quartiles and the mean values during the minimum and two phases of solar cycle 23 are marked by different colors. The size of the circles depends on the probability given by the sign test as: probability ≤0.15 (small), 0.1 (medium), and 0.05 (large), i.e., when the null hypothesis (H0: V trail V lead 0 $ \widetilde{V_{\mathrm{trail}}} - \widetilde{V_{\mathrm{lead}}} \leq 0 $) is rejected at the 15, 10, and 5% level based in the sign test. The events considered in each subset with probability ≤0.05 are also shown in the lower right panel.

thumbnail Fig. 11.

Same as Fig. 10 but for solar cycle 24 and for the entire time range considered in this work (blue color).

The spatially-related CMEs (marked with circles in Figs. 10 and 11) show the largest speed differences. Overall, the behavior of the Vdiff is clearly seen in the bottom panels of Figs. 10 and 11, where the upper, meaning the high velocity quartile (left), and the V diff ¯ $ \overline{V_{\mathrm{diff}}} $ (right) differences, are shown. The speed difference has higher values for WTs ≤ 2 h and decreases in a quasi-exponential way up to WTs ∼ 7 − 8 h. Then, it slowly decreases towards zero. Finally, two important facts that maintain the CME speed bounded: (i) the speed difference reaches high values when we consider spatially-related CMEs, whereas the values remain low for all CMEs; and (ii) the speed difference tends to zero for large WTs during the minimum activity period.

9. Discussion

The statistical characterization of the properties of consecutive CMEs is important because it sheds light on the generation processes and early stages of the CME phenomena. In particular, the WT statistical distribution has been addressed by several authors (see, Telloni et al. 2014, and references therein). Although, as far as we know, the associated PA and speed differences between successive CMEs have not been explored until now. A major subject of study is the nature of the CME initiation process: Is this a “pure” stochastic process? Or is there some kind of relationship between consecutive events? This possible relationship is called “memory” of the system or “long-term dependence,” and there are different tools to investigate if the events of a given system have this dependence. If the process and its time series are stationary, the most common approaches used to look for the existence of long-range dependence are the auto-correlation function and the power spectrum behavior at the origin. The CME process is not stationary, so we approached the problem through the shape of the statistical distribution tail and the DFA (see, Pisarenko 2010; Beran et al. 2013, for extensive reviews of this subject). In this way, a “pure” stochastic process has an exponential distribution and is considered as a memory-less system. This is the case for the waiting times of a Poisson process (Lee & Fapojuwo 2005), for example. In contrast, processes where the distribution tail decreases more slowly than the exponential, the so called heavy-tailed distributions, such as the Pareto family and the generalized student T-distributions, are often associated with long-range dependence or memory processes (Pisarenko 2010; Beran et al. 2013).

The main question in this statistical study can be formulated as: Is the CME process “purely” stochastic? Or does the CME process have a long-range dependence or “memory” in such a way that an event is influenced by the occurrence of past events?

To address this question, Wheatland (2003) analyzed a set of CMEs observed by LASCO during the ascending phase and part of the maximum of cycle 23 (1996–2001). The author fit the tail of the WT distribution (WT >  10 h) by a power law function (with index γ ≈ −2.36). The author also applied a Bayesian statistical analysis and concluded that the observed WT distribution corresponds to a time-dependent Poisson process. Performing a similar analysis but using LASCO observations in a restricted time range around the maximum of cycle 23 (1999–2001), Moon et al. (2003) reached the same conclusion that the WT distribution of 3187 CMEs may be represented by a nonstationary Poisson process.

A power law distribution, using only 113 CMEs observed by the Solar Maximum Mission, during the minimum phase of the cycle (1984–1989), was also found by Méndez Berhondo et al. (2006), with γ = 1.41. In this case, the authors associated the WT with a characteristic length Λ, (between two consecutive CMEs) as Λ ∝ VAWT, where VA is the Alfvén speed (this implies causality), and conclude that the WT power law distribution is consistent with their picture of characteristic length.

It is clear that a power law may be used to fit the tail of the WT distribution, but as mentioned earlier, 83% of the total events have a WT under 10 h. Therefore, the tail after this WT time is not necessarily representative of the whole process (only the extreme events fall in the extended tail of the distribution). In fact, a power law distribution is heavy-tailed and may represent a scaling system, which, in turn, may have self-organized criticality or related to processes with long-range dependence or “memory” (Pisarenko 2010; Beran et al. 2013). For instance, Takahashi et al. (2016) found a scaling relationship between CMEs and proton events. Also, Telloni et al. (2014) found that the WT distribution may be fit by a Weibull distribution (a heavy-tailed distribution), implying some degree of correlation between CMEs in opposition to a purely stochastic process.

Our findings are in agreement with Jiggens & Gabriel (2009; who found that SEP events appear to have some memory indicating that events are not completely random) and with the scaling relationship between CMEs and SEP found by Takahashi et al. (2016).

The possible correlation between consecutive CMEs is supported by the remarkable peak at short PAdiff, which shows that the source region of a large number of CMEs is located within close distance (Fig. 5). This result confirms the finding of Moon et al. (2003) of an excess of events within an angular difference of 10°, these authors interpreted it as evidence of CMEs that occurred in the same active region. This result is also in agreement with the hypothesis that many CMEs (and flares) are produced in the so-called active longitudes (e.g., Gyenge et al. 2017 and references therein).

The WT distribution is telling us that, on average, the reconfiguration of the magnetic field after the launch of the leading CME, in order to replenish the free energy and launch the trailing CME, takes place in ∼3.7 h during the maximum phase of cycle 24 and ∼8.5 h during the minimum phase. This is an important fact from the point of view of the CME energetics, in particular in the case where two or more consecutive CMEs originated in the same source region. It is worth to note that the maximum fraction of free energy exhausted by one CME may be 25% (Gopalswamy et al. 2005; Emslie et al. 2012; Feng et al. 2013).

Finally, the speed distributions of all leading/trailing CMEs show a small but systematic tendency to have higher mean and median velocities for shorter WTs. The mean and median speed of the trailing CME are a few km s−1 higher than that of the leading CMEs, for WTs shorter than 7 h, after a WT of 10 h the differences remain constant and negligible. Nevertheless, these differences are higher (tens of km s−1) when we only take into account spatially-related CMEs (see Fig. 8).

As the mean is not the best location indicator for highly skewed distributions, we also analyzed the quartiles (which are better suited to characterizing these distributions), the results are plotted in Fig. 9. To assess the accuracy of the analysis, we applied the sign test to show that the differences between the medians of trailing and leading CMEs are statistically significant. The phases of the solar cycle where the sign-test probability (with the null hypothesis H0: V trail V lead 0 $ \widetilde{V_{\mathrm{trail}}} - \widetilde{V_{\mathrm{lead}}} \leq 0 $) is less than 0.1 are shown in Figs. 10 and 11.

These observations may be explained by the drag force models. It is clear that the magnetic field provides the necessary energy to launch and first acceleration of CMEs. Although this energy does not play an important role in the CME dynamics after few solar radii, at least for fast CMEs (such as those that form the tail of the distributions). Recently, Sachdeva et al. (2017) investigated CME dynamics using coronographic observations, and they concluded that for fast CMEs, the Lorentz force has its maximum effect at heights lower than 2.5 R, and then, at heights of 3.5–4 R becomes negligible compared with the drag force. In the case of slow CMEs, the Lorentz force may be significant at larger distances (up to 50 R, but slow CMEs do not contribute to the tail of the distributions). The drag force F D = 1 2 C D ρ A ( V CME V SW ) 2 $ F_{\mathrm{D}} = - \frac{1}{2} C_{\mathrm{D}} \rho A (V_{\mathrm{CME}} - V_{\mathrm{SW}})^2 $ is proportional to the drag coefficient (CD), the transverse area of the CME (A), the speed difference between the CME (VCME) and the solar wind (VSW), and the ambient density (ρ) (Cargill 2004; Borgazzi et al. 2009; Vršnak et al. 2010).

In a statistical analysis of a large number of events, the main values of A and VSW are similar for both leading and trailing CMEs. On the other hand, the leading CME may sweep the ambient medium, and the trailing CME encounters a less dense medium and therefore the drag force decreases, resulting in a faster CME. In this case, the shorter the WT, the higher the CME speed. This effect is clearly seen in Figs. 10 and 11.

10. Conclusions

We performed a statistical study of 27 761 CMEs observed by SOHO/LASCO from the ascending phase of cycle 23 to the declining phase of cycle 24, looking for long-range or memory signatures in the genesis and the first stages of evolution of CMEs. The main findings of this work are:

  • The time elapsed between consecutive CMEs or WT for 98% (61%) of all events is less than 25 (5) h.

  • We found that a Pareto Type IV probability density function fits the observed WT distribution during the entire period (1996–2018), and during the different phases of the solar cycle.

  • As expected, the WT varies with the solar cycle, with large mean values (∼10 h) and high rates of change over time at low activity phases, and short mean (∼4 h) and ∼0 rates of change over time when the activity is high. These mean values are consistent with recurrence times in large active regions reported before (see Fig. 3 in Gopalswamy et al. 2005).

  • The distribution of the PA difference of consecutive CMEs (as a proxy to the CME source region) shows four components: three wide Gaussian distributions centered at 0° and ±180° related to the spatially uncorrelated CMEs; the fourth is a narrow Gaussian centered at 0° related to CMEs with very close source regions.

  • The difference between the trailing and leading CME speed follows a generalized student T-distribution.

  • The heavy-tailed distributions as well as the DFA show that there is a long-range dependence or memory in the WT and speed difference of consecutive CMEs.

  • The long-range dependence of the CME speed is confirmed by the strong statistical evidence that the speed of the trailing CMEs is slightly higher than that of leading CMEs where the WTs are < 10 h.

  • The trailing or leading speed difference is higher for the CMEs with close source-region association.

These findings point towards the fact that leading CMEs modify (sweep up) the ambient medium, causing the trailing CMEs to encounter less opposition to their movement. This difference is clear for CMEs produced by individual active regions (see Fig. 9 of Gopalswamy et al. 2004).

In summary, by analyzing a large number of consecutive CMEs, we find that their WT distribution follows a Pareto Type IV distribution, and the associated speed difference follows a generalized student T-distribution. These heavy-tailed distributions along the DFA suggest long-term dependence in the CME process. The position-angle difference distribution shows a large Gaussian peak centered at zero, indicating a close spatial relationship between consecutive CMEs. Furthermore, there is a small but statistically significant difference between the speed of consecutive CMEs indicating that the trailing CME is a few km s−1 faster than the leading CME when the WT is < 10 h. All these findings suggest a physical connection between a considerable population of consecutive CMEs.


2

Heavy-tailed distributions are characterized by the slow decreasing tail as compared with the exponential distribution.

3

The CDAW catalog also includes the measurement position angle (MPA), where the height-time measurements were taken (Yashiro et al. 2004): we do not use this angle in the present study.

4

This presentation is similar to the “box-plot”, an often-used graphical method to illustrate the differences between distributions where the quartiles (and the minimum and maximum values or “whiskers”) of the distributions that one wants to compare, are plotted side by side. The difference is that in this case, we do not plot the whiskers.

Acknowledgments

The CME catalog is generated and maintained at the CDAW Data Center by NASA and The Catholic University of America in cooperation with the Naval Research Laboratory. SOHO is a project of international cooperation between ESA and NASA. This work was partially supported by CONACyT (179588) and UNAM, PASPA-DGAPA program. N. Gopalswamy is supported by NASA’s Heliophysics LWS and GI programs. We thank the anonymous referee for his/her useful and constructive comments.

References

  1. Aschwanden, M. J., Caspi, A., Cohen, C. M. S., et al. 2017, ApJ, 836, 17 [NASA ADS] [CrossRef] [Google Scholar]
  2. Barabási, A.-L. 2005, Nature, 435, 207 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  3. Beran, J., Feng, Y., Ghosh, S., & Kulik, R. 2013, Long-Memory Processes (Berlin, Heidelberg: Springer) [CrossRef] [Google Scholar]
  4. Borgazzi, A., Lara, A., Echer, E., & Alves, M. V. 2009, A&A, 498, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, Sol. Phys., 162, 357 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cargill, P. J. 2004, Sol. Phys., 221, 135 [NASA ADS] [CrossRef] [Google Scholar]
  7. Emslie, A. G., Dennis, B. R., Shih, A. Y., et al. 2012, ApJ, 759, 71 [NASA ADS] [CrossRef] [Google Scholar]
  8. Feng, L., Wiegelmann, T., Su, Y., et al. 2013, ApJ, 765, 37 [NASA ADS] [CrossRef] [Google Scholar]
  9. Forbes, T. G. 2000, J. Geophys. Res., 105, 23153 [NASA ADS] [CrossRef] [Google Scholar]
  10. Forbes, C. 2010, Statistical Distributions, 4th edn. (Hoboken, NJ: John Wiley & Sons) [CrossRef] [Google Scholar]
  11. Gopalswamy, N. 2004, in The Sun and the Heliosphere as an Integrated System, eds. G. Poletto, & S. T. Suess, Astrophys. Space Sci. Lib., 317, 201 [NASA ADS] [CrossRef] [Google Scholar]
  12. Gopalswamy, N. 2016, Geosci. Lett., 3, 8 [NASA ADS] [CrossRef] [Google Scholar]
  13. Gopalswamy, N. 2017, ArXiv e-prints [arXiv:1709.03165] [Google Scholar]
  14. Gopalswamy, N., Yashiro, S., Kaiser, M. L., Howard, R. A., & Bougeret, J.-L. 2001, ApJ, 548, L91 [Google Scholar]
  15. Gopalswamy, N., Lara, A., Yashiro, S., & Howard, R. A. 2003, ApJ, 598, L63 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gopalswamy, N., Yashiro, S., Krucker, S., Stenborg, G., & Howard, R. A. 2004, J. Geophys. Res. (Space Phys.), 109, A12105 [NASA ADS] [CrossRef] [Google Scholar]
  17. Gopalswamy, N., Yashiro, S., Liu, Y., et al. 2005, J. Geophys. Res. (Space Phys.), 110, A09S15 [Google Scholar]
  18. Gopalswamy, N., Yashiro, S., Michalek, G., et al. 2009, Earth Moon and Planets, 104, 295 [Google Scholar]
  19. Gopalswamy, N., Mäkelä, P., Xie, H., & Yashiro, S. 2013, Space Weather, 11, 661 [NASA ADS] [CrossRef] [Google Scholar]
  20. Gyenge, N., Singh, T., Kiss, T. S., Srivastava, A. K., & Erdélyi, R. 2017, ApJ, 838, 18 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hess, P., & Colaninno, R. C. 2017, ApJ, 836, 134 [NASA ADS] [CrossRef] [Google Scholar]
  22. Jiggens, P. T. A., & Gabriel, S. B. 2009, J. Geophys. Res. (Space Phys.), 114, A10105 [NASA ADS] [CrossRef] [Google Scholar]
  23. Kantelhardt, J. W., Koscielny-Bunde, E., Rego, H. H. A., Havlin, S., & Bunde, A. 2001, Phys. A Stat. Mech. Appl., 295, 441 [CrossRef] [Google Scholar]
  24. Krista, L. D., & Reinard, A. A. 2017, ApJ, 839, 50 [NASA ADS] [CrossRef] [Google Scholar]
  25. Lara, A. 2008, ApJ, 688, 647 [NASA ADS] [CrossRef] [Google Scholar]
  26. Lara, A., Gopalswamy, N., Xie, H., et al. 2006, J. Geophys. Res. (Space Phys.), 111, A06107 [NASA ADS] [Google Scholar]
  27. Lee, I. W., & Fapojuwo, A. O. 2005, Comput. Commun., 29, 1 [NASA ADS] [CrossRef] [Google Scholar]
  28. Liggett, M., & Zirin, H. 1985, Sol. Phys., 97, 51 [NASA ADS] [CrossRef] [Google Scholar]
  29. López, F. M., Hebe Cremades, M., Nuevo, F. A., Balmaceda, L. A., & Vásquez, A. M. 2017, Sol. Phys., 292, 6 [NASA ADS] [CrossRef] [Google Scholar]
  30. Lugaz, N., Temmer, M., Wang, Y., & Farrugia, C. J. 2017, Sol. Phys., 292, 64 [Google Scholar]
  31. Marshall, A. 2007, Life Distributions: Structure of Nonparametric, Semiparametric, and Parametric Families (New York: Springer) [Google Scholar]
  32. Méndez Berhondo, A. L., Rodríguez Taboada, R. E., & Larralde, L. A. 2006, Ap&SS, 302, 213 [NASA ADS] [CrossRef] [Google Scholar]
  33. Michalek, G., & Yashiro, S. 2013, Adv. Space Res., 52, 521 [NASA ADS] [CrossRef] [Google Scholar]
  34. Moon, Y.-J., Choe, G. S., Wang, H., & Park, Y. D. 2003, ApJ, 588, 1176 [NASA ADS] [CrossRef] [Google Scholar]
  35. Murray, S. A., Guerra, J. A., Zucca, P., et al. 2018, Sol. Phys., 293, 60 [NASA ADS] [CrossRef] [Google Scholar]
  36. Niembro, T., Cantó, J., Lara, A., & González, R. F. 2015, ApJ, 811, 69 [NASA ADS] [CrossRef] [Google Scholar]
  37. Peng, C.-K., Buldyrev, S. V., Havlin, S., et al. 1994, Phys. Rev. E, 49, 1685 [NASA ADS] [CrossRef] [Google Scholar]
  38. Petrie, G. J. D. 2015, ApJ, 812, 74 [NASA ADS] [CrossRef] [Google Scholar]
  39. Pisarenko, V. 2010, Heavy-Tailed Distributions in Disaster Analysis. Advances in Natural and Technological Hazards Research (Dordrecht: Springer), 30 [Google Scholar]
  40. Sachdeva, N., Subramanian, P., Vourlidas, A., & Bothmer, V. 2017, Sol. Phys., 292, 118 [NASA ADS] [CrossRef] [Google Scholar]
  41. Samorodnitsky, G. 2007, Foundations and Trends® in Stochastic Systems, 1, 163 [CrossRef] [Google Scholar]
  42. Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
  43. Sheng, H., Chen, Y. Q., & Qiu, T. 2011, Fluct. Noise Lett., 10, 93 [CrossRef] [Google Scholar]
  44. Takahashi, T., Mizuno, Y., & Shibata, K. 2016, ApJ, 833, L8 [NASA ADS] [CrossRef] [Google Scholar]
  45. Telloni, D., Carbone, V., Lepreti, F., & Antonucci, E. 2014, ApJ, 781, L1 [NASA ADS] [CrossRef] [Google Scholar]
  46. Thompson, B. J., Cliver, E. W., Nitta, N., Delannée, C., & Delaboudinière, J.-P. 2000, Geophys. Res. Lett., 27, 1431 [Google Scholar]
  47. Tousey, R., Bartoe, J. D. F., Bohlin, J. D., et al. 1973, Sol. Phys., 33, 265 [NASA ADS] [Google Scholar]
  48. van Driel-Gesztelyi, L., & Green, L. M. 2015, Liv. Rev. Sol. Phys., 12, 1 [NASA ADS] [CrossRef] [Google Scholar]
  49. Vršnak, B., Žic, T., Falkenberg, T. V., et al. 2010, A&A, 512, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Wang, Y.-M., & Colaninno, R. 2014, ApJ, 784, L27 [NASA ADS] [CrossRef] [Google Scholar]
  51. Webb, D. F., & Howard, T. A. 2012, Liv. Rev. Sol. Phys., 9, 3 [Google Scholar]
  52. Wheatland, M. S. 2003, Sol. Phys., 214, 361 [NASA ADS] [CrossRef] [Google Scholar]
  53. Yashiro, S., Gopalswamy, N., Michalek, G., et al. 2004, J. Geophys. Res. (Space Phys.), 109, A07105 [NASA ADS] [CrossRef] [Google Scholar]
  54. Yashiro, S., Michalek, G., & Gopalswamy, N. 2008, Annal. Geophys., 26, 3103 [Google Scholar]

All Tables

Table 1.

WT distribution parameters.

Table 2.

Parameters of the Vdiff distributions.

All Figures

thumbnail Fig. 1.

Delayed time of first observation of a CME with constant speed, occulted by a 2 R disk as a function of the helio-longitude of the source region. Continuous and dashed lines correspond to CME width of 30° and 60°, respectively.

In the text
thumbnail Fig. 2.

Relative frequency (divided by the bin size of 2 h) of observed WT distribution (black circles), along the fit exponential (blue) and Pareto (red) distributions. The inner plot is the so-called “Q–Q plot”, which graphically illustrates that the Pareto distribution (red line) is closer to the “true” distribution (represented by the gray line) than the exponential distribution.

In the text
thumbnail Fig. 3.

WT between every two consecutive CMEs (gray plus symbols) as a function of time during the analyzed period. The black dots correspond to a smoothed running average of the WT, whereas the cyan line corresponds to the Sunspot number. The different phases of the solar cycle (upper right) are marked with colored plus symbols, which correspond to the 30-day WT average.

In the text
thumbnail Fig. 4.

Similar to Fig. 2 but for different phases of solar cycles 23 and 24, marked by colors. The continuous and dash lines correspond to Pareto Type IV and exponential distributions, respectively.

In the text
thumbnail Fig. 5.

Distribution of the difference of the position angle between any two consecutive CMEs with a bin size of 10°. The green, pink, red, and blue continuous lines represent the Gaussian distributions added up (black curve) to fit the observations (black circles). The orange diamonds represent the observed distribution but for narrow (W <  50°) CMEs.

In the text
thumbnail Fig. 6.

Observed (circles) and fit generalized student T- (continuous lines) and Laplace (dotted lines) distributions of the Vdiff during different phases of the solar cycle and the entire period of study (teal color).

In the text
thumbnail Fig. 7.

Detrended fluctuation as a function of the window size or scale (s) for the entire time period (upper left) of the WT (black), Vdiff (blue) and PAdiff (green). The filled circles and continuous lines correspond to the actual data, whereas the empty circles and dashed lines correspond to the same, but randomly re-arranged (shuffled), data. The same analysis but for each phase of the cycle applied to the PAdiff (upper right), Vdiff (lower left) and WT (lower right). The αs value during each phase of the cycle is marked with the corresponding color. Open circles and dashed lines correspond to the DFA applied to a random vector created for comparison purposes.

In the text
thumbnail Fig. 8.

Differences of trailing minus leading CME mean (left) and median (right) speed of AllCMEs (blue) and spatially-related (red) subsets for all the considered WTs.

In the text
thumbnail Fig. 9.

Quartiles of speed distributions of the leading (light-blue) and trailing (pink) AllCMEs consecutive pairs, as well as spatially-related leading (blue) and trailing (red) consecutive CME pairs, as a function of the phase of the solar cycle for four different WTs (as marked in each panel).

In the text
thumbnail Fig. 10.

Difference between location parameters of the Vtrail − Vlead distributions as a function of the WT. The three quartiles and the mean values during the minimum and two phases of solar cycle 23 are marked by different colors. The size of the circles depends on the probability given by the sign test as: probability ≤0.15 (small), 0.1 (medium), and 0.05 (large), i.e., when the null hypothesis (H0: V trail V lead 0 $ \widetilde{V_{\mathrm{trail}}} - \widetilde{V_{\mathrm{lead}}} \leq 0 $) is rejected at the 15, 10, and 5% level based in the sign test. The events considered in each subset with probability ≤0.05 are also shown in the lower right panel.

In the text
thumbnail Fig. 11.

Same as Fig. 10 but for solar cycle 24 and for the entire time range considered in this work (blue color).

In the text

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.