Search

Numerous recent X-ray observations of coronal loops in both active regions (ARs) and solar flares have shown clearly that elemental abundances vary with time. Over the course of a flare, they have been found to move from coronal values towards photospheric values near the flare peak, before slowly returning to coronal values during the gradual phase. Coronal loop models typically assume that the elemental abundances are fixed, however. In this work, we introduce a time-variable abundance factor into the 0D ebtel++ code that models the changes due to chromospheric evaporation in order to understand how this affects coronal loop cooling. We find that for strong heating events ($\gtrsim$ 1 erg s$^{-1}$ cm$^{-3}$), the abundances quickly tend towards photospheric values. For smaller heating rates, the abundances fall somewhere between coronal and photospheric values, causing the loop to cool more quickly than the time-fixed photospheric cases (typical flare simulations) and more slowly than time-fixed coronal cases (typical AR simulations). This suggests heating rates in quiescent AR loops no larger than $\approx 0.1$ erg s$^{-1}$ cm$^{-3}$ to be consistent with recent measurements of abundance factors $f \gtrsim 2$., Comment: Accepted to ApJL. Comments and criticisms are more than welcome! The paper is fully reproducible, with scripts to run the simulations and generate the figures available on this Github repository: https://github.com/jwreep/ebtel_abundances

An expansion of cross-sectional area directly impacts the mass flow along a coronal loop, and significantly alters the radiative and hydrodynamic evolution of that loop as a result. Previous studies have found that an area expansion from chromosphere to corona significantly lengthens the cooling time of the corona, and appears to suppress draining from the corona. In this work, we examine the fluid dynamics to understand how the mass flow rate, the energy balance, and the cooling and draining timescales are affected by a non-uniform area. We find that in loops with moderate or large expansion (cross-sectional area expansion factors of 2, 3, 10, 30, 100 from photosphere to apex), impulsive heating, for either direct thermal heating or electron beam heating, induces a steady flow into the corona, so that the coronal density continues to rise during the cooling phase, whereas a uniform loop drains during the cooling phase. The induced upflow carries energy into the corona, balancing the losses from thermal conduction, and continues until thermal conduction weakens enough so that it can no longer support the radiative losses of the transition region (TR). As a result, the plasma cools primarily radiatively until the onset of catastrophic collapse. The speed and duration of the induced upflow both increase in proportion to the rate of area expansion. We argue that observations of blue-shifted spectral lines, therefore, could place a constraint on a loop's area expansion., Comment: Resubmitted to ApJ. Comments and criticisms are welcomed!

A significant impediment to solving the coronal heating problem is that we currently only observe active region (AR) loops in their cooling phase. Previous studies showed that the evolution of cooling loop densities and apex temperatures are insensitive to the magnitude, duration, and location of energy deposition. Still, potential clues to how energy is released are encoded in the cooling phase properties. The appearance of coronal rain, one of the most spectacular phenomena of the cooling phase, occurs when plasma has cooled below 1MK, which sets constraints on the heating frequency, for example. Most observations of coronal rain have been made by imaging instruments. Here we report rare Hinode/EUV Imaging Spectrometer (EIS) observations of a loop arcade where coronal rain forms following an X2.1 limb flare. A bifurcation in plasma composition measurements between photospheric at 1.5MK and coronal at 3.5MK suggests that we are observing post-flare driven coronal rain. Increases in non-thermal velocities and densities with decreasing temperature (2.7MK to 0.6MK) suggest that we are observing the formation and subsequent evolution of the condensations. Doppler velocity measurements imply that a 10% correction of apparent flows in imaging data is reasonable. Emission measure analysis at 0.7MK shows narrow temperature distributions, indicating coherent behaviour reminiscent of that observed in coronal loops. The space-time resolution limitations of EIS suggest that we are observing the largest features or rain showers. These observations provide insights into the heating rate, source, turbulence, and collective behaviour of coronal rain from observations of the loop cooling phase., Comment: To be published in The Astrophysical Journal. Figure 1 animation exceeds size limits but will be available in the online journal version

