Spectral Evidence for Heating at Large Column Mass in Umbral Solar Flare Kernels I: IRIS NUV Spectra of the X1 Solar Flare of 2014 Oct 25

Spectral Evidence for Heating at Large Column Mass in Umbral Solar Flare Kernels I: IRIS NUV Spectra of the X1 Solar Flare of 2014 Oct 25

Adam F. Kowalski    Elizabeth Butler    Adrian N. Daw    Lyndsay Fletcher    Joel C. Allred    Bart De Pontieu    Graham S. Kerr    Gianna Cauzzi
\move@AU\move@AF\@affiliation

Department of Astrophysical and Planetary Sciences, University of Colorado Boulder, 2000 Colorado Ave, Boulder, CO 80305, USA. \move@AU\move@AF\@affiliationNational Solar Observatory, University of Colorado Boulder, 3665 Discovery Drive, Boulder, CO 80303, USA. \move@AU\move@AF\@affiliationLaboratory for Atmospheric and Space Physics, University of Colorado Boulder, 3665 Discovery Drive, Boulder, CO 80303, USA. adam.f.kowalski@colorado.edu \move@AU\move@AF\@affiliationDepartment of Astrophysical and Planetary Sciences, University of Colorado Boulder, 2000 Colorado Ave, Boulder, CO 80305, USA. \move@AU\move@AF\@affiliationNASA Goddard Space Flight Center, Heliophysics Sciences Division, Code 671, 8800 Greenbelt Rd., Greenbelt, MD 20771, USA. \move@AU\move@AF\@affiliationSUPA School of Physics & Astronomy, University of Glasgow, Glasgow, G12 8QQ, UK. \move@AU\move@AF\@affiliationRosseland Centre for Solar Physics, University of Oslo, P.O.Box 1029 Blindern, NO-0315 Oslo, Norway. \move@AU\move@AF\@affiliationNASA Goddard Space Flight Center, Heliophysics Sciences Division, Code 671, 8800 Greenbelt Rd., Greenbelt, MD 20771, USA. \move@AU\move@AF\@affiliationLockheed Martin Solar and Astrophysics Laboratory, Org. A021S, Bldg. 252, 3251 Hanover St., Palo Alto, CA 94304, USA. \move@AU\move@AF\@affiliationRosseland Centre for Solar Physics, University of Oslo, P.O. Box 1029 Blindern, NO-0315 Oslo, Norway. \move@AU\move@AF\@affiliationInstitute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029, Blindern, Oslo, Norway. \move@AU\move@AF\@affiliationNASA Goddard Space Flight Center, Heliophysics Sciences Division, Code 671, 8800 Greenbelt Rd., Greenbelt, MD 20771, USA. \move@AU\move@AF\@affiliationNational Solar Observatory, University of Colorado Boulder, 3665 Discovery Drive, Boulder, CO 80303, USA.

