arxiv: v1 [astro-ph.sr] 16 Dec 2015

Transcription

1 Astronomy & Astrophysics manuscript no. aia iris c ESO 215 December 17, 215 Observations of solar flares with IRIS and SDO D. Li 1,2,3, D.E. Innes 2, and Z. J. Ning 1 1 Key Laboratory of Dark Matter and Space Science, Purple Mountain Observatory, CAS, 28 Nanjing, China ( lidong@pmo.ac.cn) 2 Max-Planck Institute for Solar System Research, 3777 Göttingen, Germany 3 University of Chinese Academy of Sciences, 49 Beijing, China arxiv: v1 [astro-ph.sr] 16 Dec 215 Received; accepted ABSTRACT Context. Flare kernels brighten simultaneously in all Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) channels making it difficult to determine their temperature structure. The Interface Region Imaging Spectrograph (IRIS) is able to spectrally resolve Fe xxi emission from cold chromospheric brightenings, so can be used to infer the amount of Fe xxi emission in the 131 Å AIA channel. Aims. We use observations of two small solar flares seen by IRIS and SDO to compare the emission measures (EMs) deduced from the IRIS Fexxi line and the AIA 131 Å channel to determine the fraction of Fexxi emission in flare kernels in the 131 Å channel of AIA. Methods. Cotemporal and cospatial pseudo-raster AIA images are compared with the IRIS results. We use multi-gaussian line fitting to separate the blending chromospheric emission so as to derive Fexxi intensities and Doppler shifts in IRIS spectra. Results. We define loop and kernel regions based on the brightness of the 131 and 16 Å intensities. In the loop regions the Fexxi EMs are typically 8% of the 131 Å ones, and range from 67% to 92%. Much of the scatter is due to small misalignments but the largest site with low Fexxi contributions was probably affected by a recent injection of cool plasma into the loop. In flare kernels the contribution of Fexxi increases from less than % at the low intensity 131 Å sites to 4-8% in the brighter kernels. Here the Fexxi is superimposed on bright chromospheric emission and the Fe xxi line shows blue shifts, sometimes extending up to the edge of the spectral window, 2 km s 1. Conclusions. The AIA 131 Å emission in flare loops is due to Fe xxi emission with a -2% contribution from continuum, Fe xxiii, and cooler background plasma emission. In bright flare kernels up to 52% of the 131 Å is from cooler plasma. The wide range seen in the kernels is caused by significant structure in the kernels which is seen as sharp gradients in Fe xxi EM at sites of molecular and transition region emission. Key words. Sun: flares Sun: UV radiation Line: profiles Techniques: spectroscopic 1. Introduction The flare impulsive phase is characterized by a sudden increase in chromospheric and hard X-ray emission (Fletcher & Hudson 21; Fletcher et al. 211). The emission is concentrated in small bright kernels and along ribbons, coinciding with magnetic field concentrations, where rapid chromospheric heating drives hot plasma upward into the corona (Mason et al. 1986; Teriaca et al. 26; Benz 28; Ning & Cao 211). The kernels are thought to be the chromospheric signature of magnetic reconnection in the corona (Qiu et al. 22; Fletcher et al. 24; Gan et al. 28). Spectroscopic observations have revealed hightemperature, high-velocity blue and cooler red-shift emission (Antonucci et al. 1982; Milligan et al. 26; Teriaca et al. 26; Young et al. 213) compatible with models of chromospheric evaporation. Observations of the kernels show simultaneous brightening in all AIA extreme ultraviolet (EUV) channels (Brosius & Holman 212; Fletcher et al. 213; Young et al. 213). The analysis by Brosius & Holman (212) that compared AIA images with CDS spectra, showed that in a small GOES B4.8 flare a significant fraction of the hot channel, 94 Å and 131 Å, emission could be attributed to brightening of transition region and lower coronal lines. On the other hand Fletcher et al. (213) attributed all the 131 Å brightening in an M1. flare to plasma with temperatures greater than MK. In an M1.1 flare, Young et al. (213) s analysis of EIS spectra found almost equal emission measures (EMs) across the observed temperature range from.1 to MK but did not discuss the contributions to the AIA channel images. The main contribution to the 131 Å channel in flares is the Fe xxi Å line (O Dwyer et al. 2; Milligan & McElroy 213). The same ion, Fexxi produces the strong flare line at Å in IRIS spectra. Therefore by comparing simultaneous IRIS and AIA observations it may be possible to determine the contribution of Fe xxi to the 131 channel and hence resolve the the question of the cool plasma emission contribution to the 131 channel. Using the same ion, has the advantage that the line ratios are independent of ionization and only weakly dependent on temperature. The forbidden line of Fexxi at Å has been used in several spectroscopic studies to investigate hot plasma flow during flares (Doschek et al. 1975; Cheng et al. 1979; Mason et al. 1986; Feldman et al. 2; Innes 21; Kliem et al. 22; Wang et al. 23). More recently Young et al. (215) described 1

