We perform an analysis of the line-of-sight (LOS) observables of the SDO/HMI for models of the solar atmosphere heated by precipitating high-energy electrons during solar flares. The radiative hydrodynamic (RADYN) flare models are obtained from the F-CHROMA database. The Stokes profiles for the Fe 6173 A line observed by SDO/HMI are calculated using the radiative transfer code RH1.5D assuming statistical equilibrium for atomic level populations and imposing uniform background vertical magnetic field of various strength. The SDO/HMI observing sequence and LOS data processing pipeline algorithm are applied to derive the observables (continuum intensity, line depth, Doppler velocity, LOS magnetic field). Our results reveal that the strongest deviations of the observables from the actual spectroscopic line parameters are found for the model with a total energy deposited of $E_{total}=1.0\times{}10^{12}$ erg cm$^{-2}$, injected with power-law spectral index of $\delta=3$ above a low-energy cutoff of $E_{c}=25$ keV. The magnitudes of the velocity and magnetic field deviations depend on the imposed magnetic field, and can reach 0.35 km/s for LOS velocities, 90 G for LOS magnetic field, and 3\% for continuum enhancement for the 1000 G imposed LOS magnetic field setup. For $E_{total}\geq{}3.0\times{}10^{11}$ erg cm$^{-2}$ models, the velocity and magnetic field deviations are most strongly correlated with the energy flux carried by $\sim$50 keV electrons, and the continuum enhancement is correlated with the synthesized $\sim$55-60 keV hard X-ray photon flux. The relatively low magnitudes of perturbations of the observables and absence of magnetic field sign reversals suggest that the considered radiative hydrodynamic beam heating models augmented with the uniform vertical magnetic field setups cannot explain strong transient changes found in the SDO/HMI observations.

To extract the information that the Mg II NUV spectra (observed by the Interface Region Imaging Spectrograph; IRIS), carries about the chromosphere during solar flares, and to validate models of energy transport via model-data comparison, forward modelling is required. The assumption of statistical equilibrium is typically used to obtain the atomic level populations from snapshots of flare atmospheres, due to computational necessity. However it is possible that relying on statistical equilibrium could lead to spurious results. We compare solving the atomic level populations via statistical equilibrium versus a non-equilibrium time-dependent approach. This was achieved using flare simulations from RADYN alongside the minority species version, MS_RADYN, from which the time-dependent Mg II atomic level populations and radiation transfer were computed in complete frequency redistribution. The impacts on the emergent profiles, lightcurves, line ratios, and formation heights are discussed. In summary we note that non-equilibrium effects during flares are typically important only in the initial stages and for a short period following the cessation of the energy injection. An analysis of the timescales of ionisation equilibrium reveals that for most of the duration of the flare, when the temperatures and densities are sufficiently enhanced, the relaxation timescales are short ($\tau_{\mathrm{relax}}<0.1$ s), so that the equilibrium solution is an adequate approximation. These effects vary with the size of the flare, however. In weaker flares effects can be more pronounced. We recommend that non-equilibrium effects be considered when possible, but that statistical equilibrium is sufficient at most stages of the flare.