Abstract

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, % 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  intensity. This continuum-to-line ratio is a new spectral diagnostic of significant heating at high column mass (log [g cm) 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 (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.

1 Introduction

Flares result from a sudden magnetic reconfiguration in the atmosphere of a star, producing electrons and protons that stream at near the speed of light along the directions of the reconfigured magnetic fields. The flare is the burst of electromagnetic radiation and is thought to be the result of the impact of the (mildly) relativistic electrons (“beams”) with the lower, dense stellar atmosphere; the protons and ions produce nonthermal gamma ray emission (Murphy et al., 1997; Hurford et al., 2006), but their relative contributions to the multi-thermal response has not yet been established. The pan-chromatic continuum radiation, excluding the X-rays, extreme-ultraviolet, gamma rays, and radio emission, is collectively known as the white-light continuum radiation because it appears in broadband222Various loose definitions of a “white-light flare” exist, including a flare that could be detected by the eye (which is broadband). In dMe stars, -band flares are certainly considered white-light flares even though our eye is not sensitive to these wavelengths, and flares detected in a narrow passband of SDO/HMI are often called white-light flares, assuming that continuum radiation is the source of the HMI increase. Generally, a white-light flare is a flare that produces a change in continuum radiation that could be detected in the Johnson and/or bands.. The white-light is typically one of the most impulsive signatures in a solar flare (e.g., Hudson et al., 2006; Fletcher et al., 2007; Namekata et al., 2017; Watanabe et al., 2017) and is spatially and temporally correlated with the hard X-rays from the nonthermal electrons on the Sun (Kane et al., 1985; Hudson et al., 1992; Martínez Oliveros et al., 2012).

A common assumption is that the white-light originates from increased photospheric (or upper photospheric) radiation. However, chromospheric condensations (Livshits et al., 1981; Fisher, 1989; Kowalski et al., 2015b) with low continuum optical depth also produce broadband continuum radiation with a large jump in flux near the hydrogen Balmer limit (Gan et al., 1992; Kowalski et al., 2017a). The capability to test model predictions of the Balmer jump in solar flares has largely disappeared, which is unfortunate because spectra of solar flares in the 1980s (primarily from the Universal Spectrograph) exhibit a variety of characteristics in the Balmer jump spectral region (Hiei, 1982; Acampa et al., 1982; Neidig, 1983; Donati-Falchi et al., 1984, 1985; Boyer et al., 1985; Kowalski et al., 2015a; Procházka et al., 2017), and we now have methods to model opacity from blended lines and dissolved levels (Uitenbroek, 2001; Kowalski et al., 2015b). Therefore, we must employ other spectral diagnostics that critically test these models and constrain whether the white-light results from photospheric heating, as suggested by a recent off-limb measurement of the emission height (Martínez Oliveros et al., 2012, but see Battaglia & Kontar (2011) and Krucker et al. (2015)), primarily from chromospheric heating, or significant heating throughout several layers of the lower atmosphere (Neidig et al., 1993a, b; Kleint et al., 2016). Direct measurements of the white-light emission heights in limb flares may be partially occulted (H. Hudson, priv. communication), and these events are generally difficult to compare to 1D plane-parallel flare models (but see Heinzel et al., 2017). The specific continuum intensity over a narrow wavelength range has also been used to test electron beam heating in the lower atmosphere (Heinzel & Kleint, 2014; Kowalski et al., 2017a). However, the interpretation of the measured intensity can be degenerate since an optically thin continuum source with K produces a radiation temperature of K in the optical (Kowalski & Allred, 2018) and nearly 7000 K in the near-UV (Kleint et al., 2016). Furthermore, the white-light sources may be unresolved even with current high-spatial resolution capabilities (Krucker et al., 2011; Sharykin & Kosovichev, 2014; Kowalski et al., 2015a).

Is the photosphere heated to produce an observable amount of continuum radiation that overpowers the chromospheric flare radiation? Evidence for very deep heating exists in dMe flares from the Balmer continuum radiation that is produced in absorption (Kowalski et al., 2013) or in emission333The terminology “in emission” or “in absorption” means that the Balmer continuum flux is, respectively, greater than or less than a blackbody curve extrapolation from a fit to continuum windows at Å. with a relatively small jump in continuum flux in the Balmer jump spectral region (Hawley & Pettersen, 1991; Hawley & Fisher, 1992; Kowalski et al., 2013, 2016, 2019). The mass column density444Often used to indicate the atmospheric depth. (hereafter, “column mass”) that is heated to K to reproduce these spectral properties is far smaller than the dMe photospheric column mass of log /[g cm] (Cram & Woods, 1982; Kowalski et al., 2017b). In solar and dMe electron beam simulations that produce a flare chromosphere that is optically thin to Balmer and Paschen continuum radiation, radiative backwarming heats the (upper) photosphere but the temperature increase in these deep layers is only K (Allred et al., 2005, 2006; Cheng et al., 2010), and is not enough to explain the 10,000 K blackbody-like radiation in dMe flares (see also Kowalski et al., 2018). Furthermore, it takes some time ( s; appendix of Kowalski et al. (2017a)) for the solar photosphere to heat up in these models, and it is not clear if beam heating persists in a given flare loop for such long timescales (Aschwanden et al., 1998a; Nishizuka et al., 2009; Penn et al., 2016). Nonthermal electrons exhibit a power-law with most electrons at keV; if produced in the corona (as implied by hard X-ray, energy-dependent timing differences; Aschwanden et al., 1995; Aschwanden & Schwartz, 1995; Aschwanden, 1996; Aschwanden et al., 1996a, b, c, 1998b), most of the electron beam energy is expected to be lost in the mid-to-upper chromosphere (log to ; Vernazza et al., 1981). If electron beam heating is shown insufficient to explain impulsive, large temperature increases in the photosphere or low chromosphere, then there must be additional, important sources of heating in flares, such as Alfvén waves (Russell & Fletcher, 2013; Reep & Russell, 2016; Kerr et al., 2016; Reep et al., 2018) or proton beams (Zharkova & Zharkov, 2007; Procházka et al., 2018). Instead of Balmer jump ratio measurements, we must consider other spectral signatures that indicate significant heating at large column mass in solar flares. Additional heating mechanisms at large column mass would be transformative improvements to the standard solar flare model and would have important implications for the heating sources in superflares in other stars and the young Sun (e.g., Maehara et al., 2012; Osten et al., 2016).

To determine if low-energy electron beams are sufficient to explain the heating in the low atmosphere, we present a new measurement that can be obtained from spectra with limited wavelength coverage, as is often the case with solar imaging spectrometers. This new measurement is the ratio of the NUV continuum to the Fe II line intensity. The NUV continuum-to-line ratio will be compared directly to radiative-hydrodynamic modeling in Paper II of this series to infer the largest (deepest) column mass that is heated by energetic electrons (or other energy transport mechanisms) in a flare. The flare we study is the X1 flare on 2014 Oct 25, which produced ribbons that crossed the slit of the Interface Region Imaging Spectrograph (IRIS) during a hard X-ray event, which occurred during the decay phase of the GOES X1 event.

The paper is outlined as follows. In Section 2 we describe the IRIS and the Fermi/GBM data for the X1 flare. In Section 3, we describe the properties of the NUV continuum enhancements, the Fe II  line profiles, and the ratios between the NUV continuum intensity and Fe II  line-integrated intensity. In Section 4, we describe other characteristics of this flare that are important for radiative-hydrodynamic modeling: a line of He I that has not yet been seen in solar flare spectra and the nonthermal electron power and spectral index inferred from hard X-rays. In Section 5, we summarize our findings and compare the NUV spectral properties to the modeling and IRIS observations of another well-studied X-class flare. In Section 6, we conclude with several general implications from these unique solar flare spectra. Appendix A describes the continuum-to-line ratio dependence on the temperature and density in optically thin slab model approximations, and appendix B contains identifications of observed flare lines in the IRIS NUV and discussion of the observational constraints they provide for models.

2 Data

The 2014-Oct-25 X1 solar flare was observed with a unique dataset that included a custom IRIS observing mode, which is described in Section 2.1. Hard X-ray observations were provided by Fermi/GBM throughout the entire flare as well (Section 2.2). This comprehensive dataset allows us to test RHD modeling of electron beam heating in new ways (in Paper II).

2.1 IRIS Ultraviolet Data

The Interface Region Imaging Spectrograph (IRIS; De Pontieu et al., 2014) observed the GOES X1 flare SOL2014-10-25T17:08:00 (, , ) from National Oceanic and Atmospheric (NOAA) Active Region (AR) 12192. The X-class flares from AR 12192 have been studied extensively due to the absence of associated coronal mass ejections (Thalmann et al., 2015; Inoue et al., 2016; Bamba et al., 2017; Amari et al., 2018). The IRIS spectra of the 2014-Oct-25 X1 flare were obtained in a “sit-and-stare” mode with an exposure time of 4 s, a temporal cadence of 5.4 s, and the slit oriented in the E-W direction. The slit jaw images (SJI) in the Mg II wing (hereafter, SJI 2832), C II, and Mg II were obtained with a cadence of 16 s. All spectra and slit jaw images were calibrated from DN s to specific intensity (hereafter, just “intensity”) following Kleint et al. (2016) and Kowalski et al. (2017a) by using the time-dependent effective area curves (Wülser et al., 2018). By comparing the fiducial marks in the FUV and NUV data, we found that a (downward) shift by one pixel in the FUV spectra was necessary to align them with the NUV spectra. We found that a shift of Å was required to calibrate the velocity (using the Ni 2799.474 Å line, as suggested by IRIS Technical Note 26) in the level 2 data that we used for our analysis. We note that more recently calibrated data does not require such a shift, since the whole IRIS archive has been reprocessed in 2017/2018 including an improved wavelength calibration. In general, it is always best to check the calibration of the wavelength using the instructions provided in IRIS Technical Note 26.

Usually, to strike the best balance between cadence and data rate, only narrow wavelength regions around specific lines of interest are read out from IRIS. Instead, we obtained the observations with a custom observing mode that downlinked the full spectral range of IRIS with two-pixel binning in the dispersion direction ( Å and Å at 0.0256 Å pixel in the FUV and Å at 0.051 Å pixel in the NUV) and binning by two in the spatial direction (thus, 0.33 binned-pixel corresponding to 250 km binned-pixel at the distance of the Sun). The full IRIS spectral range was used for a robust identification and characterization of the FUV and NUV continua in flares; with limited wavelength ranges of typical IRIS readout modes (“linelists”), it is not always possible to accurately determine the continuum intensity due to multiple broad, asymmetric flares lines of Mg II, Fe II and other species. The NUV spectral range of IRIS (hereafter, IRIS NUV) has a relatively line-free continuum region at Å that has been useful as a proxy of the Balmer continuum component of the white-light radiation for comparisons to the predictions of electron beam heating models (Heinzel & Kleint, 2014; Kleint et al., 2016; Kowalski et al., 2017a). In this paper, we calculate this quantity and compare to the continuum level over the entire spectral range of IRIS.

2.2 Fermi/GBM X-ray Data

The Fermi/Gamma-Ray Burst Monitor (GBM; Meegan et al., 2009) provides high-time resolution (4s, or 1s in burst mode) X-ray data of solar flares from 8 keV to 40 MeV. RHESSI observed the main peak of the flare, but the imagery indicates that the loop-top source dominated the 25-50 keV emission, and RHESSI also suffered from significant pileup (Kleint et al., 2017). We use the Fermi/GBM data late in the flare when pileup was not significant in several of the least sunward facing detectors.

The Fermi/GBM has two types of scintillation detectors: twelve NaI detectors that are sensitive at keV and two bismuth germanate (BGO) detectors that are sensitive at keV. For spectral analysis (Section 4.2), we use data from detector NaI_n0 (hereafter, “n0”), which is the third most sunward-facing detector and does not suffer from pile-up effects555NaI_n05 is the most sunward facing detector and suffers from obvious pileup at the peak of the hard X-ray event. over the time of the decay phase of this flare. For light curve analysis, we also use data from detector NaI_n4 and NaI_n5 (hereafter, “n4” and “n5”, respectively), which are more sunward-facing in the early rise phase of the flare. The Solar Soft (SSW) IDL OSPEX is used to retrieve and analyze the Fermi/GBM data. First, we retrieve a detector response (rsp2) file for the time interval of 16:35-17:36. The background is chosen before and after the X1 event (during Fermi night) in the time-intervals of 16:15-16:35 and 17:50 - 18:05. Since the detector angles relative to the Sun change over the first few minutes of this flare666The NaI_n0 detector changes from the least sunward facing orientation to the third most sunward facing orientation from 17:01 to 17:05, which affects the rise phase of the hard X-rays., we divide the OSPEX count rate by the cosine of the angle between the detector normal and the Sun.

2.3 Overview of the Flare

The Fermi/GBM data from detector n4 at hard X-ray energies, keV, are shown in Figure 2.3(a), and the light curve exhibits four main hard X-ray peaks at 16:59:17, 17:00:23, 17:03:31, and 17:17:27 with comparable amplitudes777The main peak at 17:04 may suffer from some pileup, thus making the 58-72 keV light curve an upper limit on the flux during this peak.. Also shown are the soft X-rays ( Å) from GOES888Bamba et al. (2017) claim there is an error in the time-registry for these GOES data, but we do not see evidence to support this..

The SOL2014-10-25T17:08:00 flare was a three-ribbon flare with two ribbons in the sunspot plage and a third ribbon crossing two umbrae within a triple-umbral complex (Bamba et al., 2017). The X1 flare was the largest and the last of three flares that were triggered by an intruding positive polarity to the east of the main polarity inversion line, in accord with the triggering mechanism of Kusano et al. (2012). Following the nomenclature for the ribbons used in Bamba et al. (2017), the XR3 ribbon develops into the umbrae and crosses the IRIS slit during the decay phase of the fourth Fermi/GBM hard X-ray peak in Figure 2.3(a). The XR3 ribbon is seen in Figure 2.3, where we show three context images at 17:19: panel (a) shows the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA; Lemen et al., 2012) 1700 image with the IRIS field-of-view (FOV) indicated, panel (b) shows an IRIS SJI C II image, and panel (c) is the IRIS SJI 2832 image. The XR3 ribbon is very faint in SJI 2832, but the contrast against the umbral intensity is large. Li et al. (2018) found Fe XXI emission blueshifted by 143 km s in the IRIS/FUV spectra of the XR3 ribbon. The high temperature, low-density plasma are not analyzed further here. A comparison of flare intensity in the IRIS NUV spectra and SJI 2832 images for the X1 flare has been discussed in Kleint et al. (2017), who found that the relative contribution of Balmer continuum emission to the SJI 2832 radiation is nearly 70% during the times when the XR3 ribbon crosses the slit (see the middle panels of their Figure 6). In this paper, we present a detailed analysis of the IRIS spectra of the umbral flare ribbon XR3.

3 A New Diagnostic for Deep Atmospheric Heating

3.1 Identification of Near-UV Continuum Flare Kernels

In this section, we present the C II ribbon development and describe how we identify the largest NUV continuum enhancements, which occur in the late phase when the XR3 ribbon develops into an umbra. The large contrast of the flare ribbons in umbrae (which are rather unusual environments for flare ribbon development) provide more robust comparisons to RHD models than flare spectra with low intensity contrast in the plage regions, which vary significantly due to the granulation evolution.

The flare ribbons crossed the slit at many times and locations in the 2014-Oct-25 X1 flare. The intensity integrated over each C II line from to km/s is shown in Figure 2.3(b) and readily indicates the times when, and locations where, the XR1, XR2, and XR3 flare ribbons cross the IRIS slit. However, the detailed formation of C II in flares has yet to be investigated for a range of heating levels outside of quiet-Sun conditions (Rathore & Carlsson, 2015). Since C II has previously been analyzed extensively in IRIS data of solar flares (e.g., Tian et al., 2015; Sadykov et al., 2016), and they provide a comprehensive, high-contrast overview of the flaring ribbon evolution, we show their line-integrated evolution here but do not return to a detailed analysis of these lines in this paper.

We identify candidate NUV continuum enhancements (and by proxy, white-light-emitting kernels) with a running difference of the line center intensity of Fe ii at Å, which we expect to brighten with the NUV continuum due to their similar formation conditions in flares (Kowalski et al., 2017a). Detailed justification for using Fe II to find candidate continuum enhancements is given in Appendix A. The time-evolution of the NUV continuum region from Å to 2825.9 Å is shown in Figure 2.3(c) over a limited time range. Following Kowalski et al. (2017a), we calculate the average intensity in this wavelength window as C2826999The value of C2826 was defined as the average intensity from Å in the 2014-Mar-29 spectra in Kowalski et al. (2017a). In this study, we extend the blue end of the wavelength window for C2826 to 2824.5 Å, which is line-free in the 2014-Oct-25 flare.; the XR3 flare ribbon in C2826 crosses the IRIS slit and sunspot umbra from 17:17-17:19 UT in Figure 2.3(c) and is further evident in SJI 2832 in Figure 2.3(c). Subtracting the preflare intensity gives the “excess C2826”, hereafter C2826 (following the “prime” terminology for flare-only flux in Kowalski et al. (2012) and Kowalski et al. (2019)) and dividing by the pre-flare values gives the “C2826 enhancement” (following the terminology in Kowalski et al., 2015a).

The largest values of the C2826 (intensity) and the C2826 enhancement occur in the sunspot umbra at ()(364.8, -318.4) at 17:18:24 and 17:19:02 UT during the decay of the fourth hard X-ray event that peaks at 17:17:30. Hereafter, we refer to the brightest flaring pixels in the umbra during this hard X-ray event as the “umbral flare brightenings” (UFBs). The time of UFB-3 is indicated in Figure 2.3(a) for reference. The peak of UFB-1 occurs at 17:17:31 at the spatial location of (). The peaks of UFB-2 at 17:18:24 and UFB-3 at 17:19:02 UT occur in nearly the same location at (). However, UFB-2 appears to be less than one pixel to the west (in the positive direction) of UFB-3, and the apparent motion of the ribbon eastward along the slit is evident from the time of UFB-1 to the time of UFB-3 (see Figure 2.3(c)). Hereafter, we focus our analysis on the spectra of UFB-2 and UFB-3, which produced the largest and fastest NUV continuum enhancements in the IRIS spectra.

3.2 Continuum Properties of the Umbral Flare Brightenings

The time-evolution of C2826 averaged over three spatial pixels (the spatially averaged intensity, ) centered at ()(364.8, -318.4) is shown in Figure 3.2. The UFB-2 and UFB-3 events correspond to the first and second C2826 peaks, respectively. The C2826 enhancement remarkably exceeds the pre-flare umbral intensity by a factor of ten. The keV X-ray light curve is overplotted in Figure 3.2, and UFB-3 occurs when the ribbons cross the slit in the decay phase of this hard X-ray peak. From the range of calculated FWHM values of the C2826 light curve of UFB-3, we constrain upper limits on the duration of heating to  s over a solid angle of 0.33 (slit width) (0.40; NUV spatial resolution), or an area of cm. However, the data do not exclude shorter heating durations over smaller areas. Furthermore, the rise phase of UFB-3 is unresolved since it consists of one temporal point. However, we note  s is generally the timescale range for soft X-ray derivative bursts characterized in another flare (Rubio da Costa et al., 2016). Several heating episodes may occur (at least partially) in the same location (e.g., UFB-2 and UFB-3), but this cannot be quantified in any detail even with the impressive spatial resolution of IRIS since the ribbons move along the slit from the time of UFB-1 to the time of UFB-3.

The IRIS NUV spectrum at the peak of UFB-3 is shown in Figure 3.2(a). This spectrum exhibits the most robust characterization of the IRIS NUV continuum to date during a solar flare. The flare spectrum shows that the NUV continuum intensity level extrapolated (over the horizontal red line) from C2826 extends over a spectral range of Å. The horizontal line extrapolation provides an estimate for the relative, wavelength-integrated continuum energy in the NUV spectral range of IRIS. We find that only % of the UFB-3 peak flare kernel brightness in the IRIS NUV spectral range is attributed to the continuum radiation, while the majority of the wavelength-integrated kernel brightness (nearly 70%) is due to the Mg II lines.

The C2826 intensity from IRIS during other flares has been used to compare directly to radiative-hydrodynamic flare models (Heinzel & Kleint, 2014; Kleint et al., 2016; Kowalski et al., 2017a). But are the C2826 ribbons in the 2014 Oct 25 X1 flare spatially resolved for a direct comparison of the value of C2826 ( erg s cm sr Å; Figure 3.2) to model snapshots? We use the iris_sg_deconvolve.pro routine to deconvolve the UFB-1 and UFB-3 peak spectra by the IRIS PSF as described in Courrier et al. (2018). We find that after iterations of the deconvolution procedure, the NUV spectra increase in brightness by factors of . The IRIS NUV resolution is 0.4 or 1.2 (binned) pixels, and the spatial full-widths-at-half-maxima of UFB-1, UFB-2, and UFB-3101010UFB-3 occurs adjacent to a fiducial mark in the spectra (Figure 2.3b), and we are unable to quantify the full spatial extent of the kernel. are 3-6 pixels, corresponding to km at the Sun; in Figure 2.3(c), it is apparent that the NUV continuum radiation (C2826) from the UFBs is bright over a few pixels only. In higher resolution H data of umbral flare ribbons, very narrow widths of km are observed (Sharykin & Kosovichev, 2014), which suggests there may be significant unresolved structure in our IRIS observations. Higher-quality observations (i.e., without saturation) of the widths and motions of plage and umbral kernels will be important to obtain with the Daniel K. Inouye Solar Telescope at very high time resolution.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref C2826 light curve of the umbral flare brightenings UFB-2 and UFB-3 at ()(364.8, -318.4) arcsec during the fourth major hard X-ray peak, over the same time range as indicated by the vertical dotted lines in Figure 2.3(a,b). The peaks of UFB-2 and UFB-3 are delayed with respect to the peak of the hard X-rays in Fermi/GBM (olive) because the Fermi/GBM is not spatially resolved at the Sun and the ribbons cross the IRIS slit at this spatial location at a later time. The continuum-to-line ratios (C2826/Fe II ; asterisks, right axis) exhibit very large values of in the impulsive phases (from 17:18:03 to 17:19:34) of UFB-2 and UFB-3. The light curve of C2826 at the spatial location of UFB-1 is shown in a lighter gray tone; UFB-1 reaches its maximum closer in time to the HXR peak.

3.3 The Fe II 2814.45 Å line in umbral flare brightenings

The Fe II  emission line has been recently studied with radiative-hydrodynamic models of flares (Kowalski et al., 2017a) and has revealed two emission line components: a line component centered near the rest wavelength and a spectrally resolved, broad emission component to the red of the rest wavelength. Compared to other chromospheric flare lines, the Fe II  line exhibits one of the lowest opacities, and the optical depth does not build up as fast in that line as in the Balmer lines or in Mg II. Thus, the properties of Fe II  constrain the dynamics in deep layers of the flaring chromospheric condensation and stationary flare layers below (Kowalski et al., 2017a), where the white-light and IRIS NUV continuum radiation may originate. The similar temperature and density dependencies of Fe II  and hydrogen Balmer bound-free radiation are shown in Appendix A.

In Figure 3.2(b), we show an inset of several flare emission lines around the region of the Fe II  line in UFB-2 and UFB-3. We measure the bisector at 30% maximum for the Fe II  Å emission line to be km s. This is just greater than a 1 redshift detection111111Adding a systematic uncertainty of 1.8 km s and a statistical uncertainty of 2.3 km s in quadrature.. However, the line profile shape exhibits asymmetric brightening at such that there are more spectral bins with an intensity greater than or equal to the line half-maximum on the red side than on the blue side. We show vertical dashed lines at and vertical dotted lines at km s to indicate the red extent of this profile asymmetry. The Fe II at Å line exhibits a similar asymmetric profile and redshift ( km s) to the Fe II at Å line, while the Fe I Å line exhibits no evidence of a redshift or an asymmetric profile. Furthermore, UFB-2 and other flare spectra show similarly redshifted and asymmetric Fe II  profiles to UFB-3. These three lines are of particular utility for studying flares, as they fall within a standard window common to the IRIS NUV line lists, and thus provide the opportunity to compare with a large archive of observations. With full spectral readout, there are numerous additional Fe II and Cr II lines with similar profiles available, as discussed in Section 4.3 and Appendix B. The asymmetric, redshifted profiles of these Fe II lines and other emission lines can be synthesized from RHD models to constrain the chromospheric velocity field in white-light emitting layers at high-time resolution.

3.4 The continuum-to-line ratio in IRIS spectra

The umbral flare brightening spectra reveal that the ratio of the NUV continuum radiation to the Fe II line intensity is an interesting diagnostic in flares. This continuum-to-line ratio can be readily obtained from spectra with narrow wavelength coverage; like an equivalent width, it is independent of uncertainties (e.g., Wülser et al., 2018; Courrier et al., 2018) in the absolute intensity calibration. We define the continuum-to-line ratio as the excess C2826 (C2826) divided by the wavelength-integrated, continuum-subtracted, and preflare-subtracted intensity of the Fe II  line. The continuum-to-line ratio is hereafter denoted as C2826/Fe II  (with units of Å). The evolution of this quantity varies between over the peaks of UFB-2 and UFB-3 and is shown in Figure 3.2. At the peak of UFB-1 (not shown), the continuum-to-line ratio is also large near . Subtracting the last decay phase spectrum of UFB-2 from the peak of UFB-3 at the spatial location of the light curve in Figure 3.2 gives the newly-brightened flare radiation assuming that UFB-2 and UFB-3 occur at different locations within the spatial resolution. For this spectrum, the value of C2826/Fe II  is even larger (). A comparison of the C2826/Fe II  ratios among many flares will be presented in Butler et al. (2019, in prep.).

The Fe II  line in the IRIS NUV spectrum is an important diagnostic in flares because the line is produced (assuming LTE) over a similar temperature range of K as hydrogen bound-free radiation (Kowalski et al., 2017a, and shown in detail in Appendix A), which contributes to the C2826 intensity in flares (Heinzel & Kleint, 2014). Thus, the ratio of these two quantities reflects the amount of significant heating (T K) in deep layers of the atmosphere that exhibit for continuum radiation and for Fe II  radiation. The ratios from optically thin, uniform slabs with low to moderately high densities ( g cm) at K are (Appendix A), which are far less than the observed values of in the UFBs. Radiative-hydrodynamic modeling and the detailed interpretation and analysis of large values of C2826/Fe II  in the UFBs are outside the scope of this paper and will be presented in Paper II.

4 Other Observational Properties of the Flare

In order to model the UFBs and the continnum-to-line ratios with electron beam heating simulations in Paper II, several comprehensive input parameters (electron beam energy and power-law index) and other constraints from the emission lines are presented in this section. In Section 4.1, we calculate the flare area from high spatial resolution imagery of the flare ribbons. This flare area is combined with the nonthermal electron power inferred from Fermi/GBM data in Section 4.2 to give the energy flux density of beams assumed to heat the lower atmosphere. In Section 4.3 and appendix B, we discuss properties of the many other Fe II emission lines in the UFBs, as well as species such as Cr II and Fe I, and the constraints they provide on atmospheric models. In Section 4.4, we discuss the properties of a He I emission line in the IRIS NUV that will complement the constraints from the Fe II lines.

4.1 Evolution of the Flare Ribbons

High spatial resolution flare area estimates are critical for obtaining an accurate heating rate for modeling the RHD response (Krucker et al., 2011). However, the continuum contrast outside the umbra is low, which causes flare area measurements to be difficult in solar flares for even large flare energies. Generally, the non-flaring photospheric intensity in the NUV Mg II h wing is 50x brighter than the pre-flare umbral NUV intensity, and the large amplitude, gradual variation from photospheric convection outside of the sunspot makes the identification of bona-fide flare continuum enhancements difficult (we refer the reader to Kleint et al., 2017, for some properties of NUV enhancements away from the sunspots).

We present an algorithm that identifies pixels that brighten impulsively in flares. This algorithm excludes gradual evolution of the solar granulation that can be falsely identified as flare kernels, and the method is applicable to other data sets with high-time and high-spatial resolution. Using the intensity-calibrated SJI 2832 images (with a cadence of 16 s), we set a threshold excess intensity (T.E.I.) value for impulsive pixels in SJI 2832 to erg cm s sr Å, which corresponds to % of the maximum excess intensity observed over the entire observation sequence. We consider here only the umbral region indicated in the SJI 2832 image by the box in Figure 2.3(c). A pixel at time in image is counted as an impulsive flare pixel if the following criteria are met: 1) the excess intensity in this pixel at is greater than T.E.I.; 2) the excess intensity in this pixel at  s () or  s () is less than 0.5 x T.E.I; and 3) the excess intensity at  s () is less than 0.5 x T.E.I. These criteria identify newly flaring, impulsive pixels in SJI 2832 by efficiently excluding a significant number of pixels that exhibit a much more gradual evolution with sustained excess values. It is not possible to tell whether these gradual brightenings are due to changes in the granulation or are related to the flare, as in the so-called “type II white-light flare” variations (Matthews et al., 2003; Procházka et al., 2017). Hereafter, we refer to the impulsive flare areas using this algorithm as “newly-brightened” flare areas.