2 Fig. 1. GOES Å flux during the observing period. The two flares analysed are labelled. Red lines are drawn at the times of the peak flux, while the dashed lines show the start and end times of the of IRIS raster 6 (between 2:2 UT and 2:35 UT) and raster 9 (from 21:43 UT to 22:16 UT). high spatial and spectral resolution IRIS observations of Fe xxi from hot flare kernels and loops, with temperatures about MK from an X1 class flare. Their results support the chromospheric evaporation model. The Mason et al. (1986) observations indicated that the Fexxi line is often blended with Ci. The higher resolution spectra analysed by Young et al. (215) show that the Fe xxi kernel emission may, in addition, be blended with other chromospheric and possibly molecular lines. Molecular hydrogen was identified at the footpoints of X-ray loops with SUMER spectra (Innes 28). To obtain estimates of the Fexxi EM, the strength of these chromospheric lines needs to be taken into account. For the analysis in this paper, we have investigated a number of flare kernel spectra and developed an algorithm for simultaneously fitting the Fexxi, blending lines, and continuum. The EM obtained from the IRIS 1354 Å line was compared with that derived from simultaneous SDO/AIA 131 Å images to determine the contribution of cooler line emission to the 131 Å channel at the site of flare kernels. We find that the Fexxi observed by IRIS can account for about 8% of the SDO/AIA 131 Å emission in flare loops and 4-8% of the emission in flare kernels. Assuming that an additional 2% of the Fexxi is continuum (Milligan & McElroy 213), we find that up to 52% of the 131 Å channel emission is from cooler plasma in the flare kernels. 2. Observations The observed active region, AR 11875, produced four C-class flares during the period 16:39 UT on 24 October 213 and 2:46 UT on 25 October 213. We obtained high quality cotemporal and cospatial IRIS and AIA observations for flares at 2: and 22:5 UT (Fig. 1). During the early part of the sequence SDO was off-pointing and during the 21:9 UT flare, the IRIS spectra were badly affected by a large number of particle hits. Fig. 2 shows the AIA 131 Å and corresponding slit-jaw (SJ) 14 Å images of the two flares analysed. The GOES fluxes of the two flares are shown in Fig. 1. Flare 1 peaked at about 2: UT and ended at about 2:22 UT, while the flare 2 peaked at about 22: UT and stopped at about 22:15 UT. Fig. 2. Left: AIA 131 Å images of the active region at the times of the two flares. The red boxes are the positions of the IRIS SJ images. Right: IRIS SJ 14 Å images of the two flares analysed. The position of the spectrometer slit at the start of the flare is seen as a white dashed horizontal line. The purple pluses mark the re-spike pixels in the AIA 131 Å images with the code aia respike AIA observations In this analysis, AIA images (Lemen et al. 212) from the 131 and 16 Å channels are used to compare with the IRIS data. The AIA level 1. data were downloaded and processed to level 1.5 using the standard solarsoft (SSW) routines. Then using the code drot map.pro in SSW, the active region with a FOV of are selected from the AIA full disk images (see., Fig 2). The images from different AIA filters have alignment uncertainties of about 1 2 pixels (e.g., Young et al. 213). Therefore the AIA 131 Å images are slightly shifted (i.e.,.5 2 pixels) with respect to the 16 Å images because we obtained the best alignment of the 16 Å with the IRIS 14 Å and the 131 Å with the IRIS Fexxi. The downloaded AIA data are already de-spiked which often removes flare kernel emission and these can be put back with the aia respike routine (see., Young et al. 213). In our study, we have compared the re-spike data and the level 1 data at 131 Å as shown with the purple pluses in Fig. 2. This shows that the re-spiked pixels are far away from the IRIS slit position, and they are not within the regions of our study. This was the case for all other 131 Å images. Also during big flares the AIA 131 Å channel is often saturated. The two flares in this paper are quite small and only very few pixels in the 131 Å channel were saturated. Fortunately, those saturated pixels are far away from the IRIS raster slit. In other words, neither the de-spiking nor the saturated pixels are a problem in our study. 2