We identify 182 flares on 158 stars within 100 pc of the Sun in both the near-ultraviolet (NUV: 1750-2750 \r{A}) and far-ultraviolet (FUV: 1350-1750 \r{A}) using high-cadence light curves from the Galaxy Evolution Explorer (GALEX). Ultraviolet (UV) emission from stellar flares plays a crucial role in determining the habitability of exoplanetary systems. However, whether such UV emission promotes or threatens such life depends strongly on the energetics of these flares. Most studies assessing the effect of flares on planetary habitability assume a 9000 K blackbody spectral energy distribution that produces more NUV flux than FUV flux ($R \equiv F_{\rm FUV} / F_{\rm NUV} \approx \frac{1}{6}$). Instead, we observe the opposite with the excess FUV reaching $R \approx \frac{1}{2} - 2$, roughly $3-12$ times the expectation of a 9000 K blackbody. The ratio of FUV to NUV time-integrated flare energies is 3.0 times higher on average than would be predicted by a constant 9000 K blackbody during the flare. Finally, we find that the FUV/NUV ratio at peak tentatively correlates ($\sim 2 \sigma$ significance) both with total UV flare energy and with the G - RP color of the host star. On average, we observe higher FUV/NUV ratios at peak in $E_{\text{UV}}>10^{32}$ erg flares and in flares on fully convective stars., Comment: Submitted to MNRAS, comments welcome

Recent irradiance measurements from numerous heliophysics and astrophysics missions including SDO, GOES, Kepler, TESS, Chandra, XMM-Newton, and NICER have provided critical input in understanding the physics of the most powerful transient events on the Sun and magnetically active stars, solar and stellar flares. The light curves of flare events from the Sun and stars show remarkably similar shapes, typically with a sharp rise and protracted decay phase. The duration of solar and stellar flares has been found to be correlated with the intensity of the event in some wavelengths, such as white light, but not in other wavelengths, such as soft X-rays, but it is not evident why this is the case. In this study, we use a radiative hydrodynamics code to examine factors affecting the duration of flare emission at various wavelengths. The duration of a light curve depends on the temperature of the plasma, the height in the atmosphere at which the emission forms, and the relative importance of cooling due to radiation, thermal conduction, and enthalpy flux. We find that there is a clear distinction between emission that forms low in the atmosphere and responds directly to heating, and emission that forms in the corona, indirectly responding to heating-induced chromospheric evaporation, a facet of the Neupert effect. We discuss the implications of our results to a wide range of flare energies., Comment: Accepted to ApJ

We assess the on-orbit performance of the flare event trigger for the Hinode EUV Imaging Spectrometer. Our goal is to understand the time-delay between the occurrence of a flare, as defined by a prompt rise in soft X-ray emission, and the initiation of the response observing study. Wide (266$''$) slit patrol images in the He II 256.32A spectral line are used for flare hunting, and a reponse is triggered when a pre-defined intensity threshold is reached. We use a sample of 13 $>$ M-class flares that succesfully triggered a response, and compare the timings with soft X-ray data from GOES, and hard X-ray data from RHESSI and Fermi. Excluding complex events that are difficult to interpret, the mean on orbit response time for our sample is 2 min 10 s, with an uncertainty of 84 s. These results may be useful for planning autonomous operations for future missions, and give some guidance as to how improvements could be made to capture the important impulsive phase of flares., Comment: To be published as a Brief Report in Frontiers in Astronomy and Space Sciences