The number of newly-brightened flare pixels (where each pixel corresponds to an area of cm at the Sun) in the umbral region of SJI 2832 is shown as a function of time in Figure 4.1 and is compared to the Fermi/GBM keV and keV light curves. The newly-brightened flare area and hard X-ray emissions follow each other closely (see Graham et al., in prep), which suggests that the heating in these pixels (over a small region in the flare) is related to the non-thermal electrons that produce the fourth hard X-ray peak121212This also shows that the fourth hard X-ray peak in Fermi/GBM is a result of the X1 flare ribbons and not due to a particle event with a non-solar origin.. In the decay phase of the fourth hard X-ray peak, the newly-brightened ribbon areas have decreased as they progress through the umbra. Thus, there are impulsive flaring pixels in the gradual decay of the hard X-ray peak (which occurs in the gradual decay of a much larger soft X-ray peak). From Figure 4.1, we obtain a range of newly-brightened flare areas to be cm (17:19; at the peak UFB-3) to cm (17:17:30; at the peak of the fourth HXR peak). This information is combined with constraints from the Fermi/GBM spectrum in Section 4.2 to infer the energy flux in electron beams.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref The newly-brightened (impulsive) flare area in SJI 2832 and the Fermi/GBM hard X-ray light curve during the fourth hard X-ray peak. The conversion from pixels to area for these observations (which employ on-board binning) is cm / pixel. The UFB-3 occurs during the late phase of the fourth hard X-ray peak, but the evolution of the flare area in SJI 2832 follows the hard X-ray light curves. Note, the newly-brightened flare area in SJI 2832 is calculated from the region indicated by the box indicated in panel (c) of Figure 2.3.