3 Fig. 3. Coalignment of the SDO and IRIS images. Left: the AIA image at 16 Å. Right: the SJ image at 14 Å. The white dashed line is the IRIS slit position, and the blue contours are the SJ intensity at the level of 8 DN IRIS observations IRIS is a NASA Small Explorer Mission launched in June 213, and its main science is an investigation of the dynamics of the Sun s chromosphere and transition region (Mclntosh et al. 213). As described in detail by De Pontieu et al. (214), IRIS is a high resolution spectrograph and slit-jaw (SJ) camera which obtains spatial resolution of arcsec, and spectral resolution of 26 må (or 52 må in second order). There are three IRIS wavelength bands: (i) Å which includes the strong Cii doublet, Fexxi, and a number of strong Ci lines, (ii) Å which includes the Si iv doublet, and density sensitive Oiv lines, (iii) Å which includes Mgii h and k lines. Here we analyse lines in the Å range. IRIS also obtains SJ images centered at 133, 14, 2796, and 2832 Å. The observations described here were designed with the aim of capturing the structure and spectra of flares by making rapid rasters across a flaring active region for as long a period as telemetry constraints allowed. The IRIS rasters were obtained by taking steps across a region 174 x 63 with the roll angle 9 (i.e. the slit was oriented E-W). Simultaneous SJ images at 14 Å with a field-of-view 174 x166 and cadence 32 s were also obtained and used mainly for coalignment with 16 Å AIA images. Because continuum emission from the temperature minimum is dominant in many of the bright features seen in both sets of images, coalignment is possible to within the AIA pixel size,.6 (Fig. 3). The step cadence was 31.6 s and the exposure time 3 s, thus the raster cadence was 33 m 42 s. 18 rasters were obtained in the roughly hour period. The pixel size along the slit was at the highest resolution.167 but to save telemetry four times spectral binning and a restricted number of spectral windows were obtained. In this case we used the flare list of lines which consisted of the Cii, 1343, Fexii, Oi, Siiv FUV windows. Here we only discuss spectra from windows in the short wavelength FUV band. The spectral resolution was 4*12.72 må/pixel, equivalent to.88 km s 1 /pixel. IRIS level 2 data were downloaded. These data have been calibrated and corrected for image distortions. Several of the images were badly affected by particle hits and hot or dead pixels. These were removed using a despiking routine that detected persistent hot/dead pixels, and sudden changes in intensity since these can seriously distort the line fitting results. Care is taken to ensure that only isolated bright pixels or small Fig. 4. IRIS raster images before (left) and after despiking (right). The middle panel is an image after despiking with the bad pixels marked with green plus ( + ) signs. pixel groups are removed, and to leave intensity changes due to sudden line broadening. Fig. 4 shows the effect of despiking. In addition, small spectral shifts caused by thermal drifts and the spacecraft orbital velocity were corrected with the routine iris orbitvar corr 12.pro (Mclntosh et al. 213; Tian et al. 214; Cheng et al. 215) in the SSW package Flare ribbon spectra The ribbons are characterized by many narrow, bright emission lines. The Fe xxi emission blends with both known and unknown lines from neutral and singly ionized species, as well as molecular fluorescence lines. To determine the Fexxi intensity we need to extract the blending chromospheric emission. We have identified the main blending lines by looking at kernel spectra from other IRIS datasets at higher spectral resolution, as well as many spectra in this dataset to find out which other observed lines behave similarly to the chromospheric lines in the Fe xxi window. An example of a flare footpoint spectrum taken with this sequence is shown in Fig. 5. For comparison, the chromospheric lines seen during a C1.5 flare from AR11861 (-133, -285 ) at :46 UT on 213 October 12, taken with the full spectra resolution, are shown in Fig. 6. The main chromospheric lines in the O i window, apart from the Ci Å, are the Feii lines at , , and Å, the Siii lines at and Å, and the unidentified lines at and Å. Actually, other papers (e.g., Graham & Cauzzi 215; Li et al. 215; Tian et al. 215) have attempted multi-gaussian fits to the Fexxi line from IRIS spectral observations. In this paper, to obtain the Fexxi intensities, we have fixed or constrained the positions and widths of these lines and set their intensities to have a specified ratio to well-resolved lines from similar species. In total we fit 17 Gaussian lines superimposed on a linear background fitted across the wavelength region (i.e., Å) covered by the four IRIS spectral windows. Table 1 lists the lines used in the fitting procedure. Lines with fixed positions are indicated by a superscript 1. The two Siii lines, indicated with a superscript 2, are constrained to match the shift of the unblended Siii at Å. The widths are fixed or constrained as given in column 4. For example, the Siii line has a maximum width of 26 må; while the line width of Feii is fixed 3

4 H-H H-H,S I C II Ni II C II H-H Fe XII Si II Fe XXI C I Fe II H-H Fig. 5. Flare footpoint spectra taken on 213 October 24. The overplotted line spectra are the average across the small bright region marked by a red line on the lefthand side of each image. The main lines are labelled and indicated by light blue vertical ticks just above the line spectra H-H H-H,S I C II Ni II C II Fe XII Si II H-H Fe XXI C I Fe II H-H Fig. 6. Flare footpoint spectra taken on 213 October 12 at full resolution. The overplotted line spectra are the average across the small bright region marked by a red line on the lefthand side of each image, and, in red, the same spectra binned by a factor 4 which is the resolution used on 213 October 24. at 41 må. The peak intensities of the blending lines are forced to be in a fixed ratio (column 6) with the lines that they are tied to (column 5). For example, the peak intensity of the blending line Siii is tied to the emission line Siii in the Fe xii window, and the intensity ratio is.49. The positions and intensity ratios of the lines were determined by fitting Gaussians to all the narrow (identified and unidentified) lines in 17 spectra in the dataset with bright molecular and chromospheric lines but Fig. 7. The spectra (black) and the multi-gaussian fits (red). The blue dashed lines mark the positions of cold lines, the blue solid line marks the Ci line, while the purple line marks the Fexxi line. no Fexxi emission. For each line, the median position from the 17 spectral fits was taken as the line s fixed position. Then for each narrow line blending with the Fexxi we calculated its intensity ratio with respect to all unblended narrow lines. We used the variation in the intensity ratios to identify which line each blended lines should be tied to, and the median of the ratio was used to determine the strength of the blended line from the intensity of its tied line. The lines used in the final fitting procedure to obtain the Fexxi intensities are marked and the regions of background are indicated by the red-dashed lines in Fig. 7. At most positions the fits appear to be successful. 3. The Emission Measures of Fexxi The observed intensity is given by I=E f G(T, n e )n 2 edz (1) where n e is the electron density, T is electron temperature, z is the line-of-sight coordinate, G(T, n e ) is the contribution function of the line including the atomic abundance, and E f is the instrument effective area or response function. The units of I, in the equation above, depend on the units of the E f. As discussed below, IRIS is a spectrometer and the intensities are given in DN s 1. AIA is an imager and the intensity unit is DN px 1 s 1. To convert intensity to emission measure, we need to divide by the respective response functions in the appropriate units. Since the contribution functions of the IRIS and AIA Fexxi lines are sharply peaked in temperature and are formed over the same heights, the contribution functions can be removed from the integral. The emission measure in the Fexxi plasma can then be calculated from EM= n 2 e dz= I/(E fg(t peak )) cm 5. (2) 4

