Differential Emission Measure Plasma Diagnostics of a Long-Lived Coronal Hole

We use Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) data to reconstruct the plasma properties from differential emission measure (DEM) analysis for a previously studied long-lived, low-latitude coronal hole (CH) over its lifetime of ten solar rotations. We initially obtain a non-isothermal DEM distribution with a dominant component centered around 0.9 MK and a secondary smaller component at 1.5 – 2.0 MK. We find that deconvolving the data with the instrument point spread function (PSF) to account for long-range scattered light reduces the secondary hot component. Using the 2012 Venus transit and a 2013 lunar eclipse to test the efficiency of this deconvolution, significant amounts of residual stray light are found for the occulted areas. Accounting for this stray light in the error budget of the different AIA filters further reduces the secondary hot emission, yielding CH DEM distributions that are close to isothermal with the main contribution centered around 0.9 MK. Based on these DEMs, we analyze the evolution of the emission measure (EM), density, and averaged temperature during the CH’s lifetime. We find that once the CH is clearly observed in EUV images, the bulk of the CH plasma reveals a quite constant state, i.e. temperature and density reveal no major changes, whereas the total CH area and the photospheric magnetic fine structure inside the CH show a distinct evolutionary pattern. These findings suggest that CH plasma properties are mostly “set” at the CH formation or/and that all CHs have similar plasma properties.


Introduction
Coronal holes are regions of reduced plasma density and temperature in the solar corona. Thus they appear darker than the surrounding quiet Sun in extreme-ultraviolet (EUV) and X-ray images. Magnetically, CHs are characterized by the prevailing dominant polarity nature of the underlying magnetic field where magnetic-field lines of the dominant polarity are "open" into interplanetary space (Wilcox, 1968;Altschuler, Trotter, and Orrall, 1972;Hundhausen, 1972). Along these open magnetic-field lines, plasma is accelerated, escaping the Sun's gravitational field forming interplanetary high-speed solar-wind streams. The solarwind parameters speed and proton temperature measured at 1 AU are correlated with the CH area (Krieger, Timothy, and Roelof, 1973;Vršnak, Temmer, and Veronig, 2007;Rotter et al., 2012).
These open flux tubes are believed to be rooted in concentrations of magnetic flux, known as magnetic elements (Hassler et al., 1999;Tu et al., 2005). Based on a study of 288 lowlatitude CHs, Hofmeister et al. (2017) found that these magnetic elements are also the main origin of the unbalanced magnetic flux of CHs. The total area covered by long-lived magnetic elements with lifetimes > 4 days is strongly correlated with the unbalanced magnetic flux of CHs (cc = 0.99). These magnetic elements account for ≈ 68% of the unbalanced "open" magnetic flux of CHs at typical coverages of about 3% of the total CH area Heinemann et al., 2019). Heinemann et al. (2018a,b) investigated a long-lived low-latitude CH observed by the Solar Dynamics Observatory (SDO) and Solar Terrestrial Relations Observatory (STEREO) spacecraft over ten solar rotations from February 2012 to October 2012 (see Figure 7). Based on the CH area, they identified a distinct three-phase evolution characterized by a growing, maximum, and decaying phase. For this CH, a strong correlation between CH area and the mean magnetic-field strength inside the CH was found. The signed magnetic-field strength increased over the course of the growing phase and reached its maximum at the time of maximal CH area followed by a steep drop in field strength marking the beginning of the CH's decaying phase. During the maximum phase, a larger fraction of the overall magnetic flux is concentrated in magnetic elements than in the growing and decaying phases, and the number of magnetic elements per area for elements with field strengths > 50 G increases nearly 100% compared to either the beginning or end of the CH's life. This maximum phase further exhibits the strongest correlation between CH area and peak speed of the associated high-speed solar-wind stream.
As the magnetic-field configuration is assumed to be responsible both for the coronal heating (e.g. Abramenko, Pevtsov, and Romano, 2006;Klimchuk, 2006) and for the acceleration of the solar-wind plasma (Wang, 2010), the clear trend in the magnetic structure and the distinct correlations with the fast solar-wind characteristics of this CH identified by Heinemann et al. (2018a,b) make it an ideal candidate for probing the effect of changing magnetic fine structure on CH plasma. The low-latitude position of the CH under study minimizes projection effects, thus making it also a suitable candidate for differential emission measure (DEM) analysis. This allows us to deduce the temperature-dependent emission of the optically thin coronal plasma by studying multi-wavelength EUV measurements in combination with elemental abundances and the wavelength dependent instrument response.
In this article we apply DEM diagnostics using the inversion code developed by Hannah and Kontar (2012) using SDO/Atmospheric Imaging Assembly (AIA) data to reconstruct the DEM distribution of a CH and infer the density and temperature evolution over its lifetime.