Whereas the Fermi/GBM hard X-rays originate from the entire flaring area, the SJI 2832 data used to calculate the flare area are rather limited in field of view (see Figure 2.3(c)). Thus, we use our algorithm131313To calculate newly-brightened emission in SDO/AIA 1700 at time we used running differences of  s,  s, and  s, and we added the number of pixels that were below 50% of a high threshold value, which we took to be 30% of the saturation ( counts). to calculate the newly-brightened flare area in SDO/1700, which provides a complete field-of-view characterization of the flare region but with lower spatial resolution. The newly-brightened flare area from SDO/1700 is cm at the time of the peak of the fourth hard X-ray light curve, and it falls to cm at the times of the (spectral) UFBs. These areas are 10-40x larger than obtained from the umbral region in SJI 2832. Note that several bright kernels in the SDO images are saturated during the flare, making the areas from SDO/1700 over-estimates of the true flaring area at some times.

We also investigate whether the bright flare sources that cross the IRIS slit (UFB-1, UFB-2, or UFB-3) are obviously brighter than other newly-brightened flare pixels in the SJI 2832 images within the box indicated in Figure 2.3(c). Many of the brightest SJI 2832 pixels occur near the maximum of the HXR light curve in Figure 4.1. To determine how the UFBs from the spectra compare to the brightest flare excess pixels that do not cross the slit, we fold the IRIS spectra at the peaks of UFB-2 and UFB-3 with the effective area curve of SJI 2832 (Figure 3.2) to produce the synthetic values of SJI 2832; the synthetic excess intensity values141414The synthetic SJI 2832 intensity values from the UFB spectra are not averaged over three spatial pixels in this analysis. are and erg cm s sr Å for UFB-2 and UFB-3, respectively, while the brightest SJI 2832 excess value is erg cm s sr Å. There are 117 pixels (or 23% of the pixels brighter than the T.E.I.) in the SJI 2832 images that become as bright as or brighter than the UFB-2 spectrum, and there are 35 (or 7% of the pixels brighter than the T.E.I.) greater than or equal to the brightness of UFB-3. Clearly, a significant number of pixels that do not intersect the IRIS slit are brighter than the synthesized SJI 2832 values from the UFB-2 and UFB-3 spectra, suggesting that the continuum brightness of the umbral flare brightenings (UFB-2 and UFB-3) that cross the IRIS slit are not uniquely large.

4.2 Fermi/GBM Spectral Analysis

The hard X-rays from Fermi/GBM indicate the presence of nonthermal electrons in the flare, which are thought to produce the heating responsible for the white-light (e.g., C2826) increase. To determine if nonthermal electrons (as constrained from Fermi) are sufficient to explain the IRIS spectra, we infer the properties of a nonthermal electron beam using the standard thick-target formulae. The nonthermal electron distribution is parameterized by a power-law index (), an energy flux density (given by erg cm s; hereafter, abbreviated by F), and a low-energy cutoff () of the distribution, which is almost always an upper limit. Following the thick-target modeling of Milligan et al. (2014), we fit the Fermi/GBM spectrum at keV with a sum of thick target bremsstrahlung emission, an albedo correction from Compton scattering off the photosphere (Kontar et al., 2006), and a multi-thermal power-law spectrum. The fits were done using the SolarSoft OSPEX software in IDL and are shown in Figure 4.2 for the time-interval 17:17 - 17:20 (indicated by vertical dotted lines in Figure 3.2). Varying the fit interval to two minutes within this window gives similar results. We fit all energy intervals at E keV using the same background times. Because it is known that there is an iodine -edge (at 33 keV) in the detector (Meegan et al., 2009) that causes spectral fitting from keV to be problematic, we also fit only keV with a thick target model but find a similar power law. We conclude that the hard X-rays are consistent with a thick-target electron spectrum that has a power-law index of and an upper limit to the low-energy cutoff of keV. Note, the Compton scattering and photo-electric absorption of outgoing hard X-ray photons from the heated chromosphere has not been taken into account in the X-ray spectral fits, whereas a limb-observation has shown that the HXR source (in another flare) is low enough to make the Compton opacity significant (Martínez Oliveros et al., 2012). The extinction of hard X-rays in the evolved atmosphere will be considered in the radiative-hydrodynamic flare model atmospheres in Paper II.

At 60 keV there is clearly a break in the spectrum in Figure 4.2. While there is a transient burst of hard X-ray and gamma ray photons at keV to MeV that corresponds to the evolution of the C2826 light curve in Figure 3.2, this gamma-ray event cannot be constrained to originate from the Sun since Fermi was approaching the South Atlantic Anomaly at these times. Thus, the break is likely due to variable background radiation.

The last input parameter for RHD modeling is an estimate of the energy flux density (F) of accelerated electrons. Using the fits in Figure 4.2, we obtain a nonthermal power of erg s, which is an average over the time interval. A fit to the times of the fourth HXR peak at 17:17:20-17:17:40 gives a power of erg s (with similar cutoff and power-law indices). The hard X-ray count rate at the time of UFB-3 is 1/10 of the max count rate in the fourth peak (see Figure 2.3(a) and Figure 4.1), which allows us to estimate that the nonthermal electron power is erg s at the time of the UFBs in the IRIS spectra. The area estimated in Section 4.1 from SDO/1700 gives a flux density estimate of erg cm s, while the area from SJI 2832 gives erg cm s. We thus obtain a range of 2F102F11 for the nonthermal electron energy flux density for input to RADYN modeling.

The power-law index of the hard X-ray emission is very soft at these late times in the flare. The analysis of RHESSI data for this flare also indicates a large power-law index at earlier times (Kleint et al., 2017). The RHESSI imagery suggest the hard X-ray (25-50 keV) source is coronal, which may be due to a thick-target nonthermal coronal source or due to pileup of low-energy photons from a thermal coronal source. Because Fermi/GBM is spatially unresolved, we cannot rule out that the keV emission originates from a thick target coronal source (Veronig & Brown, 2004), but we have ensured that we use data from the detectors of Fermi that do not suffer from pileup of lower energy photons. We note that Thalmann et al. (2015) used RHESSI data of other flares from AR 12192 and find very soft electron spectra of in the decay phase and very large nonthermal electron energies as well. There are many other examples of soft power-law indices in flares in Solar Cycle 24 (Milligan et al., 2014; Thalmann et al., 2015; Kerr et al., 2015; Warmuth & Mann, 2016).

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref A multi-component fit to the keV Fermi/GBM (NaI n0 detector) spectrum from 17:17-17:20 with therm_multi_pow albedo thick2_vnorm (“total”). Parameters as returned by OSPEX for each component are given in the legend. Fitting during times and using other methods of background subtraction give similar parameters. Note, the keV energy range is excluded from the fit due to known calibration issues with Fermi/GBM.

4.3 Other Emission Line Properties of the Umbral Flare Brightenings

The full spectral readout employed for these data allow a unique and comprehensive identification of flare emission lines in the IRIS NUV wavelength range. There are many lines of Fe II in the flare spectrum in Figure 3.2, as well as for other species, and multiple lines from the same species provide powerful capability to disambiguate blends. A number of line identifications in the flaring NUV spectrum are shown in Figures 4.3 and 4.3, and a full list of identifications and intensities () at the peak of UFB-3 are given in Appendix B. While the strongest lines are those of Mg II, there are numerous strong lines of Fe II, Cr II, Fe I, and other species, for which ratios of intensity among emission lines (e.g., Fe II /2832.39) can also provide valuable constraints on models.

Another valuable capability provided by the full spectral range of these data is to characterize the contributions to SJI 2832 images in detail. The spectral range from Figure 3.2 at the peak of UFB-3 is enlarged in Figure 4.3 to show the line and continuum contributions. The brightest emission lines in the SJI 2832 bandpass are: Fe II 2826.58, 2826.86, 2828.26, 2828.73, 2829.46, 2831.75, and 2832.39; He I 2829.91 (see Section 4.4); Cr II 2831.30 and 2833.29; Fe I 2826.39; and Ti II 2833.015. Accurate SJI 2832 predictions from RHD models of flares need to include the contribution from these lines in addition to the continuum radiation. We calculate the filter-weighted specific intensity (Sirianni et al., 2005) from the spectrum with the SJI effective area curve. This calculation indicates that % of the flaring SJI 2832 count rate is due to the NUV flare continuum radiation, which is consistent with the findings of Kleint et al. (2017) for these data. In Figures 4.3 and 4.3, we show the estimated C2826 continuum level lowered by 5%, which seems to slightly better represent the continuum intensity between emission lines at Å.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref Enlarged view of the UFB-3 peak flare spectrum over the wavelength range of SJI 2832. Many Fe II and Cr II lines exhibit similar profile asymmetries, and the He I 2829.91 Å line is redshifted and broad. We adjust the value of C2826 by 5% to show that it better represents the excess NUV continuum intensity at Å.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref Some line identifications of the brighter lines in the UFB-3 peak flare spectrum at shorter NUV wavelengths around Mg II. We adjust the value of C2826 by 5% to show that it better represents the excess NUV continuum intensity at Å.