5 Table 1. The 17 emission lines from the four IRIS spectral windows used in the multi-gaussian fit. Li et al.: Observations of solar flares with IRIS and SDO IRIS window Wavelength(Å) Ion Width (må) Intensity tied to Ratio H Cii H 2, Si H Fexii Si ii Siii 26 Siii Feii 41 Feii H H Oi Siii 26 Siii Feii 41 Feii Fe xxi Ci Feii H 2 31 H : The emission lines with fixed positions in the multi-gaussian fit. 2: The emission lines with constrained positions in the multi- Gaussian fit. Fig. 8. The variation of the Fexxi / ratio as a function of density, which is from CHIANTI V (Landi et al. 213). The temperature response, E f G(T peak ), for the AIA 131 Å channel and the IRIS 1354 Å line are shown in Fig. 9. The AIA function has been computed with the procedure aia get response.pro with the calibration appropriate to the time of the observation at the default pressure in CHANTI of 15 cm 3 K (Boerner et al. 212). The keyword chiantifix is also set to account for emission that is not included in the CHIANTI database (Testa et al. 212; DeRosa & Slater 213). The response of AIA to high and low temperature flare emission is believed to be accurate to 25% with these empirical corrections (Boerner et al. 212). There are two peaks: one due to Feviii formed around 6 5 K and the other due to Fexxi formed around 7 K. The main variation in the contribution function around MK comes from the dependence of ionization on temperature. The IRIS Fe xxi line is a fine structure transition within the ground state so its excitation rate is independent of temperature but at high densities (n e > 12 cm 3 ) it suffers from collisional de-excitation. Fig. 8 shows the variation of the / line ratio, computed with CHIANTI (Landi et al. 213), as a function of density for three temperatures close to the peak of the contribution function. The temperature variation at low densities varies by less than 7% over the range shown. Only at temperatures below 4 MK is there a significant temperature dependence on the ratio. We assumed the low density limit and temperature 11 MK and computed the Fexxi response function as the line emissivity divided by the radiometric conversion coefficient, 296 (De Pontieu et al. 214), and the spectral scale (.128 Å). Again, the electron pressure was set to 15 cm 3 K. The radiometric conversion for the IRIS spectra is based on the International Ultraviolet Explorer (IUE) spectral radiances, which has an uncertainty of 15% (one σ). Both the AIA and IRIS emissivities are computed at the same pressure, elemental abundance and ionization. In the analysis of the Fexxi contribution to the 131 channel, we assume that the ratio of the intensities is constant. Fe xxi emission measure maps can be obtained by dividing the observed intensities (DN s 1 ) by the value of the peak response, at 11 MK, of the respective functions. To compare the IRIS and AIA results pixel-by-pixel pseudo rasters were constructed from a series of AIA images by selecting, for each IRIS Fig. 9. Temperature response of the AIA 131 Å channel (top) and the IRIS Fe xxi 1354 Å line (bottom). The red lines are drawn at the peak of the response curves, 11 MK. raster position, the overlapping AIA pixels from the closest in time AIA image. Along the slit, each AIA pixel (.6 ) is about 3.6 IRIS pixels (.167 ). We map the IRIS pixels onto the AIA grid and average the IRIS intensity in each AIA pixel. The IRIS raster step size (1.1 ) is larger than the AIA pixels so in the raster direction the single overlapping IRIS pixel is used. Thus, the IRIS data have been averaged across 3.6 pixels in the slit direction, so that the IRIS data are re-binned to the scale of the AIA data. In the raster direction, the AIA pseudo rasters are taken from the closest time AIA data. 4. Results 4.1. Images and spectra There were two well-observed solar flares, as shown in Fig. 2. The IRIS spectra of the two flares were fitted, as described, with 17 Gaussians superimposed on a linear background. The primary 5

