Peter J. Cargill

The Enthalpy Based Thermal Evolution of Loops approximate model for static and dynamic coronal loops is developed to include the effect of a loop cross-sectional area which increases from the base of the transition region (TR) to the corona. The TR is defined as the part of a loop between the top of the chromosphere and the location where thermal conduction changes from an energy loss to an energy gain. There are significant differences from constant area loops due to the manner in which the reduced volume of the TR responds to conductive and enthalpy fluxes from the corona. For static loops with modest area variation the standard picture of loop energy balance is retained, with the corona and TR being primarily a balance between heating and conductive losses in the corona, and downward conduction and radiation to space in the TR. As the area at the loop apex increases, the TR becomes thicker and the density in TR and corona larger. For large apex areas, the coronal energy balance changes to one primarily between heating and radiation, with conduction playing an increasingly unimportant role, and the TR thickness becoming a significant fraction of the loop length. Approximate scaling laws are derived that give agreement with full numerical solutions for the density, but not the temperature. For non-uniform areas, dynamic loops have a higher peak temperature and are denser in the radiative cooling phase by of order 50 per cent than the constant area case for the examples considered. They also show a final rapid cooling and draining once the temperature approaches 1 MK. Although the magnitude of the emission measure will be enhanced in the radiative phase, there is little change in the important observational diagnostic of its temperature dependence.

Despite its prediction over two decades ago, the detection of faint, high-temperature (“hot”) emission due to nanoflare heating in non-flaring active region cores has proved challenging. Using an efficient two-fluid hydrodynamic model, this paper investigates the properties of the emission expected from repeating nanoflares (a nanoflare train) of varying frequency as well as the separate heating of electrons and ions. If the emission measure distribution (EM(T)) peaks at T = T _m , we find that EM(T _m ) is independent of details of the nanoflare train, and EM(T) above and below T _m reflects different aspects of the heating. Below T _m , the main influence is the relationship of the waiting time between successive nanoflares to the nanoflare energy. Above T _m , power-law nanoflare distributions lead to an extensive plasma population not present in a mono-energetic train. Furthermore, in some cases, characteristic features are present in EM(T). Such details may be detectable given adequate spectral resolution and a good knowledge of the relevant atomic physics. In the absence of such resolution we propose some metrics that can be used to infer the presence of “hot” plasma.

It has long been established that gradients in the Alfvén speed, and in particular the plasma density, are an essential part of the damping of waves in the magnetically closed solar corona by mechanisms such as resonant absorption and phase mixing. While models of wave damping often assume a fixed density gradient, in this paper the self-consistency of such calculations is assessed by examining the temporal evolution of the coronal density. It is shown conceptually that for some coronal structures, density gradients can evolve in a way that the wave-damping processes are inhibited. For the case of phase mixing we argue that (a) wave heating cannot sustain the assumed density structure and (b) inclusion of feedback of the heating on the density gradient can lead to a highly structured density, although on long timescales. In addition, transport coefficients well in excess of classical are required to maintain the observed coronal density. Hence, the heating of closed coronal structures by global oscillations may face problems arising from the assumption of a fixed density gradient, and the rapid damping of oscillations may have to be accompanied by a separate (non-wave-based) heating mechanism to sustain the required density structuring.

Recent observations from the Hinode and Solar Dynamics Observatory spacecraft have provided major advances in understanding the heating of solar active regions (ARs). For ARs comprising many magnetic strands or sub-loops heated by small, impulsive events (nanoflares), it is suggested that (i) the time between individual nanoflares in a magnetic strand is 500–2000 s, (ii) a weak ‘hot’ component (more than 10 ^6.6 K) is present, and (iii) nanoflare energies may be as low as a few 10 ²³ ergs. These imply small heating events in a stressed coronal magnetic field, where the time between individual nanoflares on a strand is of order the cooling time. Modelling suggests that the observed properties are incompatible with nanoflare models that require long energy build-up (over 10 s of thousands of seconds) and with steady heating.

We examine the radiative cooling of coronal loops and demonstrate that the recently identified catastrophic cooling (Reale and Landi, 2012) is due to the inability of a loop to sustain radiative / enthalpy cooling below a critical temperature, which can be > 1 MK in flares, 0.5 - 1 MK in active regions and 0.1 MK in long tenuous loops. Catastrophic cooling is characterised by a rapid fall in coronal temperature while the coronal density changes by a small amount. Analytic expressions for the critical temperature are derived and show good agreement with numerical results. This effect limits very considerably the lifetime of coronal plasmas below the critical temperature.

The effect of the numerical spatial resolution in models of the solar corona and corona / chromosphere interface is examined for impulsive heating over a range of magnitudes using one dimensional hydrodynamic simulations. It is demonstrated that the principle effect of inadequate resolution is on the coronal density. An underresolved loop typically has a peak density of at least a factor of two lower than a resolved loop subject to the same heating, with larger discrepencies in the decay phase. The temperature for under-resolved loops is also lower indicating that lack of resolution does not "bottle up" the heat flux in the corona. Energy is conserved in the models to under 1% in all cases, indicating that this is not responsible for the low density. Instead, we argue that in under-resolved loops the heat flux "jumps across" the transition region to the dense chromosphere from which it is radiated rather than heating and ablating transition region plasma. This emphasises the point that the interaction between corona and chromosphere occurs only through the medium of the transition region. Implications for three dimensional magnetohydrodynamic coronal models are discussed.