Figure 1
Temperature response functions for the utilized AIA EUV filters generated using CHIANTI 9 and the effective areas provided by aia_get_response.pro.

Data
We made use of six of the seven EUV channels (94, 131, 171, 193, 211, and 335 Å) of the AIA instrument (Lemen et al., 2012) onboard the SDO (Pesnell, Thompson, and Chamberlin, 2012) spacecraft, sensitive over a temperature range from 10 5 K to 10 7 K. The 304 Å filter was omitted because it images optically thick plasma making it unsuitable for DEM analysis, which relies on optically thin emission. With the exception of the He II imaging 304 Å channel, AIA's EUV filters primarily observe different Fe ions (with ionization states from Fe VIII to Fe XXIV).
The AIA temperature response functions of the filters utilized are given in Figure 1. The AIA pixel scale is 0.6 arcsec (Lemen et al., 2012). To enhance the signal-to-noise ratio, the data were binned by 8 × 8 pixels while conserving flux, resulting in an effective spatial resolution of 4.8 arcsec per pixel. The AIA data were deconvolved with the instrument point spread function (PSF) in order to reduce the instrument stray light and its effect on the reconstructed DEMs (Section 3). The data were further processed using the standard SSWIDL data reduction routine aia_prep.pro.
We then perform the DEM analysis at a cadence of one hour during the ± 6-hour time window where the CH is nearest to the disk center to reduce line-of-sight (LOS) integration effects of the emission of close-by structures. We applied the same CH boundaries as in Heinemann et al. (2018a,b) which were generated using an intensity-based threshold method applied to the AIA 193 Å images (Rotter et al., 2012;Hofmeister et al., 2017).