The Neutron star Interior Composition ExploreR (NICER) X-ray observatory observed two powerful X-ray flares equivalent to superflares from the nearby young solar-like star, kappa1 Ceti, in 2019. NICER follows each flare from the onset through the early decay, collecting over 30 cts s-1 near the peak, enabling a detailed spectral variation study of the flare rise. The flare in September varies quickly in ~800 sec, while the flare in December has a few times longer timescale. In both flares, the hard band (2-4 keV) light curves show typical stellar X-ray flare variations with a rapid rise and slow decay, while the soft X-ray light curves, especially of the September flare, have prolonged flat peaks. The time-resolved spectra require two temperature plasma components at kT ~0.3-1 keV and ~2-4 keV. Both components vary similarly, but the cool component lags by ~200 sec with a 4-6 times smaller emission measure (EM) compared to the hot component. A comparison with hydrodynamic flare loop simulations indicates that the cool component originates from X-ray plasma near the magnetic loop footpoints, which mainly cools via thermal conduction. The time lag represents the travel time of the evaporated gas through the entire flare loop. The cool component has several times smaller EM than its simulated counterpart, suggesting a suppression of conductive cooling possibly by the expansion of the loop cross-sectional area or turbulent fluctuations. The cool component's time lag and small EM ratio provide important constraints on the flare loop geometry., Comment: 21 pages, 12 figures, to be published in the Astrophysical Journal

In coronal loop modeling, it is commonly assumed that the loops are semi-circular with a uniform cross-sectional area. However, observed loops are rarely semi-circular, and extrapolations of the magnetic field show that the field strength decreases with height, implying that the cross-sectional area expands with height. We examine these two assumptions directly to understand how they affect the hydrodynamic and radiative response of short, hot loops to strong, impulsive electron beam heating events. Both the magnitude and rate of area expansion impact the dynamics directly, and an expanding cross-section significantly lengthens the time for a loop to cool and drain, increases upflow durations, and suppresses sound waves. The standard $T \sim n^{2}$ relation for radiative cooling does not hold with expanding loops, which cool with relatively little draining. An increase in the eccentricity of loops, on the other hand, only increases the draining timescale, and is a minor effect in general. Spectral line intensities are also strongly impacted by the variation in the cross-sectional area since they depend on both the volume of the emitting region as well as the density and ionization state. With a larger expansion, the density is reduced, so the lines at all heights are relatively reduced in intensity and, because of the increase of cooling times, the hottest lines remain bright for significantly longer. Area expansion is critical to accurate modeling of the hydrodynamics and radiation, and observations are needed to constrain the magnitude, rate, and location of the expansion or lack thereof., Comment: Accepted to ApJ

We examine SDO/EVE data to better understand solar flare irradiance, and how that irradiance may vary for large events. We measure scaling laws relating GOES flare classes to irradiance in 21 lines measured with SDO/EVE, formed across a wide range of temperatures, and find that this scaling depends on the line formation temperature. We extrapolate these irradiance values to large events, exceeding X10. In order to create full spectra, however, we need a physical model of the irradiance. We present the first results of a new physical model of solar flare irradiance, NRLFLARE, that sums together a series of flare loops to calculate the spectral irradiance ranging from the X-rays through the far ultraviolet (~ 0 to 1250 Angstroms), constrained by GOES/XRS observations. We test this model against SDO/EVE data. The model spectra and time evolution compares well in high temperature emission, but cooler lines show large discrepancies. We speculate that the discrepancies are likely due to both a non-uniform cross section of the flaring loops as well as opacity effects. We then show that allowing the cross-sectional area to vary with height significantly improves agreement with observations, and is therefore a crucial parameter needed to accurately model the intensity of spectral lines, particularly in the transition region from 4.7 < log T < 6.0., Comment: Accepted to ApJ

The solar X-ray irradiance is significantly heightened during the course of a solar flare, which can cause radio blackouts due to ionization of the atoms in the ionosphere. As the duration of a solar flare is not related to the size of that flare, it is not directly clear how long those blackouts can persist. Using a random forest regression model trained on data taken from X-ray light curves, we have developed a direct forecasting method that predicts how long the event will remain above background levels. We test this on a large collection of flares observed with GOES-15, and show that it generally outperforms simple linear regression. This forecast is computationally light enough to be performed in real time, allowing for the prediction to be made during the course of a flare., Comment: Resubmitted to Space Weather. Comments welcome!

