Statistical Approach on Differential Emission Measure of Coronal Holes using the CATCH Catalog

Coronal holes are large-scale structures in the solar atmosphere that feature a reduced temperature and density in comparison to the surrounding quiet Sun and are usually associated with open magnetic fields. We perform a differential emission measure analysis on the 707 non-polar coronal holes in the Collection of Analysis Tools for Coronal Holes (CATCH) catalog to derive and statistically analyze their plasma properties (i.e. temperature, electron density, and emission measure). We use intensity filtergrams of the six coronal EUV filters from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory, which cover a temperature range from ≈105.5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\approx10^{5.5}$\end{document} to 107.5K\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$10^{7.5}~\mbox{K}$\end{document}. Correcting the data for stray and scattered light, we find that all coronal holes have very similar plasma properties with an average temperature of 0.94±0.18MK\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$0.94 \pm0.18~\mbox{MK}$\end{document}, a mean electron density of (2.4±0.7)×108cm−3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(2.4 \pm0.7) \times10^{8}~\mbox{cm}^{-3}$\end{document}, and a mean emission measure of (2.8±1.6)×1026cm−5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$(2.8 \pm1.6) \times10^{26}~\mbox{cm}^{-5}$\end{document}. The temperature distribution within the coronal holes was found to be largely uniform, whereas the electron density shows a 30 to 40% linear decrease from the boundary towards the inside of the coronal hole. At distances greater than 20″ (≈15Mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\approx15~\mbox{Mm}$\end{document}) from the nearest coronal hole boundary, the density also becomes statistically uniform. The coronal hole temperature may show a weak solar-cycle dependency, but no statistically significant correlation of plasma properties with solar-cycle variations could be determined throughout the observed period between 2010 and 2019.


Introduction
Coronal holes are large-scale magnetic structures that extend from the solar photosphere into interplanetary space and are characterized by their open-to-interplanetary-space magneticfield configuration. This distinct magnetic topology enables plasma to be accelerated to high speeds of up to 780 km s −1 and escape along the open field lines (Schwenn, 2006). Coronal holes are defined by a lower density and temperature in comparison to the surrounding corona and are thus observed as large-scale regions of reduced emission in X-ray and extreme ultraviolet (EUV) wavelengths (see reviews by Cranmer, 2002Cranmer, , 2009, and the references therein).
Studies using differential emission measure (DEM) techniques on spectroscopic data from the EUV imaging spectrometer onboard Hinode (Hinode/EIS: Hahn, Landi, and Savin, 2011) and the Solar Ultraviolet Measurements of Emitted Radiation on the Solar and Heliospheric Observatory (SOHO/SUMER: Landi, 2008) as well as on Atmospheric Imaging Assembly (AIA: Lemen et al., 2012) onboard the Solar Dynamics Observatory (SDO: Pesnell, Thompson, and Chamberlin, 2012) EUV filtergrams (Saqri et al., 2020) revealed that coronal holes have a peak in the emission at a temperature of T ≈ 0.9 MK. Coronal Diagnostic Spectrometer (SOHO/CDS) spectroscopy also suggests this (Fludra, Del Zanna, and Bromage, 1999). This is significantly lower than the temperature of the surrounding quiet corona, which is often estimated at temperatures around 1.4 -1.6 MK (Landi and Chiuderi Drago, 2008;Hahn, Landi, and Savin, 2011;Del Zanna, 2013;Mackovjak, Dzifčáková, and Dudík, 2014;Hahn and Savin, 2014). Wendeln and Landi (2018) and Saqri et al. (2020) used Hinode/EIS and SDO/AIA data to derive DEM profiles and revealed that besides the dominant contribution at around 0.9 MK a secondary peak around ≈ 1.4 -1.6 MK was present within the DEM of coronal holes. It was suggested that it is mostly due to the presence of stray light. Electron densities in coronal holes have been estimated ranging from 1.0 -2.5 × 10 8 cm −3 (Fludra, Del Zanna, and Bromage, 1999;Warren and Hassler, 1999;Hahn, Landi, and Savin, 2011;Saqri et al., 2020). Using coronagraphic white-light images, Guhathakurta and Holzer (1994) showed that the density in coronal holes varies in height above the solar surface but not over latitude.
The plasma properties (i.e. emission measure, temperature, and density) of coronal holes have been analyzed using different observations and methods; however, it has been done only in the scope of observational campaigns or case studies but not in the scope of a large statistical approach. By using long-term EUV observations from 2010 to 2019 covering nearly the full Solar Cycle 24, we are able to statistically investigate the distributions of the plasma properties for a large variety of coronal holes of different sizes, and whether these properties change over the solar cycle.
In this study, we performed DEM analysis using data from SDO/AIA on coronal holes of the extensive Collection of Analysis Tools for Coronal Holes (CATCH) catalog 1 (Heinemann et al., 2019). We derive the distribution of plasma properties and their dependence on the solar cycle, and we relate them to the primary coronal hole parameters such as area and magnetic-field density.