6 Fig.. Images of the first flare: (top) IRIS Fexxi Å and Ci Å intensity; (middle) AIA 131 Å and 16 Å pseudo-raster intensity; (bottom) IRIS Fe xxi Doppler velocity and continuum intensity. The green rectangle marked with + s indicates the region and sites of the spectra shown in Fig. 11. The arrow indicates the time direction, and the duration of these images is 11 m 4 s. The duration of the raster is marked in Fig. 1. lines of interest are the Fexxi Å and, to a lesser extent, Ci Å. Figs. and 12 compare IRIS Fexxi and Ci intensity images with AIA 131 and 16 Å pseudo-rasters of flare 1 and flare 2. The regions of enhanced Ci identify the flare kernels. The Fexxi Doppler velocity is shown in the bottom left panel. From the velocity images we can see that the Fexxi is blue shifted near the start of the flare near the flare kernels. As shown, bright Ci overlaps very well with the bright 16 Å emission. The highest velocity measured is about 2 km s 1. To look closely at the emission from the kernels, we extract spectra from almost neighboring positions along the slit. The points are marked by a plus ( + ) in Figs and 12. In Figs 11 (flare 1) and 13 (flare 2), we show spectral images and spectra from the first raster step with Fe xxi kernel emission, as well as spectra taken at the same positions from the earlier and later raster steps. The resultant fits appear to be very good. In both flares, the Fexxi appears as a broad line below the narrow chromospheric lines. The main sources of uncertainty in the Fe xxi intensity come from the background fit, the edge of the spectral window, and the possibly non-gaussian Fe xxi profile. The background is fitted across all spectral windows so is fairly well constrained, and the window edge only affects the few positions where Fexxi is blue-shifted to by more than 15 km s 1. As well as assessing the fits by eyes we also computed the χ 2 for the spectral region around the Fexxi line. The values were generally less than one. At one or two positions, theχ 2 was large (> ). Here the peak of the Ci line could not be fit correctly by a Gaussian possibly due to optical depth effects. Occasionally individual narrow lines blending with Fexxi causedχ 2 of about 3 at the position of the Fexxi line. Since these were much narrower than the Fexxi they did not affect the Fexxi intensity. Fig. 11. Flare-kernel spectra from the region inside the green rectangle on Fig.: (left-middle) the spectrogram of the region marked with a green rectangle in Fig., with (left-top) before and (left-bottom) after spectrograms of it. The right panels show the spectra (solid) and their fits (dashed) at the positions marked by the green lines in the left panel. The vertical solid lines mark the rest wavelength of Fe xxi and the vertical dashed-dot lines mark the rest wavelength of Ci line. The green solid line is the background, and the blue dashed line is the background+ Fexxi line fit, the red dashed line is the background+ Ci line fit The EM of AIA 131 Å and Fexxi From equation (2), we can estimate the EMs of the AIA 131 and IRIS Fexxi emitting plasma. Figs. 14 shows the result for the first flare. The upper panels are the images of the EMs with contours of the IRIS Fexxi EM. The images of the EMs have a logarithmic scaling, and show good qualitative agreement. The lower panel shows the comparison of the two EMs. To look closer at the overlap, we select the smaller region, surrounded by the white box, and compare row-by-row the EMs of AIA 131 and Fexxi in the panel underneath. The structures are well matched. The AIA 131 EMs are more diffuse and greater than Fexxi EMs except around positions 1 and 2. As a check on the Gaussian line fit, we show the spectra and fits at the numbered positions in Fig. 15. All fits appear to be very good and there is no special feature of position 1 or 2 to indicate why the computed Fexxi EM is more than the 131 EM. We checked whether the higher Fexxi than 131 EMs were due to temporal mismatch between IRIS and AIA. AIA produced two 131 Å exposures during each IRIS exposure (3 s) so we compared with both of the overlapping AIA images separately for these and all other points and found only very small differences. We therefore conclude that it is not a temporal effect. It could be an effect of the two different spatial resolutions. The IRIS slit only samples about half of the AIA pixel so if the features are very concentrated, then the filling factor may be less in AIA than in IRIS. This would result in a lower AIA than IRIS EM. Similar plots for the second flare are shown in Figs 16 and 17. The region of bright Fexxi is slightly larger in the second flare so we split the region in two, labeled R1 and R2, in Fig. 16. Again selected spectra are shown (see, Fig. 17). The overall behavior is similar to that seen in the first flare: the 131 EMs are higher and more diffuse than the Fexxi EMs. At posi- 6