The nanoflare paradigm of coronal heating has proven extremely promising for explaining the presence of hot, multi-million degree loops in the solar corona. In this paradigm, localized heating events supply enough energy to heat the solar atmosphere to its observed temperatures. Rigorously modeling this process, however, has proven difficult, since it requires an accurate treatment of both the magnetic field dynamics and reconnection as well as the plasma's response to magnetic perturbations. In this paper, we combine fully 3D magnetohydrodynamic (MHD) simulations of coronal active region plasma driven by photospheric motions with spatially-averaged, time-dependent hydrodynamic (HD) modeling of coronal loops to obtain physically motivated observables that can be quantitatively compared with observational measurements of active region cores. We take the behavior of reconnected field lines from the MHD simulation and use them to populate the HD model to obtain the thermodynamic evolution of the plasma and subsequently the emission measure distribution. We find the that the photospheric driving of the MHD model produces only very low-frequency nanoflare heating which cannot account for the full range of active region core observations as measured by the low-temperature emission measure slope. Additionally, we calculate the spatial and temporal distributions of field lines exhibiting collective behavior, and argue that loops occur due to random energization occurring on clusters of adjacent field lines.

The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute.

We have investigated an M1.3 limb flare, which develops as a magnetic loop/arch that fans out from an X-ray jet. Using Hinode/EIS, we found that the temperature increases with height to a value of over 10$^{7}$ K at the loop-top during the flare. The measured Doppler velocity (redshifts of 100$-$500 km s$^{-1}$) and the non-thermal velocity ($\geq$100 km s$^{-1}$) from Fe XXIV also increase with loop height. The electron density increases from $0.3\times10^{9}$ cm$^{-3}$ early in the flare rise to $1.3\times10^{9}$ cm$^{-3}$ after the flare peak. The 3-D structure of the loop derived with STEREO/EUVI indicates that the strong redshift in the loop-top region is due to upflowing plasma originating from the jet. Both hard X-ray and soft X-ray emission from RHESSI were only seen as footpoint brightenings during the impulsive phase of the flare, then, soft X-ray emission moves to the loop-top in the decay phase. Based on the temperature and density measurements and theoretical cooling models, the temperature evolution of the flare arch is consistent with impulsive heating during the jet eruption followed by conductive cooling via evaporation and minor prolonged heating in the top of the fan loop. Investigating the magnetic field topology and squashing factor map from SDO/HMI, we conclude that the observed magnetic-fan flaring arch is mostly heated from low atmospheric reconnection accompanying the jet ejection, instead of from reconnection above the arch as expected in the standard flare model., Comment: 55 pages, 21 Figures, accepted for publication in ApJ

Understanding how energy is released in flares is one of the central problems of solar and stellar astrophysics. Observations of high temperature flare plasma hold many potential clues as to the nature of this energy release. It is clear, however, that flares are not composed of a few impulsively heated loops, but are the result of heating on many small-scale threads that are energized over time, making it difficult to compare observations and numerical simulations in detail. Several previous studies have shown that it is possible to reproduce some aspects of the observed emission by considering the flare as a sequence of independently heated loops, but these studies generally focus on small-scale features while ignoring the global features of the flare. In this paper, we develop a multithreaded model that encompasses the time-varying geometry and heating rate for a series of successively-heated loops comprising an arcade. To validate, we compare with spectral observations of five flares made with the MinXSS CubeSat as well as light curves measured with GOES/XRS and SDO/AIA. We show that this model can successfully reproduce the light curves and quasi-periodic pulsations in GOES/XRS, the soft X-ray spectra seen with MinXSS, and the light curves in various AIA passbands. The AIA light curves are most consistent with long duration heating, but elemental abundances cannot be constrained with the model. Finally, we show how this model can be used to extrapolate to spectra of extreme events that can predict irradiance across a wide wavelength range including unobserved wavelengths., Comment: Accepted to ApJ. v2: Mistake in temperature plots in Figure 1 corrected