A number of lines in the UFB-3 peak excess intensity spectrum are identified and discussed further in Appendix B. Briefly, rest wavelengths for Fe II, Cr II, and Fe I are adopted primarily from Nave & Johansson (2013), and Sansonetti & Nave (2014), and Nave et al. (1994), respectively, with additional data on wavelengths and -values obtained from the NIST (Kramida et al., 2018), Kurucz (2018), and R. L. Kelly (https://www.cfa.harvard.edu/ampcgi/kelly.pl) databases and references therein. Identifications were verified by: (1) comparing line positions and profiles to those from the same species, and (2) comparing intensities from LTE calculations at different temperatures to the observed intensity ratios for each species, with opacity considered as discussed below. In Appendix B, we give measured wavelengths in this flare spectrum using 30% bisectors. Most bisector values for Fe II, Cr II, and Mg II are systematically redshifted by an amount greater than or approximately equal to the Fe II  line discussed in Section 4.3.

Nineteen Fe II line profiles, comprised by 24 of the 31 Fe II lines listed in Appendix B, are shown for the UFB-3 peak in Figure 4.3. These are reasonably strong, isolated features with either no evidence of significant blends from other species, or minor blends that can be accounted for, since multiple other lines are observed from those species to constrain any impact the blends would have. This provides many ratios between Fe II line intensities that can be used to constrain model atmospheres. Compared to the ratios in the flare, the LTE intensity ratios are systematically over-predicted for the strong Fe II lines. We speculate that optical depth in the brighter Fe II prevents a relatively larger amount of emission from escaping from the stationary flare layers (below the chromospheric condensation; see Kowalski et al., 2017a) in these lines. This speculation is supported by the fact that for the two brightest Fe II lines, 2784.512 and 2832.394, the rest component is less pronounced relative to the red wing than it is for the weaker Fe II lines. These two Fe II lines also exhibit the largest 30% bisectors (7-8 km s) among all Fe II line bisectors ( to km s; see Appendix B) that are reliable measured. This evidence for optical depth attenuation from the Fe II line ratios and profiles is consistent with the evidence from the continuum-to-line ratio (Appendix A, Section 3.4, Paper II). Detailed examination of this hypothesis can be done by simulating all the observed Fe II lines and lines of other species as discussed further in Appendix B.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHrefLine profiles in the NUV spectrum of UFB-3 for all definitively identified Fe II lines. All wavelength scales are from -0.2 to 0.4 Å from the rest wavelength, where the stronger Fe II line is used when two are present. Rest wavelengths of Fe II lines and other blends are indicated by dashed and dotted lines, respectively. Many of the Fe II profiles exhibit a red-wing asymmetry with more spectral bins with an intensity greater than or equal to the line half-maximum on the red side than on the blue side .

4.4 A New Flare Emission Line: He I 2829.91

We noticed an emission line in the SJI 2832 bandpass that is broader than other lines in this wavelength range. This “line” is the closely-spaced He I 2829.91 Å multiplet, which has been previously identified in the laboratory (Kramida et al., 2018; Drake, 2006) and in absorption in the spectra of hot (O- and B-type) stars (Dufton & McKeith, 1980; Fahey, 1984). This transition results from a high-lying upper level in He I compared to several other important He I flare emission lines (Figure 4.4). The peak of the profile in the UFB-3 spectrum is redshifted, as for the Mg II lines. The He I profile has comparable broadening to Mg II 2791.6, but it exhibits a larger 30% bisector velocity ( km s) than Mg II 2791.6 ( km s) and Fe II  ( km s). The width of the profile is much larger than the instrumental resolution (De Pontieu et al., 2014), which is readily apparent by comparing to the narrower Fe I and Fe II lines. Adding the instrumental and thermal widths in quadrature, we find that the width of He I 2829.91 equates to a thermal broadening corresponding to K, which is much larger than expected for a neutral helium line. The large FWHM (0.28 Å) of the He I 2829.91 Å line and the widths of the Mg II lines are likely due to a nonthermal broadening mechanism, such as a large optical depth and/or a spectrally unresolved redshifted component. The larger bisector velocity than Mg II 2791.6 Å and Fe II  could also be due to a larger optical depth in a chromospheric condensation: in this case, the redshifted component to the emergent intensity would be relatively brighter compared to the component of the intensity originating from the stationary flare layers below. There may also be a different brightness evolution in He I due to the variation in temperature sensitivity among different species formed in a condensation cooling from high to low temperature. Ostensibly, He I lines probe higher temperatures than Fe II, Mg II, and the NUV continuum radiation but detailed modeling is necessary to account for non-equilibrium ionization of helium (Allred et al., 2005; Golding et al., 2014; Allred et al., 2015; Simões et al., 2016). Modeling of this line will be included in Paper II.

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHrefPartial term diagram showing several important flare lines in the triplet system of He I, including the multiplet in the IRIS SJI 2832 bandpass. For reference, the ionization energy of He I is 198,310.6664 cm (Kandula et al., 2010).

5 Summary & Discussion

In this paper, we present a comprehensive analysis of umbral flare ribbon spectra in the NUV. These observations were obtained with a custom observing mode and high-time resolution, thus allowing us to characterize the flare continuum enhancements and constrain the time-evolution of heating in a single flaring location in the Sun. We present a new spectral ratio measurement, a continuum-to-line ratio, that can be obtained with NUV solar flare spectra with very limited wavelength coverage or large uncertainties in the intensity calibration. Because the Fe II  line and the NUV continuum radiation originate from (roughly) similar temperatures and exhibit (relatively) low optical depths in current RHD flare models (e.g., Kowalski et al., 2017a), the NUV continuum-to-line ratio (C2826/Fe II ) indicates the relative amounts of heating to K at high column mass (log /[g cm] ) compared at moderate-to-low column mass (log /[g cm] ). This interpretation will be described in detail in Paper II with new RHD models. Further, the ratio is generally independent of the spatial resolution assuming that the Fe II  line and the C2826 continuum radiation are produced at the same locations and do not exhibit different spatial structures below the instrumental resolution. Thus, the ratio facilitates a robust comparison of 1D models to observations that are spatially unresolved (as in dMe flares using the ratio of H line flux to the blue 4170 Å continuum; Kowalski et al., 2013; Silverberg et al., 2016; Kowalski et al., 2019). or have low spatial resolution, in which case the value of the continuum intensity is not accurately inferred due to spatial smoothing.

In the umbral flare brightening spectra of the 2014-Oct-25 X1 flare, C2826/Fe II , which is a factor of seven to eight larger than in the HXR impulsive phase spectra of the 2014-Mar-29 X1 flare: in the “BFP2” flare spectrum of the 2014-Mar-29 flare (Kowalski et al., 2017a), the C2826/Fe II  ratio is 1.0 but the value of C2826 is nearly a factor of 20 larger than at the peak of UFB-3. A 5F11 RHD simulation was used to model the red-wing asymmetry and bright continuum radiation in the 2014-Mar-29 flare in Kowalski et al. (2017a). This model predicts a continuum-to-line ratio of 1.1 (with a microturbulence parameter) and 1.8 (without a microturbulence parameter) and thus does not explain the large values in the UFBs in the 2014-Oct-25 flare. In the 2014-Mar-29 flare, the Fe II  line exhibits a very bright and broad red wing component (see Figure 3.2(b)). A lack of a bright, spectrally resolved red wing component in the 2014-Oct-25 flare contributes to a larger C2826/Fe II , but this is not nearly enough to explain the large values. These shortcomings of the 5F11 RHD model motivate using the Fermi/GBM data to model the nonthermal electron parameters in the 2014-Oct-25 X1 flare. The Fermi X-ray data (Section 2.2) suggest a similar energy flux as for the 2014-Mar-29 flare (Kleint et al., 2016) but a much steeper power law index ( compared to 4).

The observed range of red-wing properties in X-class flares suggest that the SJI 2832 contributions from emission lines may be significantly different in other flares that produce a spectrally resolved, bright red-wing asymmetry in the Fe II, Cr II, and Helium I lines, which all contribute to the SJI 2832 bandpass. In addition to providing detailed information on the contributions to the IRIS SJI 2832 during flares, we used the full readout data to characterize a He I flare line in the NUV. The helium emission will be used to provide constraints on the heating at higher temperatures than where Fe II and NUV continuum form, and it will help understand the origin of the non-thermal broadening in NUV flare lines. To our knowledge, this is the first report of this He I line in a solar flare.

6 Conclusions

We detect bona-fide NUV continuum radiation in IRIS flare spectra with a contrast of 1000% in a solar umbra. The spectra were obtained during the fourth hard X-ray peak at keV as a ribbon (that was part of a larger, three-ribbon X-class solar flare) developed into a sunspot. The main result of this analysis is the characterization of the ratio of NUV continuum radiation to the Fe II  line-integrated intensity, obtained from the IRIS NUV spectra. This continuum-to-line ratio is a new diagnostic of the relative heating rates at high and low column mass in the flare chromosphere because they are formed over similar temperatures (Appendix A) with moderately different optical depths (see Kowalski et al., 2017a, and Paper II). The ratios vary over values of in the umbral flare brightenings, attaining values of over the peak times. New RHD models are required to explain these large values, since previous high-beam flux models predict much lower values and highly redshifted Fe II emission line components (Kowalski et al., 2017a).

With the full spectral range of IRIS, we establish, for this flare, that the intensity in the narrow continuum window from Å (C2826) adequately represents the continuum level throughout the full NUV range of IRIS. We also identified the flare emission line landscape in the IRIS NUV and IRIS SJI 2832 images and found that the He I Å line becomes bright and broad in solar flares. This line is in the SJI 2832 wavelength range and requires full spectral readout of the IRIS NUV, as it is not included in standard IRIS line lists. Full spectral readout of the IRIS NUV also provides multiple lines of Fe II, Fe I, Cr II, and other species, yielding sets of profiles to constrain properties of the flaring solar atmosphere and disambiguate blends. The emission lines exhibit red profile asymmetries, but these are much less redshifted and broad than in another well-studied X-class flare with IRIS data.

The hard X-rays from Fermi/GBM combined with an algorithm (from Graham et al., in prep) to estimate the flare area provide starting-point inputs into RHD flare models of these intriguing umbral flare brightenings, in order to determine if the IRIS NUV flare spectra can be explained by electron beam heating alone and whether significant photospheric heating is required to produce large continuum-to-line ratios. In Paper II, we will use these two constraints (the continuum-to-line ratio and the Fe II line profile asymmetries) in addition to the constraints at higher temperatures from He I to determine the relative heating in the photosphere, in the chromosphere at high column mass, and in the chromosphere at low column mass in solar flares.


We thank an anonymous referee for comments that significantly improved the manuscript and the presentation of the results. We gratefully acknowledge the IRIS observation planners and conversations at Dr. Paola Testa’s workshop on solar microflares at the International Space Science Institute in Bern, Switzerland. AFK and EB acknowledge support from NASA Helio GI Grant NNX17AD62G. LF acknowledges support from the UK’s Science and Technology Facilities Council under grant ST/P000533/1. IRIS is a NASA small explorer mission developed and operated by LMSAL with mission operations executed at NASA Ames Research center and major contributions to downlink communications funded by ESA and the Norwegian Space Center.

References

  • Acampa et al. (1982) Acampa, E., Smaldone, L. A., Sambuco, A. M., & Falciani, R. 1982, A&AS, 47, 485
  • Aller (1963) Aller, L. H. 1963, Astrophysics. The atmospheres of the sun and stars
  • Allred et al. (2005) Allred, J. C., Hawley, S. L., Abbett, W. P., & Carlsson, M. 2005, ApJ, 630, 573
  • Allred et al. (2006) —. 2006, ApJ, 644, 484
  • Allred et al. (2015) Allred, J. C., Kowalski, A. F., & Carlsson, M. 2015, ApJ, 809, 104
  • Amari et al. (2018) Amari, T., Canou, A., Aly, J.-J., Delyon, F., & Alauzet, F. 2018, Nature, 554, 211
  • Aschwanden (1996) Aschwanden, M. J. 1996, in American Institute of Physics Conference Series, Vol. 374, American Institute of Physics Conference Series, ed. R. Ramaty, N. Mandzhavidze, & X.-M. Hua, 300–310
  • Aschwanden et al. (1996a) Aschwanden, M. J., Hudson, H., Kosugi, T., & Schwartz, R. A. 1996a, ApJ, 464, 985
  • Aschwanden et al. (1998a) Aschwanden, M. J., Kliem, B., Schwarz, U., et al. 1998a, ApJ, 505, 941
  • Aschwanden et al. (1996b) Aschwanden, M. J., Kosugi, T., Hudson, H. S., Wills, M. J., & Schwartz, R. A. 1996b, ApJ, 470, 1198
  • Aschwanden & Schwartz (1995) Aschwanden, M. J., & Schwartz, R. A. 1995, ApJ, 455, 699
  • Aschwanden et al. (1995) Aschwanden, M. J., Schwartz, R. A., & Alt, D. M. 1995, ApJ, 447, 923
  • Aschwanden et al. (1998b) Aschwanden, M. J., Schwartz, R. A., & Dennis, B. R. 1998b, ApJ, 502, 468
  • Aschwanden et al. (1996c) Aschwanden, M. J., Wills, M. J., Hudson, H. S., Kosugi, T., & Schwartz, R. A. 1996c, ApJ, 468, 398
  • Bamba et al. (2017) Bamba, Y., Inoue, S., Kusano, K., & Shiota, D. 2017, ApJ, 838, 134
  • Battaglia & Kontar (2011) Battaglia, M., & Kontar, E. P. 2011, A&A, 533, L2
  • Berry et al. (1972) Berry, H. G., Desesquelles, J., & Dufay, M. 1972, Phys. Rev. A, 6, 600
  • Boyer et al. (1985) Boyer, R., Sotirovsky, P., Machado, M. E., & Rust, D. M. 1985, SoPh, 98, 255
  • Cheng et al. (2010) Cheng, J. X., Ding, M. D., & Carlsson, M. 2010, ApJ, 711, 185
  • Courrier et al. (2018) Courrier, H., Kankelborg, C., De Pontieu, B., & Wülser, J.-P. 2018, SoPh, 293, 125
  • Cram & Woods (1982) Cram, L. E., & Woods, D. T. 1982, ApJ, 257, 269
  • De Pontieu et al. (2014) De Pontieu, B., Title, A. M., Lemen, J. R., et al. 2014, SoPh, 289, 2733
  • Donati-Falchi et al. (1985) Donati-Falchi, A., Falciani, R., & Smaldone, L. A. 1985, A&A, 152, 165
  • Donati-Falchi et al. (1984) Donati-Falchi, A., Smaldone, L. A., & Falciani, R. 1984, A&A, 131, 256
  • Drake (2006) Drake, G. 2006, High Precision Calculations for Helium, ed. G. Drake (New York, NY: Springer New York), 199–219
  • Dufton & McKeith (1980) Dufton, P. L., & McKeith, C. D. 1980, A&A, 81, 8
  • Fahey (1984) Fahey, R. P. 1984, ApJS, 55, 507
  • Fisher (1989) Fisher, G. H. 1989, ApJ, 346, 1019
  • Fletcher et al. (2007) Fletcher, L., Hannah, I. G., Hudson, H. S., & Metcalf, T. R. 2007, ApJ, 656, 1187
  • Fuhr & Wiese (2006) Fuhr, J. R., & Wiese, W. L. 2006, Journal of Physical and Chemical Reference Data, 35, 1669
  • Gan et al. (1992) Gan, W. Q., Rieger, E., Zhang, H. Q., & Fang, C. 1992, ApJ, 397, 694
  • Golding et al. (2014) Golding, T. P., Carlsson, M., & Leenaarts, J. 2014, ApJ, 784, 30
  • Halenka & Grabowski (1984) Halenka, J., & Grabowski, B. 1984, A&AS, 57, 43
  • Hawley & Fisher (1992) Hawley, S. L., & Fisher, G. H. 1992, ApJS, 78, 565
  • Hawley & Pettersen (1991) Hawley, S. L., & Pettersen, B. R. 1991, ApJ, 378, 725
  • Heinzel & Kleint (2014) Heinzel, P., & Kleint, L. 2014, ApJL, 794, L23
  • Heinzel et al. (2017) Heinzel, P., Kleint, L., Kašparová, J., & Krucker, S. 2017, ApJ, 847, 48
  • Hiei (1982) Hiei, E. 1982, SoPh, 80, 113
  • Hudson et al. (1992) Hudson, H. S., Acton, L. W., Hirayama, T., & Uchida, Y. 1992, PASJ, 44, L77
  • Hudson et al. (2006) Hudson, H. S., Wolfson, C. J., & Metcalf, T. R. 2006, SoPh, 234, 79
  • Hurford et al. (2006) Hurford, G. J., Krucker, S., Lin, R. P., et al. 2006, ApJL, 644, L93
  • Inoue et al. (2016) Inoue, S., Hayashi, K., & Kusano, K. 2016, ApJ, 818, 168
  • Kandula et al. (2010) Kandula, D. Z., Gohle, C., Pinkert, T. J., Ubachs, W., & Eikema, K. S. E. 2010, Physical Review Letters, 105, 063001
  • Kane et al. (1985) Kane, S. R., Love, J. J., Neidig, D. F., & Cliver, E. W. 1985, ApJL, 290, L45
  • Kerr et al. (2016) Kerr, G. S., Fletcher, L., Russell, A. J. B., & Allred, J. C. 2016, ApJ, 827, 101
  • Kerr et al. (2015) Kerr, G. S., Simões, P. J. A., Qiu, J., & Fletcher, L. 2015, A&A, 582, A50
  • Kleint et al. (2016) Kleint, L., Heinzel, P., Judge, P., & Krucker, S. 2016, ApJ, 816, 88
  • Kleint et al. (2017) Kleint, L., Heinzel, P., & Krucker, S. 2017, ApJ, 837, 160
  • Kontar et al. (2006) Kontar, E. P., MacKinnon, A. L., Schwartz, R. A., & Brown, J. C. 2006, A&A, 446, 1157
  • Kowalski & Allred (2018) Kowalski, A. F., & Allred, J. C. 2018, ApJ, 852, 61
  • Kowalski et al. (2017a) Kowalski, A. F., Allred, J. C., Daw, A., Cauzzi, G., & Carlsson, M. 2017a, ApJ, 836, 12
  • Kowalski et al. (2015a) Kowalski, A. F., Cauzzi, G., & Fletcher, L. 2015a, ApJ, 798, 107
  • Kowalski et al. (2015b) Kowalski, A. F., Hawley, S. L., Carlsson, M., et al. 2015b, SoPh, 290, 3487
  • Kowalski et al. (2012) Kowalski, A. F., Hawley, S. L., Holtzman, J. A., Wisniewski, J. P., & Hilton, E. J. 2012, SoPh, 277, 21
  • Kowalski et al. (2013) Kowalski, A. F., Hawley, S. L., Wisniewski, J. P., et al. 2013, ApJS, 207, 15
  • Kowalski et al. (2018) Kowalski, A. F., Mathioudakis, M., & Hawley, S. L. 2018, arXiv e-prints, arXiv:1810.07226
  • Kowalski et al. (2016) Kowalski, A. F., Mathioudakis, M., Hawley, S. L., et al. 2016, ApJ, 820, 95
  • Kowalski et al. (2017b) Kowalski, A. F., Allred, J. C., Uitenbroek, H., et al. 2017b, ApJ, 837, 125
  • Kowalski et al. (2019) Kowalski, A. F., Wisniewski, J. P., Hawley, S. L., et al. 2019, ApJ, 871, 167
  • Kramida et al. (2018) Kramida, A., Yu. Ralchenko, Reader, J., & and NIST ASD Team. 2018, NIST Atomic Spectra Database (ver. 5.6.1), [Online]. Available: https://physics.nist.gov/asd [2019, February 10]. National Institute of Standards and Technology, Gaithersburg, MD., ,
  • Krucker et al. (2011) Krucker, S., Hudson, H. S., Jeffrey, N. L. S., et al. 2011, ApJ, 739, 96
  • Krucker et al. (2015) Krucker, S., Saint-Hilaire, P., Hudson, H. S., et al. 2015, ApJ, 802, 19
  • Kurucz (2018) Kurucz, R. L. 2018, in Astronomical Society of the Pacific Conference Series, Vol. 515, Workshop on Astrophysical Opacities, 47
  • Kusano et al. (2012) Kusano, K., Bamba, Y., Yamamoto, T. T., et al. 2012, ApJ, 760, 31
  • Lemen et al. (2012) Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, SoPh, 275, 17
  • Li et al. (2018) Li, X., Zhang, J., Yang, S., & Hou, Y. 2018, PASJ, arXiv:1811.00281
  • Livshits et al. (1981) Livshits, M. A., Badalian, O. G., Kosovichev, A. G., & Katsova, M. M. 1981, SoPh, 73, 269
  • Maehara et al. (2012) Maehara, H., Shibayama, T., Notsu, S., et al. 2012, Nature, 485, 478
  • Martínez Oliveros et al. (2012) Martínez Oliveros, J.-C., Hudson, H. S., Hurford, G. J., et al. 2012, ApJL, 753, L26
  • Matthews et al. (2003) Matthews, S. A., van Driel-Gesztelyi, L., Hudson, H. S., & Nitta, N. V. 2003, A&A, 409, 1107
  • Meegan et al. (2009) Meegan, C., Lichti, G., Bhat, P. N., et al. 2009, ApJ, 702, 791
  • Mihalas (1978) Mihalas, D. 1978, Stellar atmospheres /2nd edition/
  • Milligan et al. (2014) Milligan, R. O., Kerr, G. S., Dennis, B. R., et al. 2014, ApJ, 793, 70
  • Murphy et al. (1997) Murphy, R. J., Share, G. H., Grove, J. E., et al. 1997, ApJ, 490, 883
  • Namekata et al. (2017) Namekata, K., Sakaue, T., Watanabe, K., et al. 2017, ApJ, 851, 91
  • Nave & Johansson (2013) Nave, G., & Johansson, S. 2013, The Astrophysical Journal Supplement Series, 204, 1
  • Nave et al. (1994) Nave, G., Johansson, S., Learner, R. C. M., Thorne, A. P., & Brault, J. W. 1994, ApJS, 94, 221
  • Neidig (1983) Neidig, D. F. 1983, SoPh, 85, 285
  • Neidig et al. (1993a) Neidig, D. F., Kiplinger, A. L., Cohl, H. S., & Wiborg, P. H. 1993a, ApJ, 406, 306
  • Neidig et al. (1993b) Neidig, D. F., Wiborg, P. H., & Gilliam, L. B. 1993b, SoPh, 144, 169
  • Nishizuka et al. (2009) Nishizuka, N., Asai, A., Takasaki, H., Kurokawa, H., & Shibata, K. 2009, ApJL, 694, L74
  • Osten et al. (2016) Osten, R. A., Kowalski, A., Drake, S. A., et al. 2016, ApJ, 832, 174
  • Penn et al. (2016) Penn, M., Krucker, S., Hudson, H., et al. 2016, ApJL, 819, L30
  • Procházka et al. (2017) Procházka, O., Milligan, R. O., Allred, J. C., et al. 2017, ApJ, 837, 46
  • Procházka et al. (2018) Procházka, O., Reid, A., Milligan, R. O., et al. 2018, ApJ, 862, 76
  • Rathore & Carlsson (2015) Rathore, B., & Carlsson, M. 2015, ApJ, 811, 80
  • Reep & Russell (2016) Reep, J. W., & Russell, A. J. B. 2016, ApJL, 818, L20
  • Reep et al. (2018) Reep, J. W., Russell, A. J. B., Tarr, L. A., & Leake, J. E. 2018, ApJ, 853, 101
  • Rubio da Costa et al. (2016) Rubio da Costa, F., Kleint, L., Petrosian, V., Liu, W., & Allred, J. C. 2016, ApJ, 827, 38
  • Russell & Fletcher (2013) Russell, A. J. B., & Fletcher, L. 2013, ApJ, 765, 81
  • Rutten (2003) Rutten, R. J. 2003, Radiative Transfer in Stellar Atmospheres
  • Sadykov et al. (2016) Sadykov, V. M., Kosovichev, A. G., Sharykin, I. N., Zimovets, I. V., & Vargas Dominguez, S. 2016, ApJ, 828, 4
  • Sansonetti & Nave (2014) Sansonetti, C. J., & Nave, G. 2014, The Astrophysical Journal Supplement Series, 213, 28
  • Seaton (1960) Seaton, M. J. 1960, Reports on Progress in Physics, 23, 313
  • Sharykin & Kosovichev (2014) Sharykin, I. N., & Kosovichev, A. G. 2014, ApJL, 788, L18
  • Silverberg et al. (2016) Silverberg, S. M., Kowalski, A. F., Davenport, J. R. A., et al. 2016, ApJ, 829, 129
  • Simões et al. (2016) Simões, P. J. A., Fletcher, L., Labrosse, N., & Kerr, G. S. 2016, in Astronomical Society of the Pacific Conference Series, Vol. 504, Coimbra Solar Physics Meeting: Ground-based Solar Observations in the Space Instrumentation Era, ed. I. Dorotovic, C. E. Fischer, & M. Temmer, 197
  • Sirianni et al. (2005) Sirianni, M., Jee, M. J., Benítez, N., et al. 2005, PASP, 117, 1049
  • Tayal & Zatsarinny (2018) Tayal, S. S., & Zatsarinny, O. 2018, PhRvA, 98, 012706
  • Thalmann et al. (2015) Thalmann, J. K., Su, Y., Temmer, M., & Veronig, A. M. 2015, ApJL, 801, L23
  • Tian et al. (2015) Tian, H., Young, P. R., Reeves, K. K., et al. 2015, ApJ, 811, 139
  • Uitenbroek (2001) Uitenbroek, H. 2001, ApJ, 557, 389
  • Vernazza et al. (1981) Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635
  • Veronig & Brown (2004) Veronig, A. M., & Brown, J. C. 2004, ApJL, 603, L117
  • Warmuth & Mann (2016) Warmuth, A., & Mann, G. 2016, A&A, 588, A116
  • Watanabe et al. (2017) Watanabe, K., Kitagawa, J., & Masuda, S. 2017, ApJ, 850, 204
  • Wülser et al. (2018) Wülser, J.-P., Jaeggli, S., De Pontieu, B., et al. 2018, SoPh, 293, 149
  • Zharkova & Zharkov (2007) Zharkova, V. V., & Zharkov, S. I. 2007, ApJ, 664, 573
\onecolumngrid

APPENDIX

Appendix A The formation of Fe Ii  and continuum radiation at Å in LTE

In this appendix, we present the temperature and density sensitivities for hydrogen Balmer recombination radiation, the line intensity in Fe II , and their ratios. The similar temperature sensitivities (assuming LTE) justify using the Fe II  to constrain the velocity field in the layers where the NUV continuum radiation is formed, and the ratios from optically thin uniform slabs suggest that a more sophisticated approach to the modeling (e.g., with RHD models) is necessary.

We calculate the LTE, Å hydrogen bound-free continuum emissivity and the Fe II  line-integrated emissivity for a range of gas densities () and temperatures (). We use the standard equations from Rutten (2003), Aller (1963), and Mihalas (1978) (Eq. 7-4) for the spontaneous thermal, LTE line and continuum emissivities. For hydrogen () recombination to a given principle quantum number , the LTE continuum emissivity, , follows from (where is the opacity corrected for stimulated emission) and varies as (see also Eq. 3 of Kowalski et al., 2015b). Including the dependence, the continuum emissivity reduces to the following:

(A1)

in units of [erg s cm sr Å], and Å. The values of the hydrogen bound-free gaunt factor () are taken from Seaton (1960). We calculate the recombination to ; recombination to and free-free emissivity contribute a moderate amount (20%) to the continuum emissivity only at the highest temperatures that we consider ( K). At the lowest temperatures, H recombination may contribute at the highest densities (but is not included here).

For the spontaneous thermal line emissivity, we use Eq. 2.69 of Rutten (2003) with calculated from LTE using the partition function of Halenka & Grabowski (1984) and the ionization potential lowered by 0.1 eV for Fe I and Fe II and 0.25 eV for Fe III. We integrate over the wavelength of the line (thus, giving Eq. 2.70 of Rutten (2003)).

\H@refstepcounter

figure \hyper@makecurrentfigure

Figure 0. \Hy@raisedlink\hyper@@anchor\@currentHref The C2826/Fe II  emissivity (also, optically thin LTE emergent intensity) ratios as a function of and . In RHD model atmospheres, the ratios of the emergent intensity differ from these values due to non-uniform and as a function of height and time, possible non-equilibrium effects (e.g., non-LTE ionization of Fe II ), and differences in the optical depth between the emergent C2826 continuum radiation and the Fe II  emission line. These slab models do not reasonably explain the observed range of ratios (; Figure 3.2) in the umbral flare brightenings.

The emissivity grids are shown in Figure A. For the same density, the bound-free emissivity peaks at slightly higher ( K) temperature and has a relatively brighter tail at higher temperatures. However, the peaks in the continuum emissivity curves are rather broad: 90% of the maximum for each value falls within the peak of the Fe II  emissivity, and at 50%, the emissivity of Fe II  spans the temperatures of K.

In Figure A, we show the ratios of the emissivities. This figure also gives the emergent intensity ratio from isothermal, isobaric optically thin, slab models in LTE. The ratios are typically much less than unity except at very high temperatures at 16,000 K. In RHD models (Kowalski et al., 2017a), the Fe II  line and NUV continuum radiation are formed at significantly lower temperatures than K because most of the flare chromospheric mass is at lower temperatures. For example, about 95% of the emergent NUV continuum intensity originates from K in the evolved chromospheric condensation in the 5F11 electron beam-heated model atmosphere in Kowalski et al. (2017a). Therefore, optically thin slab predictions over a reasonable temperature range do not explain the large continuum-to-line ratios in the umbral flare brightenings. As we will show in Paper II, large optical depths from significant heating at high column mass can produce consistent continuum-to-line ratios.

Appendix B Flare Line Identifications in the IRIS NUV

A number of lines in the UFB-3 peak excess intensity spectrum are identified in Table B1 and shown in Figure B. The line-integrated, continuum-subtracted line intensities were estimated by fitting the local continuum using a first-order polynomial and subtracting across the emission line(s). The 30% bisectors are calculated from the intensity after subtracting a local continuum (without subtracting a pre-flare spectrum) to indicate the measured wavelengths. In the last two rows, the excess continuum intensity, C2826, and the IRIS NUV bandpass-integrated brightness are also given. Rest wavelengths for Fe II, Cr II, and Fe I are adopted primarily from Nave & Johansson (2013), Sansonetti & Nave (2014), and Nave et al. (1994), respectively, and additional data on wavelengths and A-values were obtained from the NIST (Kramida et al., 2018), Kurucz (2018), and R. L. Kelly (https://www.cfa.harvard.edu/ampcgi/kelly.pl) databases and references therein. Most observed species have more than one line in the IRIS NUV, and many have lines from the same multiplet. Consistency for each species was established by comparing line positions and profiles to those from the same species, and by comparing the observed intensity ratios for each species to those from Boltzmann-distributed populations for a range of temperatures of formation ( 7,000 to 22,000 K for Fe II; see Fig. A), that is, optically thin LTE calculations, with opacity considered in the comparisons as discussed below.

Of the 31 Fe II lines listed in the UFB-3 peak spectrum, the majority are reasonably strong, isolated features with no evidence of significant blends, providing many ratios between Fe II line intensities that can be used to constrain model atmospheres. As discussed in Section 4.3, compared to the ratios in the flare, the LTE intensity ratios are systematically over-predicted for the bright Fe II lines, and we speculate that optical depth in these brighter Fe II lines prevents a relatively larger amount of emission from escaping from the stationary flare layers below the chromospheric condensation. This speculation is supported by the fact that for the two brightest Fe II lines, 2784.512 and 2832.394, the rest component is less pronounced relative to the red wing than it is for the weaker Fe II lines (Figure 4.3). Recent electron impact excitation calculations by Tayal & Zatsarinny (2018) are now available for 25 of the 31 observed Fe II lines, which can be used to facilitate accurate modeling of these lines to test this hypothesis.

Two observed Fe II lines are not in fact listed in the most recent work by Nave & Johansson (2013): Fe II 2820.173 and 2828.724, although these most recent energy levels are certainly consistent with the observed wavelengths. In the case of Fe II 2820.173, the lack of an identification in the most recent laboratory work is possibly due to blending with the Fe II 2820.166 line, which is readily excited in laboratory spectra, but comes from a higher level of excitation than would be expected for solar plasmas. As with the other flare line identifications for UFB-3, evidence is provided by the similarity of the Fe II 2820.173 and 2828.724 profiles to other Fe II lines, and of intensity levels to preliminary LTE calculations. In this case, both these lines and Fe II 2824.159 are from the same multiplet as well, so we are fairly confident of these identifications. There is an Fe I line at 2828.724 Å to be aware of, but it is expected to be an order of magnitude weaker than the other nearby observed lines of Fe I. The -values reported by Kurucz (2018) and Fuhr & Wiese (2006) are in reasonable agreement for all of the IRIS NUV Fe II lines that are in both databases, although with differences in some cases of up to a factor of two. -values for six of the lines, including 2828.724, are reported only by Kurucz, so in Table B1 the -value for those lines are from Kurucz, while the balance are from Fuhr & Wiese.

While the most extensive set of profiles is provided by Fe II, there are quite a few Cr II lines with profiles having red wings similar to those of Fe II, which provide another set of lines to test the hypothesis above by comparison to future modeling. Multiple lines of Fe I, Ti II, Ni I, Ni II, and Mn I are observed as well. Note that, like the He I 2829.91 Å line discussed in section 4.4 and Mg II h & k, the Al II 2817.014 Å line peaks to the red in UFB-3. There is a line of Ti II in the red wing of Al II, but the intensities of the other observed lines of Ti II do not indicate that Ti II is a significant contribution to this Al II profile. Finally, note that O V is listed not as a definitive identification, but as a possible contribution to the feature around 2788 Å. If present, the O V 2787.814 Å line appears blended with other lines, and it would be accompanied by the weaker line from the same multiplet at 2790.669 Å, for which there is a blend with a weak Fe I line.

Interestingly, there is a line from a doubly-excited state of He I that has been observed at 2819.2 Å using beam-foil spectroscopy (Berry et al., 1972), which could conceivably contribute to the observed feature at that wavelength in the UFB-3 flare spectrum. In solar plasmas, such doubly-excited states would be populated through electron capture by non-thermal (accelerated) helium ions upon hitting a thick target, but not in thermal plasmas. Unfortunately, while there are stronger such lines outside the IRIS NUV band pass, this one has a blend with Cr II 2819.18, which is the dominant contribution to the observed feature. It is nevertheless possible upon further evaluation of all the Cr II profiles that some evidence of a contribution from doubly-excited helium in that profile would provide a signature of accelerated helium ion beams.

In summary, a number of IRIS NUV lines are identified in this appendix that will be useful to model and compare to the data in order to constrain the physical structure of the flaring atmosphere observed by the IRIS spectrograph.

Intensity Species Bis Multiplet A-value
[] [Å] [km s] [cm] [cm] [s
112 Fe II 2784.512 7.5 H bH H zG 26170.181 62083.118 1.06e08
4 (Fe II) 2785.847   aF D zD 3117.488 39013.216 7.57e03
19 Fe II 2786.014 4.4 D zF D eD 41968.070 77861.650 1.53e08
8 Cr II 2786.514 6.7 G bG G F 33618.936 69506.065 2.09e08
24 (O V) 2787.814 8.3 S P 546972.700 582843.100 1.41e08
6 Cr II 2788.440 6.0 P bP P yP 30864.433 66726.782 1.50e08
19 Fe I 2788.927 0.7 F aF HP yG 6928.268 42784.349 6.30e07
2 (Fe I) 2790.624   G aG DP tF 21715.731 57550.006 2.36e07
  (O V) 2790.669   S P 546972.700 582806.400 1.43e08
1175 Mg II 2791.600 8.2 P D 35669.310 71491.063 4.08e08
8 Cr II 2792.978 3.6 G bG G F 33694.142 69498.214 2.30e08
53 Fe II 2794.711 3.2 G aG H zI 25805.327 61587.205 1.30e07
1 Mn I 2795.641   aS SP yP 0.000 35769.970 3.62e08
7036 Mg II 2796.352 9.9 S P 0.000 35760.880 2.68e08
2 (Fe II) 2797.868   aF D zD 3117.488 38858.970 2.02e03
  Mg II 2798.754   P D 35760.880 71491.063 8.09e07
1683 Mg II 2798.823 6.0 P D 35760.880 71490.190 4.81e08
1 Ni I 2799.474   D D FP D 879.813 36600.805 5.77e06
20 Fe II 2800.120 3.8 H bH F yF 26352.767 62065.528 1.55e07
3 (Fe II) 2800.548 5.4 G aG H zI 25805.327 61512.630 5.00e05
10 Cr II 2801.591 7.3 G bG G yH 33694.142 69388.151 2.20e08
7 Mn I 2801.907 -0.4 aS SP yP 0.000 35689.980 3.69e08
6169 Mg II 2803.531 10.8 S P 0.000 35669.310 2.62e08
23 Fe II 2804.846 -0.4 F aF F xD 27620.403 63272.981 1.60e06
7 Fe I 2805.347 -0.5 F aF HP yG 7376.764 43022.982 1.05e07
1 (Fe II) 2805.826   F yG F eF 65109.691 100749.825 1.39e08
1 (Fe II) 2806.145   D bD F yD 31364.455 67000.530 2.50e06
1 Ni II 2806.491   G) G D) F 32499.530 68131.210 1.30e07
5 Fe II 2806.614 -0.9 F aF F xD 27314.918 62945.045 3.20e06
9 Fe I 2807.811 0.0 F aF HP zH 7376.764 42991.694 1.15e07
4 Fe II 2810.610 5.9 D zP D eD 42658.244 78237.709 3.10e07
3 Ti II 2811.061 2.6 F zG F eH 29734.540 65308.300 5.09e08
1 Ti II 2811.133   aP P yP 9975.920 45548.760 1.24e07
8 Fe II 2812.097 0.8 G aG H zH 25428.789 60989.444 1.20e06
9 Cr II 2812.828 7.9 G bG G yH 33618.936 69170.353 2.05e08
6 Fe II 2813.322 4.8 P bP P yP 25787.582 61332.753 2.90e06
11 Fe I 2814.115 0.1 F aF HP yG 7376.764 42911.914 3.42e07
15 Fe II 2814.445 3.7 G aG H zI 25981.645 61512.630 3.40e06
22 Al II 2817.014 10.9 P S 59852.020 95350.600 3.93e08
1 Cr II 2817.670   P bP P yP 30864.433 66354.757 1.04e08
1 (Fe II) 2817.916   D zP D eD 43238.607 78725.822 3.40e07
7 Cr II 2819.184 7.4 G bG G yH 33521.090 68992.347 2.21e08
  Fe II 2820.166   D xD G D 74606.864 110065.766 8.25e05
8 Fe II 2820.173 2.7 G aG H zH 25428.784 60887.598 1.00e06
2 Ni I 2822.120 0.6 D D FP F 204.786 35639.148 4.87e06
6 Cr II 2822.842 7.2 G bG G yH 33417.981 68843.273 2.29e08
28 Cr II 2823.199 7.6 H aH H zI 30391.831 65812.649 2.28e08
  Fe I 2824.107   F aF HP yG 7728.059 43137.484 1.51e07
30 Fe II 2824.159 0.2 G aG H zH 25428.789 60837.560 2.10e06
3 Ni II 2826.062 -0.3 P) P F) F 23108.280 58493.210 2.65e06
14 Fe I 2826.387 1.3 F aF HP zH 7728.059 43108.914 1.32e07
18 Fe II 2826.579 0.5 G aG H zG 25428.789 60807.239 1.40e06
16 Fe II 2826.859 1.6 F aF P yD 27314.918 62689.874 4.50e06
25 Fe II 2828.260 2.3 H bH H zI 26170.181 61527.610 2.40e06
  Fe I 2828.724   aD F zG 415.933 35767.562 1.48e05
7 Fe II 2828.734 5.2 G aG H zH 25805.328 61156.835 1.51e06
36 Fe II 2829.459 7.9 H bH H zI 26170.181 61512.630 6.90e06
  Fe II 2829.510   F aF P yD 27620.403 62962.215 9.00e06
  Fe I 2829.640   F aF HP zH 7985.784 43325.961 1.87e06
  He I 2829.911   S P 159855.974 195192.777 1.94e06
81 He I 2829.913 10.9 S P 159855.974 195192.746 1.94e06
  He I 2829.914   S P 159855.974 195192.743 1.94e06
27 Cr II 2831.299 9.3 H aH H zI 30298.468 65617.946 2.54e08
  Cr II 2831.449   H aH F zG 30391.831 65709.442 3.71e07
3 Fe II 2831.754   G bG G xG 30764.474 66078.272 7.00e05
  Fe II 2831.917   P bP F yF 26932.735 62244.515 5.43e05
71 Fe II 2832.394 7.3 P bP P zD 25787.582 61093.406 7.60e07
6 Ti II 2833.015 2.6 F a