7 Fig. 12. Images of the second flare. (top) IRIS Fe xxi Å and Ci Å intensity; (middle) AIA 131 Å and 16 Å pseudo-raster intensity; (bottom) IRIS Fe xxi Doppler velocity and continuum intensity. The green rectangle marked with + s indicates the region and sites of the spectra shown in Fig. 13. The arrow indicates the time direction, and the duration of these images is 11 m 4 s. The duration of the raster is marked in Fig. 1. Fig. 13. Flare-kernel spectra from the region inside the green rectangle on Fig.12: (left-middle) the spectrogram at the time of the rectangle region in Fig. 12, with (top) before and (bottom) after spectrogram. Right: spectra and their fits (dashed) at the positions marked by the green lines in the left panel. The vertical solid lines indicate the rest wavelength of Fe xxi, and the vertical dashed-dot lines indicate the rest wavelength of C i. The green solid line is the background, and the blue dashed line is the background+fexxi line fit, the red dashed line is the background+ Ci line fit. tion 4, the IRIS EM is larger than the 131. As shown in Fig. 17, this point is also on the edge of a loop, so it may again be indicating a slight coalignment mismatch. We also notice that at position 2 which is near the peak in the 131 emission the IRIS EM is only about two thirds of the 131 EM. As seen in Fig. 17, the IRIS spectra at this position show a blue shift of about.5 Å or km s 1. Blue shifts in flares are usually linked to chromospheric evaporation which implies that one may expect a range of ionization states, including Fe viii at this position. Further support for this conclusion comes from the time evolution that showed brightening at this site in all AIA channels at the time of the IRIS observations. To obtain quantitative values we separate the two flare regions into loops and kernels and compare the IRIS-to-AIA EM ratios as a function of Fexxi or 131 intensity. We plot the ratio of 131 EM to the sum of 131 and Fexxi EMs to avoid going off the scale when the Fexxi EM is very low. Thus when the 131 EM is due to Fexxi alone, the ratio is.5. A contribution to the 131 EM from cooler lines, continuum, or Fe xxiii would produce a value closer to one, and where there is no Fexxi contribution, the value is one. The results for the two flares are shown in Figs. 18 and 19. The selection of loop and kernel regions is based on the AIA 131 and 16 Å intensities. Pixels with 131 Å intensities greater than twice 131 Å average and 16 Å intensities less than twice 16 Å average are selected as loops and those with 16 Å intensities greater than twice 16 Å average are selected as kernels. For loops the average value of the ratio is.54 in flare 1 and.55 in flare 2. The ratio is almost constant across the intensity range. This implies that typically about 8-85% of the Fig. 14. Upper: The EMs computed from AIA 131 Å and IRIS Fexxi at 11 MK for the first flare. The white box is the region over which we compare EMs in the panels below. Bottom: Pixelby-pixel and row-by-row EMs across the white box region, starting from the upper left. The black line is the Fexxi EM, and the red line is the 131 Å EM. The plus ( + ) indicates those pixels where the Fexxi emission is higher, and marks the selected pixels, shown in Fig Å loop emission is from Fexxi which is consistent with Milligan & McElroy (213) who estimate a 2% continuum contribution to the 131 Å channel during flares. For flare 1, the ratio ranges from.2 to.6. The low values can be attributed to alignment or filling factor issues. For example the single point in the highest intensity bin in flare 1 is from 7