Coronal rain is ubiquitous in flare loops, forming shortly after the onset of the solar flare. Rain is thought to be caused by a thermal instability, a localized runaway cooling of material in the corona. The models that demonstrate this require extremely long duration heating on the order of the radiative cooling time, localized near the footpoints of the loops. In flares, electron beams are thought to be the primary energy transport mechanism, driving strong footpoint heating during the impulsive phase that causes evaporation, filling and heating flare loops. Electron beams, however, do not act for a long period of time, and even supposing that they did, their heating would not remain localized at the footpoints. With a series of numerical experiments, we show directly that these two issues mean that electron beams are incapable of causing the formation of rain in flare loops. This result suggests that either there is another mechanism acting in flare loops responsible for rain, or that the modeling of the cooling of flare loops is somehow deficient. To adequately describe flares, the standard model must address this issue to account for the presence of coronal rain., Comment: Published in ApJ

Solar flares are observed and classified according to their intensity measured with the GOES X-ray Sensors. We show that the duration of a flare, as measured by the full width at half maximum (FWHM) in GOES is not related to the size of the flare as measured by GOES intensity. The durations of X-class flares range from a few minutes to a few hours, and the same is true of M- and C-class flares. In this work, we therefore examine the statistical relationships between the basic properties of flares -- temperature, emission measure, energy, etc. -- in comparison to both their size and duration. We find that the size of the flare is directly related to all of these basic properties, as previously found by many authors. The duration is not so clear. When examining the whole data set, the duration appears to be independent of all of these properties. In larger flares, however, there are direct correlations between the GOES FWHM and magnetic reconnection flux and ribbon area. We discuss the possible explanations, finding that this discrepancy may be due to large uncertainties in small flares, though we cannot rule out the possibility that the driving physical processes are different in smaller flares than larger ones. We discuss the implications of this result and how it relates to the magnetic reconnection process that releases energy in flares., Comment: ApJ, in press

Understanding the dynamics of the solar chromosphere is crucial to understanding energy transport across the solar atmosphere. The chromosphere is optically thick at many wavelengths and described by non-local thermodynamic equilibrium (NLTE), making it difficult to interpret observations. Furthermore, there is considerable evidence that the atmosphere is filamented, and that current instruments do not sufficiently resolve small scale features. In flares, it is likely that multithreaded models are required to describe the heating. The combination of NLTE effects and multithreaded modeling requires computationally demanding calculations, which has motivated the development of a model that can efficiently treat both. We describe the implementation of a solver in a hydrodynamic code for the hydrogen level populations that approximates the NLTE solutions. We derive an accurate electron density across the atmosphere, that includes the effects of non-equilibrium ionization for helium and metals. We show the effects on hydrodynamic simulations, which are used to synthesize light curves using a post-processing radiative transfer code. We demonstrate the utility of this model on IRIS observations of a small flare. We show that the Doppler shifts in Mg II, C II, and O I can be explained with a multithreaded model of loops subjected to electron beam heating, so long as NLTE effects are treated. The intensities, however, do not match observed values very well, which is due to assumptions about the initial atmosphere. We briefly show how altering the initial atmosphere can drastically alter line profiles, and therefore derived quantities, and suggest that it should be tuned to pre-flare observations., Comment: Now accepted to ApJ. Comments and criticisms encouraged!

Solar flares form and release energy across a large number of magnetic loops. The global parameters of flares, such as the total energy released, duration, physical size, etc., are routinely measured, and the hydrodynamics of a coronal loop subjected to intense heating have been extensively studied. It is not clear, however, how many loops comprise a flare, nor how the total energy is partitioned between them. In this work, we employ a hydrodynamic model to better understand the energy partition by synthesizing Si IV and Fe XXI line emission and comparing to observations of these lines with IRIS. We find that the observed temporal evolution of the Doppler shifts holds important information on the heating duration. To demonstrate this we first examine a single loop model, and find that the properties of chromospheric evaporation seen in Fe XXI can be reproduced by loops heated for long durations, while persistent red-shifts seen in Si IV cannot be reproduced by any single loop model. We then examine a multi-threaded model, assuming both a fixed heating duration on all loops, and a distribution of heating durations. For a fixed heating duration, we find that durations of 100 -- 200 s do a fair job of reproducing both the red- and blue-shifts, while a distribution of durations, with a median of about 50 -- 100 s, does a better job. Finally, we compare our simulations directly to observations of an M-class flare seen by IRIS, and find good agreement between the modeled and observed values given these constraints., Comment: Accepted to ApJ

