This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

Magnetic reconnection, a fundamentally important process in many aspects of astrophysics, is believed to be initiated by the tearing instability of an electric current sheet, a region where magnetic field abruptly changes direction and electric currents build up. Recent studies have suggested that the amount of magnetic shear in these structures is a critical parameter for the switch-on nature of magnetic reconnection in the solar atmosphere, at fluid spatial scales much larger than kinetic scales. We present results of simulations of reconnection in 3D current sheets with conditions appropriate to the solar corona. Using high-fidelity simulations, we follow the evolution of the linear and non-linear 3D tearing instability, leading to reconnection. We find that, depending on the parameter space, magnetic shear can play a vital role in the onset of significant energy release and heating via non-linear tearing. Two regimes in our study exist, dependent on whether the current sheet is longer or shorter than the wavelength of the fastest growing parallel mode (in the corresponding infinite system), thus determining whether sub-harmonics are present in the actual system. In one regime, where the fastest growing parallel mode has sub-harmonics, the non-linear interaction of these sub-harmonics and the coalescence of 3D plasmoids dominates the non-linear evolution, with magnetic shear playing only a weak role in the amount of energy released. In the second regime, where the fastest growing parallel mode has no-sub-harmonics, then only strongly sheared current sheets, where oblique mode are strong enough to compete with the dominant parallel mode, show any significant energy release. We expect both regimes to exist on the Sun, and so our results have important consequences for the the question of reconnection onset in different solar physics applications.

For some forms of steady heating, coronal loops are in a state of thermal nonequilibrium and evolve in a manner that includes accelerated cooling, often resulting in the formation of a cold condensation. This is frequently confused with thermal instability, but the two are in fact fundamentally different. We explain the distinction and discuss situations where they may be interconnected. Large-amplitude perturbations, perhaps associated with MHD waves, likely play a role in explaining phenomena that have been attributed to thermal nonequilibrium but also seem to require cross-field communication.

