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.

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.

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.

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.

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.

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.

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.

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.

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.