Issue |
A&A
Volume 511, February 2010
|
|
---|---|---|
Article Number | L2 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200913784 | |
Published online | 26 February 2010 |
LETTER TO THE EDITOR
Time-energy correlations in solar flare occurrence
E. Lippiello1 - L. de Arcangelis2,3 - C. Godano1
1 - Department of Environmental Sciences and CNISM, Second University of Naples, 81100 Caserta, Italy
2 - Department of Information Engineering and CNISM, Second University of Naples, 81031 Aversa (CE), Italy
3 - Institute Computational Physics for Engineering Materials, ETH, Zürich, Switzerland
Received 1 December 2009 / Accepted 5 February 2010
Abstract
Context. The existence of time-energy correlations in flare occurrence is still an open and much debated problem.
Aims. This study addresses the question whether statistically
significant correlations are present between energies of successive
flares as well as energies and waiting times.
Methods. We analyze the GOES catalog with a statistical approach
based on the comparison of the real catalog with a reshuffled one where
energies are decorrelated. This analysis reduces the effect of
background activity and is able to reveal the role of obscuration.
Results. We show the existence of non-trivial correlations
between waiting times and energies, as well as between energies of
subsequent flares. More precisely, we find that flares close in time
tend to have the second event with large energy. Moreover, after large
flares the flaring rate significantly increases, together with the
probability of other large flares.
Conclusions. Results suggest that correlations between energies
and waiting times are a physical property and not an effect of
obscuration. These findings could give important information on the
mechanisms for energy storage and release in the solar corona.
Key words: Sun: flares - methods: data analysis
1 Introduction
Solar flares are violent explosions of magnetic energy in the solar
corona. The theory of magnetic reconnection (see Priest & Forbes 2002,
for a review) represents the most plausible
and widely accepted explanation for flare occurrence, although a
definite and clear explanation of the mechanisms at the basis of flare
triggering is still lacking. A well-established property of flare
occurrence is the power law decay of the peak-flux energy
distribution (Lee et al. 1993; Aschwanden et al. 1998; Crosby et al. 1998), a property shared by
other stochastic physical phenomena like earthquakes (de Arcangelis et al. 2006).
Several mechanisms have been proposed to reproduce the above
experimental findings. In the Rosner and Vaiana (RV) model (Rosner
& Vaiana 1978)
it is assumed that flare occurrence is an uncorrelated Poisson process
where the energy grows at a rate proportional to the internal energy of
the system. Given this assumption, the energy storage follows an
exponential temporal growth, interrupted at random when the system
releases all the stored energy in a flare. This mechanism accounts for
the power law in the size distribution and predicts correlations
between the storage time and the released energy. More precisely,
the later a flare will occur the higher the energy will be.
Also the avalanche model (Lu & Hamilton 1991; Hamon et al. 2002),
which describes solar flares as energy relaxation events in a system
driven at a constant rate, correctly predicts scale-free behavior for
the peak-flux energy distribution. Within this approach, flare
occurrence is again a Poisson process, as in the RV model,
but occurrence times and energies are uncorrelated like in standard
self-organized models (Jensen 1998).
Finally, power law behavior for the size distribution is also
consistent with models assuming correlations between the waiting
times
of successive bursts, as in a shell model of magneto-hydrodynamic turbulence (Boffetta et al. 1999). This model reproduces the experimental power law decay
of the waiting time distribution
which can also be obtained however by means of an uncorrelated piecewise Poisson process (Wheatland et al. 2002; Norman et al. 2001).
The investigation of correlations between flare energies and waiting times is a useful tool to distinguish among different triggering mechanisms. This problem has been addressed in a series of papers (Crosby et al. 1998; Wheatland et al. 1998; Wheatland 2000), providing no clear evidence for time-energy correlations. In particular, Wheatland et al. (1998) detected small correlations at short time and interpreted them as a spurious effect of obscuration, i.e. long-lasting events may hide subsequent small events. By properly taking into account obscuration effects, Wheatland concluded that no significant correlations are present. In this paper we present a novel analysis of experimental catalogs providing evidence for time-energy correlations in flare occurrence.
2 Data analysis and methods
In our analysis we used data from the GOES (Geostationary Operational Environmental Satellite) catalog (ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_FLARES)
for the solar cycles 21-23 in the years 1977 to 2007. We
include only C1.4 class flares in the analysis (peak flux energy
larger than E0=1.4
10-6
). Each flare in the catalog is characterized by its starting time ti and its peak-flux energy Ei, in the following called ``energy'' to simplify notation and measured in unit of
,
counting 38 575 flares with
.
In order to enhance the eventual presence of temporal and energy correlations
among flares, we performed a recently proposed (Lippiello et al. 2008a) statistical analysis of the catalog. We computed the probability to have a given quantity xi greater than X conditioned to
being smaller than T,
.
Here,
is the waiting time and the quantity xi represents, depending on cases, the energy of the ith flare xi=Ei, the energy of the previous flare
xi=Ei-1 or the the energy ratio between successive flares,
xi = Ei/Ei-1. The quantity
is given by
,
where N(X,T) is the number of couples of subsequent events fulfilling both xi > X and
,
whereas N(T) is the number of couples satisfying only the condition
.
In order to detect eventual correlations in the catalog, we introduced the probability difference
,
where P(xi>X) is the unconditional probability to have xi>X. A positive (negative)
indicates a higher (lower) probability to find flares with xi>X if one restricts the analysis to couples with
.
Hence, values of
different from zero suggest the existence of correlations between xi and
.
Evidently, because of statistical fluctuations,
is never exactly equal to zero even in catalogs where xi and
are uncorrelated. In order to explicitly take into account the role of statistical fluctuations, we computed the quantity
,
defined as P,
but in a catalog where flare energies were randomly reshuffled. The
comparison with thereshuffled catalog is the basis of the surrogate
data technique (Theiller et al. 1992; Schreiber & Schmitz 2000), recently applied to investigate coherent structures in space plasma (Sahraoui 2008). In the reshuffled catalog xi was by construction uncorrelated to
and
fluctuated around its average value P(xi>X). The amplitude of these fluctuations defined the significance level,
,
which allowed us to distiguish between the presence and absence of
correlations. The method is schematically presented in Fig. 1 for xi=Ei, X=E=3 E0 and T=1 h.
takes different values for each realization of the reshuffled catalog. We produced 105 independent realizations of the catalog with reshuffled energies and observed that
is Gaussian-distributed with a mean P(Ei>E) and a standard deviation
.
Similar results were obtained for other values of E and T and for the other definitions of xi. Therefore,
indicates that energies in the real catalog follow a significantly
different organization than events in the reshuffled catalog.
3 Results
In Fig. 2a we plot
and in Fig. 2b
for different values of T. The standard deviation for each data point
is represented as the error bar. In Fig. 2a we notice that
takes always positive values beyond error bars. This implies that in
the real catalog the number of couples fulfilling both conditions is
greater than in catalogs where energies and intertimes are
uncorrelated. More precisely we find that for each given value of E,
decreases by increasing T. This implies that the probability to find flare couples with the second flare energy higher than E decreases if one includes events with larger
in the analysis. This result is disagrees with what is expected
according to the RV model, which predicts larger flares after
longer waiting times. On the other hand, obscuration can represent a
possible explanation of the above result. Indeed, according to the
selection procedure, for flares close in time, the second event is
recorded in the catalog only if it produces an increase in the flux of
at least
of the level of the previous flare. Moreover, the light curve of large
flares may hide smaller flares close in time. Obscuration effects then
reduce the
probability to find a small flare after a short waiting time and may
introduce spurious correlations between waiting time and the successive
flare energy.
![]() |
Figure 1:
Distribution of
|
Open with DEXTER |
![]() |
Figure 2:
Probability difference
|
Open with DEXTER |
We then explicitly explored the role of obscuration on
.
As a first analysis we considered only flare couples with a temporal distance greater than k times the duration of the first flare. Spurious effects due to obscuration should disappear by increasing k. We observed (Fig. 3a)
that this extra condition did not reduce correlations between waiting
times and energies, since they are still present also for the highest
values of k. To further prove that correlations cannot
be attributed to spurious effects related to obscuration, we performed
the same analysis as in Fig. 1a, setting different lower energy thresholds
.
Indeed, by increasing the value of the minimum energy required for
flares to be included in the analysis, the percentage of flares
hidden by obscuration was reduced and obscuration effects should
progressively decrease. In Fig. 3b we plot the results for different energy thresholds
as a function of
.
No significant dependence on
is detected, indicating that correlations cannot be attributed to
obscuration effects. We conclude that for couples of events close in
time it is more probable to have the second event with a high energy
than for events distant in time.
![]() |
Figure 3:
a) Probability difference
|
Open with DEXTER |
![]() |
Figure 4:
Probability difference
|
Open with DEXTER |
In Fig. 2b we performed the same analysis as in Fig. 2a for the correlations between the waiting time and the previous flare energy. In contrast to Fig. 2a the quantity
is not always positive now. In particular for T=1 h,
is positive at small E and negative at large E. In order to verify if the behavior at small T is affected by obscuration, we performed the same analysis as in Fig. 3b,
imposing the condition on the lower energy threshold
.
In Fig. 4a we present
for T=1 h and different energy thresholds. We found that data strongly depend on
,
in particular for
,
is comparable with the statistical fluctuations. This analysis suggests that the correlations observed at small
are related to obscuration effects. Conversely, data for T=10 h are not significantly affected by
and
always takes positive values beyond error bars (Fig. 4b).
This result indicates that correlations detected for events with
temporal distances between 1 h and 10 h are not significantly
affected by obscuration. We conclude that couples of events distant in
time between 1 h and 10 h tend to have the first event
with a higher energy than the one expected
for a process where time and energy are uncorellated.
To obtain more detailed information on correlations between Ei-1 and ,
we considered
the quantity
,
i.e. the difference between conditional and unconditional probability densities to have Ei-1=E if
.
This quantity can be obtained deriving the cumulated probability
with respect to E. Being a derivative, it is affected by higher fluctuations. For every choice of
,
the value of
was chosen to have the same number of events, about 9000, for all curves. We noticed that for
h the difference in the probability density dp was very close to zero, indicating the absence of correlations.
![]() |
Figure 5:
Probability density difference
|
Open with DEXTER |
We now discuss the results of Fig. 5 at fixed E/E0 and for varying .
We first observed that for small of E, dp is always negative for all values of
.
The negative region is restricted to very small values of
for
and increases to
for larger
.
Within this interval dp is non monotonic in
,
in particular,
it reaches the minimum value for
h and then tends to zero. For
in a given range, negative values of dp imply that it is less probable to have a flare with an energy Ei=E compared to the case where Ei and
are uncorrelated. Results in Fig. 5 then imply that it is improbable to have small
after very small flares and that these anticorrelations vanish for large
.
In particular the highest negative values of dp are obtained for
after an
flare.
If these anticorrelations were a spurious effect of obscuration,
we should rather observe a different behavior. Indeed, energy
reshuffling produces more events shortly after a large flare compared
to the real catalog. Consequently, we should observe an excess, and not
a deficit, of events after a small flare in the real catalog
compared to the reshuffled one because the number of events is
conserved. We now discuss higher values of
.
In this case dp moves from negative values for
and
h to positive values for
h and
h and then tends to zero for higher values of
.
This indicates that there is a deficit of events with
h after energetic flares, whereas for
the number of observed flares is larger than the one expected for an
uncorrelated process. The above results suggest that after a large
flare a recovery time of about 2 h takes place, when only a
small number of flares is observed. After this time the number of
observed flares reaches a maximum value for
h and then decreases to the background level for larger
.
According to the dependence on
observed in Fig. 4a, the recovery time can be a spurious effect related to obscuration.
The number of events n(t) occurring at the time t after an energetic flare can be explicitly evaluated following the method of Lippiello (2008b). More precisely, we define as a ``main'' flare an event with the energy
.
Indicating with ti the occurrence times of main-flares we compute the quantity
,
where
is the Heaviside step function, and the sum extends
over all main-flares. Here n(t-ti) is the number of flares with an energy lower than 10E0
occurring during the time t-ti after the ith main flare, and the sum extends over all main-flares in the catalog. Assuming that n(t-ti) is time translationally invariant we obtain
which represents the average number of flares occurring in a period of a duration t after a mainflare. The data (Fig. 6a) agree with the previous results, with a maximum number of flares at