8 Fig. 15. Selected IRIS spectra and AIA intensities for the first flare. The vertical dashed lines in the AIA images are the positions of the IRIS slits. The numbers on the IRIS images show the positions of the selected points in Fig. 14. The + and in the AIA images mark the positions of these points in the AIA images. Underneath are the IRIS spectra (solid black line) and their fits (black dashed line) for the selected positions. The green solid line is the background, the blue dashed line is background + Fexxi line, while the vertical line is the rest wavelengths of the Fe xxi line. Fig. 16. The EMs from AIA 131 Å and IRIS Fexxi at 11 MK for the second flare. The two white boxes are the selected regions (R1 and R2), respectively. In the lower panels, the black line is the Fexxi EM, and the red line is the 131 Å EM. The pluses ( + ) mark pixels where the Fexxi EM is higher, and the marks selected pixels, shown in Fig. 17 with higher 131 EM. a bright point on the edge of the kernels. The two points in the next highest intensity bin with a ratio below.5 are also from the edge of the kernel region. The other points with ratios less than.5 are from the bright western edge of the 131 Å region where both the Fexxi and 131 Å emission rapidly decreases. The higher ratio values are mostly from the lowest intensity bin on the edge of the loop regions. These may be due to cooler material on the edge of the loop or background 131 Å emission. Therefore the typical loop value in flare 1 is the mean of the low intensity bins,.55 or 81% Fexxi. The flare 2 loops have more scatter in their ratio. In the highest intensity bin, there are four points with a ratio around.6 (67% Fexxi). These are from the bright region at the center of R1. As mentioned, these points have blueshifted Fe xxi and at the same time brightened in all AIA channels, suggesting it is a site of chromospheric evaporation, not a typical loop. The scatter in the low intensity bins is mostly caused by coalignment issues in the R2 region as can be seen in Fig. 16. The last row, around position 7, produced many of the high ratio points. These points are along the edge of the loop and may, similar to the edges of the loops in flare 1, maybe due to cooler material on the edge of the loop or background 131 Å emission. In flare 2, the average value.55 (81% Fe xxi is typical for the 131 Å loops. The ratio in the kernels covers the full range from % Fexxi. We have plotted the ratio as a function of both Fexxi and 131 intensity. As the intensities increase the contribution of Fexxi to the 131 Å channel increases up to the value seen in the loops. In flare 1, as noted above and shown in Fig. 18, the Fig. 17. Selected IRIS spectra and AIA intensities for the second flare. The vertical lines in the AIA images mark the positions of the IRIS slits. The numbers on the IRIS images show the positions of the selected points in Fig. 16. The + and in the AIA images mark the positions of these points in the AIA images. Underneath are the IRIS spectra (solid black line) and their fits (black dashed line). The green solid line is the background, the blue dashed line is background + Fe xxi line, while the vertical line is the rest wavelength of Fexxi line. 8

9 Fig. 18. Quantitative analysis of the Fexxi and 131 Å EMs for flare1. The top row from left to right shows Fexxi, 16 Å and 131 Å intensity. The middle panels show the mean and median (red line) ratios EM(131)/(EM(131)+EM(Fexxi)) vs Fe xxi intensity in loops (bin=4 s 1 ) and kernels (bin=4 s 1 ), and on the right the ratio EM(131)/(EM(131)+EM(Fexxi)) vs intensity of 131 Å (bin=4 s 1 ). The bottom panels show the same ratios vs intensity but as scatter plots. In the loops, the mean, standard deviation and median ratios are.54±.7, and.55. Fexxi EM is greater than the 131 EM in the brightest pixels. We attribute this to imperfect coalignment or different filling factors for the two instruments, so the ratio in the two highest intensity bins of flare 1 are not included in our final conclusions. Below a 131 Å intensity of 6 DN s 1, a significant fraction of pixels have no Fexxi emission. Thus there is structure within the kernels on the scale of the AIA pixel size,.6. Considering the 131 pixels in the middle of the intensity range, one sees from Figs. 18 and 19 that the ratio is between.55 and.75 which implies the 131 kernel EM is between 2-6% greater than the Fexxi EM. This broad range again implies high resolution structure in the kernels not seen in the 131 Å images. If the continuum emission is about 2% of the Fexxi as found by Milligan & McElroy (213), then our result implies from zero to 52% of the flare kernel 131 Å emission is due to cooler plasma emission. These results assume that both instruments are well calibrated and the data are coaligned. The fact that we obtain a consistent value of 8% Fexxi contribution to the 131 channel in the loops suggests that the calibration is good. A 2% error in the intercalibration would lead to a % change in the Fexxi contribution. 5. Conclusions and Discussions Using the high resolution IRIS spectral data, together with the SDO data, we study two solar flares in AR that occurred on 24 October 213. We obtain the intensity of Fexxi from IRIS data and compare it with AIA data at 131 Å. Coalignment to an accuracy of one AIA pixel (.6 ) is achieved by coaligning Fig. 19. Quantitative analysis of the Fexxi and 131 EMs for flare 2. The top row from left to right shows Fexxi, 16 Å and 131 Å intensity. The middle panels show the mean and median (red line) ratios EM(131)/(EM(131)+EM(Fexxi)) vs Fexxi intensity in loops (bin=3 DN s 1 ) and kernels (bin=3 DN s 1 ), and on the right the ratio EM(131)/(EM(131)+EM(Fexxi)) vs intensity of 131 Å (bin= 3 DN s 1 ). The bottom panels show the same ratios vs intensity but as scatter plots. In the loops, the mean, standard deviation and median ratios are.55±.6, and.54. the IRIS Ci to AIA 16 Å (see Fig. 3) and the IRIS Fexxi to AIA 131 Å images. Fexxi is a hot and broad line which, in flare kernels, is blended with many cold lines. Previous studies (Doschek et al. 1975; Cheng et al. 1979; Mason et al. 1986) at lower resolution only consider blending of Fexxi with the Ci line. However, as IRIS reveals, there are several additional cold chromospheric lines blending with the Fe xxi (see Fig. 7). After studying other IRIS flare kernel spectra, we were able to identify lines in different parts of the IRIS spectrum that can be used to constrain the intensity, widths and positions of the blending cold lines (see Table 1). Using 17 Gaussian lines and a linear continuum we are able to fit the spectra and obtain reliable Fexxi intensities at all positions. The spectra show regions at the onset of the flares where of Fe xxi is coincident with cold-line emission (see Fig. and 12). We have carefully checked the fits to make sure that the procedure has worked at these positions (see Fig. 11 and 13). In these regions, the Fexxi line shows high blue shifts up to 2 km s 1, which is the edge of the spectral window. This velocity is also consistent with other hot spectral lines seen in flares, such as Fe XIX (Czaykowska et al. 21; Teriaca et al. 23; Brosius & Phillips 24), Fe XXIII and Fe XXIV (Milligan & Dennis 29). In our study, the EMs have the magnitude of about 28 cm 5, some pixels can reach 29 cm 5 (see Fig. 14 and 16), this also agrees with others (Graham et al. 213; Fletcher et al. 213). It is well known that the 131 Å channel in AIA has contributions from the Fexxi and the Feviii line. Brosius & Holman 9

10 (212) list a number of other transition region lines that may also be present. We compare the EMs of AIA 131 Å and the IRIS Fe xxi line pixel-by-pixel, assuming all emission is from hot plasma. For loop regions the AIA 131 Å EMs are about 2% greater than Fexxi EMs (see Fig. 18 and 19). This is consistent with the analysis of Milligan & McElroy (213) that 2% of the 131 Å channel emission is due to continuum in flares. We also identified a loop region with about 33% greater 131 than Fe xxi EM. Since, this region simultaneously brightened in all AIA channels and was associated with Fexxi blue shifts, we attribute the lower than average Fexxi to additional Feviii due to a recent injection of cool plasma into the loop. In the brighter 131 flare kernels (> 6 DN s 1 ), the 131 EM is about 2-6% more than the IRIS Fexxi suggesting that there is an up to 52% contribution from cool plasma emission (Brosius & Holman 212). Small-scale kernel structure results in a broad range of 131/Fexxi EM ratios and sharp gradients in IRIS Fexxi emission at sites of molecular and transition region emission. At the start of the flare, consecutive raster steps showed a large increase in Fe xxi across the kernels indicating that the kernel structure is both temporal and spatial. Future work will focus on the analysis of flares obtained at higher spectral resolution with broader wavelength windows and higher temporal cadence to resolve the temporal/spatial ambiguities. Ning, Z. & Cao, W. 211, Sol. Phys., 269, 283 O Dwyer, B., Del Zanna, G., Mason, H. E., Weber, M. A., & Tripathi, D. 2, A&A, 521, A21 Qiu, J., Lee, J., Gary, D. E., & Wang, H. 22, ApJ, 565, 1335 Teriaca, L., Falchi, A., Cauzzi, G., et al. 23, ApJ, 588, 596 Teriaca, L., Falchi, A., Falciani, R., Cauzzi, G., & Maltagliati, L. 26, A&A, 455, 1123 Testa, P., Drake, J. J., & Landi, E. 212, ApJ, 745, 111 Tian, H., DeLuca, E., & Reeves, K. K., et al. 214, ApJ, 786, 137 Tian, H., Young, P. R., Reeves, K. K., et al. 215, ApJ, 811, 139 Wang, T. J., Solanki, S. K., Curdt, W., et al. 23, A&A, 46, 15 Young, P. R., Doschek, G. A., Warren, H. P., & Hara, H. 213, ApJ, 766, 127 Young, P. R., Tian, H., & Jaeggli, S. 215, ApJ, 799, 218 Acknowledgements. The authors would like to thank the anonymous referee for valuable comments to improve the manuscript. 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 down-link communications funded by the Norwegian Space Center (NSC, Norway) through an ESA PRODEX contract. We are indebted to the IRIS and SDO teams for providing the high resolution data. We would also thank the people at MPS, in particular to Professor Hardi Peter, Doctors Li, Leping, Chen, Feng and Guo, Lijia. References Antonucci, E., Gabriel, A. H., Acton, L. W., et al. 1982, Sol. Phys., 78, 7 Benz, A. O. 28, Living Reviews in Solar Physics, 5, 1 Boerner, P., Edwards, C., Lemen, J., et al. 212, Sol. Phys., 275, 41 Brosius, J. W. & Holman, G. D. 212, A&A, 54, A24 Brosius, J. W. & Phillips, K. J. H. 24, ApJ, 613, 58 Cheng, C.-C., Feldman, U., & Doschek, G. A. 1979, ApJ, 233, 736 Cheng, X., Ding, M. D., & Fang, C. 215, ApJ, 84, 82 Czaykowska, A., Alexander, D., & De Pontieu, B. 21, ApJ, 552, 849 De Pontieu, B., Title, A. M., Lemen, J. R., et al. 214, Sol. Phys., 289, 2733 DeRosa, M. & Slater, G. 213, Guide to SDO Data Analysis edited on February 19, 213 Doschek, G. A., Dere, K. P., Sandlin, G. D., et al. 1975, ApJ, 196, L83 Feldman, U., Curdt, W., Landi, E., & Wilhelm, K. 2, ApJ, 544, 58 Fletcher, L., Dennis, B. R., Hudson, H. S., et al. 211, Space Sci. Rev., 159, 19 Fletcher, L., Hannah, I. G., Hudson, H. S., & Innes, D. E. 213, ApJ, 771, 4 Fletcher, L. & Hudson, H. 21, Sol. Phys., 24, 69 Fletcher, L., Pollock, J. A., & Potts, H. E. 24, Sol. Phys., 222, 279 Gan, W. Q., Li, Y. P., & Miroshnichenko, L. I. 28, Advances in Space Research, 41, 98 Graham, D. R., Hannah, I. G., Fletcher, L., & Milligan, R. O. 213, ApJ, 767, 83 Graham, D. R., & Cauzzi, G. 215, ApJ, 87, L22 Innes, D. E. 21, A&A, 378, 67 Innes, D. E. 28, A&A, 481, L41 Kliem, B., Dammasch, I. E., Curdt, W., & Wilhelm, K. 22, ApJ, 568, L61 Landi, E., Young, P. R., Dere, K. P., Del Zanna, G., & Mason, H. E. 213, ApJ, 763, 86 Lemen, J. R., Title, A. M., Akin, D. J., et al. 212, Sol. Phys., 275, 17 Li, D., Ning, Z. J., & Zhang, Q. M. 215, ApJ, 87, 72 Mason, H. E., Shine, R. A., Gurman, J. B., & Harrison, R. A. 1986, ApJ, 39, 435 Mclntosh, S. W., De Pontieu, B., Hansteen, V., & Boerner, P. 213, A User s Guide To IRIS Data Retrieval, Reduction and Analysis, 213, October, 3 Milligan, R. O. & Dennis, B. R. 29, ApJ, 699, 968 Milligan, R. O., Gallagher, P. T., Mathioudakis, M., et al. 26, ApJ, 638, L117 Milligan, R. O. & McElroy, S. A. 213, ApJ, 777, 12