A comprehensive understanding of the structure of Doppler motions in transition region including the center-to-limb variation and its relationship with the magnetic field structure is vital for the understanding of mass and energy transfer in the solar atmosphere. In this paper, we have performed such a study in an active region using the Si IV 1394~{\AA} emission line recorded by the Interface Region Imaging Spectrograph (IRIS) and the line-of-sight photospheric magnetic field obtained by the Helioseismic and Magnetic Imager (HMI) on-board the Solar Dynamics Observatory (SDO). The active region has two opposite polarity strong field regions separated by a weak field corridor, which widened as the active region evolved. On average the strong field regions (corridor) show(s) redshifts of 5{--}10 (3{--}9)~km~s$^{-1}$ (depending on the date of observation). There is, however, a narrow lane in the middle of the corridor with near-zero Doppler shifts at all disk positions, suggesting that any flows there are very slow. The Doppler velocity distributions in the corridor seem to have two components---a low velocity component centered near 0 km/s and a high velocity component centered near 10~km~s$^{-1}$. The high velocity component is similar to the velocity distributions in the strong field regions, which have just one component. Both exhibit a small center-to limb variation and seem to come from the same population of flows. To explain these results, we suggest that the emission from the lower transition region comes primarily from warm type II spicules, and we introduce the idea of a `chromospheric wall'---associated with classical cold spicules---to account for a diminished center-to-limb variation.

One of the greatest challenges in solar physics is understanding the heating of the Sun's corona. Most theories for coronal heating postulate that free energy in the form of magnetic twist/stress is injected by the photosphere into the corona where the free energy is converted into heat either through reconnection or wave dissipation. The magnetic helicity associated with the twist/stress, however, is expected to be conserved and appear in the corona. In previous work we showed that helicity associated with the small-scale twists undergoes an inverse cascade via stochastic reconnection in the corona, and ends up as the observed large-scale shear of filament channels. Our ``helicity condensation'' model accounts for both the formation of filament channels and the observed smooth, laminar structure of coronal loops. In this paper, we demonstrate, using helicity- and energy-conserving numerical simulations of a coronal system driven by photospheric motions, that the model also provides a natural mechanism for heating the corona. We show that the heat generated by the reconnection responsible for the helicity condensation process is sufficient to account for the observed coronal heating. We study the role that helicity injection plays in determining coronal heating and find that, crucially, the heating rate is only weakly dependent on the net helicity preference of the photospheric driving. Our calculations demonstrate that motions with 100\% helicity preference are least efficient at heating the corona; those with 0\% preference are most efficient. We discuss the physical origins of this result and its implications for the observed corona.

Thermal non-equilibrium (TNE) is a fascinating situation that occurs in coronal magnetic flux tubes (loops) for which no solution to the steady-state fluid equations exists. The plasma is constantly evolving even though the heating that produces the hot temperatures does not. This is a promising explanation for isolated phenomena such as prominences, coronal rain, and long-period pulsating loops, but it may also have much broader relevance. As known for some time, TNE requires that the heating be both (quasi) steady and concentrated at low coronal altitudes. Recent studies indicate that asymmetries are also important, with large enough asymmetries in the heating and/or cross-sectional area resulting in steady flow rather than TNE. Using reasonable approximations, we have derived two formulae for quantifying the conditions necessary for TNE. As a rough rule of thumb, the ratio of apex to footpoint heating rates must be less than about 0.1, and asymmetries must be less than about a factor of 3. The precise values are case dependent. We have tested our formulae with 1D hydrodynamic loop simulations and find a very acceptable agreement. These results are important for developing physical insight about TNE and assessing how widespread it may be on the Sun.

Much evidence suggests that the solar corona is heated impulsively, meaning that nanoflares may be ubiquitous in quiet and active regions (ARs). Hard X-ray (HXR) observations with unprecedented sensitivity $>$3~keV are now enabled by focusing instruments. We analyzed data from the \textit{Focusing Optics X-ray Solar Imager (FOXSI)} rocket and the \textit{Nuclear Spectroscopic Telescope Array (NuSTAR)} spacecraft to constrain properties of AR nanoflares simulated by the EBTEL field-line-averaged hydrodynamics code. We generated model X-ray spectra by computing differential emission measures for homogeneous nanoflare sequences with heating amplitudes $H_0$, durations $\tau$, delay times between events $t_N$, and filling factors $f$. The single quiescent AR observed by \textit{FOXSI-2} on 2014 December 11 is well fit by nanoflare sequences with heating amplitudes 0.02 erg cm$^{-3}$ s$^{-1}$ $<$ $H_0$ $<$ 13 erg cm$^{-3}$ s$^{-1}$ and a wide range of delay times and durations. We exclude delays between events shorter than $\sim$900 s at the 90\% confidence level for this region. Three of five regions observed by {\nustar} on 2014 November 1 are well fit by homogeneous nanoflare models, while two regions with higher fluxes are not. Generally, the {\nustar} count spectra are well fit by nanoflare sequences with smaller heating amplitudes, shorter delays, and shorter durations than the allowed \textit{FOXSI-2} models. These apparent discrepancies are likely due to differences in spectral coverage between the two instruments and intrinsic differences among the regions. Steady heating ($t_N$ = $\tau$) was ruled out with $>$99\% confidence for all regions observed by either instrument.

We report on high-resolution imaging and spectral observations of eruptions of a spiral structure in the transition region, which were taken with the Interface Region Imaging Spectrometer (IRIS), the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI). The eruption coincided with the appearance of two series of jets, with velocities comparable to the Alfv\'en speeds in their footpoints. Several pieces of evidence of magnetic braiding in the eruption are revealed, including localized bright knots, multiple well-separated jet threads, transition region explosive events and the fact that all these three are falling into the same locations within the eruptive structures. Through analysis of the extrapolated three-dimensional magnetic field in the region, we found that the eruptive spiral structure corresponded well to locations of twisted magnetic flux tubes with varying curl values along their lengths. The eruption occurred where strong parallel currents, high squashing factors, and large twist numbers were obtained. The electron number density of the eruptive structure is found to be $\sim3\times10^{12}$ cm$^{-3}$, indicating that significant amount of mass could be pumped into the corona by the jets. Following the eruption, the extrapolations revealed a set of seemingly relaxed loops, which were visible in the AIA 94 \AA\ channel indicating temperatures of around 6.3 MK. With these observations, we suggest that magnetic braiding could be part of the mechanisms explaining the formation of solar eruption and the mass and energy supplement to the corona.

Numerous observations have revealed that power-law distributions are ubiquitous in energetic solar processes. Hard X-rays, soft X-rays, extreme ultraviolet radiation, and radio waves all display power-law frequency distributions. Since magnetic reconnection is the driving mechanism for many energetic solar phenomena, it is likely that reconnection events themselves display such power-law distributions. In this work, we perform numerical simulations of the solar corona driven by simple convective motions at the photospheric level. Using temperature changes, current distributions, and Poynting fluxes as proxies for heating, we demonstrate that energetic events occurring in our simulation display power-law frequency distributions, with slopes in good agreement with observations. We suggest that the braiding-associated reconnection in the corona can be understood in terms of a self-organized criticality model driven by convective rotational motions similar to those observed at the photosphere.

This is a review of the observational and theoretical evidence for nanoflare heating of the magnetically-closed corona.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

Even in the absence of resolved flares, the corona is heated to several million degrees. However, despite its importance for the structure, dynamics, and evolution of the solar atmosphere, the origin of this heating remains poorly understood. Several observational and theoretical considerations suggest that the heating is driven by small, impulsive energy bursts which could be Parker-style "nanoflares" (Parker 1988) that arise via reconnection within the tangled and twisted coronal magnetic field. The classical "smoking gun" (Klimchuk 2009; Cargill et al. 2013) for impulsive heating is the direct detection of widespread hot plasma (T > 6 MK) with a low emission measure. In recent years there has been great progress in the development of Transition Edge Sensor (TES) X-ray microcalorimeters that make them more ideal for studying the Sun. When combined with grazing-incidence focusing optics, they provide direct spectroscopic imaging over a broad energy band (0.5 to 10 keV) combined with extremely impressive energy resolution in small pixels, as low as 0.7 eV (FWHM) at 1.5 keV (Lee 2015), and 1.56 eV (FWHM) at 6 keV (Smith 2012), two orders of magnitude better than the current best traditional solid state photon-counting spectrometers. Decisive observations of the hot plasma associated with nanoflare models of coronal heating can be provided by new solar microcalorimeters. These measurements will cover the most important part of the coronal spectrum for searching for the nanoflare-related hot plasma and will characterize how much nanoflares can heat the corona both in active regions and the quiet Sun. Finally, microcalorimeters will enable to study all of this as a function of time and space in each pixel simultaneously a capability never before available.

Coronal loops, constituting the basic building blocks of the active Sun, serve as primary targets to help understand the mechanisms responsible for maintaining multi-million Kelvin temperatures in the solar and stellar coronae. Despite significant advances in observations and theory, our knowledge on the fundamental properties of these structures is limited. Here, we present unprecedented observations of accelerating slow magnetoacoustic waves along a coronal loop that show differential propagation speeds in two distinct temperature channels, revealing the multi-stranded and multi-thermal nature of the loop. Utilizing the observed speeds and employing nonlinear force-free magnetic field extrapolations, we derive the actual temperature variation along the loop in both channels, and thus are able to resolve two individual components of the multi-thermal loop for the first time. The obtained positive temperature gradients indicate uniform heating along the loop, rather than isolated footpoint heating.

We use the cellular automaton model described in L\'opez Fuentes \& Klimchuk (2015, ApJ, 799, 128) to study the evolution of coronal loop plasmas. The model, based on the idea of a critical misalignment angle in tangled magnetic fields, produces nanoflares of varying frequency with respect to the plasma cooling time. We compare the results of the model with active region (AR) observations obtained with the Hinode/XRT and SDO/AIA instruments. The comparison is based on the statistical properties of synthetic and observed loop lightcurves. Our results show that the model reproduces the main observational characteristics of the evolution of the plasma in AR coronal loops. The typical intensity fluctuations have an amplitude of 10 to 15\% both for the model and the observations. The sign of the skewness of the intensity distributions indicates the presence of cooling plasma in the loops. We also study the emission measure (EM) distribution predicted by the model and obtain slopes in log(EM) versus log(T) between 2.7 and 4.3, in agreement with published observational values.

We study a 2D cellular automaton (CA) model for the evolution of coronal loop plasmas. The model is based on the idea that coronal loops are made of elementary magnetic strands that are tangled and stressed by the displacement of their footpoints by photospheric motions. The magnetic stress accumulated between neighbor strands is released in sudden reconnection events or nanoflares that heat the plasma. We combine the CA model with the Enthalpy Based Thermal Evolution of Loops (EBTEL) model to compute the response of the plasma to the heating events. Using the known response of the XRT telescope on board Hinode we also obtain synthetic data. The model obeys easy to understand scaling laws relating the output (nanoflare energy, temperature, density, intensity) to the input parameters (field strength, strand length, critical misalignment angle). The nanoflares have a power-law distribution with a universal slope of -2.5, independent of the input parameters. The repetition frequency of nanoflares, expressed in terms of the plasma cooling time, increases with strand length. We discuss the implications of our results for the problem of heating and evolution of active region coronal plasmas.

Among the many ways of investigating coronal heating, the time lag method of Viall & Klimchuk (2012) is becoming increasingly prevalent as an analysis technique complementary to those traditionally used. The time lag method cross correlates light curves at a given spatial location obtained in spectral bands that sample different temperature plasmas. It has been used most extensively with data from the Atmospheric Imaging Assembly on the Solar Dynamics Observatory. We have previously applied the time lag method to entire active regions and surrounding quiet Sun and create maps of the results (Viall & Klimchuk 2012; Viall & Klimchuk 2015). We find that the majority of time lags are consistent with the cooling of coronal plasma that has been impulsively heated. Additionally, a significant fraction of the map area has a time lag of zero. This does not indicate a lack of variability. Rather, strong variability must be present, and it must occur in phase in the different channels. We have shown previously that these zero time lags are consistent with the transition region response to coronal nanoflares (Viall & Klimchuk 2015; Bradshaw & Viall 2016), but other explanations are possible. A common misconception is that the zero time lag indicates steady emission resulting from steady heating. Using simulated and observed light curves, we demonstrate here that highly correlated light curves at zero time lag are not compatible with equilibrium solutions. Such light curves can only be created by evolution.

The million degree plasma of the solar corona must be supplied by the underlying layers of the atmosphere. The mechanism and location of energy release, and the precise source of coronal plasma, remain unresolved. In earlier work we pursued the idea that warm plasma is supplied to the corona via direct heating of the chromosphere by nanoflares, contrary to the prevailing belief that the corona is heated in-situ and the chromosphere is subsequently energized and ablated by thermal conduction. We found that single (low-frequency) chromospheric nanoflares could not explain the observed intensities, Doppler-shifts, and red/blue asymmetries in Fe XII and XIV emission lines. In the present work we follow up on another suggestion that the corona could be powered by chromospheric nanoflares that repeat on a timescale substantially shorter than the cooling/draining timescale. That is, a single magnetic strand is re-supplied with coronal plasma before the existing plasma has time to cool and drain. We perform a series of hydrodynamic experiments and predict the Fe XII and XIV line intensities, Doppler-shifts, and red/blue asymmetries. We find that our predicted quantities disagree dramatically with observations and fully developed loop structures cannot be created by intermediate- or high-frequency chromospheric nanoflares. We conclude that the mechanism ultimately responsible for producing coronal plasma operates above the chromosphere, but this does not preclude the possibility of a similar mechanism powering the chromosphere; extreme examples of which may be responsible for heating chromospheric plasma to transition region temperatures (e.g. type II spicules).

Aspects of solar flare dynamics, such as chromospheric evaporation and flare light-curves, have long been studied using one-dimensional models of plasma dynamics inside a static flare loop, subjected to some energy input. While extremely successful at explaining the observed characteristics of flares, all such models so far have specified energy input ad hoc, rather than deriving it self-consistently. There is broad consensus that flares are powered by magnetic energy released through reconnection. Recent work has generalized Petschek's basic reconnection scenario, topological change followed by field line retraction and shock heating, to permit its inclusion into a one-dimensional flare loop model. Here we compare the gas dynamics driven by retraction and shocking to those from more conventional static loop models energized by ad hoc source terms. We find significant differences during the first minute, when retraction leads to larger kinetic energies and produces higher densities at the loop top, while ad hoc heating tends to rarify the loop top. The loop-top density concentration is related to the slow magnetosonic shock, characteristic of Petschek's model, but persists beyond the retraction phase occurring in the outflow jet. This offers an explanation for observed loop-top sources of X-ray and EUV emission, with advantages over that provided by ad hoc heating scenarios. The cooling phases of the two models are, however, notably similar to one another, suggesting observations at that stage will yield little information on the nature of energy input.

The detailed shapes of spectral line profiles provide valuable information about the emitting plasma, especially when the plasma contains an unresolved mixture of velocities, temperatures, and densities. As a result of finite spectral resolution, the intensity measured by a spectrometer is the average intensity across a wavelength bin of non-zero size. It is assigned to the wavelength position at the center of the bin. However, the actual intensity at that discrete position will be different if the profile is curved, as it invariably is. Standard fitting routines (spline, Gaussian, etc.) do not account for this difference, and this can result in significant errors when making sensitive measurements. Detection of asymmetries in solar coronal emission lines is one example. Removal of line blends is another. We have developed an iterative procedure that corrects for this effect. It can be used with any fitting function, but we employ a cubic spline in a new analysis routine called Intensity Conserving Spline Interpolation (ICSI). As the name implies, it conserves the observed intensity within each wavelength bin, which ordinary fits do not. Given the rapid convergence, speed of computation, and ease of use, we suggest that ICSI be made a standard component of the processing pipeline for spectroscopic data.

We highlight ten key aspects of coronal heating that must be understood before we can consider the problem to be solved. (1) All coronal heating is impulsive. (2) The details of coronal heating matter. (3) The corona is filled with elemental magnetic stands. (4) The corona is densely populated with current sheets. (5) The strands must reconnect to prevent an infinite buildup of stress. (6) Nanoflares repeat with different frequencies. (7) What is the characteristic magnitude of energy release? (8) What causes the collective behavior responsible for loops? (9) What are the onset conditions for energy release? (10) Chromospheric nanoflares are not a primary source of coronal plasma. Significant progress in solving the coronal heating problem will require a coordination of approaches: observational studies, field-aligned hydrodynamic simulations, large-scale and localized 3D MHD simulations, and possibly also kinetic simulations. There is a unique value to each of these approaches, and the community must strive to coordinate better.

Our knowledge of the diffuse emission that encompasses active regions is very limited. In the present paper we investigate two off-limb active regions, namely AR10939 and AR10961, to probe the underlying heating mechanisms. For this purpose we have used spectral observations from Hinode/EIS and employed the emission measure (EM) technique to obtain the thermal structure of these diffuse regions. Our results show that the characteristic EM distributions of the diffuse emission regions peak at log T = 6.25 and the cool-ward slopes are in the range 1.4 - 3.3. This suggests that both low as well as high frequency nanoflare heating events are at work. Our results provide additional constraints on the properties of these diffuse emission regions and their contribution to the background/foreground when active region cores are observed on-disk.

It has been suggested that the hot plasma of the solar corona comes primarily from impulsive heating events, or nanoflares, that occur in the lower atmosphere, either in the upper part of the ordinary chromosphere or at the tips of type II spicules. We test this idea with a series of hydrodynamic simulations. We find that synthetic Fe XII (195) and Fe XIV (274) line profiles generated from the simulations disagree dramatically with actual observations. The integrated line intensities are much too faint; the blue shifts are much too fast; the blue-red asymmetries are much too large; and the emission is confined to low altitudes. We conclude that chromospheric nanoflares are not a primary source of hot coronal plasma. Such events may play an important role in producing the chromosphere and powering its intense radiation, but they do not, in general, raise the temperature of the plasma to coronal values. Those cases where coronal temperatures are reached must be relatively uncommon. The observed profiles of Fe XII and Fe XIV come primarily from plasma that is heated in the corona itself, either by coronal nanoflares or a quasi-steady coronal heating process. Chromospheric nanoflares might play a role in generating waves that provide this coronal heating.

Observations reveal a correspondence between chromospheric type II spicules and bright upwardly moving fronts in the corona observed in the EUV band. However, theoretical considerations suggest that these flows are unlikely to be the main source of heating in coronal magnetic loops. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes, and the possible production of emission in the EUV band. We simulate the propagation of a dense $10^4$ K chromospheric jet upwards along a coronal loop, by means of a 2-D cylindrical MHD model, including gravity, radiative losses, thermal conduction and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool ($\sim 0.8$ MK) corona, linked through a steep transition region. In our reference model, the jet's initial speed is 70 km/s, its initial density is $10^{11}$ cm$^{-3}$, and the ambient uniform magnetic field is 10 G. We explore also other values of jet speed and density in 1-D, and of magnetic field in 2-D, and the jet propagation in a hotter ($\sim 1.5$ MK) background loop. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about $10^6$ K. As a result, a bright moving front becomes visible in the 171 \AA\ channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. We study the area response with a time-dependent 2D MHD loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 MK. We find that the area can change substantially with the quasi-steady heating rate, e.g. by ~40% at 0.5 MK as the loop temperature varies between 1 and 4 MK, and, therefore, affects the interpretation of DEM(T) curves.

Recent solar spectroscopic observations have shown that coronal spectral lines can exhibit asymmetric profiles, with enhanced emissions at their blue wings. These asymmetries correspond to rapidly upflowing plasmas at speeds exceeding ~ 50 km/s. Here, we perform a study of the density of the rapidly upflowing material and compare it to that of the line core which corresponds to the bulk of the plasma. For this task we use spectroscopic observations of several active regions taken by the Extreme Ultraviolet Imaging Spectrometer of the Hinode mission. The density sensitive ratio of the Fe XIV lines at 264.78 and 274.20 A is used to determine wing and core densities. We compute the ratio of the blue wing density to the core density and find that most values are of order unity. This is consistent with the predictions for coronal nanoflares if most of the observed coronal mass is supplied by chromospheric evaporation driven by the nanoflares. However, much larger blue wing-to-core density ratios are predicted if most of the coronal mass is supplied by heated material ejected with type II spicules. Our measurements do not rule out a spicule origin for the blue wing emission, but they argue against spicules being a primary source of the hot plasma in the corona. We note that only about 40% of the pixels where line blends could be safely ignored have blue wing asymmetries in both Fe XIV lines. Anticipated sub-arcsecond spatial resolution spectroscopic observations in future missions could shed more light on the origin of blue, red, and mixed asymmetries.

It has previously been argued that 1. spicules do not provide enough pre-heated plasma to fill the corona, and 2. even if they did, additional heating would be required to keep the plasma hot as it expands upward. We here address the question of whether spicules play an important role by injecting plasma at cooler temperatures ($< 2$ MK), which then gets heated to coronal values at higher altitudes. We measure red-blue asymmetries in line profiles formed over a wide range of temperatures in the bright moss areas of two active regions. We derive emission measure distributions from the excess wing emission. We find that the asymmetries and emission measures are small and conclude that spicules do not inject an important (dominant) mass flux into the cores of active regions at temperatures $> 0.6$ MK ($\log T > 5.8$). These conclusions apply not only to spicules, but to any process that suddenly heats and accelerates chromospheric plasma (e.g., a chromospheric nanoflare). The traditional picture of coronal heating and chromospheric evaporation appears to remain the most likely explanation of the active region corona.

We will use new data from the High-resolution Coronal Imager (Hi-C) with unprecedented spatial resolution of the solar corona to investigate the structure of coronal loops down to 0.2 arcsec. During a rocket flight Hi-C provided images of the solar corona in a wavelength band around 193 A that is dominated by emission from Fe XII showing plasma at temperatures around 1.5 MK. We analyze part of the Hi-C field-of-view to study the smallest coronal loops observed so far and search for the a possible sub-structuring of larger loops. We find tiny 1.5 MK loop-like structures that we interpret as miniature coronal loops. These have length of the coronal segment above the chromosphere of only about 1 Mm and a thickness of less than 200 km. They could be interpreted as the coronal signature of small flux tubes breaking through the photosphere with a footpoint distance corresponding to the diameter of a cell of granulation. We find loops that are longer than 50 Mm to have a diameter of about 2 arcsec or 1.5 Mm, consistent with previous observations. However, Hi-C really resolves these loops with some 20 pixels across the loop. Even at this greatly improved spatial resolution the large loops seem to have no visible sub-structure. Instead they show a smooth variation in cross-section. The fact that the large coronal loops do not show a sub-structure at the spatial scale of 0.1 arcsec per pixel implies that either the densities and temperatures are smoothly varying across these loops or poses an upper limit on the diameter of strands the loops might be composed of. We estimate that strands that compose the 2 arcsec thick loop would have to be thinner than 15 km. The miniature loops we find for the first time pose a challenge to be properly understood in terms of modeling.

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.

A solar flare is composed of impulsive energy release events by magnetic reconnection, which forms and heats flare loops. Recent studies have revealed a two-phase evolution pattern of UV 1600\AA\ emission at the feet of these loops: a rapid pulse lasting for a few seconds to a few minutes, followed by a gradual decay on timescales of a few tens of minutes. Multiple band EUV observations by AIA further reveal very similar signatures. These two phases represent different but related signatures of an impulsive energy release in the corona. The rapid pulse is an immediate response of the lower atmosphere to an intense thermal conduction flux resulting from the sudden heating of the corona to high temperatures (we rule out energetic particles due to a lack of significant hard X-ray emission). The gradual phase is associated with the cooling of hot plasma that has been evaporated into the corona. The observed footpoint emission is again powered by thermal conduction (and enthalpy), but now during a period when approximate steady state conditions are established in the loop. UV and EUV light curves of individual pixels may therefore be separated into contributions from two distinct physical mechanisms to shed light on the nature of energy transport in a flare. We demonstrate this technique using coordinated, spatially resolved observations of UV and EUV emission from the footpoints of a C3.2 thermal flare.

One of the outstanding problems in all of space science is uncovering how the solar corona is heated to temperatures greater than 1 MK. Though studied for decades, one of the major difficulties in solving this problem has been unraveling the line-of-sight (LOS) effects in the observations. The corona is optically thin, so a single pixel measures counts from an indeterminate number (perhaps tens of thousands) of independently heated flux tubes, all along that pixel's LOS. In this paper we model the emission in individual pixels imaging the active region corona in the Extreme Ultraviolet. If LOS effects are not properly taken into account, erroneous conclusions regarding both coronal heating and coronal dynamics may be reached. We model the corona as a LOS integration of many thousands of completely independently heated flux tubes. We demonstrate that despite the superposition of randomly heated flux tubes, nanoflares leave distinct signatures in light curves observed with multi-wavelength and high time cadence data, such as those data taken with the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. These signatures are readily detected with the time-lag analysis technique of Viall & Klimchuk (2012). Steady coronal heating leaves a different and equally distinct signature that is also revealed by the technique.

The time-dependence of heating in solar active regions can be studied by analyzing the slope of the emission measure distribution cool-ward of the peak. In a previous study we showed that low-frequency heating can account for 0% to 77% of active region core emission measures. We now turn our attention to heating by a finite succession of impulsive events for which the timescale between events on a single magnetic strand is shorter than the cooling timescale. We refer to this scenario as a "nanoflare train" and explore a parameter space of heating and coronal loop properties with a hydrodynamic model. Our conclusions are: (1) nanoflare trains are consistent with 86% to 100% of observed active region cores when uncertainties in the atomic data are properly accounted for; (2) steeper slopes are found for larger values of the ratio of the train duration $\Delta_H$ to the post-train cooling and draining timescale $\Delta_C$, where $\Delta_H$ depends on the number of heating events, the event duration and the time interval between successive events ($\tau_C$); (3) $\tau_C$ may be diagnosed from the width of the hot component of the emission measure provided that the temperature bins are much smaller than 0.1 dex; (4) the slope of the emission measure alone is not sufficient to provide information about any timescale associated with heating - the length and density of the heated structure must be measured for $\Delta_H$ to be uniquely extracted from the ratio $\Delta_H/\Delta_C$.

We examine the suggestion that most of the hot plasma in the Sun's corona comes from type II spicule material that is heated as it is ejected from the chromosphere. This contrasts with the traditional view that the corona is filled via chromospheric evaporation that results from coronal heating. We explore the observational consequences of a hypothetical spicule dominated corona and conclude from the large discrepancy between predicted and actual observations that only a small fraction of the hot plasma can be supplied by spicules (<2% in active regions, <5% in the quiet Sun, and <8% in coronal holes). The red-blue asymmetries of EUV spectral lines and the ratio of lower transition region (LTR; T<0.1 MK) to coronal emission measures are both predicted to be 2 orders of magnitude larger than observed. Furthermore, hot spicule material would cool dramatically by adiabatic expansion as it rises into the corona, so substantial coronal heating would be needed to maintain the high temperatures that are seen at all altitudes. We suggest that the corona contains a mixture of thin strands, some of which are populated by spicule injections, but most of which are not. A majority of the observed hot emission originates in non-spicule strands and is explained by traditional coronal heating models. However, since these models predict far too little emission from the LTR, most of this emission comes from the bulk of the spicule material that is only weakly heated and visible in He II (304 A) as it falls back to the surface.

Observational measurements of active region emission measures contain clues to the time-dependence of the underlying heating mechanism. A strongly non-linear scaling of the emission measure with temperature indicates a large amount of hot plasma relative to warm plasma. A weakly non-linear (or linear) scaling of the emission measure indicates a relatively large amount of warm plasma, suggesting that the hot active region plasma is allowed to cool and so the heating is impulsive with a long repeat time. This case is called {\it low-frequency} nanoflare heating and we investigate its feasibility as an active region heating scenario here. We explore a parameter space of heating and coronal loop properties with a hydrodynamic model. For each model run, we calculate the slope $\alpha$ of the emission measure distribution $EM(T) \propto T^\alpha$. Our conclusions are: (1) low-frequency nanoflare heating is consistent with about 36% of observed active region cores when uncertainties in the atomic data are not accounted for; (2) proper consideration of uncertainties yields a range in which as many as 77% of observed active regions are consistent with low-frequency nanoflare heating and as few as zero; (3) low-frequency nanoflare heating cannot explain observed slopes greater than 3; (4) the upper limit to the volumetric energy release is in the region of 50 erg cm$^{-3}$ to avoid unphysical magnetic field strengths; (5) the heating timescale may be short for loops of total length less than 40 Mm to be consistent with the observed range of slopes; (6) predicted slopes are consistently steeper for longer loops.

A well known behavior of EUV light curves of discrete coronal loops is that the peak intensities of cooler channels or spectral lines are reached at progressively later times than hotter channels. This time lag is understood to be the result of hot coronal loop plasma cooling through these lower respective temperatures. However, loops typically comprise only a minority of the total emission in active regions. Is this cooling pattern a common property of active region coronal plasma, or does it only occur in unique circumstances, locations, and times? The new SDO/AIA data provide a wonderful opportunity to answer this question systematically for an entire active region. We measure the time lag between pairs of SDO/AIA EUV channels using 24 hours of images of AR 11082 observed on 19 June 2010. We find that there is a time-lag signal consistent with cooling plasma, just as is usually found for loops, throughout the active region including the diffuse emission between loops for the entire 24 hour duration. The pattern persists consistently for all channel pairs and choice of window length within the 24 hour time period, giving us confidence that the plasma is cooling from temperatures of greater than 3 MK, and sometimes exceeding 7 MK, down to temperatures lower than ~ 0.8 MK. This suggests that the bulk of the emitting coronal plasma in this active region is not steady; rather, it is dynamic and constantly evolving. These measurements provide crucial constraints on any model which seeks to describe coronal heating.

Studying the Doppler shifts and the temperature dependence of Doppler shifts in moss regions can help us understand the heating processes in the core of the active regions. In this paper we have used an active region observation recorded by the Extreme-ultraviolet Imaging Spectrometer (EIS) onboard Hinode on 12-Dec-2007 to measure the Doppler shifts in the moss regions. We have distinguished the moss regions from the rest of the active region by defining a low density cut-off as derived by Tripathi et al. (2010). We have carried out a very careful analysis of the EIS wavelength calibration based on the method described in Young et al. (2012). For spectral lines having maximum sensitivity between log T = 5.85 and log T = 6.25 K, we find that the velocity distribution peaks at around 0 km/s with an estimated error of 4-5 km/s. The width of the distribution decreases with temperature. The mean of the distribution shows a blue shift which increases with increasing temperature and the distribution also shows asymmetries towards blue-shift. Comparing these results with observables predicted from different coronal heating models, we find that these results are consistent with both steady and impulsive heating scenarios. However, the fact that there are a significant number of pixels showing velocity amplitudes that exceed the uncertainty of 5 km s$^{-1}$ is suggestive of impulsive heating. Clearly, further observational constraints are needed to distinguish between these two heating scenarios.

This paper develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating, especially impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. The new features in EBTEL are important for accurate tracking of, in particular, the density. The 0D results are compared to a 1D hydro code (Hydrad) with generally good agreement. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function, (b) extensive parameter surveys and (c) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Since the launch of NASA's Solar Dynamics Observatory on 2010 February 11, the Extreme ultraviolet Variability Experiment (EVE) has observed numerous flares. One interesting feature observed by EVE is that a subset of flares exhibit an additional enhancement of the 2-3 million K emission several hours after the flare's soft X-ray emission. From the Atmospheric Imaging Assembly (AIA) images, we observe that this secondary emission, dubbed the EUV late phase, occurs in the same active region as the flare but not in the same coronal loops. Here, we examine the C8.8 flare that occurred on 2010 May 5 as a case study of EUV late phase flares. In addition to presenting detailed observations from both AIA and EVE, we develop a physical model of this flare and test it using the Enthalpy Based Thermal Evolution of Loops (EBTEL) model.

The heating of the solar corona is one of the big questions in astrophysics. Rapid pulses called nanoflares are among the best candidate mechanisms. The analysis of the time variability of coronal X-ray emission is potentially a very useful tool to detect impulsive events. We analyze the small-scale variability of a solar active region in a high cadence Hinode/XRT observation. The dataset allows us to detect very small deviations of emission fluctuations from the distribution expected for a constant rate. We discuss the deviations in the light of the pulsed-heating scenario.

Using data from the Extreme-ultraviolet Imaging Spectrometer aboard Hinode, we have studied the coronal plasma in the core of two active regions. Concentrating on the area between opposite polarity moss, we found emission measure distributions having an approximate power-law form EM$\propto T^{2.4}$ from $\log\,T = 5.5$ up to a peak at $\log\,T = 6.55$. We show that the observations compare very favorably with a simple model of nanoflare-heated loop strands. They also appear to be consistent with more sophisticated nanoflare models. However, in the absence of additional constraints, steady heating is also a viable explanation.

It is largely agreed that many coronal loops---those observed at a temperature of about 1 MK--- are bundles of unresolved strands that are heated by storms of impulsive nanoflares. The nature of coronal heating in hotter loops and in the very important but largely ignored diffuse component of active regions is much less clear. Are these regions also heated impulsively, or is the heating quasi steady? The spectacular new data from the Atmospheric Imaging Assembly (AIA) telescopes on the Solar Dynamics Observatory (SDO) offer an excellent opportunity to address this question. We analyze the light curves of coronal loops and the diffuse corona in 6 different AIA channels and compare them with the predicted light curves from theoretical models. Light curves in the different AIA channels reach their peak intensities with predictable orderings as a function the nanoflare storm properties. We show that while some sets of light curves exhibit clear evidence of cooling after nanoflare storms, other cases are less straightforward to interpret. Complications arise because of line-of-sight integration through many different structures, the broadband nature of the AIA channels, and because physical properties can change substantially depending on the magnitude of the energy release. Nevertheless, the light curves exhibit predictable and understandable patterns consistent with impulsive nanoflare heating.

It is generally agreed that small impulsive energy bursts called nanoflares are responsible for at least some of the Sun's hot corona, but whether they are the explanation for most of the multi-million degree plasma has been a matter of ongoing debate. We here present evidence that nanoflares are widespread in an active region observed by the X-Ray Telescope on-board the Hinode mission. The distributions of intensity fluctuations have small but important asymmetries, whether taken from individual pixels, multi-pixel subregions, or the entire active region. Negative fluctuations (corresponding to reduced intensity) are greater in number but weaker in amplitude, so that the median fluctuation is negative compared to a mean of zero. Using MonteCarlo simulations, we show that only part of this asymmetry can be explained by Poisson photon statistics. The remainder is explainable with a tendency for exponentially decreasing intensity, such as would be expected from a cooling plasma produced from a nanoflare. We suggest that nanoflares are a universal heating process within active regions.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.

Using a full spectral scan of an active region from the Extreme-Ultraviolet Imaging Spectrometer (EIS) we have obtained Emission Measure EM$(T)$ distributions in two different moss regions within the same active region. We have compared these with theoretical transition region EMs derived for three limiting cases, namely \textit{static equilibrium}, \textit{strong condensation} and \textit{strong evaporation} from \cite{ebtel}. The EM distributions in both the moss regions are strikingly similar and show a monotonically increasing trend from $\log T[\mathrm{K}]=5.15 -6.3$. Using photospheric abundances we obtain a consistent EM distribution for all ions. Comparing the observed and theoretical EM distributions, we find that the observed EM distribution is best explained by the \textit{strong condensation} case (EM$_{con}$), suggesting that a downward enthalpy flux plays an important and possibly dominant role in powering the transition region moss emission. The downflows could be due to unresolved coronal plasma that is cooling and draining after having been impulsively heated. This supports the idea that the hot loops (with temperatures of 3{-}5 MK) seen in the core of active regions are heated by nanoflares.

We develop and analyze a simple cellular automaton (CA) model that reproduces the main properties of the evolution of soft X-ray coronal loops. We are motivated by the observation that these loops evolve in three distinguishable phases that suggest the development, maintainance, and decay of a self-organized system. The model is based on the idea that loops are made of elemental strands that are heated by the relaxation of magnetic stress in the form of nanoflares. In this vision, usually called "the Parker conjecture" (Parker 1988), the origin of stress is the displacement of the strand footpoints due to photospheric convective motions. Modeling the response and evolution of the plasma we obtain synthetic light curves that have the same characteristic properties (intensity, fluctuations, and timescales) as the observed cases. We study the dependence of these properties on the model parameters and find scaling laws that can be used as observational predictions of the model. We discuss the implications of our results for the interpretation of recent loop observations in different wavelengths.

Any successful model of coronal loops must explain a number of observed properties. For warm (~ 1 MK) loops, these include: 1. excess density, 2. flat temperature profile, 3. super-hydrostatic scale height, 4. unstructured intensity profile, and 5. 1000--5000 s lifetime. We examine whether thermal nonequilibrium can reproduce the observations by performing hydrodynamic simulations based on steady coronal heating that decreases exponentially with height. We consider both monolithic and multi-stranded loops. The simulations successfully reproduce certain aspects of the observations, including the excess density, but each of them fails in at least one critical way. Monolithic models have far too much intensity structure, while multi-strand models are either too structured or too long-lived. Our results appear to rule out the widespread existence of heating that is both highly concentrated low in the corona and steady or quasi-steady (slowly varying or impulsive with a rapid cadence). Active regions would have a very different appearance if the dominant heating mechanism had these properties. Thermal nonequilibrium may nonetheless play an important role in prominences and catastrophic cooling events (e.g., coronal rain) that occupy a small fraction of the coronal volume. However, apparent inconsistencies between the models and observations of cooling events have yet to be understood.

It was once thought that all coronal loops are in static equilibrium, but observational and modeling developments over the past decade have shown that this is clearly not the case. It is now established that warm (~1 MK) loops observed in the EUV are explainable as bundles of unresolved strands that are heated impulsively by storms of nanoflares. A raging debate concerning the multi-thermal versus isothermal nature of the loops can be reconciled in terms of the duration of the storm. We show that short and long storms produce narrow and broad thermal distributions, respectively. We also examine the possibility that warm loops can be explained with thermal nonequilibrium, a process by which steady heating produces dynamic behavior whenever the heating is highly concentrated near the loop footpoints. We conclude that this is not a viable explanation for monolithic loops under the conditions we have considered, but that it may have application to multi-stranded loops. Serious questions remain, however.

Nanoflares, short and intense heat pulses within spatially unresolved magnetic strands, are now considered a leading candidate to solve the coronal heating problem. However, the frequent occurrence of nanoflares requires that flare-hot plasma be present in the corona at all times. Its detection has proved elusive until now, in part because the intensities are predicted to be very faint. Here we report on the analysis of an active region observed with five filters by Hinode/XRT in November 2006. We have used the filter ratio method to derive maps of temperature and emission measure both in soft and hard ratios. These maps are approximate in that the plasma is assumed to be isothermal along each line-of-sight. Nonetheless, the hardest available ratio reveals the clear presence of plasma around 10 MK. To obtain more detailed information about the plasma properties, we have performed Monte Carlo simulations assuming a variety of non-isothermal emission measure distributions along the lines-of-sight. We find that the observed filter ratios imply bi-modal distributions consisting of a strong cool (log T ~ 6.3-6.5) component and a weaker (few percent) and hotter (6.6 < log T < 7.2) component. The data are consistent with bi-modal distributions along all lines of sight, i.e., throughout the active region. We also find that the isothermal temperature inferred from a filter ratio depends sensitively on the precise temperature of the cool component. A slight shift of this component can cause the hot component to be obscured in a hard ratio measurement. Consequently, temperature maps made in hard and soft ratios tend to be anti-correlated. We conclude that this observation supports the presence of widespread nanoflaring activity in the active region.

EUV observations of warm coronal loops suggest that they are bundles of unresolved strands that are heated impulsively to high temperatures by nanoflares. The plasma would then have the observed properties (e.g., excess density compared to static equilibrium) when it cools into the 1-2 MK range. If this interpretation is correct, then very hot emission should be present outside of proper flares. It is predicted to be vey faint, however. A critical element for proving or refuting this hypothesis is the existence of hot, very faint plasmas which should be at amounts predicted by impulsive heating. We report on the first comprehensive spectroscopic study of hot plasmas in active regions. Data from the EIS spectrometer on Hinode were used to construct emission measure distributions in quiescent active regions in the 1-5 MK temperature range. The distributions are flat or slowly increasing up to approximately 3 MK and then fall off rapidly at higher temperatures. We show that active region models based on impulsive heating can reproduce the observed EM distributions relatively well. Our results provide strong new evidence that coronal heating is impulsive in nature.

The X-Ray Telescope (XRT) on the Japanese/USA/UK {\it Hinode (Solar-B)} spacecraft has detected emission from a quiescent active region core that is consistent with nanoflare heating. The fluxes from 10 broadband X-ray filters and filter combinations were used to constructed Differential Emission Measure (DEM) curves. In addition to the expected active region peak at Log T = 6.3-6.5, we find a high-temperature component with significant emission measure at Log T $>$ 7.0. This emission measure is weak compared to the main peak -- the DEM is down by almost three orders of magnitude -- which accounts of the fact that it has not been observed with earlier instruments. It is also consistent with spectra of quiescent active regions: no Fe XIX lines are observed in a CHIANTI synthetic spectrum generated using the XRT DEM distribution. The DEM result is successfully reproduced with a simple two-component nanoflare model.

This paper investigates the temporal evolution of temperature, emission measure, energy loss and velocity in a C-class solar flare from both an observational and theoretical perspective. The properties of the flare were derived by following the systematic cooling of the plasma through the response functions of a number of instruments -- RHESSI (>5 MK), GOES-12 (5-30 MK), TRACE 171 A (1 MK) and SOHO/CDS (~0.03-8 MK). These measurements were studied in combination with simulations from the 0-D EBTEL model. At the flare on-set, upflows of ~90 km s-1 and low level emission were observed in Fe XIX, consistent with pre-flare heating and gentle chromospheric evaporation. During the impulsive phase, upflows of ~80 km s-1 in Fe XIX and simultaneous downflows of 20 km s-1 in He I and O V were observed, indicating explosive chromospheric evaporation. The plasma was subsequently found to reach a peak temperature of ~13 MK in approximately 10 minutes. Using EBTEL, conduction was found to be the dominant loss mechanism during the initial ~300s of the decay phase. It was also found to be responsible for driving gentle chromospheric evaporation during this period. As the temperature fell below ~8 MK, and for the next ~4,000s, radiative losses were determined to dominate over conductive losses. The radiative loss phase was accompanied by significant downflows of <40 km s-1 in O V. This is the first extensive study of the evolution of a canonical solar flare using both spectroscopic and broad-band instruments in conjunction with a hydrodynamic model. While our results are in broad agreement with the standard flare model, the simulations suggest that both conductive and non-thermal beam heating play important roles in heating the flare plasma during the impulsive phase of at least this event.

Commission 10 deals with solar activity in all of its forms, ranging from the smallest nanoflares to the largest coronal mass ejections. This report reviews scientific progress over the roughly two-year period ending in the middle of 2008. This has been an exciting time in solar physics, highlighted by the launches of the Hinode and STEREO missions late in 2006. The report is reasonably comprehensive, though it is far from exhaustive. Limited space prevents the inclusion of many significant results. The report is divided into the following sections: Photosphere and Chromosphere; Transition Region; Corona and Coronal Heating; Coronal Jets; Flares; Coronal Mass Ejection Initiation; Global Coronal Waves and Shocks; Coronal Dimming; The Link Between Low Coronal CME Signatures and Magnetic Clouds; Coronal Mass Ejections in the Heliosphere; and Coronal Mass Ejections and Space Weather. Primary authorship is indicated at the beginning of each section.

The physical modeling of active regions (ARs) and of the global coronal is receiving increasing interest lately. Recent attempts to model ARs using static equilibrium models were quite successful in reproducing AR images of hot soft X-ray (SXR) loops. They however failed to predict the bright EUV warm loops permeating ARs: the synthetic images were dominated by intense footpoint emission. We demonstrate that this failure is due to the very weak dependence of loop temperature on loop length which cannot simultaneously account for both hot and warm loops in the same AR. We then consider time-dependent AR models based on nanoflare heating. We demonstrate that such models can simultaneously reproduce EUV and SXR loops in ARs. Moreover, they predict radial intensity variations consistent with the localized core and extended emissions in SXR and EUV AR observations respectively. We finally show how the AR morphology can be used as a gauge of the properties (duration, energy, spatial dependence, repetition time) of the impulsive heating.

Observational and theoretical evidence suggests that coronal heating is impulsive and occurs on very small cross-field spatial scales. A single coronal loop could contain a hundred or more individual strands that are heated quasi-independently by nanoflares. It is therefore an enormous undertaking to model an entire active region or the global corona. Three-dimensional MHD codes have inadequate spatial resolution, and 1D hydro codes are too slow to simulate the many thousands of elemental strands that must be treated in a reasonable representation. Fortunately, thermal conduction and flows tend to smooth out plasma gradients along the magnetic field, so "0D models" are an acceptable alternative. We have developed a highly efficient model called Enthalpy-Based Thermal Evolution of Loops (EBTEL) that accurately describes the evolution of the average temperature, pressure, and density along a coronal strand. It improves significantly upon earlier models of this type--in accuracy, flexibility, and capability. It treats both slowly varying and highly impulsive coronal heating; it provides the differential emission measure distribution, DEM(T), at the transition region footpoints; and there are options for heat flux saturation and nonthermal electron beam heating. EBTEL gives excellent agreement with far more sophisticated 1D hydro simulations despite using four orders of magnitude less computing time. It promises to be a powerful new tool for solar and stellar studies.

We study the effect of the coronal background in the determination of the diameter of EUV loops, and we analyze the suitability of the procedure followed in a previous paper (L\'opez Fuentes, Klimchuk & D\'emoulin 2006) for characterizing their expansion properties. For the analysis we create different synthetic loops and we place them on real backgrounds from data obtained with the Transition Region and Coronal Explorer (\textit{TRACE}). We apply to these loops the same procedure followed in our previous works, and we compare the results with real loop observations. We demonstrate that the procedure allows us to distinguish constant width loops from loops that expand appreciably with height, as predicted by simple force-free field models. This holds even for loops near the resolution limit. The procedure can easily determine when loops are below resolution limit and therefore not reliably measured. We find that small-scale variations in the measured loop width are likely due to imperfections in the background subtraction. The greatest errors occur in especially narrow loops and in places where the background is especially bright relative to the loop. We stress, however, that these effects do not impact the ability to measure large-scale variations. The result that observed loops do not expand systematically with height is robust.

We study the temporal evolution of coronal loops using data from the Solar X-ray Imager (SXI) on board of GOES-12. This instrument allows us to follow in detail the full lifetime of coronal loops. The observed light curves suggest three somewhat distinct evolutionary phases: rise, main, and decay. The durations and characteristic timescales of these phases are much longer than a cooling time and indicate that the loop-averaged heating rate increases slowly, reaches a maintenance level, and then decreases slowly. This suggests that a single heating mechanism operates for the entire lifetime of the loop. For monolithic loops, the loop-averaged heating rate is the intrinsic energy release rate of the heating mechanism. For loops that are bundles of impulsively heated strands, it is an indication of the frequency of occurrence of individual heating events, or nanoflares. We show that the timescale of the loop-averaged heating rate is proportional to the timescale of the observed intensity variation. The ratios of the radiative to conductive cooling times in the loops are somewhat less than 1, putting them intermediate between the values measured previously for hotter and cooler loops. Our results provide further support for the existence of a trend suggesting that all loops are heated by the same mechanism, or that different mechanisms have fundamental similarities (e.g., are all impulsive or are all steady with similar rates of heating).

This article assesses the current state of understanding of coronal heating, outlines the key elements of a comprehensive strategy for solving the problem, and warns of obstacles that must be overcome along the way.

Previous studies have found that coronal loops have a nearly uniform thickness, which seems to disagree with the characteristic expansion of active region magnetic fields. This is one of the most intriguing enigmas in solar physics. We here report on the first comprehensive one-to-one comparison of observed loops with corresponding magnetic flux tubes obtained from cotemporal magnetic field extrapolation models. We use EUV images from TRACE, magnetograms from the MDI instrument on SOHO, and linear force-free field extrapolations. For each loop, we find the particular value of the force-free parameter (alpha) that best matches the observed loop axis and then construct flux tubes using different assumed cross sections at one footpoint (circle and ellipses with different orientations). We find that the flux tubes expand with height by typically twice as much as the corresponding loops. We also find that many flux tubes are much wider at one footpoint than the other, whereas the corresponding loops are far more symmetric. It is clear that the actual coronal magnetic field is more complex than the models we have considered. We suggest that the observed symmetry of loops is related to the tangling of elemental magnetic flux strands produced by photospheric convection.

We use a combination of analytical theory, numerical simulation, and data analysis to study the propagation of acoustic waves along coronal loops. We show that the intensity perturbation of a wave depends on a number of factors, including dissipation of the wave energy, pressure and temperature gradients in the loop atmosphere, work action between the wave and a flow, and the sensitivity properties of the observing instrument. In particular, the scale length of the intensity perturbation varies directly with the dissipation scale length (i.e., damping length) and the scale lengths of pressure, temperature, and velocity. We simulate wave propagation in three different equilibrium loop models and find that dissipation and pressure and temperature stratification are the most important effects in the low corona where the waves are most easily detected. Velocity effects are small, and cross-sectional area variations play no direct role for lines-of-sight that are normal to the loop axis. The intensity perturbation scale lengths in our simulations agree very well with the scale lengths we measure in a sample of loops observed by TRACE. The median observed value is 4.35x10^9 cm. In some cases the intensity perturbation increases with height, which is likely an indication of a temperature inversion in the loop (i.e., temperature that decreases with height). Our most important conclusion is that thermal conduction, the primary damping mechanism, is accurately described by classical transport theory. There is no need to invoke anomalous processes to explain the observations.

We propose a new model for the initiation of a solar coronal mass ejection (CME). The model agrees with two properties of CMEs and eruptive flares that have proved to be very difficult to explain with previous models. a) Very low-lying magnetic field lines, down to the photospheric neutral line, can open toward infinity during an eruption. b) The eruption is driven solely by magnetic free energy stored in a closed, sheared arcade; consequently, the magnetic energy of the closed state is well above that of the post-eruption open state. The key new feature of our model is that CMEs occur in multi-polar topologies, in which reconnection between a sheared arcade and neighboring flux systems triggers the eruption. In this ``magnetic breakout'' model, reconnection removes the unsheared field above the low-lying, sheared core flux near the neutral line, thereby allowing this core flux to burst open. We present numerical simulations which demonstrate that our model can account for the energy requirements for CMEs. We discuss the implication of the model for CME/flare prediction.