Alfv\'enic waves have been proposed as an important energy transport mechanism in coronal loops, capable of delivering energy to both the corona and chromosphere and giving rise to many observed features, of flaring and quiescent regions. In previous work, we established that resistive dissipation of waves (ambipolar diffusion) can drive strong chromospheric heating and evaporation, capable of producing flaring signatures. However, that model was based on a simplified assumption that the waves propagate instantly to the chromosphere, an assumption which the current work removes. Via a ray tracing method, we have implemented traveling waves in a field-aligned hydrodynamic simulation that dissipate locally as they propagate along the field line. We compare this method to and validate against the magnetohydrodynamics code Lare3D. We then examine the importance of travel times to the dynamics of the loop evolution, finding that (1) the ionization level of the plasma plays a critical role in determining the location and rate at which waves dissipate; (2) long duration waves effectively bore a hole into the chromosphere, allowing subsequent waves to penetrate deeper than previously expected, unlike an electron beam whose energy deposition rises in height as evaporation reduces the mean-free paths of the electrons; (3) the dissipation of these waves drives a pressure front that propagates to deeper depths, unlike energy deposition by an electron beam., Comment: Accepted to ApJ

We report on the structure and evolution of a current sheet that formed in the wake of an eruptive X8.3 flare observed at the west limb of the Sun on September 10, 2017. Using observations from the EUV Imaging Spectrometer (EIS) on Hinode and the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO), we find that plasma in the current sheet reaches temperatures of about 20 MK and that the range of temperatures is relatively narrow. The highest temperatures occur at the base of the current sheet, in the region near the top of the post-flare loop arcade. The broadest high temperature line profiles, in contrast, occur at the largest observed heights. Further, line broadening is strong very early in the flare and diminishes over time. The current sheet can be observed in the AIA 211 and 171 channels, which have a considerable contribution from thermal bremsstrahlung at flare temperatures. Comparisons of the emission measure in these channels with other EIS wavelengths and AIA channels dominated by Fe line emission indicate a coronal composition and suggest that the current sheet is formed by the heating of plasma already in the corona. Taken together, these observations suggest that some flare heating occurs in the current sheet while additional energy is released as newly reconnected field lines relax and become more dipolar., Comment: Submitted to ApJ, but comments are welcome!

GOES soft X-ray light curves are used to measure the timing and duration of solar flare emission. The timing and duration of the magnetic reconnection and subsequent energy release which drives solar flares are unknown, though the light curves are presumably related. It is therefore critical to understand the physics which connects the two: how does the time scale of reconnection produce an observed GOES light curve? In this work, we model the formation and expansion of an arcade of loops with a hydrodynamic model, which we then use to synthesize GOES light curves. We calculate the FWHM and the e-folding decay time of the light curves and compare them to the separation of the centroids of the two ribbons which the arcade spans, which is representative of the size scale of the loops. We reproduce a linear relation between the two, as found observationally in previous work. We show that this demonstrates a direct connection between the duration of energy release and the evolution of these light curves. We also show that the cooling processes of individual loops comprising the flare arcade directly affect the measured time scales. From the clear consistency between the observed and modeled linearity, we conclude that the primary factors that control the flare time scales are the duration of reconnection and the loop lengths., Comment: Accepted to ApJ

Solar flares are driven by the release of magnetic energy from reconnection events in the solar corona, whereafter energy is transported to the chromosphere, heating the plasma and causing the characteristic radiative losses. In the collisional thick-target model, electrons accelerated to energies exceeding 10 keV traverse the corona and impact the chromosphere, where they deposit their energy through collisions with the much denser plasma in the lower atmosphere. While there are undoubtedly high energy non-thermal electrons accelerated in flares, it is unclear whether these electron beams are the sole mechanism of energy transport, or whether they only dominate in certain phases of the flare's evolution. Alfv\'enic waves are generated during the post-reconnection relaxation of magnetic field lines, so it is important to examine their role in energy transport., Comment: Next Generation Solar Physics Mission White paper, 4 pages, 2 figures

