<figure><img src="image_1.png" /></figure>

# Transient Inverse-FIP Plasma Composition Evolution within a Solar Flare

Deborah Baker<fnref n="1" /> ID, Lidia van Driel-Gesztelyi<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID, David H. Brooks<fnref n="4" /> ID, Gherardo Valori<fnref n="1" /> ID, Alexander W. James<fnref n="1" /> ID, J. Martin Laming<fnref n="5" /> ID, David M. Long<fnref n="1" /> ID, Pascal Démoulin<fnref n="2" /> ID, Lucie M. Green<fnref n="1" /> ID, Sarah A. Matthews<fnref n="1" /> ID, Katalin Oláh<fnref n="6" /> ID, and Zsolt Kővári<fnref n="6" /> ID

Received 2018 December 11; revised 2019 February 7; accepted 2019 February 15; published 2019 April 10

## Abstract

Understanding elemental abundance variations in the solar corona provides an insight into how matter and energy flow from the chromosphere into the heliosphere. Observed variations depend on the first ionization potential (FIP) of the main elements of the Sun's atmosphere. High-FIP elements (>10 eV) maintain photospheric abundances in the corona, whereas low-FIP elements have enhanced abundances. Conversely, inverse FIP (IFIP) refers to the enhancement of high-FIP or depletion of low-FIP elements. We use spatially resolved spectroscopic observations, specifically the Ar XIV/Ca XIV intensity ratio, from *Hinode's* Extreme-ultraviolet Imaging Spectrometer to investigate the distribution and evolution of plasma composition within two confined flares in a newly emerging, highly sheared active region. During the decay phase of the first flare, patches above the flare ribbons evolve from the FIP to the IFIP effect, while the flaring loop tops show a stronger FIP effect. The patch and loop compositions then evolve toward the preflare basal state. We propose an explanation of how flaring in strands of highly sheared emerging magnetic fields can lead to flare-modulated IFIP plasma composition over coalescing umbrae which are crossed by flare ribbons. Subsurface reconnection between the coalescing umbrae leads to the depletion of low-FIP elements as a result of an increased wave flux from below. This material is evaporated when the flare ribbons cross the umbrae. Our results are consistent with the ponderomotive fractionation model for the creation of IFIP-biased plasma.

**Key words:** Sun: abundances – Sun: corona – Sun: magnetic fields

**Supporting material:** animations

## 1. Introduction

Elemental abundance variations are tracers of physical processes throughout the universe, with the cosmic reference standard being the solar chemical composition. Understanding how the Sun's chemical composition varies in time and space provides an insight into how mass and energy flow from the Sun's chromosphere into the heliosphere and in turn, from the chromospheres of solar-like stars into their astropheres (e.g., Testa 2010; Laming 2015; Testa et al. 2015).

From over 50 yr of spectroscopic observations we know that the solar corona has different elemental composition from that of the photosphere (Pottasch 1963; Meyer 1985a, 1985b). In the corona, elements of low first ionization potential (FIP; ≤10 eV) such as Fe, Si, Mg, and Ca are enhanced by a factor of two to four compared to their photospheric abundances. This is known as the FIP effect and it is typically expressed in terms of the FIP bias parameter which is the ratio of an element's coronal and photospheric abundances. The enhancement of low-FIP elements varies depending on the coronal structure, e.g., coronal holes show little plasma fractionation (FIP bias ~1 or photospheric composition; Feldman & Widing 1993; Doschek et al.

1998; Feldman et al. 1998; Brooks & Warren 2011; Baker et al. 2013), quiet-Sun regions typically have an FIP bias of 1.5–2 (Feldman & Widing 1993; Doschek et al. 1998; Warren 1999; Baker et al. 2013; Ko et al. 2016; Baker et al. 2018), whereas active regions (ARs) consist of highly fractionated plasma with FIP bias of 3–4 (Feldman 1992; Widing & Feldman 1995; Sheeley 1996; Brooks & Warren 2011, 2012; Baker et al. 2013; Del Zanna & Mason 2014; Baker et al. 2015; Brooks et al. 2015). (See reviews by e.g., Feldman & Widing 2003; Schmelz et al. 2012; Laming 2015 and Chapter 14 of Del Zanna & Mason 2018.)

Unresolved Sun-as-a-star observations using full disk integrated spectra show the solar FIP effect of the active Sun. For temperatures greater than ~1 MK, Laming et al. (1995) found low-FIP elements were enhanced by a factor of 3–4 which is in line with the FIP bias value for the element Fe deduced by Schonfeld et al. (2015). The FIP effect is reduced at lower temperatures (Laming et al. 1995). Recent work of Brooks et al. (2017, 2018) based on Sun-as-a-star spectra obtained by the *Solar Dynamics Observatory* (SDO) Extreme-ultraviolet Variability Experiment (EVE) demonstrated that the variation of coronal composition is highly correlated with the solar cycle.

The coronae of solar-like stars exhibit a varying degree of the solar-like FIP effect in their X-ray spectra whereas more active cool dwarf stars show an inverse FIP (IFIP) effect. In

<figure><img src="image_2.png" /></figure>