![]() |
Figure 6:
a) Average number of flares n(t) as a function of the time t occurring after a main flare with an energy
|
Open with DEXTER |
We now discuss the existence of correlations between energies of
subsequent flares. In order to do so, we computed the energy
distribution
of the first m flares occurring after a flare with a larger energy than a reference value Er. Since the analysis is restricted to all flares with
,
obviously
coincides with the flare energy distribution of the whole catalog
,
for Er=E0. In Fig. 6b we plot
for m=20 and different values of Er. Similar results are obtained for other values of m. Deviations from
become more and more evident for increasing values of Er. In particular, the higher Er is,
the larger is the probability to have subsequent flares with higher
energy and the lower is the probability to have small flares.
In order to get further insights in the correlations between flare energies, we computed the quantity
,
i.e. the difference between the conditional and the unconditional
probability density to have the energy of the subsequent flare
times the previous one. We found (Fig. 7a) that for T=1 h,
is significantly different from zero for all
values.
This indicates that the energies of two subsequent flares are
correlated. Moreover we observed that these correlations depend on the
time separation between the two flares and are practically zero for
h (Fig. 7a). We verified that this result is not affected by obscuration performing the same analysis for T=1 h and different lower energy thresholds
(Fig. 7b). We found that curves for different
coincide within statistical fluctuations. This indicates that energy
correlations and their dependence on time separation are a physical
property, not a spurious effect due to obscuration. Curves present
a maximum for
,
indicating that it is more probable to find the next flare with an
energy close to but slightly higher than the previous one.
![]() |
Figure 7:
Probability density difference
|
Open with DEXTER |
4 Conclusions
In conclusion, we presented a statistical analysis of the GOES catalog indicating the existence of time-energy correlations between successive events not to be attributed to obscuration effects. More precisely, we observed that for couples of events close in time (T<1 h), the second event tends to have a high energy. Moreover couples of events distant in time between 1 h and 10 h have the first event with a high energy. The analysis of the rate decay after large flares shows evidence that the largest number of events is detected about 4 h after the occurrence of the main event. Finally the distribution of flare energies confirms that the higher the flare energy, the larger the number of subsequent events with high energy. The existence of time-energy correlations suggests the possibility of scaling laws relating time with the energy released in a flare. This is a still open question, which could provide interesting insights in the energy storage and release mechanisms at the origin of solar flare occurrence.
References
- Aschwanden, M. J., Dennis, B. R., & Benz, A. O. 1998, ApJ, 497, 972 [NASA ADS] [CrossRef] [Google Scholar]
- Boffetta, G., Carbone, V., Giuliani, P., Veltri, P., & Vulpiani, A. 1999, Phys. Rev. Lett., 83, 4662 [NASA ADS] [CrossRef] [Google Scholar]
- Crosby, N., Vilmer, N., Lund, N., & Sunyaev, R. 1998, A&A, 334, 299 [NASA ADS] [Google Scholar]
- de Arcangelis, L., Godano, C., Lippiello, E., & Nicodemi, M. 2006, Phys. Rev. Lett., 96, 051102 [Google Scholar]
- Hamon, D., Nicodemi, M., & Jensen, H. J. 2002, A&A, 387, 326 [Google Scholar]
- Jensen, H. J. 1998, Self Organized Criticality (Cambridge University Press) [Google Scholar]
- Lee, T. T., Petrosian, V., & McTiernan, J. M. 1993, ApJ, 412, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Lippiello, E., de Arcangelis, L., & Godano, C. 2008a, Phys. Rev. Lett., 100, 038501 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Lippiello, E., de Arcangelis, L., & Godano, C. 2008b, A&A, 488, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lu, E. T., & Hamilton, R. J. 1991, ApJ, 380, L89 [NASA ADS] [CrossRef] [Google Scholar]
- Norman, J. P., Charbonneau, P., McIntosh, S. W., & Liu, H. 2001, ApJ, 557, 891 [NASA ADS] [CrossRef] [Google Scholar]
- Omori, F. 1894, J. Coll. Sci. Imp. Univ. Tokyo, 7, 111 [Google Scholar]
- Priest, E. R., & Forbes, T. G. 2002, A&ARv, 10, 313 [NASA ADS] [CrossRef] [Google Scholar]
- Rosner, R., & Vaiana, G. S. 1978, ApJ, 222, 1104 [NASA ADS] [CrossRef] [Google Scholar]
- Sahraoui, F. 2008, Phys. Rev. E, 78, 026402 [NASA ADS] [CrossRef] [Google Scholar]
- Schreiber, T., & Schmitz, A. 2000, Physica D, 142, 346 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Theiler, J., Eubank, S., Longtin, A., Galdrikian, B., & Doyne Farmer, J. 1992, Physica D, 58, 77 [NASA ADS] [CrossRef] [Google Scholar]
- Wheatland, M. S. 2000, Sol. Phys., 191, 381 [Google Scholar]
- Wheatland, M. S., & Litvinenko, Y. E. 2002, Sol. Phys., 211, 255 [NASA ADS] [CrossRef] [Google Scholar]
- Wheatland, M. S., Sturrock, P. A., & McTiernan, J. M. 1998, ApJ, 509, 448 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
![]() |
Figure 1:
Distribution of
|
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Probability difference
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
a) Probability difference
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Probability difference
|
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Probability density difference
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
a) Average number of flares n(t) as a function of the time t occurring after a main flare with an energy
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
Probability density difference
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.