The Interface Region Imaging Spectrograph (IRIS) has routinely observed the flaring Mg II NUV spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward modelled both the general properties of these lines and specific flare observations. Generally these have forward modelled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward modelling MgII in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II: (1) PRD is still required during flare simulations despite the increased densities, (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot, (3) including Mg I in NLTE results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum, (4) only hydrogen and Mg II need to be included in NLTE, (5) ideally the non-equilibrium hydrogen populations, with non-thermal collisional rates, should be used rather than the statistical equilibrium populations, (6) an atom consisting of only the ground state, h & k upper levels, and continuum level is insufficient to model the resonance lines, and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of transition region and coronal irradiation in flares.

The GOES X1 flare SOL2014-10-25T17:08:00 was a three-ribbon solar flare observed with IRIS in the near and far ultraviolet. One of the flare ribbons crossed a sunspot umbra, producing a dramatic, $\sim1000$% increase in the near-ultraviolet (NUV) continuum radiation. We comprehensively analyze the ultraviolet spectral data of the umbral flare brightenings, which provide new challenges for radiative-hydrodynamic modeling of the chromospheric velocity field and the white-light continuum radiation. The emission line profiles in the umbral flare brightenings exhibit redshifts and profile asymmetries, but these are significantly smaller than in another, well-studied X-class solar flare. We present a ratio of the NUV continuum intensity to the Fe II 2814.45 Ang intensity. This continuum-to-line ratio is a new spectral diagnostic of significant heating at high column mass (log $m/$[g cm$^{-2}] >-2$) during solar flares because the continuum and emission line radiation originate from relatively similar temperatures but moderately different optical depths. The full spectral readout of these IRIS data also allow for a comprehensive survey of the flaring NUV landscape: in addition to many lines of Fe II and Cr II, we identify a new solar flare emission line, He I $\lambda2829.91$ (as previously identified in laboratory and early-type stellar spectra). The Fermi/GBM hard X-ray data provide inputs to radiative-hydrodynamic models (which will be presented in Paper II) in order to better understand the large continuum-to-line ratios, the origin of the white-light continuum radiation, and the role of electron beam heating in the low atmosphere.

We present a statistical study of chromospheric evaporation in solar flares using simultaneous observations by the RHESSI X-ray telescope and the IRIS UV spectrograph. The results are compared with radiation hydrodynamic flare models from the F-CHROMA RADYN database. For each event, we study spatially-resolved Doppler shifts of spectral lines formed in the transition region (C\,II\,1334.5\,{\AA}) and hot coronal plasma (Fe\,XXI\,1354.1\,{\AA}) to investigate the dynamics of the solar atmosphere during the flare impulsive phase. We estimate the energy fluxes deposited by high-energy electrons using X-ray imaging spectroscopy and assuming the standard thick-target model. Using the RADYN flare models, the RH 1.5D radiative transfer code and the Chianti atomic line database, we calculate C\,II and Fe\,XXI line profiles and compare with the observations. While the RADYN models predict a correlation between the Doppler shifts and deposited energy flux for both lines, this was only observed in the C II data. Several quantitative discrepancies are found between the observations and models: the Fe\,XXI Doppler shifts are substantially stronger in the models than in the data; the C\,II mean blueshifts are absent in the observations but predicted by the models. The transition energy between `gentle' and `explosive' evaporation regimes estimated from the observations ($2-8\times{}10^{9}$~erg~cm$^{-2}$~s$^{-1}$) and derived from the models ($2.2-10.1\times{}10^{9}$~erg~cm$^{-2}$~s$^{-1}$) are comparable with each other. The results illustrate relationships among the processes of chromospheric evaporation, response of the colder layers, and the flare energy flux deposited by high-energy electrons, although demonstrating discrepancy between analyzed observations and RADYN models.

The Interface Region Imaging Spectrograph (IRIS) routinely observes the Si IV resonance lines. When analyzing observations of these lines it has typically been assumed they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron beam driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres we solved the radiation transfer equations to obtain the non-LTE, non-equilibrium populations, line profiles, and opacities for a model Silicon atom, including charge exchange processes. This was achieved using the `minority species' version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionisation equilibrium. All simulations with an injected energy flux $F>5\times10^{10}$ erg cm$^{-2}$ s$^{-1}$ resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to $F\approx5\times10^{9}$ erg cm$^{-2}$ s$^{-1}$) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of $5\times10^{-6}$ g cm$^{-2}$ over the temperature range $40$ to $100$ kK, and the Si IV formation temperatures were between $30$ and $60$ kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

The hydrogen Lyman lines ( 91.2 nm < lambda < 121.6 nm) are significant contributors to the radiative losses of the solar chromosphere, and are enhanced during flares. We have shown previously that the Lyman lines observed by the Extreme Ultraviolet Variability instrument onboard the Solar Dynamics Observatory exhibit Doppler motions equivalent to speeds on the order of 30 km/s. But contrary to expectation, no dominant flow direction was observed, with both redshifts and blueshifts present. To understand the formation of the Lyman lines, particularly their Doppler motions, we have used the radiative hydrodynamic code, RADYN, and the radiative transfer code, RH, to simulate the evolution of the flaring chromosphere and the response of the Lyman lines during solar flares. We find that upflows in the simulated atmospheres lead to blueshifts in the line cores, which exhibit central reversals. We then model the effects of the instrument on the profiles using the EVE instrument's properties. What may be interpreted as downflows (redshifted emission) in the lines after they have been convolved with the instrumental line profile may not necessarily correspond to actual downflows. Dynamic features in the atmosphere can introduce complex features in the line profiles which will not be detected by instruments with the spectral resolution of EVE, but which leave more of a signature at the resolution of the Spectral Investigation of the Coronal Environment (SPICE) instrument on Solar Orbiter.

Observations of the Sun with the Atacama Large Millimeter Array have now started, and the thermal infrared will regularly be accessible from the NSF's Daniel K. Inouye Solar Telescope. Motivated by the prospect of these new data, and by recent flare observations in the mid infrared, we set out here to model and understand the source of the infrared continuum in flares, and to explore its diagnostic capability for the physical conditions in the flare atmosphere. We use the 1D radiation hydrodynamics code RADYN to calculate mid-infrared continuum emission from model atmospheres undergoing sudden deposition of energy by non-thermal electrons. We identify and characterise the main continuum thermal emission processes relevant to flare intensity enhancement in the mid- to far-infrared (2-200 $\mu$m) spectral range as free-free emission on neutrals and ions. We find that the infrared intensity evolution tracks the energy input to within a second, albeit with a lingering intensity enhancement, and provides a very direct indication of the evolution of the atmospheric ionization. The prediction of highly impulsive emission means that, on these timescales, the atmospheric hydrodynamics need not be considered in analysing the mid-IR signatures.

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.

We use radiation hydrodynamic simulations to examine two models of solar flare chromospheric heating: Alfv\'en wave dissipation and electron beam collisional losses. Both mechanisms are capable of strong chromospheric heating, and we show that the distinctive atmospheric evolution in the mid-to-upper chromosphere results in Mg II k-line emission that should be observably different between wave-heated and beam-heated simulations. We also present Ca II 8542A profiles which are formed slightly deeper in the chromosphere. The Mg II k-line profiles from our wave-heated simulation are quite different from those from a beam-heated model and are more consistent with IRIS observations. The predicted differences between the Ca II 8542A in the two models are small. We conclude that careful observational and theoretical study of lines formed in the mid-to-upper chromosphere holds genuine promise for distinguishing between competing models for chromospheric heating in flares.

We explore the response of the He II 304 {\AA} and He I 584 {\AA} line intensities to electron beam heating in solar flares using radiative hydrodynamic simulations. Comparing different electron beams parameters, we found that the intensities of both He lines are very sensitive to the energy flux deposited in the chromosphere, or more specifically to the heating rate, with He II 304 {\AA} being more sensitive to the heating than He I 584 {\AA}. Therefore, the He line ratio increases for larger heating rates in the chromosphere. A similar trend is found in observations, using SDO/EVE He irradiance ratios and estimates of the electron beam energy rate obtained from hard X-ray data. From the simulations, we also found that spectral index of the electrons can affect the He ratio but a similar effect was not found in the observations.

The bulk of the radiative output of a solar flare is emitted from the chromosphere, which produces enhancements in the optical and UV continuum, and in many lines, both optically thick and thin. We have, until very recently, lacked observations of two of the strongest of these lines: the Mg II h & k resonance lines. We present a detailed study of the response of these lines to a solar flare. The spatial and temporal behaviour of the integrated intensities, k/h line ratios, line of sight velocities, line widths and line asymmetries were investigated during an M class flare (SOL2014-02-13T01:40). Very intense, spatially localised energy input at the outer edge of the ribbon is observed, resulting in redshifts equivalent to velocities of ~15-26km/s, line broadenings, and a blue asymmetry in the most intense sources. The characteristic central reversal feature that is ubiquitous in quiet Sun observations is absent in flaring profiles, indicating that the source function increases with height during the flare. Despite the absence of the central reversal feature, the k/h line ratio indicates that the lines remain optically thick during the flare. Subordinate lines in the Mg II passband are observed to be in emission in flaring sources, brightening and cooling with similar timescales to the resonance lines. This work represents a first analysis of potential diagnostic information of the flaring atmosphere using these lines, and provides observations to which synthetic spectra from advanced radiative transfer codes can be compared.

This paper presents measurements of the energy radiated by the lower solar atmosphere, at optical, UV, and EUV wavelengths, during an X-class solar flare (SOL2011-02-15T01:56) in response to an injection of energy assumed to be in the form of nonthermal electrons. Hard X-ray observations from RHESSI were used to track the evolution of the parameters of the nonthermal electron distribution to reveal the total power contained in flare accelerated electrons. By integrating over the duration of the impulsive phase, the total energy contained in the nonthermal electrons was found to be $>2\times10^{31}$ erg. The response of the lower solar atmosphere was measured in the free-bound EUV continua of H I (Lyman), He I, and He II, plus the emission lines of He II at 304\AA\ and H I (Ly$\alpha$) at 1216\AA\ by SDO/EVE, the UV continua at 1600\AA\ and 1700\AA\ by SDO/AIA, and the WL continuum at 4504\AA, 5550\AA, and 6684\AA, along with the Ca II H line at 3968\AA\ using Hinode/SOT. The summed energy detected by these instruments amounted to $\sim3\times10^{30}$ erg; about 15% of the total nonthermal energy. The Ly$\alpha$ line was found to dominate the measured radiative losses. Parameters of both the driving electron distribution and the resulting chromospheric response are presented in detail to encourage the numerical modelling of flare heating for this event, to determine the depth of the solar atmosphere at which these line and continuum processes originate, and the mechanism(s) responsible for their generation.

White light flares (WLFs) are observational rarities, making them understudied events. However, optical emission is a significant contribution to flare energy budgets and the emission mechanisms responsible could have important implications for flare models. Using Hinode SOT optical continuum data taken in broadband red, green and blue filters, we investigate white-light emission from the X2.2 flare SOL2011-02-15T01:56:00. We develop a technique to robustly identify enhanced flare pixels and, using a knowledge of the RGB filter transmissions, determined the source color temperature and effective temperature. We investigated two idealized models of WL emission - an optically thick photospheric source, and an optically thin chromospheric slab. Under the optically thick assumption, the color temperature and effective temperature of flare sources in sunspot umbra and penumbra were determined as a function of time and position. Values in the range of 5000-6000K were found, corresponding to a blackbody temperature increase of a few hundred kelvin. The power emitted in the optical was estimated at $\sim 10^{26}$ergs s$^{-1}$. In some of the white-light sources the color and blackbody temperatures are the same within uncertainties, consistent with a blackbody emitter. In other regions this is not the case, suggesting that some other continuum emission process is contributing. An optically thin slab model producing hydrogen recombination radiation is also discussed as a potential source of WL emission; it requires temperatures in the range 5,500 - 25,000K, and total energies of $\sim 10^{27}$ergs s$^{-1}$.

In this paper, we present comprehensive analysis of a two-ribbon flare observed in UV 1600{\AA} by Transition Region and Coronal Explorer and in HXRs by Reuven Ramaty High Energy Solar Spectroscopic Imager. HXR (25-100 keV) imaging observations show two kernels of size (FWHM) 15?? moving along the two UV ribbons. We find the following results. (1) UV brightening is substantially enhanced wherever and whenever the compact HXR kernel is passing, and during the HXR transit across a certain region, the UV count light curve in that region is temporally correlated with the HXR total flux light curve. After the passage of the HXR kernel, the UV light curve exhibits smooth monotonical decay. (2)We measure the apparent motion speed of the HXR sources and UV ribbon fronts, and decompose the motion into parallel and perpendicular motions with respect to the magnetic polarity inversion line (PIL). It is found that HXR kernels and UV fronts exhibit similar apparent motion patterns and speeds. The parallel motion dominates during the rise of the HXR emission, and the perpendicular motion starts and dominates at the HXR peak, the apparent motion speed being 10-40 km s-1. (3) We also find that UV emission is characterized by a rapid rise correlated with HXRs, followed by a long decay on timescales of 15-30 minutes. The above analysis provides evidence that UV brightening is primarily caused by beam heating, which also produces thick-target HXR emission. The thermal origin of UV emission cannot be excluded, but would produce weaker heating by one order of magnitude. The extended UV ribbons in this event are most likely a result of sequential reconnection along the PIL, which produces individual flux tubes (post-flare loops), subsequent non-thermal energy release and heating in these flux tubes, and then the very long cooling time of the transition region at the feet of these flux tubes.