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 therm
al 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.
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.