Dataset
For the presented statistical study, we used the CATCH catalog, which contains 707 observations of well-defined non-polar coronal holes extracted from SDO/AIA 193 Å filtergrams. The catalog contains the extracted boundaries and properties of the coronal holes such as area, intensity, signed and unsigned magnetic-field strength, and magnetic flux including uncertainty estimates. The coronal holes are distributed between latitudes of ± 60 • and cover nearly the full Solar Cycle 24 from 2010 to 2019. A description of the catalog has been given by Heinemann et al. (2019).

Data Processing
For the DEM analysis, the level 1.6 data (processed with aia_prep.pro and point-spreadfunction-corrected) of the six coronal channels of AIA/SDO (Lemen et al., 2012) was used, namely 94 Å (Fe XVIII), 131 Å (Fe VIII, XXI), 171 Å (Fe IX), 193 Å (Fe XII, XXIV), 211 Å (Fe XIV), and 335 Å (Fe XVI). To enable us to pre-process all coronal holes in a reasonable amount of time and to enhance the signal-to-noise ratio for the DEM analysis, the data were rebinned by a factor of eight to a plate scale of 4.8 per pixel.

Point Spread Function Correction
From previous DEM studies of coronal holes it is known that stray light and scattered light significantly affect the analysis, since the emission in coronal holes is much lower than the surrounding quiet Sun and active regions (Wendeln and Landi, 2018;Saqri et al., 2020). Therefore, we use the point spread function (PSF) correction available in the SolarSoftware package of the Interactive Data Language (SSW IDL) aia psf (aia_calc_psf.pro written by M. Weber, SAO), to remove contributions from bright sources as well as possible.
We tested the performance of the PSF corrections using lunar-eclipse observations. At the eclipse boundary, the counts should drop to zero and all measured counts are due to stray light and noise. A well-performing PSF correction should show such a behavior. To verify, we analyzed ≈ 20 lunar eclipses that occurred between 2010 and 2019 in all six wavelengths. Figure 1 shows for each AIA filter, superposed light profiles of level 1.6 data across the boundaries of multiple solar eclipses. We find that, for all channels, some counts remain, however, without an obvious correlation to the solar cycle or the mean intensity of the solar disk. The hot channels (94 Å, 131 Å, 335 Å) show no dependence on the distance from the transition, and we assume the remaining counts to be primarily isotropic noise. This seems especially true for the 94 Å filter, where the remaining eclipse counts may be of the order of the quiet-Sun counts. The 171 Å, 193 Å, and 211 Å channels display a weak dependence on the distance, which indicates the presence of some uncorrected long-range scattered light (Wendeln and Landi, 2018). We found no working solution to correct for this; however, we can account for the remaining counts in the DEM analysis by increasing the input error of the counts (DN) in each pixel (see Saqri et al., 2020). For each filter, we derive the remaining counts by averaging over the lunar-eclipse light profiles (as shown in Figure 1) at a distance between 50 to 200 from the eclipse boundary. These distances were chosen to avoid contribution of the transition and associated boundary effects in the eclipse data (50 ) and because there are no coronal hole pixels with distances exceeding 200 . The average eclipse counts for level 1.5 and level 1.6 are listed in Table 1. We find Figure 1 Superposed intensity profiles across the solar disk during ≈ 20 lunar eclipses of level 1.6 data for different wavelengths. The eclipse border is centered on 0, the negative direction marks the visible solar disk and the positive direction the part covered by the lunar eclipse. Table 1 Mean lunar-eclipse counts (DN) of level 1.5 and level 1.6 data for all wavelengths derived from the light profiles in Figure 1, the resulting eclipse correction, as well as the mean counts of the 707 coronal holes using the level 1.6 data. that the level 1.6 data show strongly reduced remaining eclipse counts in the 171 Å, 193 Å, and 211 Å channels. The other channels are primarily dominated by isotropic noise and are comparable between level 1.5 and 1.6 data. The eclipse correction for the DEM input is given as the mean remaining eclipse counts (as stated above and shown in Table 1) plus one sigma. Note, that with this method we ignore the distance dependence and possible effects such as the roughly one-sided illumination of the eclipse in contrast to coronal holes in the center of the disk. This is done because we cannot reliably make an estimate for such effects.