Original content from this work may be used under the terms of the [Creative Commons Attribution 3.0 licence](https://creativecommons.org/licenses/by/3.0/). Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

<footnote n="1">¹ University College London, Mullard Space Science Laboratory, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK; deborah.baker@ucl.ac.uk</footnote>

<footnote n="2">² LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, F-92195 Meudon, France</footnote>

<footnote n="3">³ Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, Konkoly Thege út 15-17., H-1121, Budapest, Hungary</footnote>

<footnote n="4">⁴ College of Science, George Mason University, 4400 University Drive, Fairfax, VA 22030, USA</footnote>

<footnote n="5">⁵ Space Science Division, Naval Research Laboratory, Code 7684, Washington, DC 20375, USA</footnote>

<footnote n="6">⁶ Konkoly Observatory of the Hungarian Academy of Sciences, Budapest, Hungary</footnote>
surveys of M dwarf stars and active binaries, low-FIP elements are under-abundant relative to high-FIP elements, e.g., O, Ne, and Ar (e.g., Brinkman et al. 2001; Argiroffi et al. 2005; Robrade & Schmitt 2005, 2006; Telleschi et al. 2005; Wood & Linsky 2006; Liefke et al. 2008; Wood et al. 2012). Wood & Linsky (2010) firmly established the dependence of the (I)FIP effects on F to K spectral type stars with X-ray luminosity $<10^{29}$ erg s$^{-1}$. The solar FIP effect decreases from G to early K-type stars and becomes zero at ~K5 then reverses to the IFIP effect for later K stars and M dwarfs. Wood et al. (2018) extended the FIP bias-spectral type relationship to include stars of earlier spectral types A and F. It is worth noting that the observed composition of stellar coronae can be more complex than having either a straightforward FIP or IFIP effect. For example, Peretz et al. (2015) found that all elements were consistently depleted in the coronae of six main-sequence stars of spectral type F7–K1 compared to their respective photospheres, whether compared to solar abundances or the individual stellar abundances.

Surveys of plasma composition in solar flares show notable variability in elemental abundances. Warren (2014) measured absolute abundances for 21 M9.3 to X6.9 class flares using SDO/EVE spectra. The mean FIP bias in their sample is close to photospheric composition, which is consistent with the earlier results of Veck & Parkinson (1981), Feldman & Widing (1990), McKenzie & Feldman (1992), and more recently with Del Zanna & Woods (2013). However, in large samples of solar flares observed by multiple instruments, low-FIP element abundances are enhanced by a factor of ~2–3, depending on the emission lines and instruments used to determine the FIP bias, the atomic data used at the time of the measurements, and temperature effects (e.g., Doschek et al. 1985; Sterling et al. 1993; Bentley et al. 1997; Fludra & Schmelz 1999; Phillips & Dennis 2012; Dennis et al. 2015; Sylwester et al. 2015). Spatially resolved flare observations from Hinode/EIS provide clues to understanding the abundance variability observed in solar flares. Doschek et al. (2018) show that the FIP bias varies from coronal composition in the post-flare loops to photospheric composition in the loop footpoints for an X8.3 flare and Warren et al. (2018) find coronal abundances in the current sheet of the same limb flare. At least for the solar case, it is possible that the spatially unresolved observations (e.g., Warren 2014) show little FIP effect in flares because photospheric abundances are dominant and the coronal composition of specific features such as post-flare loops is only evident in spatially resolved images.

Studies of abundance changes in stellar flares have been mainly limited to large flares on active stars, i.e., to stars with IFIP composition. Using high-resolution X-ray spectroscopy, Nordon & Behar (2007, 2008) analyzed abundance variations during 14 large flares observed by the XMM-Newton and Chandra observatories. In 7 of 14 flares, they found a trend of enhanced low-FIP elements tending toward photospheric composition during the flares for stars with IFIP-biased quiescent coronae. The opposite trend was observed in two stars with FIP bias dominated plasma comprising their coronae and five stars showed no effect during flaring. Their results are consistent with case studies of solar-like stars (e.g., Testa et al. 2007), active M dwarf EV Lac (Laming & Hwang 2009), active M dwarf CN Leo (Liefke et al. 2010), and RS CVn binary stars (e.g., Güdel et al. 1999; Audard et al. 2001, 2003). Indeed, our knowledge is presently limited for flares in stars

because of the difficulties in determining abundance variations in stellar coronae; see the extensive discussion in Testa (2010).

Any theoretical framework for a fractionation mechanism must be able to account for all of the solar and stellar observations including the IFIP effect. Early models based on, for example, thermal diffusion and Coulomb drag could reproduce some observational aspects of the FIP effect but not those of the IFIP effect (Laming 2015, and references therein). To date, only the ponderomotive force model is able to produce IFIP composition. The Laming model invokes the ponderomotive force acting only on ions to separate ions from neutrals in the chromosphere of the Sun and other stars (Laming 2012, 2015, 2017). Fractionation takes place in the chromosphere at temperatures where low-FIP elements are mainly ionized, but high-FIP elements remain neutral. The ponderomotive force arises as Alfvén waves reflect from or refract at the high density gradient in the chromosphere of the Sun. Alfvén waves originating in the corona cause FIP fractionation at the chromospheric level at loop footpoints, which are magnetically connected to reconnection sites of flares or nanoflares (Laming 2017). Conversely, upwardly directed photospheric acoustic waves may mode convert to fast mode waves as the plasma transitions from $β > 1$ to $β < 1$, where plasma $β$ is the ratio of plasma pressure and magnetic pressure. The IFIP effect arises as these new fast mode waves refract in the chromosphere back down again, producing a downward-directed ponderomotive acceleration (Laming 2015).

Recently, Doschek et al. (2015) reported the first observations of the IFIP effect on the Sun. Highly localized regions of IFIP were captured by Hinode/EIS (Culhane et al. 2007) near large sunspots during flares (Doschek et al. 2015; Doschek & Warren 2016, 2017). In the most extreme case, the IFIP effect locally exceeded that of the integrated value of M dwarf stars of spectral type ~M5 at the extreme end of the FIP bias-spectral type relationship of Wood & Linsky (2010), Laming (2015), and Wood et al. (2018).

In this paper, we present a detailed spectroscopic study of the evolution of plasma composition in the very active AR 11429. We find that patches of IFIP bias plasma appear and disappear within the FIP bias dominated composition of the AR during a confined flare, while IFIP bias plasma was not observed in another confined flare nine hours later. We discuss how the relatively rare magnetic field configuration of the AR and its evolution relate to the spatial locations and temporal evolution of IFIP patches and we consider whether flaring plays a role in creating the anomalous plasma composition or simply reveals the already present IFIP plasma.

## 2. Overview of AR 11429

### 2.1. Coronal Activity

AR 11429 appeared at the NE limb on 2012 March 3. During its disk transit, AR 11429 was the source region of 3 X-class, 14 M-class, 32 C-class flares, and 4 coronal mass ejections (CMEs), making it one of the most flare and CME productive ARs of Solar Cycle 24. The high activity level and rare magnetic field configuration have inspired a large number of studies of the AR (e.g., Petrie 2012; Donea & Hanson 2013; Chintzoglou et al. 2015; Sun et al. 2015; Patsourakos et al. 2016; Syntelis et al. 2016; Polito et al. 2017). In this study, we focus on two confined M-class flares observed by Hinode/EIS on March 6. They are identified as FL1 and FL2 in the GOES
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Top panel: GOES soft X-ray curve from 12:00 UT on 2012 March 6 to 00:00 on March 7 with Hinode/EIS raster times indicated by the dashed red lines. Flare 1 (M2.2) peaked at ~12:37 UT and Flare 2 (M1.4) peaked at ~21.12 UT. Bottom panel: clockwise from upper left: SDO/HMI line-of-sight (LOS) magnetogram, SDO/AIA 1600 Å, 94 Å, and 171 Å maps of AR 11429 at 12:39 on 2012 March 6 overplotted with the Hinode/EIS field of view (dashed white box). The animation of the bottom panels begins at 11:44:44 and ends at 22:29:44. Its duration is 28 s. Note that the cadence of the animation is 45 s for the flare periods but is 5 minutes from 14:00 UT to 20:15 UT.
(An animation of this figure is available.)</figcaption></figure>

1–8 Å X-ray flux curve in Figure 1. The first flare starts at 12:23 UT, peaks at 12:37 UT, and decays to background flux levels by ~13:45 UT and the second flare begins at 21:04 UT, peaks at 21:12 UT, and quickly decays by ~21:25 UT. FL1 and FL2 flare classifications are M2.2 and M1.4, respectively.

## 2.2. Coronal Evolution

The preflare coronal configuration of AR 11429 was a highly sheared structure. The hot-channel passbands, e.g., 94 Å and 131 Å of the SDO's Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) have sheared bright loops evident on March 5. One day later, the bright sheared loop system has evolved into a highly twisted “corkscrew” on its NE end at the beginning of the rise phase of flare FL1. Figure 1 (bottom panel) shows a frame from its animation at 12:39 UT on March 6 in which the “corkscrew” feature is prominent in each of the SDO/AIA passbands (1600, 171, and 94 Å). Figure 1's animation spans from 11:44 UT to 22:29 UT, covering both

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Location of ribbon pairs R1–R2 and R3–R4 in SDO/AIA 1600 Å image at 12:37 UT during flare FL1.</figcaption></figure>

confined flares at a cadence of 45 s but the cadence is 5 minutes from 14:00 UT to 20:15 UT.

Syntelis et al. (2016) conducted a spectroscopic analysis of the pre-eruptive configuration of AR 11429 prior to the eruption of two CMEs that occurred early on 2012 March 7. Using the same Hinode/EIS data as in our present study, they found substantial spectroscopic evidence for the presence of a hot flux rope in the NE section of the AR that formed during flare FL1. Two distinct plasma components were identified by Syntelis et al. (2016), one at 1.6–2.5 MK ($\log_{10} T = 6.2–6.4$ K) and the other at 6.3–12.6 MK ($\log_{10} T = 6.8–7.1$ K). The hotter component contained regions of enhanced nonthermal line broadening, relatively strong Doppler upflow velocities, and lower plasma densities. Collectively, these spectral parameters are characteristic of hot flux ropes.

## 2.3. Flare Ribbons

The lower solar atmosphere evolved in parallel with the corona during confined flare FL1. Elongated bright ribbons appear early in the rise phase of FL1; however, there is evidence of preflare heating at these locations in the SDO/AIA 1600 Å and 171 Å passbands (see Figure 1's animation). At the peak of the flare, two pairs of flare ribbons, R1–R2 and R3–R4, are visible in the SDO/AIA 1600 Å passband in Figure 2. Ribbon pair R1–R2 is centered on the N–S aligned, near-vertical section of the polarity inversion line (PIL; see SDO's Helioseismic and Magnetic Imager—HMI; Scherrer et al. 2012—magnetogram in Figure 1) in the NE and the second pair, R3–R4, runs along the E–W aligned, horizontal segment of the PIL in the SW of the AR. Ribbon pair R1–R2 appears ~1 minute before R3–R4. They are brightest at the flare peak (12:37 UT in Figure 2) and gradually decrease in intensity during the decay phase.

There was continuous, low-level reconfiguration of the corona in the period between the two flares when the main activity shifted to the vicinity of the ribbon pair R3–R4 of flare FL1, along the horizontal portion of the PIL. At 20:43 UT, brightenings begin to reappear in SDO/AIA 1600 Å passband in the NE section of the AR, while the main activity continues in the SW (see Figure 1's animation). The ribbons have faded in the NE by 21:06 UT, at the same time that FL2 begins in the vicinity of the former ribbons R3–R4. We observe highly
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Maps of SDO/HMI SHARP CEA continuum (left panel) and radial magnetic field, $B_r$ (right panel). The positive/negative red/blue radial magnetic field is saturated at ±1500 G. Key flux emergence episodes are labeled using red, blue, and green text (see continuum map at 02:58 on March 6). Footpoints of magnetic flux ropes 1 and 2 (red and blue, respectively) in the continuum map at 12:34 on March 6 correspond to the NLFF extrapolation in Figure 11 of Chintzoglou et al. (2015). MFR1 is rooted in the sunspot umbrae in the north of AR 11429 whereas the footpoints of MFR2 are located further south in the quiet Sun. In the animation of this figure, the locations of IFIP-biased and photospheric plasma observed by Hinode/EIS are identified with arrows. The video begins 2012 March 3 11:58:14 and ends the same day at 23:58:14. The video duration is 9 s.</figcaption></figure>

(An animation of this figure is available.)
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Left to right: *Hinode/EIS* Ar XIV 194.4 Å and Ca XIV 193.87 Å intensity maps, Ar XIV/Ca XIV ratio maps without and with *SDO/HMI* contours of ±500 (white/purple), *SDO/AIA* 1600 Å and 94 Å maps. (*SDO/AIA* 94 Å maps are shown using reverse color table.) Top to bottom: observations are from 12:38 UT to 13:32 UT during flare FL1. The color bar scale shows the FIP effect as blue/green, photospheric composition as orange, and IFIP effect as yellow. All *Hinode/EIS* maps are coaligned to *SDO/AIA* and HMI maps at the times shown. Black vertical stripes indicate data drop out periods.</figcaption></figure>

packed ribbons over the magnetically complex horizontal PIL during the rise phase. Activity continues in the central part of the AR throughout the decay phase of FL2.

## 2.4. Magnetic Field Evolution

When AR 11429 rotated onto the disk on 2012 March 3, it was already a mature sunspot group containing several umbrae in a common penumbra. The magnetic structure of the AR was bipolar, with an anti-Hale orientation for the N hemisphere (Figure 3, continuum image at 22:10 UT on March 3). The

orientation of the bipolar fields also deviated significantly from Joy's law (Hale et al. 1919). Over the next few days, major flux emergence took place along the AR's inversion line, with the opposite-polarity new flux concentrations diverging at roughly a right angle to the line connecting the pre-existing spots. By March 6, the magnetic structure became αβγ, with two–three major bipoles still in emergence.

In Figure 3, second row from the top, the main emerging bipoles are indicated with arrows of different colors. Since the emerging flux was highly sheared, the magnetic inversion line
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Left to right: *Hinode/EIS* Ar XIV 194.4 Å and Ca XIV 193.87 Å intensity maps, Ar XIV/Ca XIV ratio maps without and with *SDO/HMI* contours of ±500 (white/purple), *SDO/AIA* 1600 Å and 94 Å maps. (*SDO/AIA* 94 Å maps are shown using reverse color table.) Top to bottom: observations are from 21:10 UT to 22:04 UT during flare FL2. The color bar scale shows the FIP effect as blue/green, photospheric composition as orange, IFIP effect as yellow. All *Hinode/EIS* maps are coaligned to *SDO/AIA* and HMI maps at the times shown.</figcaption></figure>

maintained its NE-SW orientation from March 3 throughout the emergence process. High shear in the emerging fields was evidenced in the radial magnetic field as a *yin-yang* magnetic polarity pattern (magnetic tongues; Luoni et al. 2011), indicating strong negative, left-handed helicity in the emerging flux. Early on March 6, Chintzoglou et al. (2015), using *SDO/HMI* SHARP vector magnetograms, found the mean shear angle along the PIL to be around 68°, which showed a slow increase throughout the day. Chintzoglou et al. (2015) pointed out that there were two sources of shear: highly sheared flux

emergence in the NE part of the AR, and shearing motions between emerging fields (their positive polarity) and the pre-existing negative-polarity spot in the middle of the AR, where the magnetic PIL was nearly of E-W orientation, and had a quadrupolar structure (see Figure 3, right column).

Takasao et al. (2015) carried out magnetohydrodynamic (MHD) simulations showing that the emergence of a subsurface-kinked flux rope can spontaneously form a quadrupolar non-Hale, non-Joy AR with multiple magnetic inversion lines in its midst. The complex inversion line structure forms in their
<table><caption>**Table 1**
Hinode/EIS Study Details</caption><tr><td>Study Name</td><td>FlareResponse 01</td></tr><tr><td>IFIP Emission Lines</td><td>Ca XIV 193.87 Å<br>Ar XIV 194.40 Å</td></tr><tr><td>Field of View</td><td>240" × 304"</td></tr><tr><td>Rendering</td><td>2" slit, 80 positions, 3" coarse steps</td></tr><tr><td>Exposure Time</td><td>5 s</td></tr></table>

simulation by the submergence of emerged fields. This simulation gives a good description of AR 11429 and its multiple inversion lines on March 6 (see Figure 3, right column).

Using nonlinear-force-free (NLFF) magnetic extrapolations, Chintzoglou et al. (2015) found two flux ropes on March 6 in AR 11429: magnetic flux rope 1 (MFR1) in the NE emerging fields, while flux rope 2 (MFR2) formed around the pre-existing negative spot, where strong shearing was observed around the complex E-W inversion line area. MFR1 had its footpoints in the forming new umbrae, while the latter (MFR2) flux rope had its footpoints over quiet-Sun areas. The locations of the flux rope footpoints are indicated in the continuum image at 12:34 UT on March 6 in Figure 3. Chintzoglou et al. (2015) suggested that MFR1, and probably MFR2, formed during the M2.2 flare (FL1) that we analyze. The flare ribbon pair R1-R2 is linked to the formation of MFR1, while ribbon pair R3-R4 is linked to MFR2 (see Figure 2 and Section 2.3).

For about a day prior to FL1, the highly sheared divergence of opposite polarities emerging along the PIL led to an inflow of negative magnetic flux moving toward the isolated pre-existing spot at the eastern footpoint of MFR1 (R1). Flux approached and collided with the spot, forcing the coalescence of the smaller flux fragments into a growing, strongly coherent umbra surrounded by a common penumbra. At the western footpoint (R2), the emerging positive flux fragments crashed into and coalesced with the main positive polarity in the center of the AR. Repeated flux emergence episodes drove umbral coalescence until ~13:00 UT during FL1. By ~15:00 UT, the umbra at R1 had ceased coalescing and by ~16:00 UT, the umbra at R2 started to break apart/decay. This process is evident in Figure 3's animation beginning at 11:22 UT on March 5 where arrows indicate each footpoint region in the continuum and radial field images.

### 3. Hinode/EIS Observations of AR 11429

The observations featured in Figures 4 and 5 were acquired while *Hinode/EIS* was operating in an autonomous observing mode during a Major Flare Watch campaign. When an intense brightening was detected in the He II 256.32 Å lines, a high cadence flare response study was triggered. A series of six rasters at a cadence of 9 minutes was run for each flare. Observations spanned 12:38 UT to 13:32 UT (FL1) and from 21:10 UT to 22:04 UT (FL2) on 2012 March 6. EIS captured the peak and the decay phase for both flares (see top panel of Figure 1). Some of the key details of the flare response study are listed in Table 1.

#### 3.1. Ar XIV/Ca XIV Intensity Ratio

The EIS flare response study contains two lines that are suitable for measuring plasma composition at temperatures higher than

<figure><img src="image_8.png" /><figcaption>**Figure 6.** SDO/HMI Continuum at 12:47 UT overlaid with contours of IFIP (green) and flare ribbons (orange). The contours overlap at the umbrae associated with ribbon pairs R1-R2. IFIP contours are from the *Hinode/EIS* Ar XIV/Ca XIV intensity ratio map at 12:47 UT in Figure 4. Flare ribbon contours are from the SDO/AIA 1600 image Å at 12:37 UT in Figure 2 adjusted for solar rotation.</figcaption></figure>

those expected in nonflaring ARs: low-FIP Ca XIV (FIP = 6.11 eV) at 193.87 Å and high-FIP Ar XIV (FIP = 15.76) at 194.40 Å (Feldman et al. 2009). Both lines are relatively strong, close in wavelength, and contain no strong blends. The two ions are formed in ionization equilibrium at similar temperatures of ~3.5 MK ($log_{10} T = 6.55$ K) (Feldman et al. 2009; Doschek et al. 2015), and they have very similar emissivity temperature dependences: the theoretical coronal ratio is about 0.25 ± 0.10 over the temperature range 1.6–6.3 MK ($6.2 \le log_{10} T \le 6.8$ K) in coronal conditions at an electron density of $log_{10} N = 10$. However, the ratio rises to levels exceeding 0.50 at a temperatures above ~8 MK ($log_{10} T = 6.9$ K). (See the Appendix for a plot of the ratio of Ar XIV and Ca XIV contribution functions at different densities calculated using the abundances given below.) With similar contribution functions, the intensity ratio of these lines can be used to determine (IFIP bias levels.

In line with Doschek et al. (2015) and Doschek & Warren (2016), we use $Log_{10}$ abundance values, relative to H (with $Log_{10}$ (H density) set to 12), as follows: corona–Ca = 6.93 (Feldman 1992), Ar = 6.50 and photosphere–Ca = 6.33 (Caffau et al. 2011), and Ar = 6.50 (Lodders 2008), yielding a coronal FIP bias of $10^{0.6} = 4$. The CHIANTI Atomic Database, Version 8.0 (Dere et al. 1997; Del Zanna et al. 2015) was employed to carry out the calculations of the contribution functions using these abundances. The abundances of Ar and Ca affect mostly the magnitude of the respective contribution functions but not their shapes as a function of temperature, so the intensity ratio allows us to determine the relative abundances of Ar to Ca. Typically, the FIP bias is the ratio of the low-FIP element's coronal to photospheric abundances and the high-FIP element's coronal to photospheric abundances; however, the convention for IFIP bias is to invert the ratio. We adopt the latter convention such that high-FIP Ar XIV/low-FIP Ca XIV line intensity ratio >1 indicates the IFIP effect, =1 is unfractionated photospheric plasma, and <1 is the solar FIP effect. The estimated uncertainty based on an intensity error of 20% is ±0.28.
All *Hinode/EIS* data were reduced using the eis_prep routine available in Solar Software (Freeland & Handy 1998). The routine converts the CCD signal in each pixel into calibrated intensity units and removes/flags cosmic rays, dark current, dusty, warm, and hot pixels. The original *Hinode/EIS* calibration was used instead of the newer calibrations of Del Zanna (2013) and Warren et al. (2014) because the lines are close in wavelength and within the same spectral window. The ratio of the corrected intensities using either calibration is within 2% of the original. We fit three Gaussians to the Ar XIV 194.40 Å line to remove two unidentified weak lines in its blue wing (Brown et al. 2008; Doschek et al. 2015) and to the spectral region around the Ca XIV 193.87 Å line to separate it from two nearby lines. Two sample spectra for the spectral window are shown in the Appendix, one of the FIP effect and the other of the IFIP effect. Typically, the unidentified weak lines in the blue wing of the Ar XIV are evident in the FIP effect spectra but are negligible in the IFIP effect spectra. Columns 3 and 4 of Figures 4 and 5 show the Ar XIV/Ca XIV intensity ratio maps.

In this work, we have followed the methodology of Doschek et al. (2015) and Doschek & Warren (2016, 2017) in producing the Ar XIV/Ca XIV intensity ratio maps. They give a full account of the assumptions and issues concerning this EIS composition diagnostic. We refer the reader to their extensive discussions contained in this series of papers, especially Doschek & Warren (2017). This method has the advantage of being simple, but does not fully account for the (relatively weak) temperature and density sensitivity of the ratio. To alleviate any concerns that these effects might explain the detection of patches of inverse FIP, as an independent check, we have analyzed one of the IFIP patches using the more complete method of computing the FIP bias from a differential emission measure (DEM) analysis. The DEM was inferred by measuring the electron density and using it to compute contribution functions for a set of emission lines covering a wide range of temperatures (Fe VIII-XXII, Ca XIV-XV). The ratio of the calculated to observed Ar XIV intensity gives the FIP bias. This analysis showed that accounting for the temperature and density sensitivity increased the magnitude of the inverse FIP effect. The reason is that the Ar XIV line is brighter than expected from theory, so taking the observed Ar XIV/Ca XIV ratio underestimates the inverse FIP effect. In addition, when the contribution functions are convolved with the DEM, it appears that the two lines are formed closer to 4.5 MK ($log_{10} T = 6.65$ K), which falls within the temperature range of the ratio values quoted above. They are also formed over a much narrower temperature range than expected from theory, so any larger variations in the ratio outside of the temperature range 3.2-6.3 MK ($6.5 \le log_{10} T \le 6.8$ K) are not important in this event (see plots of the ratio in the Appendix). The theoretical Ar XIV and Ca XIV contribution functions ($G(T)$) along with their contribution functions convolved with the DEM are shown in the Appendix.

## 3.2. IFIP Plasma Evolution During Flare FL1

Figure 4 is composed of *Hinode/EIS* and *SDO/AIA* images during flare FL1 as follows: columns 1 and 2 contain EIS Ar XIV and Ca XIV intensity images; the Ar XIV/Ca XIV ratio maps are in columns 3 and 4, without and with *SDO/HMI* line of sight (LOS) contours of ±500 G, and *SDO* 1600 and 94 Å intensity images are displayed in the last two columns. The saturation levels are fixed for each column of intensity images

to discern the relative changes in brightness with time for a given wavelength or passband.

The first EIS observation at 12:38 UT coincides with the flare peak. At that time, the AR contains plasma that is predominantly of coronal composition. FIP bias ranges from 0.25 to 0.40 (using the ratio of high-FIP element Ar XIV over low-FIP element Ca XIV but equivalent to 2.5–4 in the solar FIP bias convention as stated in the previous section). The highest FIP bias (~0.25) is located in the N–S directed bright loops evident in the EIS Ar XIV, Ca XIV, and AIA 94 Å intensity images. Nine minutes later at 12:47 UT, the extent of the strong FIP bias material has spread as the bright loops expand during the confined eruption. This is also when the first highly localized patches of IFIP plasma composition are visible in the EIS Ar XIV/Ca XIV ratio images. Both IFIP regions are bordered by a “moat” of photospheric composition, which in turn is surrounded by coronal plasma composition. Throughout the remainder of the decay phase of FL1, the western IFIP patch is still clearly visible from 12:56 UT to 13:23 UT; however, the eastern patch evolves from IFIP to photospheric composition (see Ar XIV/Ca XIV intensity ratio images at 13:14 UT and 13:23 UT in column 3 of Figure 4). IFIP values are in the range [1.5, 2.1]. The IFIP plasma appears within the 9 minute time cadence of the EIS rasters and then rapidly evolves toward photospheric composition within 18 minutes during the decay phase of the M2.2 flare.

Figure 4, column 4, shows that the IFIP patches are located at the footpoints of the bright loops connecting the strong magnetic field concentrations of the northern bipole. We have overlaid green contours of the strong IFIP regions observed at 12:47 UT on the *SDO/HMI* continuum image at the corresponding time in Figure 6 to emphasize the very specific position of the anomalous composition within the strong field. The IFIP plasma is located entirely within the umbrae of the northern bipole. Furthermore, the flare ribbons of FL1 also cross/end at the same places as the IFIP patches over the umbrae (see the intersection of the orange contours of the flare ribbons and the green contours at the umbrae associated with ribbon R1 in the NE and R2 along the horizontal PIL). In summary, IFIP composition is found over umbral areas crossed by flare ribbons. The footpoints of flux rope MFR1 of Chintzoglou et al. (2015) are rooted in the same umbrae of the emerging bipole where we observe the IFIP plasma.

## 3.3. Photospheric Plasma Composition of Flare FL2

Figure 5 contains the same series of images for flare FL2 as is used for FL1 in Figure 4. Once again, the EIS observations began after the flare trigger was initiated and continued throughout and beyond the short decay phase of the flare. The locations of the IFIP plasma within AR 11429 remained fixed throughout the decay phase of FL1 and appear to correspond to where photospheric-like composition is observed ~8 hr later during flare FL2 (see the orange patches in the EIS Ar XIV/Ca XIV intensity ratio maps in Figure 5). Very little composition evolution is evident during FL2. The anomalous composition is no longer found to be in distinct IFIP patches; rather it has near-photospheric composition and is more dispersed especially in the eastern region. In the western region, much of the plasma over the positive polarity has evolved back toward coronal composition. The umbra where the eastern footpoint of MFR1 is rooted remains a coherent structure during and after FL2, though it has stopped coalescing, whereas, the umbra at the other (western
footpoint has broken up and dispersed (see the continuum images at 12:34 and 23:58 of Figure 3 and see the location indicated by the right arrow in each arrow-pair of Figure 3's animation).

# 4. Discussion and Interpretation

## 4.1. Summary of Plasma Composition Evolution in the Confined Flares

In this study, we have analyzed the evolution of plasma composition in the highly active AR 11429 observed by Hinode/EIS during two confined M-class flares. During the first confined flare, FIP bias plasma comprises the AR when the flare is at peak intensity. Minutes into the flare's decay phase, IFIP bias plasma appears at the footpoints of bright flare loops within the AR. The IFIP patches are surrounded by a ring of photospheric material. Co-temporally, the spatial extent of FIP bias increases in the bright flare loops. By the time the GOES soft X-ray intensity has returned to preflare levels, there is asymmetric composition evolution in the flare loop footpoints; the plasma at the eastern footpoints evolves to photospheric composition whereas at the western footpoints, it maintains IFIP composition encircled by photospheric plasma composition. Hinode/EIS observes these extremes in plasma evolution in less than one hour, primarily within the decay phase of the flare. Approximately seven hours later the second confined flare occurs and photospheric plasma composition is present in the vicinity of the flare loop footpoints where IFIP plasma was observed during the first flare.

Distinct IFIP patches occurred near the footpoints of one of the two magnetic flux ropes identified by Chintzoglou et al. (2015) where flare ribbons cross the umbrae of the emerging bipole in the northern section of the AR. These are very particular locations within the unusually complex magnetic configuration of the AR. The fact that the IFIP plasma is only observed for a short time during the decay phase of a moderate flare in highly localized places within the AR's magnetic field raises key questions as to what roles, if any, the AR's magnetic field configuration and flaring activity have in the creation and observation of anomalous plasma composition.

## 4.2. What is the Significance of the Emergence of Highly Sheared Magnetic Field and Coalescing Sunspots?

Such field represents different strands of highly sheared field that are converging toward each other to form sunspots, and therefore meet below the photosphere/chromosphere in the location of the coalescing umbrae. The high shear in coalescing strands of the same magnetic polarity suggests that the strands with a nonzero component of anti-parallel magnetic field have the possibility for subchromospheric magnetic reconnection when brought together. In the high-$\beta$ plasma regime of the photosphere/low chromosphere, upward-moving acoustic waves are generated which can mode convert to fast mode waves at $\sim 1$ Mm above the photosphere in the chromosphere where plasma-$\beta$ is equal to unity (Bourdin 2017). Over regions of high magnetic field concentration such as sunspots, the transition to a low-$\beta$ plasma occurs at lower heights within the photosphere. The atmospheric models of Avrett et al. (2015), constrained by observations, show a plasma pressure lower by a factor of 20–70 in a sunspot compared to the quiet Sun. This quiet-Sun pressure model is close to the one used by Gary (2001) to derive the plasma-$\beta$ versus height. With a field strength of 2500 G, Gary (2001) found plasma-$\beta$ to be around or below 0.2 within the photosphere. With the new sunspot atmospheric model of Avrett et al. (2015), the plasma-$\beta$ is a

factor between 20 and 70 lower, setting the sunspot umbra well in the low plasma-$\beta$ regime. Furthermore, the umbral magnetic field expands less with height than does the penumbral field. Consequently in the decreasing density, the Alfvén and fast mode speeds increase faster in the umbra, leading to the strongest refraction of fast mode waves, and the strongest downward ponderomotive acceleration.

This is significant in the context of the ponderomotive force fractionation model (Laming 2012, 2015, 2017). When Alfvén wave flux originates in the corona, the ponderomotive force points upward, bringing low-FIP elements up from the chromosphere so that the FIP effect is observed. Fast mode waves coming from below the chromosphere means the ponderomotive force is directed downward thereby depleting low-FIP elements from chromospheric plasma. This is consistent with Brooks (2018) who found that depletion of low-FIP elements instead of enhanced abundances of high-FIP elements yields IFIP plasma in post-flare loops.

In the Laming model, Ar/Ca is fractionated between just above the plasma-$\beta$ = 1 layer, and below a height of $\sim 1$ Mm where H starts to be ionized. The plasma-$\beta$ = 1 layer with 300 G is in the photosphere. The deeper this layer lies, the more likely acoustic waves generated by sub-photospheric reconnection will mode convert to magnetoacoustic or fast mode waves when reflected/refracted at high density gradients while H is still neutral, and causes the IFIP fractionation. In the case of AR 11429, the umbral magnetic field exceeds 500 G, therefore it is plausible that the plasma-$\beta$ = 1 layer could be low down in the photosphere or even below that as discussed above, enabling the IFIP fractionation to take place.

Chromospheric dynamics generally occur on timescales much faster than those for ionization and recombination, and so the chromospheric ionization balance is almost static, although elevated from that expected in equilibrium (Carlsson & Stein 2002). This is reproduced in the Avrett & Loeser (2008), Avrett et al. (2015) models. A more important concern resulting from chromospheric dynamics would be the effect on the wave physics, especially in the low chromosphere where the IFIP fractionation occurs. We speculate that the extra wave interactions with density structures would increase the reflectivity of the chromosphere to Alfvén and fast mode waves, thus reinforcing our conclusions. Extra dynamics in the chromosphere are probably required for IFIP, since it is worth noting that waves generated from solar or stellar convection are not strong enough, by an order of magnitude in amplitude, or two orders of magnitude in energy, to cause sufficient IFIP fractionation. Typical turbulent amplitudes are of order $1 \text{ km s}^{-1}$ (Bruntt et al. 2010) and relatively constant with stellar spectral type (Chaplin et al. 2009; Baudin et al. 2011; Kjeldsen & Bedding 2011), while around $10 \text{ km s}^{-1}$ is necessary to produce sufficient ponderomotive acceleration (Laming 2015).

The inference in this paper that subsurface reconnection generates the waves responsible for the IFIP fractionation is supported by surveys of FIP and IFIP in stars. Wood et al. (2018) find FIP fractionation reducing and becoming IFIP fractionation in stellar coronae as the stellar spectral type becomes later (i.e., cooler). Interestingly, the IFIP fractionation appears for stars where the magnetic field saturates when plotted against Rossby number, which is defined as the ratio of a star's rotational period to its convective turnover time (Reiners et al. 2009). Testa et al. (2015) extend this plot to X-ray emission and FIP/IFIP fractionation against Rossby number. The saturation of magnetic
field (or equivalently X-ray emission) implies the magnetic field generated by the rotation is being quenched, presumably by reconnection, which must presumably be subsurface reconnection because the X-ray emission also saturates. And, we emphasize, this subsurface reconnection coincides with the IFIP fractionation appearing in the coronae of these stars.

In the case of the pair of flare ribbons R3–R4 linked to flux rope MFR2 identified by Chintzoglou et al. (2015), no anomalous plasma composition was observed at any time during either flare. MFR2 has footpoints in the quiet Sun (Chintzoglou et al. 2015) and there were no coalescing umbrae observed in these locations (see the location of the footpoints of flux rope 2 in the continuum, radial magnetic field images at 12:34 UT on March 6 in Figure 3 and the corresponding locations in its animation). Therefore, the conditions for subchromospheric reconnection and the resulting Alfvénic wave generation are not satisfied at least at the locations of the flare ribbons. A comparison of the conditions associated with the two flux ropes highlights that the IFIP patches are found in very specific locations within the AR.

### 4.3. Does Flaring Play a Role in Creating IFIP Bias Plasma?

During a solar flare, magnetic reconnection high up in the corona accelerates electron beams and initiates Alfvén waves downward along magnetic loops connected to the chromosphere (e.g., Fletcher et al. 2011; Laming 2017; Reep et al. 2018). In the very dense chromosphere, kinetic energy is converted to thermal energy, causing the plasma to be heated up to temperatures of $log_{10} T \sim 7.0$ K (10 MK). The overdense, heated plasma expands upward, and due to the high overpressure, fills the coronal loops driving mass flows from the chromosphere into the corona i.e., the process of chromospheric evaporation (Neupert 1968).

In most flares, evaporation leads to plasma of photospheric or near-photospheric composition as both low- and high-FIP elements are ionized and evaporated together (e.g., Warren 2014 and other references given in Section 1). The unfractionated plasma arises during flares because the upward flow speed is too fast for sufficient fractionation to occur (Laming 2009). However, in FL1 we have chromospheric plasma that is depleted of low-FIP elements located above coalescing umbrae in the highly sheared emerging field. The atomic masses of Argon and Calcium are similar, 39.948 amu and 40.078 amu, respectively, so that there is no preferential evaporation upflow of one element over another with the Ar XIV/Ca XIV intensity ratio composition diagnostic used in this study. Chromospheric evaporation brings up what is there, and in the very specific locations within AR 11429, this is low-FIP depleted plasma, consequently we observe IFIP patches in the corona.

The flare ribbons cover a variety of magnetic field strengths from quiet-Sun to umbral (see the location of flare ribbon contours in Figure 2 and radial magnetograms at approximately the same time in Figure 3's animation). At those locations away from the strong field, the Hinode/EIS observations show that the flare energy releases low-FIP-biased plasma first, creating strong low-FIP composition in flare loops as observed in FL1 and to a much lesser extent in FL2. In Figure 4, the loops are filled with low-FIP enhanced plasma at 12:38 UT, at the time of the flare peak, and continue to be filled with the low-FIP biased plasma from 12:47 UT, when we first observe IFIP-biased plasma at the strong umbral field locations, until 13:05 UT.

Over the coalescing umbrae, where the composition deeper in the chromosphere is depleted of low-FIP elements prior to

the flare, the bottom of flare loops will be filled with plasma of IFIP composition. The IFIP patches appear to last as long as either chromospheric evaporation is triggered by ionizing electron beams produced by magnetic field reconnection in the corona and/or thermal conduction of coronal plasma persists (Bagalá et al. 1995; Cheung et al. 2018). The only observations of high energy electrons of FL1 are from the Fermi Gamma-ray Space Telescope. The peak in the 25–50 keV energy bin occurs at the peak of the flare at 12:38 UT before returning to background levels by ~12:45 UT (not shown here). As Fermi then entered its night phase, we do not have observations beyond 12:45 UT, however, the fact that the electron energy levels have returned to background levels by this time suggests that chromospheric evaporation is not driven by electron beams but rather by thermal conduction during the later phase of the flare. In either case, once the process of chromospheric evaporation is finished, with radiative and conductive cooling taking place, the IFIP plasma would no longer be observable using the Ar XIV/Ca XIV diagnostic ratio. Therefore, the patches appear to decay within the gradual phase of the flare.

The EIS observations indicate that IFIP composition plasma is present in the vicinity of footpoints in certain flare loops. The flaring reconnected loops inherit the composition of their progenitor loops, i.e., FIP composition. In addition, when the top of the chromosphere is being evaporated at the start of the flare, low-FIP-element enhanced plasma created during intermittent heating episodes enters the reconnected loops. Only later, when bremsstrahlung heating reaches deeper chromospheric layers, will the IFIP-biased plasma be evaporated. As IFIP plasma is injected into the flare loops, chromospheric evaporation creates the strongest IFIP bias at loop footpoints. At greater heights, plasma mixing creates a transition from IFIP bias through photospheric to FIP bias composition. A composition gradient from photospheric to coronal (FIP effect) along flare loops was observed by Doschek et al. (2018) in an AR at the limb, but in this case there was no IFIP plasma at the loop footpoints. In a quiescent newly emerged AR, Baker et al. (2013) also observed a composition gradient from the AR's loop footpoints to greater heights along the loops.

Our scenario for the role of flaring activity in the creation and observation of IFIP plasma on the Sun is similar to the interpretation presented in Laming (2009) for a flare on the M dwarf star EV Lac which has an IFIP quiescent corona. Element abundances during a moderate flare showed a near stellar photospheric composition. Laming (2009) argue that the downward-directed ponderomotive force in the chromosphere of the IFIP dominated star increases the heat conduction from the flare, consequently enhancing chromospheric evaporation and for active stars like EV Lac, stronger chromospheric evaporation leads to element abundance variation during flares. However, in the case of the less active Sun, comparatively reduced evaporation leads to little variation in the abundances observed during flares, at least on large scales.

Furthermore, in a multi-wavelength study of a giant flare observed on active M dwarf CN Leo with an IFIP dominated corona, Liefke et al. (2010) found more than a twofold increase in the low-FIP Fe abundance during the flare i.e., the FIP effect, compared to the preflare quiescent abundance level. The peak in Fe abundance enhancement approximately coincided with the flare's peak in soft X-ray emission. By the end of the decay phase, the Fe abundance returned to the preflare level of the quiescent corona, similar to the progression of plasma
composition we observed in the flare loops during FL1 in AR 11429. In each case, the plasma appears to approach photospheric composition and then returns to preflare quiescent coronal composition during the time period of flare decay phase.

### 4.4. Why do We Observe Photospheric Composition at the Locations of the IFIP Patches During the Second Flare?

The Ar XIV/Ca XIV ratio maps in Figure 5 contain photospheric plasma where IFIP patches had been observed during FL1. As evidenced in the animation of Figure 3 and described in Section 2.4, by the time of FL2 at ~21:00 UT, the sunspot coalescence was completed and the umbrae had entered their decay phase, halting subchromospheric reconnection. Therefore, the fast mode wave flux arriving from below the chromosphere would have ceased, together with the resulting ponderomotive force, which had depleted low-FIP elements from these locations. In effect, the IFIP mechanism was “turned off.” Once the IFIP plasma is no longer generated, plasma mixing with the surrounding field containing low-FIP bias material would begin to change what we observe in the corona in these specific locations; from IFIP to photospheric composition in flaring loops. Repeated flaring activity after FL1 may have provided a steady supply of unfractionated plasma to the vicinity of IFIP patches, accelerating the transition from IFIP to photospheric composition though this time period was not observed by Hinode/EIS. The pattern of evolution is partly consistent with what Nordon & Behar (2007, 2008) observed on much larger scales in their study; flares tend to decrease both IFIP bias in the IFIP dominated coronae of active stars and FIP bias in the coronae of solar-like stars (Testa 2010).

## 5. Conclusion

In this case study, highly localized inverse-FIP composition patches are appearing during a confined flare in an overall FIP bias dominated AR. These patches evolve and fade during the decay phase of the flare. The IFIP patches are observed in highly sheared emerging flux over coalescing umbrae crossed by flare ribbons. We propose that subsurface magnetic reconnection between coalescing umbrae led to an increased fast mode wave flux from below the fractionation height and resulted in the depletion of low-FIP elements. When these coalescing umbrae with plasma composition depleted of low-FIP elements became footpoints of flare loops (i.e., they were crossed by flare ribbons), the chromospheric evaporation of the low-FIP depleted plasma led to the appearance of IFIP patches above these umbrae. The IFIP patches were observed in the corona as long as chromospheric evaporation lasted. The insight gained from spatially resolved composition maps suggests that in Sun-as-a-star or stellar cases the composition of the chromosphere which is evaporated into flare loops may have local anomalies and a flare's effect on the overall coronal composition may depend on the filling factor of those anomalous locations. These findings are consistent with the Laming model as well as Brooks (2018). In future work we will analyze the magnetic field evolution of the eight ARs where IFIP plasma has been observed to determine how common the characteristics seen in this case study are when considering the entire sample.

We thank the referee for their constructive comments and suggestions that have improved the manuscript. The authors are very grateful to Konkoly Observatory, Budapest, Hungary, for

hosting two workshops on Elemental Composition in Solar and Stellar Atmospheres (IFIPWS-1, 2017 February 13–15 and IFIPWS-2, 2018 February 27 to 2018 March 1) and acknowledge the financial support from the Hungarian Academy of Sciences under grant NKSZ 2018_2. The workshops have fostered collaboration by exploiting synergies in solar and stellar magnetic activity studies and exchanging experience and knowledge in both research fields. Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with NAOJ as a domestic partner, and NASA and STFC (UK) as international partners. Scientific operation of Hinode is performed by the Hinode science team organized at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner countries. Support for the post-launch operation is provided by JAXA and NAOJ (Japan), STFC (UK), NASA, ESA, and NSC (Norway). SDO data were obtained courtesy of NASA/SDO and the AIA and HMI science teams. D.B. is funded under STFC consolidated grant No. ST/N000722/1. L.v. D.G. is partially funded under STFC consolidated grant No. ST/N000722/1. The work of D.H.B. was performed under contract to the Naval Research Laboratory and was funded by the NASA Hinode program. A.W.J., L.M.G., G.V., and L.v.D.G. acknowledge the support of the Leverhulme Trust Research Project Grant 2014-051 and the Royal Society. J.M.L. was supported by the NASA HSR (NNH16AC391) and LARS (NNH17AE601) programs, the Chandra GO program and by Basic Research Funds of the Chief of Naval Research. D.M.L. acknowledges support from the European Commission's H2020 Program under the following Grant Agreements: GREST (No. 653982) and Pre-EST (No. 739500) as well as support from the Leverhulme Trust for an Early-Career Fellowship (ECF-2014-792) and is grateful to the Science Technology and Facilities Council for the award of an Ernest Rutherford Fellowship (ST/R003246/1).

## Appendix

We address specific technical aspects of the Ar XIV/Ca XIV intensity ratio diagnostic described in Section 3.1. In Figure 7 we show the effect of temperature on the ratio for temperatures exceeding $log_{10} T = 6.8$ K. Figure 8 features the Gaussian

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Intensity ratios of Ar XIV and Ca XIV contribution functions vs. electron temperature for electron densities of $log_{10} N = 9.0, 10.0, 11.0$. Shaded box shows the range of the Ar XIV/Ca XIV ratio in the temperature range of [6.2, 6.8] for electron density of $log_{10} N = 10.0$ and coronal abundances as given in Section 3.1. The ratio is quite narrow for a wide temperature range but exceeds 0.50 for temperatures above $log_{10} T = 6.9$ K.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Lines in the Ca XIV spectral window for the *Hinode* /EIS raster at 12:56 UT on March 6. The Y-axis is intensity. Three Gaussians were fit to the spectral region around Ca XIV 193.87 Å line (red), one for the Fe X 193.72 Å (yellow) in the blue wing and the other for the Ni XVI–Ar XI blend (green) in the red wing. Three Gaussians were fit to the Ar XIV 194.40 Å (magenta), two of which lie in the blue wing and are unidentified lines (blue and orange; Brown et al. 2008). Top panel: FIP effect spectrum for pixel = [23, 196]. The Ca XIV 193.87 Å line is much more intense compared to the Ar XIV 194.40 Å line and FIP bias is ~0.40. Bottom panel: IFIP effect spectrum at pixel = [024, 216]. The Ar XIV line is more intense than the Ca XIV line and the IFIP bias is ~1.60. The unidentified lines in the blue wing of the Ar XIV line are visible in FIP effect spectra but are less evident in the IFIP spectra.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Log normalized contribution functions, $G(T)$, of Ar XIV and Ca XIV (smoothly peaked curves) and $G(T)$ convolved with DEM (sharply peaked curves) for the IFIP patch discussed in the last paragraph of Section 3.1.</figcaption></figure>

fittings to the Ca XIV spectral window for sample FIP and IFIP effect spectra. Finally, Figure 9 illustrates the contribution functions convolved with the DEM.

## ORCID iDs

Deborah Baker https://orcid.org/0000-0002-0665-2355

Lidia van Driel-Gesztelyi https://orcid.org/0000-0002-2943-5978

David H. Brooks https://orcid.org/0000-0002-2189-9313

Gherardo Valori https://orcid.org/0000-0001-7809-0067

Alexander W. James https://orcid.org/0000-0001-7927-9291

J. Martin Laming https://orcid.org/0000-0002-3362-7040

David M. Long https://orcid.org/0000-0003-3137-0277

Pascal Démoulin https://orcid.org/0000-0001-8215-6532

Lucie M. Green https://orcid.org/0000-0002-0053-4876

Sarah A. Matthews https://orcid.org/0000-0001-9346-8179

Katalin Oláh https://orcid.org/0000-0001-9570-3558

Zsolt Kővári https://orcid.org/0000-0001-5160-307X

## References

Argiroffi, C., Maggio, A., Peres, G., Stelzer, B., & Neuhäuser, R. 2005, A&A, 439, 1149

Audard, M., Güdel, M., & Mewe, R. 2001, A&A, 365, L318

Audard, M., Güdel, M., Sres, A., Raassen, A. J. J., & Mewe, R. 2003, A&A, 398, 1137

Avrett, E., Tian, H., Landi, E., Curdt, W., & Wülser, J.-P. 2015, ApJ, 811, 87

Avrett, E. H., & Loeser, R. 2008, ApJS, 175, 229

Bagalá, L. G., Mandrini, C. H., Rovira, M. G., Démoulin, P., & Hénoux, J. C. 1995, SoPh, 161, 103

Baker, D., Brooks, D. H., Démoulin, P., et al. 2013, ApJ, 778, 69

Baker, D., Brooks, D. H., Démoulin, P., et al. 2015, ApJ, 802, 104

Baker, D., Brooks, D. H., van Driel-Gesztelyi, L., et al. 2018, ApJ, 856, 71

Baudin, F., Barban, C., Belkacem, K., et al. 2011, A&A, 535, C1

Bentley, R. D., Sylwester, J., & Lemen, J. R. 1997, AdSpR, 20, 2275

Bourdin, P.-A. 2017, ApJL, 850, L29

Brinkman, A. C., Behar, E., Güdel, M., et al. 2001, A&A, 365, L324

Brooks, D. H. 2018, ApJ, 863, 140

Brooks, D. H., Baker, D., van Driel-Gesztelyi, L., & Warren, H. P. 2017, NatCo, 8, 183

Brooks, D. H., Baker, D., van Driel-Gesztelyi, L., & Warren, H. P. 2018, ApJ, 861, 42

Brooks, D. H., Ugarte-Urra, I., & Warren, H. P. 2015, NatCo, 6, 5947

Brooks, D. H., & Warren, H. P. 2011, ApJL, 727, L13

Brooks, D. H., & Warren, H. P. 2012, ApJL, 760, L5

Brown, C. M., Feldman, U., Seely, J. F., Korendyke, C. M., & Hara, H. 2008, ApJS, 176, 511

Bruntt, H., Bedding, T. R., Quirion, P.-O., et al. 2010, MNRAS, 405, 1907

Caffau, E., Ludwig, H.-G., Steffen, M., Freytag, B., & Bonifacio, P. 2011, SoPh, 268, 255

Carlsson, M., & Stein, R. F. 2002, ApJ, 572, 626

Chaplin, W. J., Houdek, G., Karoff, C., Elsworth, Y., & New, R. 2009, A&A, 500, L21

Cheung, M. C. M., Rempel, M., Chintzoglou, G., et al. 2018, NatAs, 3, 160

Chintzoglou, G., Patsourakos, S., & Vourlidas, A. 2015, ApJ, 809, 34

Culhane, J. L., Harra, L. K., James, A. M., et al. 2007, SoPh, 243, 19

Del Zanna, G. 2013, A&A, 555, A47

Del Zanna, G., Dere, K. P., Young, P. R., Landi, E., & Mason, H. E. 2015, A&A, 582, A56

Del Zanna, G., & Mason, H. E. 2014, A&A, 565, A14

Del Zanna, G., & Mason, H. E. 2018, LRSP, 15, 5

Del Zanna, G., & Woods, T. N. 2013, A&A, 555, A59

Dennis, B. R., Phillips, K. J. H., Schwartz, R. A., et al. 2015, ApJ, 803, 67

Dere, K. P., Landi, E., Mason, H. E., Monsignori Fossi, B. C., & Young, P. R. 1997, A&AS, 125, 149

Donea, A., & Hanson, C. 2013, Journal of Physics Conference Series, 440, 012028

Doschek, G. A., Feldman, U., & Seely, J. F. 1985, MNRAS, 217, 317

Doschek, G. A., Laming, J. M., Feldman, U., et al. 1998, ApJ, 504, 573

Doschek, G. A., & Warren, H. P. 2016, ApJ, 825, 36
Doschek, G. A., & Warren, H. P. 2017, ApJ, 844, 52

Doschek, G. A., Warren, H. P., & Feldman, U. 2015, ApJL, 808, L7

Doschek, G. A., Warren, H. P., Harra, L. K., et al. 2018, ApJ, 853, 178

Feldman, U. 1992, PhyS, 46, 202

Feldman, U., Schühle, U., Widing, K. G., & Laming, J. M. 1998, ApJ, 505, 999

Feldman, U., Warren, H. P., Brown, C. M., & Doschek, G. A. 2009, ApJ, 695, 36

Feldman, U., & Widing, K. G. 1990, ApJ, 363, 292

Feldman, U., & Widing, K. G. 1993, ApJ, 414, 381

Feldman, U., & Widing, K. G. 2003, SSRv, 107, 665

Fletcher, L., Dennis, B. R., Hudson, H. S., et al. 2011, SSRv, 159, 19

Fludra, A., & Schmelz, J. T. 1999, A&A, 348, 286

Freeland, S. L., & Handy, B. N. 1998, SoPh, 182, 497

Gary, G. A. 2001, SoPh, 203, 71

Güdel, M., Linsky, J. L., Brown, A., & Nagase, F. 1999, ApJ, 511, 405

Hale, G. E., Ellerman, F., Nicholson, S. B., & Joy, A. H. 1919, ApJ, 49, 153

Kjeldsen, H., & Bedding, T. R. 2011, A&A, 529, L8

Ko, Y.-K., Young, P. R., Muglach, K., Warren, H. P., & Ugarte-Urra, I. 2016, ApJ, 826, 126

Laming, J. M. 2009, ApJ, 695, 954

Laming, J. M. 2012, ApJ, 744, 115

Laming, J. M. 2015, LRSP, 12, 2

Laming, J. M. 2017, ApJ, 844, 153

Laming, J. M., Drake, J. J., & Widing, K. G. 1995, ApJ, 443, 416

Laming, J. M., & Hwang, U. 2009, ApJL, 707, L60

Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, SoPh, 275, 17

Liefke, C., Fuhrmeister, B., & Schmitt, J. H. M. M. 2010, A&A, 514, A94

Liefke, C., Ness, J.-U., Schmitt, J. H. M. M., & Maggio, A. 2008, A&A, 491, 859

Lodders, K. 2008, ApJ, 674, 607

Luoni, M. L., Démoulin, P., Mandrini, C. H., & van Driel-Gesztelyi, L. 2011, SoPh, 270, 45

McKenzie, D. L., & Feldman, U. 1992, ApJ, 389, 764

Meyer, J.-P. 1985a, ApJS, 57, 151

Meyer, J.-P. 1985b, ApJS, 57, 173

Neupert, W. M. 1968, ApJL, 153, L59

Nordon, R., & Behar, E. 2007, A&A, 464, 309

Nordon, R., & Behar, E. 2008, A&A, 482, 639

Patsourakos, S., Georgoulis, M. K., Vourlidas, A., et al. 2016, ApJ, 817, 14

Peretz, U., Behar, E., & Drake, S. A. 2015, A&A, 577, A93

Petrie, G. J. D. 2012, ApJ, 759, 50

Phillips, K. J. H., & Dennis, B. R. 2012, ApJ, 748, 52

Polito, V., Del Zanna, G., Valori, G., et al. 2017, A&A, 601, A39

Pottasch, S. R. 1963, ApJ, 137, 945

Reep, J. W., Russell, A. J. B., Tarr, L. A., & Leake, J. E. 2018, ApJ, 853, 101

Reiners, A., Basri, G., & Browning, M. 2009, ApJ, 692, 538

Robrade, J., & Schmitt, J. H. M. M. 2005, A&A, 435, 1073

Robrade, J., & Schmitt, J. H. M. M. 2006, A&A, 449, 737

Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, SoPh, 275, 207

Schmelz, J. T., Reames, D. V., von Steiger, R., & Basu, S. 2012, ApJ, 755, 33

Schonfeld, S. J., White, S. M., Henney, C. J., Arge, C. N., & McAteer, R. T. J. 2015, ApJ, 808, 29

Sheeley, N. R., Jr. 1996, ApJ, 469, 423

Sterling, A. C., Doschek, G. A., & Feldman, U. 1993, ApJ, 404, 394

Sun, X., Bobra, M. G., Hoeksema, J. T., et al. 2015, ApJL, 804, L28

Sylwester, B., Phillips, K. J. H., Sylwester, J., & Kępa, A. 2015, ApJ, 805, 49

Syntelis, P., Gontikakis, C., Patsourakos, S., & Tsinganos, K. 2016, A&A, 588, A16

Takasao, S., Fan, Y., Cheung, M. C. M., & Shibata, K. 2015, ApJ, 813, 112

Telleschi, A., Güdel, M., Briggs, K., et al. 2005, ApJ, 622, 653

Testa, P. 2010, SSRv, 157, 37

Testa, P., Reale, F., Garcia-Alvarez, D., & Huenemoerder, D. P. 2007, ApJ, 663, 1232

Testa, P., Saar, S. H., & Drake, J. J. 2015, RSPTA, 373, 20140259

Veck, N. J., & Parkinson, J. H. 1981, MNRAS, 197, 41

Warren, H. P. 1999, SoPh, 190, 363

Warren, H. P. 2014, ApJL, 786, L2

Warren, H. P., Brooks, D. H., Ugarte-Urra, I., et al. 2018, ApJ, 854, 122

Warren, H. P., Ugarte-Urra, I., & Landi, E. 2014, ApJS, 213, 11

Widing, K. G., & Feldman, U. 1995, ApJ, 442, 446

Wood, B. E., Laming, J. M., & Karovska, M. 2012, ApJ, 753, 76

Wood, B. E., Laming, J. M., Warren, H. P., & Poppenhaeger, K. 2018, ApJ, 862, 66

Wood, B. E., & Linsky, J. L. 2006, ApJ, 643, 444

Wood, B. E., & Linsky, J. L. 2010, ApJ, 717, 1279