Issue |
A&A
Volume 565, May 2014
|
|
---|---|---|
Article Number | A14 | |
Number of page(s) | 11 | |
Section | The Sun | |
DOI | https://doi.org/10.1051/0004-6361/201423471 | |
Published online | 21 April 2014 |
Elemental abundances and temperatures of quiescent solar active region cores from X-ray observations⋆
DAMTP, Centre for Mathematical Sciences, Wilberforce Road, Cambridge CB3 0WA, UK
e-mail: g.del-zanna@damtp.cam.ac.uk
Received: 20 January 2014
Accepted: 18 March 2014
A brief review of studies of elemental abundances and emission measures in quiescent solar active region cores is presented. Hinode EUV Imaging Spectrometer (EIS) observations of strong iron spectral lines have shown sharply peaked distributions around 3 MK. EIS observations of lines emitted by a range of elements have allowed good estimates of abundances relative to iron. However, X-ray observations are required to measure the plasma emission above 3 MK and the abundances of oxygen and neon. We revisit, using up-to-date atomic data, older X-ray observations obtained by a sounding rocket and by the Solar Maximum Mission (SMM) Flat Crystal Spectrometer (FCS). We find that the Fe/O and Fe/Ne ratios are normally increased by a factor of 3.2, compared to the photospheric values. Similar results are obtained from FCS observations of six quiescent active region cores. The FCS observations also indicate that the emission measure above 3 MK has a very steep negative slope, with very little plasma observed at 5 MK or above.
Key words: Sun: abundances / Sun: corona / Sun: X-rays, gamma rays / techniques: spectroscopic
Appendix A is available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
This paper discusses measurements of elemental abundances and temperatures in the hot active region (AR) cores obtained from remote-sensing instruments. By “AR cores” we mean the hot loops that typically have temperatures of 3 MK, hence are very bright in emission lines formed at this temperature, in particular in the X-rays. Previous measurements are briefly reviewed, and some X-ray observations revisited, using up-to-date atomic data. Measurements of elemental abundances and temperatures are closely linked, since the differential emission measure (DEM) distribution is used to both describe the temperature distribution of the plasma and obtain relative elemental abundances.
It is well known that solar coronal abundances present variations and differ from the photospheric ones. The ratio of the abundances of the low (≤10 eV) first ionisation potential (FIP) elements vs. the high-FIP ones is higher in coronal plasma, however the absolute abundances (relative to hydrogen) are still debated, as well as the amplitude of the difference. For recent reviews on the issue of solar abundances see, e.g., Asplund et al. (2009), Lodders et al. (2009), Schmelz et al. (2012).
There is a vast literature on the subject of solar coronal abundances, but comparatively little related to active region (AR) cores. In fact, most published results on chemical abundances are based on either observations with uncollimated spectrometers (often full Sun), or flare observations.
Active region abundance measurements relative to iron.
In terms of AR core temperature measurements, there is a vast literature. Rosner et al. (1978) found that Skylab X-ray imaging was consistent with the hot core loops in quiescent ARs being approximately isothermal around 3 MK, an observational aspect that was found puzzling. Early X-ray spectroscopy also found a similar result (Evans & Pounds 1968). Later spectroscopic measurements in the X-rays from the Solar Maximum Mission (SMM) Flat Crystal Spectrometer (FCS) of many ARs have often shown that little material at temperatures above 3 MK is present, and therefore most of the abundance analysis have assumed that the plasma was effectively isothermal (see, e.g. Schmelz et al. 1996).
Extreme ultraviolet (EUV) spectroscopy with the SOHO Coronal Diagnostic Spectrometer (CDS; Mason et al. 1999; Del Zanna & Mason 2003) and Hinode EUV Imaging Spectrometer (EIS; O’Dwyer et al. 2011; Tripathi et al. 2011; Winebarger et al. 2011; Warren et al. 2012; Schmelz & Pathak 2012; Del Zanna 2013b) have indicated that the plasma distribution in the 1–3 MK range is very steep, and that little material is present at higher temperatures, although this issue is still debated in the literature.
The EUV measurements such as those from Hinode/EIS constrain very well the plasma distribution up to 3 MK, but are limited at higher temperatures, because hotter lines are normally weak and/or blended. For a discussion of the Hinode EIS Fe xvii lines see Del Zanna & Ishikawa (2009), while details on the hotter Hinode EIS lines have been given e.g. in Young et al. (2007), Del Zanna (2008a, 2012, 2013b), Warren et al. (2008); Ko et al. (2009), Del Zanna et al. (2011b).
As discussed in Del Zanna (2013b), previous studies that did not consider the in-flight Hinode EIS calibration (Del Zanna 2013a) underestimated by a factor of two the intensities in lines formed around 3 MK, hence the actual steepness of the EM distribution of the active region cores was underestimated.
Assessing the temperature distribution of the AR core hot plasma above 3 MK is an important constraint to nanoflare models (see, e.g. Patsourakos & Klimchuk 2009; Bradshaw & Klimchuk 2011, and references therein). There is ample literature on measurements of this hot component in active regions, using broad-band X-ray observations from e.g. Yohkoh SXT and Hinode XRT (see, e.g. references in Reale 2010). However, these instruments are broad-band and the results are therefore uncertain.
Fe xviii emission is often observed in ARs with the Solar Dynamic Observatory (SDO) Atmospheric Imaging Assembly (AIA, see Lemen et al. 2012) 94 Å band, so it has been used to provide a further constraint on the emission measure (EM) above 3 MK. For example, Warren et al. (2012) combined Hinode EIS with AIA 94 Å observations to find the slope β of the distribution (EM = T− β) in the 3–10 MK range. The β values ranged between 6 and 10.
However, as shown in Del Zanna (2013b), a considerable part of the Fe xviii emission can come from the 3 MK plasma and the 94 Å band is highly multi-thermal, so it is very easy to overestimate the Fe xviii emission. A proper assessment of the plasma at temperatures higher than 3 MK requires observations in the X-rays.
The first observations of active region cores were carried out during total solar eclipses. The active region cores were called coronal condensations. We searched the literature of eclipse observations and found that the only clear and well-calibrated observation of a coronal condensation was done in 1952 in Kharthoum, as described by Aly et al. (1962). The coronal condensation that was observed was unusually bright and a series of calcium (including the hot Ca xv) and iron lines were observed. The continuum was also observed, so in principle these are excellent observations to constrain the absolute abundance of iron and calcium, as well as the temperature. The interpretation of the line and continuum emission is not straightforward though, and requires some modelling and some assumptions on the emitting source. This was done by Mason (1975), where the calcium abundance was found to be 6.8 dex and the iron abundance 7.85 dex (in the usual logarithmic notation where the hydrogen abundance is 12), i.e. an enhancement of about a factor of three for both elements, relative to their photospheric values. We also note that Young et al. (1997) reviewed the original measurements by Aly et al. (1962) of an Ar xiv forbidden line to find an Ar/Ca relative abundance of 0.85 ± 0.20.
With the advent of the space observations, an extensive literature on EUV measurements has appeared, however very few abundance measurements pertained to AR cores. In fact, many observations in the EUV either used emission lines formed around or below 1 MK, hence probed different AR structures, or compared emission lines formed over very different temperatures, hence likely formed in spatially different regions. One example is the paper on the OSO-6 observations by Withbroe (1975). A significant number of results pertain to either flare spectra or the legs of warm loops, i.e. only point-like features (as observed e.g. with the Skylab NRL S082A EUV spectra, see Widing & Feldman 1993).
A sample of abundance measurements of active region cores that are relevant for this paper are shown in Table 1. Monsignori Fossi et al. (1994) used SERTS EUV observations of an AR core to obtain, with a DEM analysis, relative elemental abundances for many elements. The DEM presented the usual two-peak structure, with one peak at 1 MK due to the overlying corona, and one at 3 MK due to the core of the AR. Lines from S, Ni, Ar, Fe formed around 3 MK were observed. The following abundances were obtained: Fe/S = 4, Fe/Ni = 8.3, Fe/Ar = 25.
The SOHO SUMER spectrometer carried out limited observations of ARs in the UV. Mohan et al. (2000) used the observed intensities of lines from Si xi, K xiii, S x, and Ar xii to study the relative abundances of Si, K, S, and Ar (two low-FIP and two high-FIP elements). The SUMER observations scanned an AR off-limb. The S/Ar abundance was found to be nearly constant with height at around 4.5–5, i.e. close to its “photospheric” value recommended by Asplund et al. (2009), 5.2. The Si/Ar abundance was found to decrease with height, while the K/Ar was found to increase, a puzzling result. Both Si/Ar and K/Ar relative abundances were higher than the photospheric values.
That active region cores were more consistent with a FIP enhancement of 4 rather than photospheric abundances was previously noted (see, e.g. Tripathi et al. 2011). However, the first comprehensive measurement of the abundances of several elements in an AR core was recently published in Del Zanna (2013b), using Hinode EIS observations and a recent calibration (Del Zanna 2013a). The relative abundances of several key elements were obtained, including low-FIP (e.g. Fe, Si, Ca, Ni) and high-FIP (e.g. Ar, S). A FIP enhancement of 3.2 was found. Other active region cores showed similar abundances. A consideration of the path length suggested that the iron absolute abundance must be enhanced, compared to its photospheric value, by at least a factor of three. The Ca and Ar abundances, relative to iron, were well constrained, and were 7.4 × 10-6 and 2 × 10-6. The line-to-continuum measurements of the Ca abundance by Mason (1975) was 6.3 × 10-6, in excellent agreement, while the Ar abundance of 5.4 × 10-6 by Young et al. (1997) is significantly higher.
One limitation of the EUV observations is the inability to constrain the relative abundances of oxygen and neon, two important high-FIP elements. X-ray observations are needed for this, due to the presence of the strong resonance lines from O viii 19 Å and Ne ix 13.4 Å.
An early measurement of the AR abundances of oxygen and neon using X-ray spectroscopy was obtained by Evans & Pounds (1968) with a rocket spectrum. They found the AR core to be nearly isothermal at 3 MK, and measured the relative abundances between iron, oxygen and neon. They obtained Fe/O = 0.28 and Fe/Ne = 2.5. Davis et al. (1975) also used a rocket X-ray spectrum (with a collimator) to obtain Fe/O = 0.26 and Fe/Ne = 1.7. One of the best X-ray spectra of a quiescent AR was obtained by Parkinson (1975). Parkinson (1977) used these AR observations to obtain, using essentially a line ratio method, the relative abundances of a few elements. Fixing the abundance of Si to 7.65 dex, he measured the abundances of O (8.45), Ne (7.65), Na (6.5), Mg (7.65), Fe (7.65) and Ni (6.4), i.e. found relative abundances between iron, oxygen and neon that were quite different (Fe/O = 0.16 and Fe/Ne = 1.0) from the previous measurements.
Excellent X-ray spectra of ARs were later obtained by the SMM FCS. There are several observations of the cores of quiescent ARs. The intensities of the strong lines were found by various authors (see, e.g. Waljeski et al. 1994) to be generally consistent with an FIP enhancement of a factor of 3.5 (Fe/O = 0.16 and Fe/Ne = 1.12), as suggested by the compilation of Meyer (1985).
However, large variations in the relative abundances of low-FIP vs. high-FIP elements as measured from SMM/FCS have been reported by several authors, for example Saba & Strong (1993), Schmelz et al. (1996). Schmelz et al. (1996) performed a review of several SMM/FCS quiescent AR observations, and used the EM loci method to find large variations in the Ne abundance.
McKenzie & Feldman (1992) also found a large scatter of Fe/Ne and Fe/O values, indicative of elemental abundance variations (Fe/Ne ranged between 0.3 and 1.1, while the Fe/O varied between 0.04 and 0.1). These authors studied the intensity ratios of the O viii 19 Å, Ne ix 13.4 Å, Fe xvii 15.0 Å, and Fe xviii 14.2 Å observed by the SOLEX B spectrometer on board the P78-1 satellite in active regions. Their analysis did not take into account the differences in the temperature response of these ions though.
There are various reasons to revisit the previous results obtained from the X-ray observations. First, new atomic data have been made available. The strongest Fe xvii and Fe xviii lines notoriously showed large (factors of 2–3) discrepancies between observed and predicted intensities. It turned out that the problems were largely due to the limitation of distorted wave (DW) calculations, as shown in Del Zanna (2011, 2006). Indeed, we found excellent agreement between solar observations and theory when R-matrix calculations (Liang & Badnell 2010; Witthoeft et al. 2006)1 are used. We note, however, that many previous abundance studies used the strongest Fe xvii and Fe xviii resonance lines, for which the predicted line intensities with the new atomic data changed by small amounts (typically 10–20%). The really major difference in predicted intensities comes from recent calculations of ionisation and recombination rates. The latest comprehensive set of rates produced ion fractions (made available within CHIANTI v.6, see Dere et al. 2009) that are significantly shifted in temperature compared to previous calculations such as those of Arnaud & Rothenflug (1985), which were extensively used in the previous studies. For Fe xvii, the combined changes in the atomic data produces well over a factor of two difference in the emissivity at the peak ion abundance, which is indeed significant.
In terms of measuring abundances, it is useful to revisit older observations because the standard photospheric abundances of many elements have recently changed significantly, and because most previous studies made the assumption that the AR cores are isothermal. This is a reasonable assumption, although we note that any small departures could enhance the different temperature response of the various ions, leading to incorrect results.
In the next two sections we present a discussion of the best X-ray measurements of AR cores we have available, while in the last section we draw our conclusions.
2. Parkinson (1975) AR spectrum
An instrument, built by the University of Leicester (UK), consisting of Bragg crystal spectrometers with a collimator was flown on a British Skylark sounding rocket on 30 Nov. 1971 and recorded the best X-ray solar spectrum of a quiescent active region (Parkinson 1975). The observations obtained with the KAP crystal in the 9–23 Å range are used here. The instrument had an excellent spectral resolution, was radiometrically calibrated and the whole spectral region was observed simultaneously, unlike many other X-ray instruments which scanned the wavelengths (e.g. SMM/FCS, see below). The Fe xvii lines, the He-like Ne ix 13.45 Å (1s2–1s 2p), the H-like Ne x 12.1 Å (1s–2p), the O vii, O viii were all prominent. The strongest Fe xviii line, at 14.20 Å, was barely visible. Lines formed at higher temperatures were totally absent, indicating a lack of plasma with T> 4 MK.
One simple way to estimate if the plasma is nearly isothermal is the EM loci method, by which the curves Iob/ (A × G(T)) as a function of temperature are plotted (Iob is the observed intensity, A the element abundance and G(T) the contribution function containing the atomic parameters). If the curves cross at one temperature, then the plasma can be modelled with an isothermal distribution. In general, the loci of the curves represent upper limits to the emission measure distribution. For details, see Strong (1978) and Del Zanna et al. (2002).
Del Zanna (2011) applied the EM loci method to find a near-isothermal distribution for the Parkinson (1975) data and to show that the Fe xvii lines, in conjunction with R-matrix calculations (Liang & Badnell 2010), offer an excellent diagnostic to measure the electron temperature. For the first time, excellent agreement between the observed and predicted Fe xvii intensities for the strongest lines was found, showing the potential of these lines for measuring temperatures. A value around log T[K] = 6.5, i.e. 3 MK was found. The relative Fe/Ne abundance needed to be increased by about a factor of two, compared to the recommended photospheric abundances of Asplund et al. (2009).
To obtain a proper measurement of the elemental abundances, we have performed a DEM analysis using the method described in Del Zanna (1999). The latest v.7.1 CHIANTI atomic data and ionisation equilibrium tables (Landi et al. 2013) are adopted here. We converted the Parkinson (1975) intensities assuming a collimator angular (FWHM) dimension of 3′.
We obtained the results shown in Fig. 1 and Table 2. Excellent agreement (within a relative 20%) between observed and predicted radiances is found for all the strongest lines, with the elemental abundances shown in Table 1, where the low-FIP Fe and Mg elements are enhanced by a factor of 3.2, and Ne, O maintain the photospheric abundances of Asplund et al. (2009). The points in Fig. 1 (top) are plotted at the effective temperature
which is an average temperature more indicative of where a line is formed. This is often quite different than Tmax, the temperature where the G(T) of a line has a maximum. The bottom plot of Fig. 1 shows the corresponding values at Tmax. Teff is nearly 3 MK for all the Ne ix, Fe xvii, O viii lines, as detailed in Table 2, where more information is provided. Note that the Fe xvii lines have Tmax = 5.6 × 106 K (log T = 6.8). When the DEM has a steep slope or the lines are significantly blended, Tmax can be quite different than Teff.
Fig. 1 Top: DEM for the Parkinson (1975) active region spectrum and the present elemental abundances. The ion and the wavelength (Å) is given. The numbers in parentheses are the ratios between predicted and observed intensities. The points are plotted at the effective temperature Teff, and at the theoretical vs. the observed intensity ratio multiplied by the DEM value. Bottom: same DEM, but with the points plotted at the temperature of the maximum of the G(T). |
Observed and predicted radiances for the Parkinson (1975) AR observation.
The 16.0 Å line is often a blend of O viii and Fe xviii, as discussed in the Fe xviii benchmark paper (Del Zanna 2006). However, in the observations considered here the Fe xviii contribution is small, so the observed intensity can be used to constrain very accurately the O/Fe abundance, considering that the effective temperature of this O viii line is the same as that of some of the Fe xvii strong lines, as shown in Table 2. The other O viii self-blend, at 18.97 Å, has an effective temperature slightly less. The Ne ix lines have an effective temperature in between the O viii 16.0 and 18.97 Å lines, so the Ne ix abundance can be estimated quite accurately with respect to the O abundance. The H-like Ne x has also an effective temperature close to Fe xvii, and confirms the Ne abundance obtained from Ne ix. This is an important result. The Mg xi 9.17 Å line also has a similar effective temperature. We obtain the Mg abundance shown in Table 1.
The Parkinson (1977) relative abundances, if we scale the iron value to the present one (8.0 dex), were: O (8.80), Ne (8.00), Mg (8.00), Si (8.00), and Ni (6.75). They are slightly different (by about 25%) from the values obtained here for O, Ne and Mg. The Ni value is 65% higher than the value obtained in Del Zanna (2013b). We note that the few Ni xix lines in the Parkinson (1975) spectrum are formed at high temperatures, where there are no firm constraints from iron lines. On the other hand, the Del Zanna (2013b) results are based on a few strong lines from Ni xv, Ni xvi, Ni xvii all formed around 3 MK, where good constraints from iron lines are available. So, the nickel abundances obtained by Del Zanna (2013b) are reliable, although we note that the atomic data for these nickel ions need some improvement.
3. SMM FCS
The SMM FCS database was searched for quiescent AR observations. The types of observations considered here typically started with a monochromatic image of an active region. The instrument was then pointed at the brightest region (in the AR core), where the FCS scanned the 13–20 Å region. Each scan lasted about 10 min. We obtained a long list of AR observations, summarised in Table 3. They include all those listed by previous authors, e.g. Schmelz et al. (1996), Phillips et al. (1997), and Saba et al. (1999).
List of the SMM FCS AR observations considered here.
Fig. 2 Light curves in the SMM BCS Ca xix channel (data are count rates) for the selected active regions cores analysed here. The dashed lines indicate the start and end times of the FCS observations selected for the DEM analysis. |
We searched the count rates in the SMM bent crystal spectrometer (BCS) Ca xviii/Ca xix channel for signs of variability. Since FCS scans across wavelength with time, it is in fact necessary to ensure that as little solar variability is present during each scan as possible, because we are using for the diagnostics intensities of lines observed at different wavelengths (times). We are also interested in quiescent AR core observations, i.e. not in hot loops formed after obvious small flaring events, which are quite common in active regions.
We found that many AR observations in our list were not actually taken during quiescence, and so we discarded them. Figure A.1 shows the BCS count rates for the observations which we discarded in our analysis. In addition, many other observations had a very low signal, or did not have any BCS data associated with them, and so these datasets have also been discarded. The observations selected for further analysis are listed in Table 3, and the respective BCS count rates shown in Fig. 2.
Observations of six active regions remained for further analysis. To have an idea of the type of active regions selected by this process, we searched for ground-based observations. Figure 3 shows Kitt Peak Vacuum Telescope (KPVT) magnetograms of these six ARs, obtained from the Virtual Solar Observatory. There is clearly a variety of active regions, some more dispersed, some more compact, some younger, some older.
Fig. 3 KPVT magnetograms for the selected active regions analysed here. |
Fig. 4 SMM/FCS AR spectrum of 1987 Dec. 13. |
NOAA 4729 was a fairly extended region with small sunspots, and probably at least a week old, when it was observed by SMM FCS near central meridian on 1986 May 20. NOAA 4731 was observed on 1986 May 21 about three days after its appearance at the east limb. It was a fairly large AR with many small sunspots, probably a few days old. It was observed again on May 23 and 24 when it was quiescent. NOAA 4790 was a fairly large active region with a large sunspot group. It was observed on 1987 Apr. 11. It appeared on the east limb on Apr. 6, hence it was at least a week old. It was observed again on Apr. 13 and 14 as it approached the central meridian. However, during the later observations, a lot of variability was present. NOAA 4891 appeared small and young and was at the eastern limb on 1987 Nov. 17. It was observed on Nov. 27 and 29 when it was near and at the west limb, hence when it was over 10 days old. NOAA 4901 appeared in the KPVT magnetograms of 1987 Dec. 3 near the Sun centre. It was observed with the FCS on Dec. 6–11. Unfortunately, no BCS observations for all these dates are available, except for Dec. 6 and 11. On Dec. 6, the BCS calcium channel shows signs of strong variability. On Dec. 11, the AR was near the west limb. There are two sets of FCS raster observations, but the BCS count rates are low and the FCS data are very noisy. Finally, NOAA 4906 was observed on 1987 Dec. 13 at the SE limb. It had a fairly large sunspot, and seemed young. It was observed later, during its meridian passage, on Dec. 15, 16, 18, and 20, but the BCS calcium channel showed strong variability on those dates.
Line intensities of the SMM FCS AR observations analysed.
For each SMM FCS observation selected, successive spectra taken during quiescence were averaged, to increase the signal-to-noise ratio (S/N), during the time intervals displayed in Table 3 and Fig. 2. One of the best spectra in terms of S/N was taken on 1987 Dec. 13, and is shown in Fig. 4.
The selected averaged spectra have then been rebinned in wavelength (to further enhance signal/ noise) and then flux-calibrated using the software available within Solarsoft2. This includes a correction for the in-flight degradation of the instrument, which was considerable (a factor of 10) at longer wavelengths. The correction was based on the assumption of carbon deposition on the front filter (see, e.g. Saba et al. 1999). The lines were fitted with Voigt profiles using the Solarsoft cfit package (Haugan 1997). For uncertainties in the line intensities, we have considered the sum (in quadrature) of the uncertainty derived from the fitting process, plus a reasonable 20% for calibration and atomic data errors.
Intensities have then been converted to phot cm-2 s-1 arcsec-2 by taking into account the FCS collimator angular dimensions, which are approximately 15′′ (FWHM, equivalent to 225 arcsec2), but are slightly dependent on the wavelength (K. Phillips, priv. comm.). The intensities of the main lines used for DEM analysis are shown in Table 4. There are some slight differences between our relative intensities and those published in the literature (normally within the uncertainties). We find good agreement with the Fe xvii line ratios published by Saba et al. (1999).
The Fe xvii lines are very strong in the FCS spectra. However, they are formed around 3 MK, as discussed earlier, and also shown below. One remarkable aspect of the FCS observations of quiescent ARs is the absence of hotter emission from higher ionisation stages of iron. For example, even the strongest Fe xviii line, the 14.20 Å self-blend, is often barely visible, as seen for example in Fig. 4. This confirms the general lack of hot emission above 4 MK in the cores of quiescent active regions.
The He-like Ne ix lines are normally observable, although often the best measurement is of the resonance Ne ix 13.45 Å line. The O viii 18.97 Å and 16.00 Å self-blends are often just above noise. The Mg xi line was found normally too weak to be reliable, so it has not been considered here.
We recall that Schmelz et al. (1996) used the EM loci method to find large variations in the Ne abundance, compared to O, Fe and Mg. Schmelz et al. (1996) adopted a standard set of coronal abundances close to those recommended by Meyer (1985). In particular, they adopted a Ne/H value of 3.55 × 10-5, an Fe/O standard abundance of 0.2, and an Fe/Ne abundance of 1.3 (Meyer 1985 recommended 0.16 and 1.1 respectively). These are very close to the relative abundances recommended here, 0.2 and 1.2.
In terms of atomic data, Schmelz et al. (1996) used the Arnaud & Rothenflug (1985) and Arnaud & Raymond (1992) ionisation equilibria, DW data for Fe xvii and Fe xviii, and a set of emissivities for the other ions tabulated by Mewe et al. (1985). It is not possible to reproduce exactly the Schmelz et al. (1996) results, however Fig. 5 (top) shows the EM loci curves for the SMM/FCS AR spectrum of 1987 Dec. 13 using their standard set of abundances, the Arnaud & Rothenflug (1985) ionisation equilibria for oxygen and neon, and Arnaud & Raymond (1992) for iron, with atomic data from CHIANTI v.4.2, which also had DW data for Fe xvii and Fe xviii. We can see that the curves cluster around 3 MK (log T[K] = 6.5), if one neglects the Fe xvii resonance line at 15.01 Å. We can therefore reproduce the Schmelz et al. (1996) result (their standard abundances for this dataset), although the same set of observations with the current atomic data show significant differences (see Fig. 5, bottom plot). The main differences are due to changes in the ion fractions for all the elements. Significant differences in the Fe xvii curves are also clear (see Del Zanna 2011 for detailed comparisons of DW and R-matrix calculations for this ion). The net result is a much better agreement between observations and the latest atomic data, with an isothermal temperature close to log T[K] = 6.4, but still similar relative abundances as those obtained with the older atomic data.
As in the previous case, we have set out to perform a full DEM analysis on all the datasets listed in Table 3, to improve on the strictly isothermal assumption. The results of the DEM analyses are summarised in Table 3. In all cases the Ne/O Asplund et al. (2009) photospheric abundance represents well the observations. The difference with the Schmelz et al. (1996) result (where the neon abundance was found to vary) are to be ascribed to the large differences in the ionisation equilibria mentioned previously, and the use of the DEM method, instead of the strictly isothermal assumption.
In most cases, we found the need to increase the Fe/O and Fe/Ne abundance by a factor of 3.2, compared to the photospheric values, in order to get consistent results. This is in agreement with the previous Fe/Ar and Fe/S results (Del Zanna 2013b). Only in one case (1987 Nov. 27), does the Fe abundance need to be increased by a factor of 5. The 1986 May 21 observation also indicates a need for such higher FIP bias, however the spectra are very noisy, in particular for the 16.0 Å line.
Fig. 5 EM loci curves for the SMM/FCS AR spectrum of 1987 Dec. 13. Top: with CHIANTI v.4.2, the Arnaud & Rothenflug (1985) and Arnaud & Raymond (1992) ion abundances, and the Schmelz et al. (1996) elemental abundances. Bottom: with the present atomic data and abundances. |
Fig. 6 Top: DEM distribution for the SMM/FCS AR spectrum of 1987 Dec. 13, with the Fe abundance increased by a factor of 3.2 compared to the Asplund et al. (2009) photospheric abundance (as Fig. 1). The points are displayed at the effective temperature. ul indicates an upper limit. The labels indicate the main ion and the wavelength in Å. Bottom: same curve, with the results of the XRT_DEM inversion code (black histogram) and of the Monte Carlo XRT_DEM inversion (grey histograms). |
Figure 6 (top) shows as an example the DEM distribution for the SMM/FCS AR spectrum of 1987 Dec. 13, where the Fe abundance was increased by a factor of 3.2 compared to the Asplund et al. (2009) photospheric abundance (O and Ne were kept at the Asplund et al. 2009 photospheric abundances). The SMM/FCS data provide very good constraints in the 3 MK to 5 MK range. Constraints at lower temperatures are more problematic. The O vii 18.63 Å line constraints lower temperatures but it is in fact barely visible. We have provided an intensity for this line, but it has an uncertainty of about 50% or more. The DEM distribution was obtained with the spline method, assuming that it has a decrease at temperatures lower than 3 MK, according to the EUV results mentioned in the introduction.
As shown in Del Zanna (2011), the Fe xvii lines are an excellent temperature diagnostic in the 3–5 MK region. The Fe xvii intensities are consistent with a sudden drop of the DEM from around 3 MK to 5 MK by at least 3 orders of magnitude. This is a solid and important result. In terms of emission measures, EM ≃ T-24 in the 3–5 MK range. The lack of any emission from hotter ions provides further constraints at the higher temperatures.
We have also tested various other DEM inversion methods. As we have found in previous cases, the results of automatic inversion methods are quite sensitive to the temperature range, and depend on the temperature bin size. For well-constrained cases, the results are very similar (Del Zanna et al. 2011c; O’Dwyer et al. 2014). Figure 6 (bottom) shows as an example the DEM distribution obtained from the XRT_DEM inversion code (black histogram) and of 400 Monte Carlo XRT_DEM inversions (grey histograms). The Monte Carlo simulations are obtained by randomly varying the input intensities within the estimated uncertainties. There is good agreement in the peak DEM and in the log T[K] = 6.6–6.7 range, while differences are present in the ill-constrained low- and high-temperature regions.
As in the case of the Parkinson’s spectrum, the O viii and Fe xvii have an excellent overlap in temperature, so the iron/oxygen abundance can be measured accurately. The same is true for the neon/oxygen abundance. Table 5 lists the observed and predicted radiances for the SMM FCS AR observation of 1987 Dec. 13, as in Table 2. Excellent agreement (within a relative 30%) between observed and predicted radiances is found for all the lines, with the elemental abundances shown in Table 1, where the low-FIP Fe is enhanced by a factor of 3.2, and Ne, O maintain the Asplund et al. (2009) photospheric abundances.
It is important to note the slight differences between the peak of the DEM (log T[K] = 6.47), the isothermal temperature that can be estimated from the EM loci plot (cf. Fig. 5, about log T[K] = 6.4), and the effective temperature at which most of the emission lines are formed (about log T[K] = 6.5), displayed in Table 5. The points in Fig. 6 are displayed at their effective temperature, with the corresponding uncertainties, which provide a measure on the uncertainty in the DEM at the various temperatures.
The DEM distributions for the other observations are shown in Fig. 7. They are similar, with a peak around 3 MK and a sharp drop off at higher temperatures, although the steepness at higher temperatures is not always the same. The weak Fe xviii emission that is sometimes present requires some additional emission measure around 5 MK. In general, the FCS DEM distributions are steeper than those of Parkinson (1975). We recall that the FCS measurements were made over small angular regions (15′′ FWHM) near the brightest part in the AR cores, while the Parkinson (1975) observation was over a much larger area (3′ FWHM), and probably encompassed different structures. This could result in a less sharply peaked distribution.
4. Conclusions
Active region cores are notoriously variable, especially during the early stages, with frequent small flaring activity, as has quite clearly been shown by current XUV broad-band instruments, such as Hinode XRT and SDO AIA 94, 131 Å (due to the emission from lines formed around 7–12 MK). The older SMM BCS (and, of course, Yohkoh SXT) observations of “quiescent” active regions also show frequent flaring activity. However, outside those periods of clear activity, there is very little (if any) emission above 3 MK, the peak temperature of most of the plasma. As discussed in Del Zanna (2013b), current EUV spectroscopy from Hinode EIS is limited in measuring weak emission above 3 MK. X-ray spectroscopy is needed in this regard, although as we have seen, it is limited in constraining lower temperatures below 3 MK.
We have revisited several older X-ray spectral measurements of quiescent active regions which we believe to be the best to date. The slope in the EM distribution that we found is much steeper than previous estimates, such as those of Warren et al. (2012), where values ranged between 6 and 10. For example, the strong Fe xvii lines indicate EM ≃ T-24 in the 3–5 MK range for the 1987 Dec. 13 SMM FCS observations. Often, estimating the distribution above 5 MK is not really possible from SMM FCS observations, since only upper limits can be provided, because no measurable emission is found. In the case of the Parkinson (1975) AR observation, the EM is at least as steep as ≃T-14 above 3 MK.
The fact that the slope in the EM distribution above 3 MK in quiescent active regions is extremely steep has important implications for coronal heating modelling. Models that involve impulsive heating (such as caused by nanoflares) normally require the presence of some hot plasma formed at or above 5 MK. We have found no evidence of such plasma in the quiescent cores of the active regions considered here. We should point out, however, that low densities and non-equilibrium effects can actually make such hot plasma difficult to observe (see, e.g. Reale & Orlando 2008; Bradshaw & Klimchuk 2011).
The first comprehensive measurement of the abundances of several elements in an AR core was recently published in Del Zanna (2013b), using Hinode EIS observations and a recent calibration (Del Zanna 2013a). The relative abundances of several key elements were obtained, including low-FIP (e.g. Fe, Si, Ca, Ni) and high-FIP (e.g. Ar, S). A FIP enhancement of 3.2 was found. It is interesting to see that most of the X-ray observations of AR cores we have studied here indicate the same FIP enhancement, considering the Fe/Ne and Fe/O relative abundances and adopting the photospheric abundances of Asplund et al. (2009) as reference values. We found similar results for five out of six different active regions. This is remarkable, although EUV spectroscopy from Hinode EIS of a few AR cores also showed a similar constancy (Del Zanna 2013a).
We note that our average FIP enhancement of 3.2 is slightly lower than the “average coronal value” of 4 suggested by Feldman et al. (1992), which is used extensively in the literature. We also note that different active regions structures have indicated a variability in abundances, and recommend the present abundances only for the hot core loops in quiescent active regions.
It is interesting to note that Ulysses measurements of slow solar wind (SSW) abundances have also shown on average FIP enhancements between 2 and 3 (von Steiger et al. 2000) during periods of solar activity, while the fast solar wind (FSW) streaming from coronal holes has nearly photospheric abundances (Gloeckler & Geiss 2007). Recent measurements from the ACE spacecraft during solar maximum and minimum conditions (Lepri et al. 2013) also indicate that at minimum the differences in composition between the slow and fast wind are not large, while the SSW during solar maximum has an FIP enhancement. For example, the Fe/O increases by a factor of two, from the FSW during minimum to the SSW during maximum.
The FIP enhancements observed both in-situ, when active regions are present, and with XUV spectroscopy of their cores, suggest a possible contribution of the hot core loops to the SSW, perhaps via the interchange reconnection mechanism, proposed by Del Zanna et al. (2011a) to explain the AR coronal outflows observed by Hinode EIS (see, e.g. Del Zanna 2008b; Doschek et al. 2008). Future close-by remote sensing and in-situ observations with Solar Orbiter will shed some light into this issue.
We are also eagerly awaiting new X-ray spectroscopy with instruments such as the Marshall Grazing Incidence X-ray Spectrograph (MaGIXS; Kobayashi et al. 2011), which will be initially proposed as a sounding rocket experiment. MaGIXS will have enhanced sensitivity and spatial resolution in the important 6–24 Å spectral region. This instrument, in combination with EUV observations, will be able to provide new accurate measurements of temperatures and chemical abundances in active regions.
Online material
Appendix A: SMM BCS light curves
Fig. A.1 Light curves in the SMM BCS Ca xix channel (data are count rates) for the active regions cores that were not selected for analysis. |
Work done within the UK APAP network, see http://www.apap-network.org
Acknowledgments
The present work was funded by STFC (UK) through the University of Cambridge DAMTP astrophysics grant. We thank the referee for useful suggestions.
References
- Aly, M. K., Evans, J. W., & Orrall, F. Q. 1962, ApJ, 136, 956 [NASA ADS] [CrossRef] [Google Scholar]
- Arnaud, M., & Raymond, J. 1992, ApJ, 398, 394 [NASA ADS] [CrossRef] [Google Scholar]
- Arnaud, M., & Rothenflug, R. 1985, A&AS, 60, 425 [NASA ADS] [Google Scholar]
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Bradshaw, S. J., & Klimchuk, J. A. 2011, ApJS, 194, 26 [NASA ADS] [CrossRef] [Google Scholar]
- Davis, J. M., Gerassimenko, M., Krieger, A. S., & Vaiana, G. S. 1975, Sol. Phys., 45, 393 [NASA ADS] [CrossRef] [Google Scholar]
- Del Zanna, G. 1999, Ph.D. Thesis, Univ. of Central Lancashire, UK [Google Scholar]
- Zanna, G. 2006, A&A, 459, 307 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G. 2008a, A&A, 481, L69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G. 2008b, A&A, 481, L49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G. 2011, A&A, 536, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G. 2012, A&A, 537, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G. 2013a, A&A, 555, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G. 2013b, A&A, 558, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G., & Ishikawa, Y. 2009, A&A, 508, 1517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G., & Mason, H. E. 2003, A&A, 406, 1089 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G., Landini, M., & Mason, H. E. 2002, A&A, 385, 968 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G., Aulanier, G., Klein, K.-L., & Török, T. 2011a, A&A, 526, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G., Mitra-Kraev, U., Bradshaw, S. J., Mason, H. E., & Asai, A. 2011b, A&A, 526, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Del Zanna, G., O’ Dwyer, B., & Mason, H. E. 2011c, A&A, 535, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dere, K. P., Landi, E., Young, P. R., et al. 2009, A&A, 498, 915 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Doschek, G. A., Warren, H. P., Mariska, J. T., et al. 2008, ApJ, 686, 1362 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, K., & Pounds, K. A. 1968, ApJ, 152, 319 [NASA ADS] [CrossRef] [Google Scholar]
- Feldman, U., Mandelbaum, P., Seely, J. F., Doschek, G. A., & Gursky, H. 1992, ApJS, 81, 387 [NASA ADS] [CrossRef] [Google Scholar]
- Gloeckler, G., & Geiss, J. 2007, Space Sci. Rev., 130, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Haugan, S. V. H. 1997, CDS software note, 47 [Google Scholar]
- Ko, Y., Doschek, G. A., Warren, H. P., & Young, P. R. 2009, ApJ, 697, 1956 [NASA ADS] [CrossRef] [Google Scholar]
- Kobayashi, K., Cirtain, J., Golub, L., et al. 2011, in SPIE Conf. Ser., 8147 [Google Scholar]
- Landi, E., Young, P. R., Dere, K. P., Del Zanna, G., & Mason, H. E. 2013, ApJ, 763, 86 [NASA ADS] [CrossRef] [Google Scholar]
- Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Lepri, S. T., Landi, E., & Zurbuchen, T. H. 2013, ApJ, 768, 94 [NASA ADS] [CrossRef] [Google Scholar]
- Liang, G. Y., & Badnell, N. R. 2010, A&A, 518, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lodders, K., Palme, H., & Gail, H.-P. 2009, Landolt Börnstein, 712 [Google Scholar]
- Mason, H. E. 1975, MNRAS, 171, 119 [NASA ADS] [Google Scholar]
- Mason, H. E., Landi, E., Pike, C. D., & Young, P. R. 1999, Sol. Phys., 189, 129 [NASA ADS] [CrossRef] [Google Scholar]
- McKenzie, D. L., & Feldman, U. 1992, ApJ, 389, 764 [NASA ADS] [CrossRef] [Google Scholar]
- Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 197 [NASA ADS] [Google Scholar]
- Meyer, J. 1985, ApJS, 57, 173 [NASA ADS] [CrossRef] [Google Scholar]
- Mohan, A., Landi, E., & Dwivedi, B. N. 2000, A&A, 364, 835 [NASA ADS] [Google Scholar]
- Monsignori Fossi, B. C., Landini, M., Thomas, R. J., & Neupert, W. M. 1994, Adv. Space Res., 14, 163 [NASA ADS] [CrossRef] [Google Scholar]
- O’Dwyer, B., Del Zanna, G., Mason, H. E., et al. 2011, A&A, 525, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- O’Dwyer, B., Del Zanna, G., & Mason, H. E. 2014, A&A, 561, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Parkinson, J. H. 1975, Sol. Phys., 42, 183 [NASA ADS] [CrossRef] [Google Scholar]
- Parkinson, J. H. 1977, A&A, 57, 185 [NASA ADS] [Google Scholar]
- Patsourakos, S., & Klimchuk, J. A. 2009, ApJ, 696, 760 [NASA ADS] [CrossRef] [Google Scholar]
- Phillips, K. J. H., Greer, C. J., Bhatia, A. K., et al. 1997, A&A, 324, 381 [NASA ADS] [Google Scholar]
- Reale, F. 2010, Liv. Rev. Sol. Phys., 7, 5 [Google Scholar]
- Reale, F., & Orlando, S. 2008, ApJ, 684, 715 [NASA ADS] [CrossRef] [Google Scholar]
- Rosner, R., Tucker, W. H., & Vaiana, G. S. 1978, ApJ, 220, 643 [NASA ADS] [CrossRef] [Google Scholar]
- Saba, J. L. R., & Strong, K. T. 1993, Adv. Space Res., 13, 391 [NASA ADS] [CrossRef] [Google Scholar]
- Saba, J. L. R., Schmelz, J. T., Bhatia, A. K., & Strong, K. T. 1999, ApJ, 510, 1064 [NASA ADS] [CrossRef] [Google Scholar]
- Schmelz, J. T., & Pathak, S. 2012, ApJ, 756, 126 [NASA ADS] [CrossRef] [Google Scholar]
- Schmelz, J. T., Saba, J. L. R., Ghosh, D., & Strong, K. T. 1996, ApJ, 473, 519 [NASA ADS] [CrossRef] [Google Scholar]
- Schmelz, J. T., Reames, D. V., von Steiger, R., & Basu, S. 2012, ApJ, 755, 33 [NASA ADS] [CrossRef] [Google Scholar]
- Strong, K. 1978, Ph.D. Thesis, University College London, UK [Google Scholar]
- Tripathi, D., Klimchuk, J. A., & Mason, H. E. 2011, ApJ, 740, 111 [NASA ADS] [CrossRef] [Google Scholar]
- von Steiger, R., Schwadron, N. A., Fisk, L. A., et al. 2000, J. Geophys. Res., 105, 27217 [NASA ADS] [CrossRef] [Google Scholar]
- Waljeski, K., Moses, D., Dere, K. P., et al. 1994, ApJ, 429, 909 [NASA ADS] [CrossRef] [Google Scholar]
- Warren, H. P., Ugarte-Urra, I., Doschek, G. A., Brooks, D. H., & Williams, D. R. 2008, ApJ, 686, L131 [NASA ADS] [CrossRef] [Google Scholar]
- Warren, H. P., Winebarger, A. R., & Brooks, D. H. 2012, ApJ, 759, 141 [NASA ADS] [CrossRef] [Google Scholar]
- Widing, K. G., & Feldman, U. 1993, ApJ, 416, 392 [NASA ADS] [CrossRef] [Google Scholar]
- Winebarger, A. R., Schmelz, J. T., Warren, H. P., Saar, S. H., & Kashyap, V. L. 2011, ApJ, 740, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Withbroe, G. L. 1975, Sol. Phys., 45, 301 [NASA ADS] [CrossRef] [Google Scholar]
- Witthoeft, M. C., Badnell, N. R., Del Zanna, G., Berrington, K. A., & Pelan, J. C. 2006, A&A, 446, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Young, P. R., Mason, H. E., Keenan, F. P., & Widing, K. G. 1997, A&A, 323, 243 [NASA ADS] [Google Scholar]
- Young, P. R., Del Zanna, G., Mason, H. E., et al. 2007, PASJ, 59, 857 [Google Scholar]
All Tables
All Figures
Fig. 1 Top: DEM for the Parkinson (1975) active region spectrum and the present elemental abundances. The ion and the wavelength (Å) is given. The numbers in parentheses are the ratios between predicted and observed intensities. The points are plotted at the effective temperature Teff, and at the theoretical vs. the observed intensity ratio multiplied by the DEM value. Bottom: same DEM, but with the points plotted at the temperature of the maximum of the G(T). |
|
In the text |
Fig. 2 Light curves in the SMM BCS Ca xix channel (data are count rates) for the selected active regions cores analysed here. The dashed lines indicate the start and end times of the FCS observations selected for the DEM analysis. |
|
In the text |
Fig. 3 KPVT magnetograms for the selected active regions analysed here. |
|
In the text |
Fig. 4 SMM/FCS AR spectrum of 1987 Dec. 13. |
|
In the text |
Fig. 5 EM loci curves for the SMM/FCS AR spectrum of 1987 Dec. 13. Top: with CHIANTI v.4.2, the Arnaud & Rothenflug (1985) and Arnaud & Raymond (1992) ion abundances, and the Schmelz et al. (1996) elemental abundances. Bottom: with the present atomic data and abundances. |
|
In the text |
Fig. 6 Top: DEM distribution for the SMM/FCS AR spectrum of 1987 Dec. 13, with the Fe abundance increased by a factor of 3.2 compared to the Asplund et al. (2009) photospheric abundance (as Fig. 1). The points are displayed at the effective temperature. ul indicates an upper limit. The labels indicate the main ion and the wavelength in Å. Bottom: same curve, with the results of the XRT_DEM inversion code (black histogram) and of the Monte Carlo XRT_DEM inversion (grey histograms). |
|
In the text |
Fig. 7 DEM distributions for the active region cores analysed here, as in Fig. 6. |
|
In the text |
Fig. A.1 Light curves in the SMM BCS Ca xix channel (data are count rates) for the active regions cores that were not selected for analysis. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.