Results from the Solar Maximum Mission showed a close connection between the hard X-ray and transition region emission in solar flares. Analogously, the modern combination of RHESSI and IRIS data can inform the details of heating processes in ways never before possible. We study a small event that was observed with RHESSI, IRIS, SDO, and Hinode, allowing us to strongly constrain the heating and hydrodynamical properties of the flare, with detailed observations presented in a previous paper. Long duration red-shifts of transition region lines observed in this event, as well as many other events, are fundamentally incompatible with chromospheric condensation on a single loop. We combine RHESSI and IRIS data to measure the energy partition among the many magnetic strands that comprise the flare. Using that observationally determined energy partition, we show that a proper multi-threaded model can reproduce these red-shifts in magnitude, duration, and line intensity, while simultaneously being well constrained by the observed density, temperature, and emission measure. We comment on the implications for both RHESSI and IRIS observations of flares in general, namely that: (1) a single loop model is inconsistent with long duration red-shifts, among other observables; (2) the average time between energization of strands is less than 10 seconds, which implies that for a hard X-ray burst lasting ten minutes, there were at least 60 strands within a single IRIS pixel located on the flare ribbon; (3) the majority of these strands were explosively heated with energy distribution well described by a power law of slope $\approx -1.6$; (4) the multi-stranded model reproduces the observed line profiles, peak temperatures, differential emission measure distributions, and densities., Comment: 10 pages, 8 figures. Accepted to ApJ

We exploit the high spatial resolution and high cadence of the Interface Region Imaging Spectrograph (IRIS) to investigate the response of the transition region and chromosphere to energy deposition during a small flare. Simultaneous observations from RHESSI provide constraints on the energetic electrons precipitating into the flare footpoints while observations of XRT, AIA, and EIS allow us to measure the temperatures and emission measures from the resulting flare loops. We find clear evidence for heating over an extended period on the spatial scale of a single IRIS pixel. During the impulsive phase of this event the intensities in each pixel for the Si IV 1402.770, C II 1334.535, Mg II 2796.354 and O I 1355.598 emission lines are characterized by numerous, small-scale bursts typically lasting 60s or less. Red shifts are observed in Si IV, C II, and Mg II during the impulsive phase. Mg II shows red-shifts during the bursts and stationary emission at other times. The Si IV and C II profiles, in contrast, are observed to be red-shifted at all times during the impulsive phase. These persistent red-shifts are a challenge for one-dimensional hydrodynamic models, which predict only short-duration downflows in response to impulsive heating. We conjecture that energy is being released on many small-scale filaments with a power-law distribution of heating rates., Comment: Comments are welcome!

Recent observations from the Interface Region Imaging Spectrograph (IRIS) have discovered a new class of numerous low-lying dynamic loop structures, and it has been argued that they are the long-postulated unresolved fine structures (UFS) that dominate the emission of the solar transition region. In this letter, we combine IRIS measurements of the properties of a sample of 108 UFS (intensities, lengths, widths, lifetimes) with 1-D non-equilibrium ionization simulations using the HYDRAD hydrodynamic model to examine whether the UFS are now truly spatially resolved in the sense of being individual structures rather than composed of multiple magnetic threads. We find that a simulation of an impulsively heated single strand can reproduce most of the observed properties suggesting that the UFS may be resolved, and the distribution of UFS widths implies that they are structured on a spatial scale of 133km on average. Spatial scales of a few hundred km appear to be typical for a range of chromospheric and coronal structures, and we conjecture that this could be an important clue to the coronal heating process., Comment: This is the version to be published