Differential Emission Measure
Differential emission measure analysis is the reconstruction of the plasma properties from observed intensities in different wavelengths. The DEM is defined as the emission of optically thin plasma in thermodynamic equilibrium for a specific temperature along the lineof-sight (LoS) and is given by with n e the electron number density as a function of the temperature T and h the LoS distance over which the emission observed is integrated (see Mariska, 1992, Chapter 4). According to Hannah and Kontar (2012), from the observed intensities of each filter [I λ ] the DEM can be estimated by solving the inverse problem with the emission measure EM = DEM dT and the LoS integration length [h] which we approximate with the hydrostatic scale height h = k B T mgμ . Due to the open magnetic field that does not vertically constrain the plasma, we assume the hydrostatic scale height to be valid as a first-order approximation. However, this is only valid for the regime where the field is mostly vertical and approximately uniform, i.e. it is not valid outside of coronal holes. In quiet-Sun regions, a height-dependent DEM model using a scale-height approximation modeled with an ensemble of multi-hydrostatic loops might be used (Aschwanden, 2005). We use g = 274 m s −2 , a proton mass of m = 1.67 × 10 −27 kg (Aschwanden, 2005), μ = 0.60 (Asplund et al., 2009), and the median DEM temperature (Equation 4) for the hydrostatic scale height. The input error of each pixel was calculated using aia_bp_estimate_error.pro considering shot-, dark-, read-, quantum-, compression-, and calibration noise to which the eclipse correction is added (see Table 1).
To investigate the plasma properties of the coronal holes, we applied a regularized inversion technique developed by Hannah and Kontar (2012) to reconstruct the DEM from the six optically thin EUV channels of SDO/AIA. As Equation 2 does not yield a unique solution without further constraints, the code by Hannah and Kontar (2012) gives the DEM solution with the smallest amount of plasma required to explain the observed emission (zerothorder constraint). The instrument response function [K λ (T )] was calculated assuming photospheric abundances (CHIANTI 9 database: Dere et al. 1997Dere et al. , 2019 and the AIA filterresponse function available via SSW IDL (aia_get_response.pro). We prefer photospheric over coronal abundances as it was shown that the elemental abundances in chromospheric and coronal layers of coronal holes strongly resemble the photospheric ones (Feldman, 1998;Feldman and Widing, 2003). For every pixel in each coronal hole we use 60 equally spaced temperature bins (in log 10 space) between log 10 (T ) = 5.3 and log 10 (T ) = 6.5 to calculate the DEM curves. The integral over all bins gives the total emission measure (EM). From the Figure 2 Average coronal hole DEM for 29 May, 2013 using level 1.5 data (red line), level 1.6 (PSF-corrected; blue line), and level 1.6 with eclipse correction (green line). The error bars represent the uncertainties from the DEM calculations of all pixels of the coronal hole. Note that the error bars are shown only for every third bin for better visualization.
DEM curve, the EM-weighted median temperature is calculated such that The median temperature was chosen over a mean, or EM-weighted mean temperature, because it better describes the asymmetrical DEM profile and thus better represents the dominant emission from coronal holes (also see DEM curves in Saqri et al., 2020). For each pixel i, we derive the EM-weighted median temperature [T median,i ], and calculating the mean of all pixels in a coronal hole gives the average coronal hole temperature as T CH = 1 N T median,i with N being the number of coronal hole pixels.
In Figure 2, we show the DEM solutions for a coronal hole observed on 29 May, 2013 for level 1.5, level 1.6 data, and level 1.6 plus applying the eclipse correction. The non-PSF-deconvolved solution shows a peak at high temperatures, which is lower in the PSFcorrected solution and reduces further when considering the remaining counts derived from the lunar-eclipse analysis. This finding supports that the contribution of coronal hole emission at quiet-Sun temperatures, also found by Hahn, Landi, and Savin (2011) and Saqri et al. (2020), is mainly due to stray light from regions outside the coronal hole.
When calculating the uncertainties, three components have to be considered: the error from the DEM calculation [σ DEM ], which comes from the method and the initial uncertainties in the observations, the variation of the individual pixel values [σ CH : standard deviation of pixel values], and the variation of the coronal hole mean values [σ TOTAL : standard deviation of coronal hole mean values]. The resulting uncertainty can be given as follows: withσ DEM,j being the mean DEM error and the index j running over all coronal holes.

Correlations
The correlation analysis was done using a bootstrapping method (Efron, 1979;Efron and Tibshirani, 1993) with > 10 4 repetitions to derive Pearson correlation coefficients that take The white contour is the coronal hole boundary as extracted by CATCH (Heinemann et al., 2019). The temperature is distributed around a mean of T CH = 0.94 ± 0.05 MK and thus not well visible. Note that the values outside of the coronal hole boundary should be considered with care as the DEMs were calculated using photospheric abundances, which are applicable in coronal holes but not in the quiet-Sun regions, which require coronal abundances.
into account the uncertainties of the parameters. Additionally, to show that the correlations found do not depend on the data preparation, i.e. on the choice of the PSF nor the estimated correction for the remaining counts, we calculated all correlations with two different PSF deconvolutions (aia psf and a PSF by Poduval et al. 2013) and different correction values. This is shown in Tables 2 -4 in the appendix. Although the values of the plasma parameters can change on average by up to 20% depending on the data preparation, the correlations do not change significantly, indicating that the correlations shown in the following sections are reliable estimates.

Results
We derived the plasma properties of a set of 707 coronal holes by performing a DEM analysis using the coronal EUV observations by SDO/AIA and obtained the following results. Figure 3 shows an example of the plasma properties (i.e. temperature, density, and emission measure) of a coronal hole and the surrounding quiet-Sun areas on 8 September, 2015. The temperature shows only small variations around a mean value of T CH = 0.94 ± 0.05 MK, and is statistically uniformly distributed inside the coronal hole. The density and emissionmeasure maps show a gradient from the boundary inward with regions of lesser emission and density at some distance from the boundaries. We confirm this statistically in Section 3.4. Figure 4 shows the median and the 80th and 90th percentiles of the superposed mean DEM curves of all coronal holes to present the range of values found (note that the errors in the DEM are not shown for this figure). We found that the DEM of a coronal hole is gaussian shaped around a peak temperature of roughly 0.9 MK with a tail towards the higher temperatures. We note that the curve is not meaningful for temperatures lower than 0.5 MK because the DEM is not well-constrained by the AIA filters (see Lemen et al., 2012, for filter sensitivity and response curves). The superposed DEM curve shows that the shape is very Figure 4 Median, 80th, and 90th percentiles of 707 superposed mean coronal hole DEM curves. The black line gives the median, the blue shaded area the 20th and 80th percentiles, and the gray-shaded area the 10th and 90th percentiles. similar for all coronal holes, the major variation is in the height of the peak, which can vary up to a factor of two. A contribution at quiet-Sun temperatures (≈ 1.4 -1.6 MK), which is often interpreted as stray light, is hinted at. From the average DEM curve of each coronal hole, we derived the temperature, electron density, and emission measure according to the equations in Section 2.4. Figure 5 shows the distribution of the average coronal hole plasma properties. By averaging over the 707 coronal holes, we find an average coronal hole temperature [T ] of 0.94 ± 0.18 MK, with a derived minimum of 0.86 MK and a maximum of 1.32 MK. For all coronal holes under study, the average temperature [T CH ] is significantly lower than average quiet-Sun temperatures. The electron density is distributed around (2.37 ± 0.72) × 10 8 cm −3 , with a minimum of 1.43 × 10 8 cm −3 and a maximum of 3.47 × 10 8 cm −3 . The emission measure was found in a range between 1.15 × 10 26 cm −5 and 5.76 × 10 26 cm −5 with a mean of (2.79 ± 1.57) × 10 26 cm −5 .

Plasma Properties over the Solar Cycle
Due to the large number of coronal hole observations in the CATCH catalog that span from 2010 to 2019, we can study the plasma properties over almost the full Solar Cycle 24 and investigate a possible dependence on solar activity. In Figure 6 we present the calculated plasma properties together with the international sunspot number (bottom panel) as functions of time. The coronal hole temperature shows small variations over time, which seem to follow the solar-activity cycle. When considering the uncertainties of the average coronal

Correlation of Coronal Hole Properties
To investigate how plasma properties are correlated with morphological and magnetic coronal hole properties, we plot the average coronal hole temperature [T CH ], electron density [n e ], and emission measure [EM] as a function of coronal hole area and signed mean magnetic-field density, which were obtained from the CATCH catalog. This is shown in Figure 8 Distribution of the plasma properties of individual pixels as a function of the distance to the nearest coronal hole boundary [d]. For each vertical bin (of size 5 ) the normalized distribution of pixels and their temperature, density, and emission measure (from top to bottom) was calculated. A darker shade represents a higher percentage of pixels in the according bin. Note that the distributions with d > 100 are very uncertain because of the low number of pixels available. Thus, the focus is on the region between 0 and ≈ 80 . Such a large minimum distance to the coronal hole boundary is only possible in very large coronal holes, which are rare. The negative distances represent pixel outside of the coronal hole, which should be considered with care as discussed in Section 3.4.

Spatial Distribution of Plasma Properties in Coronal Holes
In addition to the average coronal hole plasma properties, we investigated how the properties are spatially distributed within the coronal hole. To this aim, we calculated for each coronal hole pixel the distance to the closest coronal hole boundary [d in arcsec] and derived the temperature, electron density, and emission measure of each individual pixel as a function of the distance, which is shown in Figure 8. In each bin of size 5 , the pixel distribution of T median , n e , and EM is given. Each bin is normalized to reflect the probability for a pixel in a given distance bin to have a certain value. In addition to the pixels within the coronal holes, the DEM for pixels surrounding the coronal hole boundary was calculated (represented by negative distances) to show the general trend outside of the coronal hole. We note that the calculations outside are not reliable as they were calculated in the same way as for coronal hole pixels, but here the assumptions of a hydrostatic scale height and photospheric abundances are not valid. Figure 8 shows that the temperature (top panel) within the coronal holes is very uniform and does not depend on the distance from the EUV extracted boundary. Near the boundary, small variations are seen (within ≈ 20 ) and the temperature does not change strongly at the coronal hole boundary. For the electron density as well as for the emission measure we find a dependence on the distance to the closest coronal hole boundary. Within a distance of 20 from the coronal hole boundary, on average the electron density drops by a factor of ≈ 2 from roughly 2.8 to 1.5 × 10 8 cm −3 and the emission measure by a factor of ≈ 2.5 from 4.0 to 1.5 × 10 8 cm −5 . At distances > 20 the dependence ceases and the individual pixels are distributed around (1.87 ± 0.61) × 10 8 cm −3 for the density and around (2.13 ± 1.32) × 10 26 cm −5 for the emission measure (also considering the individual DEM errors). These values were derived for all pixels that are located at distances larger than 20 from the closest coronal hole boundary inside coronal holes. When comparing the average pixel values for the electron density and emission measure at distances smaller and larger than 20 from the closest coronal hole boundary, we find that the mean pixel densities close (d < 20 ) to the boundary are a factor 1.3 higher than further away (d > 20 ). For the emission measure we find a factor of 1.5. The gradient in the electron density and emission measure shows that the coronal hole boundary extracted using an intensity-threshold technique is a good tracer for an area of reduced density.

Discussion
The coronal hole properties derived from the presented statistical DEM study are in good agreement with the results of individual case studies and observation campaigns. Fludra, Del Zanna, and Bromage (1999) used the CDS on SOHO to derive coronal hole densities and temperatures as a function of height. The temperatures range from 0.75 MK at a radial distance of 1.0 R to 0.85 MK at a radial distance of 1.1 R . This approximately agrees with the results derived in this study (T = 0.94 ± 0.18 MK). The derived average coronal hole temperatures are also in good agreement with the 0.9 MK derived by Landi (2008), Hahn, Landi, and Savin (2011), Wendeln and Landi (2018), and Saqri et al. (2020). Further, it is notable that the temperature is almost uniformly distributed (σ T = 0.05 MK) within the coronal holes.
When using the DEM analysis to derive electron densities, multiple assumptions have to be made. The densities are calculated for the emission of a LoS column, which we defined as the hydrostatic scale height, which is believed to be a reasonable first-order approximation due to the open field lines in coronal holes. In the quiet Sun a hydrostatic scale height cannot be used due to the abundance of primarily closed fields. As the magnetic topology at the coronal hole boundaries is not well known, it is unclear how valid our assumptions are in this regime. Additionally, it is known that stray light in EUV images can contribute up to 40% of the derived electron density value (Shearer et al., 2012), and it is unclear, even after a correction, how much stray and scattered light still remains in the images (Poduval et al., 2013). The derived electron densities of (2.36 ± 0.72) × 10 8 cm −3 are in fair agreement with the DEM study of the evolution of one particular coronal hole performed by Saqri et al. (2020), who derived values between 1.9 and 2.4 × 10 8 cm −3 , and Warren and Hassler (1999) who derived values between 1.8 and 2.4 × 10 8 cm −3 . Hahn, Landi, and Savin (2011) used Hinode/EIS spectroscopy of a polar coronal hole and derived density values of 1.0 × 10 8 cm −3 and 1.8 × 10 8 cm −3 from the Fe VIII and Fe XIII line ratios, respectively. They suggested that the density derived from the cooler lines (Fe VIII) probes the cool coronal hole plasma and the hotter lines show a contribution from quiet-Sun plasma. This also agrees with the results by Pascoe, Smyrli, and Van Doorsselaere (2019) that estimated coronal hole electron densities are around ≈ 10 8 cm −3 . Figure 8 shows that the density decreases as a function of the distance from the coronal hole boundary up to d ≈ 20 . Under the assumption that the PSF correction used does indeed remove most of the stray light, this result suggest that the extracted coronal hole boundary does not represent a vertical separation of two magnetic regimes but rather the position where the emission drop is the strongest. Thus, the LoS column might represent a mixture of coronal hole plasma emission and quiet-Sun plasma emission, maybe due to field lines that are bent rather than radial (e.g. forming an inclined separation between coronal hole and quiet Sun). We notice that the decrease in electron density from the coronal hole boundary to d ≈ 20 within the coronal hole is approximately 30%, which is smaller than the 50% that has been found by Doschek et al. (1997).
The plasma properties show no correlation with the area and magnetic-field density of coronal holes. This suggests that the size of a coronal hole does not determine the plasma properties nor vice versa (we find this to be valid for non-polar coronal holes). A trend can be seen in the correlation of the mean magnetic-field density and the average coronal hole temperature, however this is not significant (cc T = 0.11, CI 90% = [0.05, 0.17]) when considering the large uncertainties.
The investigation of the plasma properties over the course of the solar cycle between 2010 and 2019 revealed only small variations in density and emission measure. This is especially intriguing as other coronal hole parameters indeed reveal a solar-cycle dependence. Heinemann et al. (2019) found that the mean 193Å intensity and the mean magnetic-field density show a dependence on solar activity and so does the long-term evolution of large coronal holes . Only the average coronal hole temperature does show some variation, which might be linked to the solar cycle where slightly higher spreads during solar maximum are observed than during solar minimum. However, the Pearson correlation coefficient does only show a very weak correlation of cc T = 0.13 and CI 90% = [0.07, 0.19]. Doyle et al. (2010) showed that the temperature in coronal holes varies with the solar cycle between 0.8 MK at solar minimum and 1.04 MK at solar maximum, which we find to be in fair agreement with this study when considering the average values in the correlation, without taking into account the uncertainties (cc T = 0.57, CI 90% = [0.52, 0.61]).

Summary and Conclusions
In this statistical study on the DEM analysis of coronal holes we investigated temperature, electron density, and emission measure of the 707 coronal holes from the CATCH catalog and analyzed their distribution, variability over Solar Cycle 24, and correlations to the coronal hole area and magnetic-field density. Our major findings can be summarized as follows: i) DEM: The shape of the average DEM curve for coronal holes is very stable and resembles a Gaussian profile centered around a peak temperature with an extended tail towards higher temperatures. ii) Temperature: We find that coronal holes show a EM-weighted median temperature of 0.94 ± 0.18 MK, which is in accordance with previous studies. Additionally, we found that the temperature is spatially very uniform inside the coronal hole.
iii) Electron Density: Using the DEM analysis we derive electron number densities for the coronal holes to be (2.4 ± 0.7) × 10 8 cm −3 . The density within coronal holes decreases from the boundary to 20 by ≈ 30% and approaches a constant level further inside. iv) Solar Activity: We observe only small variations in the temperature, electron number density, and emission measure during the period from 2010 to 2019 but these appear not to be correlated to changes in solar activity. Although the average coronal hole temperature may hint toward some solar-cycle dependence, due to the uncertainties this dependence is not significant.
From the statistical analysis we find that coronal holes show a strong similarity in their DEM curves and derived properties. Thus, coronal holes can be can be clearly defined by their plasma properties. This not only enables a deeper understanding of the structure of coronal holes but can also serve to constrain the input for models (e.g. solar-wind models such as Odstrčil andPizzo, 1999 or Pomoell andPoedts, 2018) and studies of coronal waves interacting with coronal holes (e.g. Podladchikova et al., 2019;Piantschitsch, Terradas, and Temmer, 2020). a Stray-light correction derived from lunar eclipse data as given in Table 1. b Stray-light correction derived from Venus transit data as given by Saqri et al. (2020). Table 3 Pearson correlation coefficients for the coronal hole electron density against coronal hole area, signed mean magnetic-field density, and solar activity approximated by the international sunspot number and calculated from different input configurations. The correlation coefficients are given as calculated with and without uncertainties. In the square brackets the CI 90 % are given. a Stray-light correction derived from lunar eclipse data as given in Table 1. b Stray-light correction derived from Venus transit data as given by Saqri et al. (2020). analysis in a reasonable amount of time. The resulting bootstrapped Pearson correlation coefficients for the average coronal hole temperature, electron density, and emission measure are presented in Tables 2, 3, and 4 respectively. We find that the data preparation, i.e. which Table 4 Pearson correlation coefficients for the coronal hole emission measure against coronal hole area, signed mean magnetic-field density, and solar activity approximated by the international sunspot number, and calculated from different input configurations. The correlation coefficients are given as calculated with and without uncertainties. In the square brackets the CI 90 % are given. PSF was used and what correction is applied, does not significantly change the correlations, and as such the results are reliable.