A&A 486, L5-L8 (2008)
LETTER TO THE EDITOR
S. J. Bradshaw
Space and Atmospheric Physics, Blackett Laboratory, Imperial College London, Prince Consort Road, SW7 2BZ, London, UK
Received 26 March 2008 / Accepted 27 May 2008
Context. New observations made by the Hinode EUV Imaging Spectrometer have revealed persistent redshifts in solar active region loops in the temperature range K. The presence of redshifts, interpreted as bulk downflows, indicates that the loops are undergoing radiative cooling rather than continuous heating. This has significant consequences for current ideas regarding the physics of the ubiquitous 1 MK loops observed by instruments such as TRACE and SoHO-EIT.
Aims. The aim of this work is to assess the validity of the radiative cooling interpretation of the physics governing these loops.
Methods. This interpretation requires an enthalpy flux to balance the radiative energy loss in order to avoid catastrophic cooling. An analytical model is developed and used to find the downflows needed to drive an enthalpy flux sufficient to balance the known radiative loss rates, which are then compared with the observationally measured redshifts.
Results. The model results are found to be in good agreement with the observed redshifts.
Conclusions. Based upon the current analysis it appears that the majority of active region loops observed in the temperature range K are in the radiative cooling phase of their lifecycle.
Key words: Sun: corona - Sun: transition region - Sun: UV radiation
In a recent paper Del Zanna (2008) showed new observations made by the Hinode EUV Imaging Spectrometer (EIS, Culhane et al. 2007) of an active region (NOAA 10926) as it crossed the centre of the solar disk. Dopplergrams of the active region revealed persistent redshifts, more pronounced in cooler emission lines (20-30 km s-1 in Fe VIII, 5-10 km s-1 in Fe XII and 5 km s-1 in Fe XV), in almost all of the loop structures that were observed in the temperature range K. Observations of active regions made within this temperature range by imaging instruments such as TRACE and SoHO-EIT, which cannot resolve individual spectral lines, have in the past been interpreted in terms of static, steady-state loop structures undergoing some form of constant heating (Aschwanden et al. 2001; Aschwanden & Schrijver 2002), whereas EIS has now revealed significant dynamic activity which calls these previous interpretations into question.
The presence of redshifts, interpreted as bulk downflows, indicates that the loops are in the radiative cooling phase of their lifecycle (Bradshaw & Cargill 2005) and not the heating/thermal conduction phase that would be associated with upflows (chromospheric evaporation). This leads to a significant re-interpretation of the nature of the energy balance and the key physical processes that operate within these ubiquitous structures. If the loops are indeed in their radiative cooling phase, the following two conditions must be satisfied: (a) the downflow must drive an enthalpy flux sufficient to balance the optically-thin radiative losses from the lower corona/transition region, otherwise the loops would be radiatively unstable and undergo catastrophic cooling in contradiction to their observed lifetimes of the order of hours (Lenz et al. 1999a; Lenz et al. 1999b); and (b) the magnitude of the downflow associated with the required enthalpy flux must be consistent with the observationally measured redshifts.
In the present work an analytical model is developed and used to assess the validity of the radiative cooling interpretation. In Sect. 2 an equation for the critical downflow velocity needed for an enthalpy flux to balance the radiative losses is derived and applied to a selection of loops with varying temperature and density structures. The results are compared with Del Zanna (2008). Conclusions are drawn and work to follow is discussed in Sect. 3.
An optically thin, radiating plasma cools on a timescale given by
Let us assume that the plasma is gravitationally stratified at and above s=s0. Thus,
In order for the lower corona/transition region radiative emission to be sustained, the condition
must be satisfied in these regions. This implies, from Eq. (8),
The critical parameter in Eq. (9) is the velocity of the plasma,
Table 1: The critical velocity that the downflow must reach in order to drive an enthalpy flux sufficient to sustain the radiative emission.
Table 1 shows values of for a range of n=n(s0) and . The value of s is held fixed at cm (a height of about 2000 km for the lower corona/transition region) and T = 105.8 K (the temperature of peak Fe VIII abundance in equilibrium). Thus, the values in Table 1 are the velocities one would expect to observe from Doppler shifts in Fe VIII emission lines if the enthalpy flux sustains the radiative emission. This parameter space was chosen to allow for direct comparison with the Dopplergrams of Del Zanna (2008). Note that depends sensitively on erg cm3 s-1, the optically-thin radiative emissivity. depends most sensitively on the values selected for the element abundances, which highlights the need for accurate spectroscopic determinations of these quantities.
During radiative cooling (Bradshaw & Cargill 2005) and so the corona cools more rapidly than it drains (for example, a factor of 2 decrease in n corresponds to a factor 4 decrease in T). Table 1 shows that is weakly dependent on (and thus on the temperature structure of the loop) and strongly dependent on n; therefore, after significant cooling and only moderate draining the value of can substantially decrease. For example, consider a coronal loop of moderate temperature K and density cm-3 rising by a factor of 10 to cm-3 in the lower corona/transition region: the radiative cooling timescale in the coronal part of the loop is about h and a coronal temperature of order 106 K and a density ratio of 0.1 implies a loop length of order 105 km (in the case of gravitational stratification) which, according to Table 1, implies a draining timescale of about h. Therefore, while the loop remains bright at EUV temperatures it should not catastrophically empty and sufficient plasma remains in the corona to act as a reservoir for the energy needed to sustain the radiative emission from lower-lying cooler and denser plasma.
In order to compare the critical velocities in Table 1 with the redshifts shown in the Dopplergrams of Del Zanna (2008) (for example, his Fig. 2), some idea of the range of densities is required. There are no density maps corresponding to the Dopplergrams in Del Zanna (2008); however, Young et al. (2007) analysed Hinode-EIS observations of an active region and used a Mg VII ( T=105.8 K) diagnostic to find a density of cm-3 at the foot-point of a loop visible in the TRACE 171 Å channel and Tripathi et al. (2008) found densities of up to 1010.5 cm-3 in the core of an active region. Assuming that these properties do not vary too markedly from active region to active region then values of a few 109 cm-3 for loop foot-point densities and a few 1010 cm-3 for active region core densities would appear to be reasonable choices. The critical velocities in Table 1 that correspond to this density range are in good agreement with the redshifts measured by Del Zanna (2008). For example, in Fig. 2 of Del Zanna (2008) the redshifts in the Fe VIII emission from the active region core are between 24 and 40 km s-1, and Table 1 predicts downflow speeds of around 30 km s-1 for a density of cm-3 (1010.5 cm-3). Moving out from the active region core to regions of lower density, the redshifts in the Fe VIII emission fall to values of 10 km s-1, as also predicted by Table 1.
These results may also explain the origin and longevity of the ubiquitous active region moss (Berger et al. 1999), now thought to be the relatively cool (1 MK) emission from the upper transition regions of hot (3-10 MK) coronal loops (Martens et al. 2000). Given typical moss values of order T=106 K and cm-3 (Berger et al. 1999) then in the range 3 MK 10 MK Eq. (11) gives 3.4 km s km s-1. Such gentle downflows could explain the apparent stability of the moss emission in the imaging observations made by TRACE while also providing an explanation for its longevity. Furthermore, Gray & Kilkenny (1980) showed that thermal conduction can be strongly attenuated from its classical value (Spitzer & Härm 1953) at high temperatures, which may explain the lack of a significant role for thermal conduction in the energy balance.
An analytical model has been developed and used to assess the validity of a radiative cooling interpretation for the physics of active region loops observed in the temperature range K. These loops have previously been interpreted in terms of a static, steady-state model, in which their properties and evolution were understood in terms of an energy balance between heating and thermal conduction. In the new model the key physical processes lead to an energy balance between optically-thin radiative emission and a mechanical transport of energy (an enthalpy flux). The critical velocities predicted by the model and needed to drive an enthalpy flux sufficient to balance the radiative losses were found to be in good agreement with the observationally measured redshifts in the Dopplergrams of Del Zanna (2008). Thus, one may conclude that the majority of these loops are in the radiative cooling phase of their lifecycle.
Moreover, the radiative cooling interpretation is in agreement with Warren et al. (2002) who suggested that coronal loops are made up of many strands heated separately to X-ray temperatures and detected at EUV temperatures as they cool through the sensitivity range of the observing instrument's wavelength filters. This process also provides a natural explanation for the observed overdensity of EUV loops compared with the density predicted by static, steady-state scaling laws (Rosner et al. 1978; Serio et al. 1981), as confirmed by the numerical work of Spadaro et al. (2003).
Pneuman & Kopp (1977, 1978) reached similar conclusions to those presented here after they showed that a downflow of spicular material can provide an enthalpy flux to the transition region that is significantly greater than the energy supplied by thermal conduction from the corona.
The observational study will be extended in forthcoming work to include more active regions and density maps will be produced to accompany Dopplergrams for a larger number of emission lines. This will facilitate an improved comparison with the analytical model and allow the generality of the results and conclusions presented here to be determined.
The author is grateful to Dr. Giulio Del Zanna for bringing these observational results to his attention and to Professor Peter Cargill for insightful discussions during the preparation of this work. Thanks also to the referee for their detailed and interesting comments. The author is supported by a Post-Doctoral Fellowship awarded by the STFC.