-
The Flare Likelihood and Region Eruption Forecasting (FLARECAST) Project: Flare forecasting in the big data & machine learning era
Authors:
M. K. Georgoulis,
D. S. Bloomfield,
M. Piana,
A. M. Massone,
M. Soldati,
P. T. Gallagher,
E. Pariat,
N. Vilmer,
E. Buchlin,
F. Baudin,
A. Csillaghy,
H. Sathiapal,
D. R. Jackson,
P. Alingery,
F. Benvenuto,
C. Campi,
K. Florios,
C. Gontikakis,
C. Guennou,
J. A. Guerra,
I. Kontogiannis,
V. Latorre,
S. A. Murray,
S. -H. Park,
S. von Stachelski
, et al. (3 additional authors not shown)
Abstract:
The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) differe…
▽ More
The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) different ML techniques, also on equal footing, to optimize the immense Big Data parameter space created by these many predictors; third, the establishment of a robust, three-pronged communication effort oriented toward policy makers, space-weather stakeholders and the wider public. FLARECAST pledged to make all its data, codes and infrastructure openly available worldwide. The combined use of 170+ properties (a total of 209 predictors are now available) in multiple ML algorithms, some of which were designed exclusively for the project, gave rise to changing sets of best-performing predictors for the forecasting of different flaring levels. At the same time, FLARECAST reaffirmed the importance of rigorous training and testing practices to avoid overly optimistic pre-operational prediction performance. In addition, the project has (a) tested new and revisited physically intuitive flare predictors and (b) provided meaningful clues toward the transition from flares to eruptive flares, namely, events associated with coronal mass ejections (CMEs). These leads, along with the FLARECAST data, algorithms and infrastructure, could help facilitate integrated space-weather forecasting efforts that take steps to avoid effort duplication. In spite of being one of the most intensive and systematic flare forecasting efforts to-date, FLARECAST has not managed to convincingly lift the barrier of stochasticity in solar flare occurrence and forecasting: solar flare prediction thus remains inherently probabilistic.
△ Less
Submitted 12 May, 2021;
originally announced May 2021.
-
Differential Emission Measure Evolution as a Precursor of Solar Flares
Authors:
C. Gontikakis,
I. Kontogiannis,
M. K. Georgoulis,
C. Guennou,
P. Syntelis,
S. H. Park,
E. Buchlin
Abstract:
We analyse the temporal evolution of the Differential Emission Measure (DEM) of solar active regions and explore its usage in solar flare prediction. The DEM maps are provided by the Gaussian Atmospheric Imaging Assembly (GAIA-DEM) archive, calculated assuming a Gaussian dependence of the DEM on the logarithmic temperature. We analyse time-series of sixteen solar active regions and a statistically…
▽ More
We analyse the temporal evolution of the Differential Emission Measure (DEM) of solar active regions and explore its usage in solar flare prediction. The DEM maps are provided by the Gaussian Atmospheric Imaging Assembly (GAIA-DEM) archive, calculated assuming a Gaussian dependence of the DEM on the logarithmic temperature. We analyse time-series of sixteen solar active regions and a statistically significant sample of 9454 point-in-time observations corresponding to hundreds of regions observed during solar cycle 24. The time-series analysis shows that the temporal derivatives of the Emission Measure dEM/dt and the maximum DEM temperature dTmax/dt frequently exhibit high positive values a few hours before M- and X-class flares, indicating that flaring regions become brighter and hotter as the flare onset approaches. From the point-in-time observations we compute the conditional probabilities of flare occurrences using the distributions of positive values of the dEM/dt, and dTmax/dt and compare them with corresponding flaring probabilities of the total unsigned magnetic flux, a conventionally used, standard flare predictor. For C-class flares, conditional probabilities have lower or similar values with the ones derived for the unsigned magnetic flux, for 24 and 12 hours forecast windows. For M- and X-class flares, these probabilities are higher than those of the unsigned flux for higher parameter values. Shorter forecast windows improve the conditional probabilities of dEM/dt, and dTmax/dt in comparison to those of the unsigned magnetic flux. We conclude that flare forerunner events such as preflare heating or small flare activity prior to major flares reflect on the temporal evolution of EM and Tmax. Of these two, the temporal derivative of the EM could conceivably be used as a credible precursor, or short-term predictor, of an imminent flare.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
Testing predictors of eruptivity using parametric flux emergence simulations
Authors:
Chloé Guennou,
Etienne Pariat,
Nicole Vilmer,
James E. Leake
Abstract:
Solar flares and coronal mass ejections (CMEs) are among the most energetic events in the solar system, impacting the near-Earth environment. Flare productivity is empirically known to be correlated with the size and complexity of active regions. Several indicators, based on magnetic-field data from active regions, have been tested for flare forecasting in recent years. None of these indicators, o…
▽ More
Solar flares and coronal mass ejections (CMEs) are among the most energetic events in the solar system, impacting the near-Earth environment. Flare productivity is empirically known to be correlated with the size and complexity of active regions. Several indicators, based on magnetic-field data from active regions, have been tested for flare forecasting in recent years. None of these indicators, or combinations thereof, have yet demonstrated an unambiguous eruption or flare criterion. Furthermore, numerical simulations have been only barely used to test the predictability of these parameters. In this context, we used the 3D parametric MHD numerical simulations of the self-consistent formation of the flux emergence of a twisted flux tube, inducing the formation of stable and unstable magnetic flux ropes of Leake (2013, 2014). We use these numerical simulations to investigate the eruptive signatures observable in various magnetic scalar parameters and provide highlights on data analysis processing. Time series of 2D photospheric-like magnetograms are used from parametric simulations of stable and unstable flux emergence, to compute a list of about 100 different indicators. This list includes parameters previously used for operational forecasting, physical parameters used for the first time, as well as new quantities specifically developed for this purpose. Our results indicate that only parameters measuring the total non-potentiality of active regions associated with magnetic inversion line properties, such as the Falconer parameters $L_{ss}$, $WL_{ss}$, $L_{sg}$ and $WL_{sg}$, as well as the new current integral $WL_{sc}$ and length $L_{sc}$ parameters, present a significant ability to distinguish the eruptive cases of the model from the non-eruptive cases, possibly indicating that they are promising flare and eruption predictors.
△ Less
Submitted 15 June, 2017;
originally announced June 2017.
-
Energetic characterisation and statistics of solar coronal brightenings
Authors:
Vincent Joulin,
Éric Buchlin,
Jacques Solomon,
Chloé Guennou
Abstract:
To explain the high temperature of the corona, much attention has been paid to the distribution of energy in dissipation events. Indeed, if the event energy distribution is steep enough, the smallest, unobservable events could be the largest contributors to the total energy dissipation in the corona. Previous observations have shown a wide distribution of energies but remain inconclusive about the…
▽ More
To explain the high temperature of the corona, much attention has been paid to the distribution of energy in dissipation events. Indeed, if the event energy distribution is steep enough, the smallest, unobservable events could be the largest contributors to the total energy dissipation in the corona. Previous observations have shown a wide distribution of energies but remain inconclusive about the precise slope. Furthermore, these results rely on a very crude estimate of the energy. On the other hand, more detailed spectroscopic studies of structures such as coronal bright points do not provide enough statistical information to derive their total contribution to heating. We aim at getting a better estimate of the distributions of the energy dissipated in coronal heating events using high-resolution, multi-channel Extreme Ultra-Violet (EUV) data. To estimate the energies corresponding to heating events and deduce their distribution, we detect brightenings in five EUV channels of the Atmospheric Imaging Assembly (AIA) on-board the Solar Dynamics Observatory (SDO). We combine the results of these detections and we use maps of temperature and emission measure derived from the same observations to compute the energies. We obtain distributions of areas, durations, intensities, and energies (thermal, radiative, and conductive) of events. These distributions are power-laws, and we find also power-law correlations between event parameters. The energy distributions indicate that the energy from a population of events like the ones we detect represents a small contribution to the total coronal heating, even when extrapolating to smaller scales. The main explanations for this are how heating events can be extracted from observational data, and the incomplete knowledge of the thermal structure and processes in the coronal plasma attainable from available observations.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Evidence for evaporation-incomplete condensation cycles in warm solar coronal loops
Authors:
Clara Froment,
Frédéric Auchère,
Karine Bocchialini,
Eric Buchlin,
Chloé Guennou,
Jacques Solomon
Abstract:
Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the plasma: thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of prominences and coronal rain. However, it is often discarded to be involved in the heating of warm loops as the models do not reproduce observations. Recent simulations have shown that these incon…
▽ More
Quasi-constant heating at the footpoints of loops leads to evaporation and condensation cycles of the plasma: thermal non-equilibrium (TNE). This phenomenon is believed to play a role in the formation of prominences and coronal rain. However, it is often discarded to be involved in the heating of warm loops as the models do not reproduce observations. Recent simulations have shown that these inconsistencies with observations may be due to oversimplifications of the geometries of the models. In addition, our recent observations reveal that long-period intensity pulsations (several hours) are common in solar coronal loops. These periods are consistent with those expected from TNE. The aim of this paper is to derive characteristic physical properties of the plasma for some of these events to test the potential role of TNE in loop heating. We analyzed three events in detail using the six EUV coronal channels of SDO/AIA. We performed both a Differential Emission Measure (DEM) and a time-lag analysis, including a new method to isolate the relevant signal from the foreground and background emission. For the three events, the DEM undergoes long-period pulsations, which is a signature of periodic heating even though the loops are captured in their cooling phase, as is the bulk of the active regions. We link long-period intensity pulsations to new signatures of loop heating with strong evidence for evaporation and condensation cycles. We thus witness simultaneously widespread cooling and TNE. Finally, we discuss the implications of our new observations for both static and impulsive heating models.
△ Less
Submitted 26 August, 2015; v1 submitted 30 April, 2015;
originally announced April 2015.
-
Relative Abundance Measurements in Plumes and Interplumes
Authors:
Chloé Guennou,
Michael Hahn,
Daniel Wolf Savin
Abstract:
We present measurements of relative elemental abundances in plumes and interplumes. Plumes are bright, narrow structures in coronal holes that extend along open magnetic field lines far out into the corona. Previous work has found that in some coronal structures the abundances of elements with a low first ionization potential (FIP) < 10 eV are enhanced relative to their photospheric abundances. Th…
▽ More
We present measurements of relative elemental abundances in plumes and interplumes. Plumes are bright, narrow structures in coronal holes that extend along open magnetic field lines far out into the corona. Previous work has found that in some coronal structures the abundances of elements with a low first ionization potential (FIP) < 10 eV are enhanced relative to their photospheric abundances. This coronal-to-photospheric abundance ratio, commonly called the FIP bias, is typically 1 for element with a high-FIP (> 10 eV). We have used EIS spectroscopic observations made on 2007 March 13 and 14 over an ~24 hour period to characterize abundance variations in plumes and interplumes. To assess their elemental composition, we have used a differential emission measure (DEM) analysis, which accounts for the thermal structure of the observed plasma. We have used lines from ions of iron, silicon, and sulfur. From these we have estimated the ratio of the iron and silicon FIP bias relative to that for sulfur. From the results, we have created FIP-bias-ratio maps. We find that the FIP-bias ratio is sometimes higher in plumes than in interplumes and that this enhancement can be time dependent. These results may help to identify whether plumes or interplumes contribute to the fast solar wind observed in situ and may also provides constraints on the formation and heating mechanisms of plumes.
△ Less
Submitted 20 March, 2015;
originally announced March 2015.
-
Can the Differential Emission Measure constrain the timescale of energy deposition
Authors:
Chloé Guennou,
Frédéric Auchère,
James A. Klimchuk,
Karine Bocchialini,
Susanna Parenti
Abstract:
In this paper, the ability of the Hinode/EIS instrument to detect radiative signatures of coronal heating is investigated. Recent observational studies of AR cores suggest that both the low and high frequency heating mechanisms are consistent with observations. Distinguishing between these possibilities is important for identifying the physical mechanism(s) of the heating. The Differential Emissio…
▽ More
In this paper, the ability of the Hinode/EIS instrument to detect radiative signatures of coronal heating is investigated. Recent observational studies of AR cores suggest that both the low and high frequency heating mechanisms are consistent with observations. Distinguishing between these possibilities is important for identifying the physical mechanism(s) of the heating. The Differential Emission Measure (DEM) tool is one diagnostic that allows to make this distinction, through the amplitude of the DEM slope coolward of the coronal peak. It is therefore crucial to understand the uncertainties associated with these measurements. Using proper estimations of the uncertainties involved in the problem of DEM inversion, we derive confidence levels on the observed DEM slope. Results show that the uncertainty in the slope reconstruction strongly depends on the number of lines constraining the slope. Typical uncertainty is estimated to be about $\pm 1.0$, in the more favorable cases.
△ Less
Submitted 13 June, 2013;
originally announced June 2013.
-
On the Accuracy of the Differential Emission Measure Diagnostics of Solar Plasmas. Application to AIA/SDO. Part I: Isothermal plasmas
Authors:
Chloé Guennou,
Frédéric Auchère,
Elie Soubrié,
Karine Bocchialini,
Susanna Parenti
Abstract:
DEM analysis is a major diagnostic tool for stellar atmospheres. But both its derivation and its interpretation are notably difficult because of random and systematic errors, and the inverse nature of the problem. We use simulations with simple thermal distributions to investigate the inversion properties of SDO/AIA observations of the solar corona. This allows a systematic exploration of the para…
▽ More
DEM analysis is a major diagnostic tool for stellar atmospheres. But both its derivation and its interpretation are notably difficult because of random and systematic errors, and the inverse nature of the problem. We use simulations with simple thermal distributions to investigate the inversion properties of SDO/AIA observations of the solar corona. This allows a systematic exploration of the parameter space and using a statistical approach, the respective probabilities of all the DEMs compatible with the uncertainties can be computed. Following this methodology, several important properties of the DEM inversion, including new limitations, can be derived and presented in a very synthetic fashion.
In this first paper, we describe the formalism and we focus on isothermal plasmas, as building blocks to understand the more complex DEMs studied in the second paper. The behavior of the inversion of AIA data being thus quantified, and we provide new tools to properly interpret the DEM. We quantify the improvement of the isothermal inversion with 6 AIA bands compared to previous EUV imagers. The maximum temperature resolution of AIA is found to be 0.03 log Te, and we derive a rigorous test to quantify the compatibility of observations with the isothermal hypothesis. However we demonstrate limitations in the ability of AIA alone to distinguish different physical conditions.
△ Less
Submitted 8 October, 2012;
originally announced October 2012.
-
On the Accuracy of the Differential Emission Measure Diagnostics of Solar Plasmas. Application to AIA/SDO. Part II: Multithermal plasmas
Authors:
Chloé Guennou,
Frédéric Auchère,
Elie Soubrié,
Karine Bocchialini,
Susanna Parenti
Abstract:
The Differential Emission Measure (DEM) analysis is one of the most used diagnostic tools for solar and stellar coronae. Being an inverse problem, it has limitations due to the presence of random and systematic errors. We present in theses series of papers an analysis of the robustness of the inversion in the case of AIA/SDO observations. We completely characterize the DEM inversion and its statis…
▽ More
The Differential Emission Measure (DEM) analysis is one of the most used diagnostic tools for solar and stellar coronae. Being an inverse problem, it has limitations due to the presence of random and systematic errors. We present in theses series of papers an analysis of the robustness of the inversion in the case of AIA/SDO observations. We completely characterize the DEM inversion and its statistical properties, providing all the solutions consistent with the data along with their associated probabilities, and a test of the suitability of the assumed DEM model.
While Paper I focused on isothermal conditions, we now consider multi-thermal plasmas and investigate both isothermal and multithermal solutions. We demonstrate how the ambiguity between noises and multi-thermality fundamentally limits the temperature resolution of the inversion. We show that if the observed plasma is multi-thermal, isothermal solutions tend to cluster on a constant temperature whatever the number of passbands or spectral lines. The multi-thermal solutions are also found to be biased toward near isothermal solutions around 1 MK. This is true even if the residuals support the chosen DEM model, possibly leading to erroneous conclusions on the observed plasma. We propose tools to identify and quantify the possible degeneracy of solutions, thus helping the interpretation of DEM inversion.
△ Less
Submitted 8 October, 2012;
originally announced October 2012.
-
Coronal Temperature Maps from Solar EUV images: a Blind Source Separation Approach
Authors:
T. Dudok de Wit,
S. Moussaoui,
C. Guennou,
F. Auchère,
G. Cessateur,
M. Kretzschmar,
L. A. Vieira,
F. F. Goryaev
Abstract:
Multi-wavelength solar images in the EUV are routinely used for analysing solar features such as coronal holes, filaments, and flares. However, images taken in different bands often look remarkably similar as each band receives contributions coming from regions with a range of different temperatures. This has motivated the search for empirical techniques that may unmix these contributions and conc…
▽ More
Multi-wavelength solar images in the EUV are routinely used for analysing solar features such as coronal holes, filaments, and flares. However, images taken in different bands often look remarkably similar as each band receives contributions coming from regions with a range of different temperatures. This has motivated the search for empirical techniques that may unmix these contributions and concentrate salient morphological features of the corona in a smaller set of less redundant source images. Blind Source Separation (BSS) precisely does this. Here we show how this novel concept also provides new insight into the physics of the solar corona, using observations made by SDO/AIA. The source images are extracted using a Bayesian positive source separation technique. We show how observations made in six spectral bands, corresponding to optically thin emissions, can be reconstructed by linear combination of three sources. These sources have a narrower temperature response and allow for considerable data reduction since the pertinent information from all six bands can be condensed in only one single composite picture. In addition, they give access to empirical temperature maps of the corona. The limitations of the BSS technique and some applications are briefly discussed.
△ Less
Submitted 1 March, 2012;
originally announced March 2012.
-
TomograPy: A Fast, Instrument-Independent, Solar Tomography Software
Authors:
Nicolas Barbey,
Chloé Guennou,
Frédéric Auchère
Abstract:
Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handl…
▽ More
Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handling of the temporal evolution of the observed structures, and the heterogeneity of the data in multi-spacecraft studies.
We present a generic software package that can perform fast tomographic inversions that scales linearly with the number of measurements, linearly with the length of the reconstruction cube (and not the number of voxels) and linearly with the number of cores and can use data from different sources and with a variety of physical models: TomograPy (http://nbarbey.github.com/TomograPy/), an open-source software freely available on the Python Package Index. For performance, TomograPy uses a parallelized-projection algorithm. It relies on the World Coordinate System standard to manage various data sources. A variety of inversion algorithms are provided to perform the tomographic-map estimation. A test suite is provided along with the code to ensure software quality. Since it makes use of the Siddon algorithm it is restricted to rectangular parallelepiped voxels but the spherical geometry of the corona can be handled through proper use of priors.
We describe the main features of the code and show three practical examples of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1 data. Static and smoothly varying temporal evolution models are presented.
△ Less
Submitted 30 March, 2011;
originally announced March 2011.