Differential Emission Measure Analysis
For optically thin emission from plasma in thermodynamic equilibrium the temperature distribution of the contributing plasma in the LOS is described by the differential emission measure defined as where n is the number density dependent on the temperature T along the LOS. The measured intensity I λ (T ) for a given AIA filter is then related to the DEM via where K λ (T ) denotes the response function of the corresponding filter, which depends on elemental abundances and the temperature of the emitting plasma as well as on the sensor sensitivity (Hannah and Kontar, 2012). We assumed photospheric element abundances, as up to the lower corona CHs show abundances close to photospheric values (Feldman, 1998;Feldman and Widing, 2003). For calculating the instrument response, abundances were taken from the CHIANTI 9 database (Dere et al., 1997(Dere et al., , 2019 and the filter response function from aia_get_response.pro in the SolarSoftware (SSWIDL) package. To infer the DEM from AIA data, we applied the regularized inversion technique developed by Hannah and Kontar (2012) to reconstruct the DEM curve for each binned pixel from the six coronal AIA EUV channels. The temperature range considered for DEM analysis was chosen from 0.2 to 5 MK. Integration of the DEM over the temperature range yields the total emission measure: Following Cheng et al. (2012) and Vanninathan et al. (2015), the mean temperature can be estimated by the emission-weighted temperature, Assuming a filling factor of unity, the plasma density can then be derived from Equation 1 as with h the column height of emitting plasma along the LOS. As an estimate for the column height of plasma contributing to the measured emission, we use the hydrostatic scale height, Using μ = 0.64 and m H = 1.67 × 10 −27 kg for a fully ionized solar plasma and assuming a typical CH plasma temperature around 0.9 MK, one arrives at a typical scale height h ≈ 42 Mm. Wendeln and Landi (2018) pointed out that any high-temperature component in the CH plasma distribution inferred by spectroscopic analysis is likely an artifact due to instrument scattered light. In order to reduce the effect of instrument stray light on the recovered DEM, we deconvolved the data with the instrument PSF provided by the aia_calc_psf.pro routine. This is in particular relevant for studying regions of low emission such as CHs. To estimate the effectiveness of the deconvolution in reducing stray light over low-emission regions, we analyze the intensity profiles over areas occulted by the Venus transit on 6 June 2012 and a lunar eclipse on 6 August 2013. Figure 2 (bottom panels) compares the intensities of data processed solely with aia_prep.pro (level 1.5) with the intensities of the deconvolved data (level 1.6) along white lines indicated in the images in the top panels. Over the areas shaded by Venus (left) or the Moon (right), deconvolution reduces the remaining intensity for most channels, but some DNs still remain owing to long-range scattered light due to the microroughness of the telescope mirrors (private communication, M. Cheung, 2019). The deconvolved data also show better contrast along the edge of the lunar eclipse.

Effect of PSF Deconvolution on Reconstructed DEMs
To quantify the significance of the remaining stray light for analyzing CHs, Table 1 shows the remaining DNs after deconvolving with the instrument PSF for typical Venus and lunar pixels compared with the mean DNs of the CH under study on 6 May 2012 for each AIA channel used. Comparing Venus and CH, for the 131 Å and 171 Å channels, the remaining intensity is small compared to typical CH intensities. For the 193 Å channel the remaining stray light is about 50% the intensity of typical CH DNs. For the 211 Å and 94 Å channels, the contribution from stray light is about as high as the measured intensity inside the CH. Since the 171 Å and 131 Å filters, which sample cooler plasma inside the CH (Figure 1), are less affected by stray light than the filters sensitive to hotter plasma (Table 1), we conclude that most of the high-temperature emission in the recovered DEM of CHs is likely due to contamination from regions outside the CH. Figure 3 shows the effect of deconvolving the data with the PSF and accounting for residual stray light on the recovered DEM of CHs. It compares the average DEM for an entire CH computed with the original, i.e. not-PSF-deconvolved level 1.5 data (red), with the DEM where the PSF was taken into account (blue), and with the DEM where we have additionally increased the presumed errors in each AIA channel due to long-range scattered stray light (green). Thereby, we have used the residual DNs of the Venus transit from Table 1 as an estimate on the long-range scattered light and add these to the error budget of the AIA channels given by the SSWIDL routine aia_bp_estimate_error.pro. We assume this to be a reasonable estimate, as the eclipse happened during the lifetime of the CH under study, so that the EUV emission from the solar disk and thus scattered light are comparable in these data sets. When taking the PSF into account, the high-temperature tail of the CH plasma is significantly reduced compared to the original data. Considering residual stray light as detailed above, the hot-temperature component even almost vanishes. To compute more accurate DEMs, we therefore account for both PSF and residual stray light.

General Features of the CH DEM Distribution
Over the ten solar rotations of the CH's lifetime, we derived the DEMs of each 8 × 8 binned pixel of the AIA images, when the CH was close to the central meridian. Figure 4 shows examples of the derived DEM distributions for CH (top panels) and quiet Sun (QS) pixels (bottom panels). The total EM for CH regions is about ten times smaller than of the QS with the shape of the DEM also differing. The DEM distribution for CH regions shows a maximum of the emission at a temperature of around 0.9 MK with a small secondary contribution from emission centered between 1.5 -2.0 MK. DEM from QS regions exhibits peak emission between 1.5 -2.0 MK with a secondary 0.9 -1.1 MK component being sometimes present. Also shown in Figure 4 are the results of fitting the DEMs with the sum of two Gaussian functions centered at those temperatures (solid red and blue lines). The finding that coronal plasma outside CHs also shows emission from a component centered around CH temperatures is consistent with the reasoning of Landi and Feldman (2008), who concluded that coronal plasma may be a superposition of CH, QS, and active-region plasma, with CHs missing the hotter components. We further investigated the dependence of the DEM on the distance from the CH boundaries. To this aim, we separated the CH and its immediate surroundings into several segments outward and inward of the CH boundary. Inside the CH, each segment has a width of 9.6 . For the first region we considered pixels that were just inside the determined CH boundary up to a distance of 9.6 . For the second we used pixels at least 9.6 and at most 19.2 from the boundary, and so on. For the regions outside the CH boundary shown in black in the upper panel of Figure 5, the widths are 19.2 , 19.2 , and 38.4 . We then took the averages of the DEM for all pixels inside those segments. The resulting averaged DEMs are shown in Figure 5. The segments show a sharp difference between regions inside (red) and outside (blue) the CH. The DEM for segments inside the determined CH boundary shows a distinct peak at T ≈ 0.9 MK and a small high-temperature tail of constant emission most prominent around T ≈ 1.5 -2.0 MK. For the segments outside the CH boundary, the total EM is strongly increasing and the contribution from T ≈ 1.5 -2.0 MK rises. However, we note that in the segment directly adjacent to the CH boundary, the T ≈ 0.9 MK peak is still dominant, and only outside > 38.4 from the CH boundary the emission peak at T ≈ 1.5 -2.0 MK contributes equally. We note that our assumption of photospheric elemental abundances is not fulfilled outside the CH region, which may affect the derived absolute values and the ratio of cool/hot emission.
The radial profiles for plasma density and DEM-weighted temperature (Equations 4, 5) derived from the averaged DEM curves are also shown in Figure 5. The plasma density decreases towards the interior of the CH and increases outside the CH boundary. The temperature inside the CH is slightly higher towards the CH center because as the total emission decreases, the constant hot tail becomes relatively more important when calculating the mean temperature from the distribution. As the temperature of the peak does not change and the tail appears to originate mostly from stray light, this increase in plasma temperature towards the CH center is likely an effect due to limitations of the available data. The temperature increases again for segments further outside the CH boundary. For reference, the density and temperature averaged over the QS region (utilizing the same photospheric abundances) indicated by the green rectangle in the top-left panel of Figure 5 are shown in green. As discussed previously, the favored interpretation of the non-vanishing hot (T ≈ 1.5 -2.0 MK) component inside the CH seems to be stray light caused by instrumental effects and/or contribution from coronal bright points inside the CH, which show a DEM profile with a main contribution of T > 1.5 MK. To reduce the influence of overlying structures and stray light from the surrounding QS, in the further analysis we limited the CH area used for the DEMs to pixels that are at least three pixels (14.4 ) inside any CH boundary, in addition to deconvolving the data with the instrument PSF and accounting for the residual counts in the error budget as described in Section 2. Using only pixels well inside the CH has the added benefit of excluding regions where the assumption of photospheric abundances may not be valid. Figure 7 shows these reduced CH masks overplotted in black alongside the initial CH masks derived from the image segmentation in white.
Next we investigate the emission measure integrated over different temperature bins to illustrate the described difference in emission between CH and QS regions. Figure 6 shows the EM maps in 0.2 MK bins over the temperature range T = 0.6 to 1.8 MK for the CH on 03 June 2012. In the cooler bins up to 1.2 MK, the CH does not appear as a distinct region of reduced emission when compared to its surroundings. In contrast, in the EM bins at T > 1.2 MK, the CH region clearly sticks out as a region of strongly reduced emission compared to the ambient corona. This finding holds for all solar rotations of the CH's lifetime (see Appendix A: Figures 11 -19).  Figure 7 shows snapshots of the CH under study in the AIA 193 Å filter for nine solar rotations from 15 February 2012 to 13 October 2012 illustrating the CH evolution with a growth of the CH's area followed by a period of maximal area and a decaying phase of diminishing area. CH boundaries used by Heinemann et al. (2018a,b) are indicated in white. Figure 8 shows maps of the mean density and mean ion temperature (Equations 4 and 5) during the CH evolution. To quantify the changes of the DEM during the CH's lifetime, we derive the mean of all of the DEM curves of pixels inside the CH at a given solar rotation, utilizing the smaller boundaries indicated in black in Figure 7. Figure 9 shows the derived mean DEMs together with the 1σ range of values for the nine solar rotations, shown in Figure 7. It is noted that the DEM in April 2012 is notably reduced. The reason for this steep reduction is unclear but it seems to be a real effect. We tested the DEM of the CH for different days during this rotation, and all of them showed the same behavior. Figure 10 shows the evolution of the total EM as well as density and averaged temperature (Equations 3 -5) together with the CH area and fraction of the CH area covered by magnetic elements > 50 G from Heinemann et al.

Discussion and Conclusion
Using the DEM code developed by Hannah and Kontar (2012), we reconstructed mean plasma temperature and density for a long-lived, low-latitude CH over its lifetime of ten solar rotations from the six coronal AIA EUV filters. We initially recovered a nonisothermal DEM distribution for CH regions (Figure 4) with the main peak centered around T = 0.9 MK and a secondary smaller component in the range 1.5 -2.0 MK. A similar finding was reported in previous studies such as Hahn, Landi, and Savin (2011) who performed DEM analysis of a polar CH using Hinode/Extreme-Ultraviolet Imaging Spectrometer (EIS) data and attributed the hot contribution from 1.5 to 2 MK to structures along the LOS, contamination from streamers or closed loops of higher temperature inside the CH itself. Wendeln and Landi (2018) studied Hinode/EIS observations of an equatorial CH ) in red together with the initial (white) and reduced (green) CH boundaries for all ten solar rotations where the CH was observed. and concluded that any emission from hot plasma in CHs may be the result of a previously underestimated contribution from instrumental stray light. A study of the off-disk thermal structure for a polar CH by Landi (2008) using spectra from the Solar and Heliospheric Observatory (SOHO)/Solar Ultraviolet Measurements of Emitted Radiation (SUMER) instrument recovered a near isothermal off-disk CH plasma with a height-dependent temperature  5.85 < log(T ) < 6.04 K, further supporting the assumption of a non-physical origin of the hot tail of the recovered CH DEMs.
We showed that deconvolving the data with the instrument PSF reduces the hot tail in the recovered DEM distribution of CHs (Figure 3). Using AIA observations of the 2012 Venus transit and a lunar eclipse, we verified the assumption that the hot tail is not actually emitted from the CH. We found that, for hot AIA channels, even after PSF deconvolution significant DNs remain for pixels in regions where the Sun is occulted (Figure 2, Table 1). We therefore used the residual counts obtained by analyzing the radial profile of the Venus transit as an additional contribution to the errors in the AIA measurements, thereby further reducing the hot tail in the derived CH DEM profiles (Figure 3).
EM maps clearly show the CH as a region of reduced intensity for temperatures > 1.2 MK, whereas for emission from cooler temperatures no such effect is noticeable (Figure 6, Appendix A). The radial profile of the density and temperature shows a steep decrease in the regions adjacent to the CH boundary, with relatively flat profiles inside the CH. While the density decreases away from the boundary, the temperature slightly increases towards the CH center due to a small, constant, high-temperature emission ( Figure 5). As the temperature of the DEM maximum remains constant, we interpret this as the result of instrumental effects that lead to the described constant higher-temperature tail.
We do not find a distinct trend linking CH area and EM/density or plasma temperature that would correspond to the reported three-phase evolution of the long-lived CH under study. But we find that once the CH is clearly observed as a region of reduced emission in the EUV images, it already has a significantly reduced density and temperature, which do not distinctly change over the CH's lifetime. Since the magnetic properties follow the same three-phase evolution as the CH area (Heinemann et al., 2018b), this indicates that the mean density and temperature evolution of the CH are not strongly governed by the changing magnetic properties and fine structure of the CH. While the magnetic field is believed to be responsible for heating coronal plasma (Abramenko, Pevtsov, and Romano, 2006), the bulk of the CH plasma appears to not be significantly altered by the changes in the photospheric magnetic fine structures of the dominant coronal-hole polarity. Our finding that the plasma properties in the CH do not evolve with the area and magnetic-field strength suggests that the CH plasma properties are mostly "set" at the CH formation and/or that all CHs have similar plasma properties. To verify these hypotheses, a statistical study covering a large set of well observed CHs is necessary.