Brown & Mallik 2008 and the Brown et al. 2010 corrigendum of it presented expressions for non-thermal recombination (NTR) in the collisionally thin- and thick-target regimes, claiming that the process could account for a substantial part of hard X-ray continuum in solar flares usually attributed entirely to thermal and non-thermal bremsstrahlung (NTB). However, we have found the thick-target expression to become unphysical for low cut-offs in the injected electron energy spectrum. We trace this to an error in the derivation, derive a corrected version which is real-valued and continuous for all photon energies and cut-offs, and show that, for thick targets, Brown et al. over-estimated NTR emission at small photon energies. The regime of small cut-offs and large spectral indices involve large (reducing) correction factors but in some other thick-target parameter regimes NTR/NTB can still be of order unity. We comment on the importance of these results to flare and to microflare modeling and spectral fitting. An empirical fit to our results shows that the peak NTR contribution comprises over half the hard X-ray signal if delta > 6 (E_0c/4 keV)^(0.4)., Comment: Accepted to ApJ. 5 pages, 3 figures

We have developed a numerical model of flare heating due to the dissipation of Alfv\'enic waves propagating from the corona to the chromosphere. With this model, we present an investigation of the key parameters of these waves on the energy transport, heating, and subsequent dynamics. For sufficiently high frequencies and perpendicular wave numbers, the waves dissipate significantly in the upper chromosphere, strongly heating it to flare temperatures. This heating can then drive strong chromospheric evaporation, bringing hot and dense plasma to the corona. We therefore find three important conclusions: (1) Alfv\'enic waves, propagating from the corona to the chromosphere, are capable of heating the upper chromosphere and the corona, (2) the atmospheric response to heating due to the dissipation of Alfv\'enic waves can be strikingly similar to heating by an electron beam, and (3) this heating can produce explosive evaporation., Comment: Accepted to ApJL

Observations of solar flares with RHESSI have shown X-ray sources traveling along flaring loops, from the corona down to the chromosphere and back up. The 28 November 2002 C1.1 flare, first observed with RHESSI by Sui et al. 2006 and quantitatively analyzed by O'Flannagain et al. 2013, very clearly shows this behavior. By employing numerical experiments, we use these observations of X-ray source height motions as a constraint to distinguish between heating due to a non-thermal electron beam and in situ energy deposition in the corona. We find that both heating scenarios can reproduce the observed light curves, but our results favor non-thermal heating. In situ heating is inconsistent with the observed X-ray source morphology and always gives a height dispersion with photon energy opposite to what is observed., Comment: Accepted to ApJ

The Extreme-ultraviolet Imaging Spectrometer (EIS) on the Hinode spacecraft observed flare footpoint regions coincident with a surge for a M3.7 flare observed on 25 September 2011 at N12 E33 in active region 11302. The flare was observed in spectral lines of O VI, Fe X, Fe XII, Fe XIV, Fe XV, Fe XVI, Fe XVII, Fe XXIII and Fe XXIV. The EIS observations were made coincident with hard X-ray bursts observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Overlays of the RHESSI images on the EIS raster images at different wavelengths show a spatial coincidence of features in the RHESSI images with the EIS upflow and downflow regions, as well as loop-top or near-loop-top regions. A complex array of phenomena was observed including multiple evaporation regions and the surge, which was also observed by the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) telescopes. The slit of the EIS spectrometer covered several flare footpoint regions from which evaporative upflows in Fe XXIII and Fe XXIV lines were observed with Doppler speeds greater than 500 km s$^{-1}$. For ions such as Fe XV both evaporative outflows (~200 km s$^{-1}$) and downflows (~30-50 km s$^{-1}$) were observed. Non-thermal motions from 120 to 300 km s$^{-1}$ were measured in flare lines. In the surge, Doppler speeds are found from about 0 to over 250 km s$^{-1}$ in lines from ions such as Fe XIV. The non-thermal motions could be due to multiple sources slightly Doppler-shifted from each other or turbulence in the evaporating plasma. We estimate the energetics of the hard X-ray burst and obtain a total flare energy in accelerated electrons of $\geq7\times10^{28}$ ergs. This is a lower limit because only an upper limit can be determined for the low energy cutoff to the electron spectrum. We find that detailed modeling of this event would require a multi-threaded model due to its complexity., Comment: 42 pages, 18 figures

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$.

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.