Observing Episodic Coronal Heating Events Rooted in Chromospheric Activity
We present results of a multi-wavelength study of episodic plasma injection into the corona of AR 10942. We exploit long-exposure images of the Hinode and Transition Region and Coronal Explorer (TRACE) spacecraft to study the properties of faint, episodic, “blobs” of plasma that are propelled upward along coronal loops that are rooted in the AR plage. We find that the source location and characteristic velocities of these episodic upflow events match those expected from recent spectroscopic observations of faint coronal upflows that are associated with upper chromospheric activity, in the form of highly dynamic spicules. The analysis presented ties together observations from coronal and chromospheric spectrographs and imagers, providing more evidence of the connection of discrete coronal mass heating and injection events with their source, dynamic spicules, in the chromosphere.
Subject headings:Sun: magnetic fields—Sun: chromosphere—Sun: transition region—Sun: corona
Observing the mechanics of the process that raises the 10,000 K chromospheric plasma to several million K has eluded the community since the inference of hot coronal plasma was made at the dawn of the rocket age (Edlén 1943). Recently, analysis of data from the Hinode spacecraft (Kosugi et al. 2007) has established the connection of chromospheric and coronal heating processes (De Pontieu et al. 2009). The key to this observational-driven leap in understanding is the discovery of a second class of “spicule” (e.g., Roberts 1945; Beckers 1968). These so-called type-II spicules were revealed in the high spatial resolution observations of the chromospheric limb provided by the Solar Optical Telescope (Tsuneta et al. 2008). De Pontieu et al. (2007b) showed that they originate in strong magnetic regions, are longer (4-8Mm), and display larger upward velocities (50-150km/s) than their classical (shock-driven) “Type-I” counterparts (2-5Mm tall with speeds 10-40km/s; De Pontieu et al. 2004; Hansteen et al. 2006; De Pontieu et al. 2007a; Rouppe van der Voort et al. 2007).
De Pontieu et al. (2009) demonstrated the spatio-temporal correlation of upper chromospheric activity, in the form of Type-II spicules, and a weak component (2-5% of the peak line emission) in the blue wings of three emission lines formed in the transition region and corona (C iv 1548Å at 150,000K, Ne viii 770Å at 600,000K, and Fe xiv 274Å at 2MK). This correlation was strongest in locations of unipolar plage, as well as the supergranular network vertices of the quiet sun and coronal holes. In each case the blue-wing component, indicative of the presence of an upward moving jet, shows a nearly identical velocity distribution to that of the dynamic spicules seen in the same unipolar regions. This led De Pontieu et al. (2009) (and McIntosh & De Pontieu 2009) to suggest that these episodic heating events rooted in chromospheric activity play a significant role in filling the Sun’s upper atmosphere with hot plasma.
In this Letter we explore a long-exposure observing sequence performed by Hinode and the Transition Region and Coronal Explorer (TRACE; Handy et al. 1999) of solar active region (AR) 10942. Observations with the Hinode X-Ray Telescope (XRT; Golub et al. 2007) revealed the presence of episodic, fast, “blobs” of coronal material (MK) that appear to emanate from the plage in the active region (Sakao et al. 2007). Our motivation comes from the fact that the reported blob velocities (qualitatively) matched those of Type-II spicules and associated weak coronal upflows (measured spectroscopically) observed in the magnetic footpoints of AR coronal loops. In the sections that follow we explore the thermal connectivity and origins of these coronal blobs using an extended observational dataset, taken the day before the observations discussed by Sakao et al. (2007).
The observations were taken on 2007 February 20 with the Hinode XRT, SOT and Extreme-ultraviolet Imaging Spectrometer (EIS; Culhane et al. 2007) instruments with support from TRACE (Handy et al. 1999). The observations were part of the joint observing campaign that focused on AR 10942 as it transited the solar disk (for details, see Sakao et al. 2007). The fields of view of the Hinode instruments used here are shown in the left panel of Fig. 1. The location of the AR plage and network is shown with intensity contours on top of the TRACE 1600Å image in the right panel.
XRT observed the AR continuously from 11:52-17:37 in the Ti-Poly filter with a mean cadence of 116s and a range of exposures (from 1 to 16s). TRACE observed the AR from 15:09 to 17:59 UT in its 171Å passband with a mean exposure and cadence of 77s and 85s, respectively. EIS rastered a 128″ 512″ region from 15:10-17:22UT acquiring data in 17 spectral windows that covered a broad range of (equilibrium) formation temperatures and a 60s exposure at each step. We select the longest exposure frames (10s) in the XRT and TRACE sequences that span the time period from 15:08-17:36UT. There are 89 and 85 frames respectively in each sequence.
SOT observed this region twice in the Ca iiH filter of BFI with relatively high cadence timeseries, once with 0.128s (14:54-16:23UT; 16s cadence) and another with 0.410s exposures (11:55-13:36; 6.4s cadence). While the latter sequence is not co-temporal with the observations from the other instruments, it has exposures deep enough to allow us to isolate the faint, highly dynamic upper chromospheric emission in the BFI Ca iiH passband (see, e.g., De Pontieu et al. 2009) from the dominant photospheric contribution on the disk.
In each case we coalign the image sequences using a Fourier cross-correlation technique, performing intra-instrument coalignment in the same fashion. The SOT images of the chromosphere can be coaligned to the XRT images by using TRACE 1600 and 171 Å images with an accuracy of about 1″ although precise coalignment is not necessary - our intent is to determine and compare the velocities of the blobs (and related signatures) across chromospheric and coronal temperatures in a region, co-spatial to the pointing tolerance of EIS.
In order to analyze the blobs visible in the TRACE and XRT sequences (Fig. 2 and associated movie) we study space-time (x-t) plots at selected locations around their origin in the plage in a fashion similar to Sakao et al. (2007). The dashed line shown in the two panels of Fig. 2 (dashed line) is a sample. The data underlying the sample “slit” is extracted from the image sequence and interpolated onto the same timescale (60s) for ease of representation using a cubic interpolation scheme. We see the strong visual correspondence between the two x-t plots shown in Fig. 3, noting that the signals of the blobs in each sequence are weak (1-10DN), at most a few percent of the coronal brightness in the remainder of the AR. The TRACE 171Å blobs show a little more contrast in the region extended eastward from the plage region. Shown in each panel are a small sample of the blob velocities (straight lines in the x-t plot) showing typical values from 100 to 140km/s. The measured velocities are lower bounds - the blobs on the periphery of the plage are likely to be (highly) inclined with respect to the horizontal and so we would expect to infer a somewhat reduced speed. These x-t plots, derived velocities and blob intensities are typical of the corona above the entire plage region.
Figure 4 shows the results of performing a profile asymmetry analysis on a selection of the available EIS spectral window rasters (De Pontieu et al. 2009). This “R-B” analysis involves several steps. First we fit a single Gaussian shape to the emission line profile at each pixel to establish the line centroid. Once determined, we sum the amount of emission in narrow (24km/s wide) spectral windows symmetrically placed about that centroid in a line profile interpolated to ten times the spectral resolution. We then subtract the red and blue wing contributions to the interpolated profile (hence R-B) to make a filtergram sampling a particular velocity range, dividing the result by the factor of ten. A positive value of R-B indicates an asymmetry in the red wing of the line, which we can interpret as the signature of excess downflowing material at that velocity while, conversely, a negative value of R-B would indicate an excess of upflowing material. The rows of the figure show a region of the EIS field-of-view in the Si vii 275Å (top row), Si x 261Å (middle row), and Fe xiv 264Å (bottom row) spectral windows which contain lines formed (in equilibrium) over a relatively broad range at 0.6, 1.3 and 2.0MK, respectively (Mazzotta et al. 1998). In the right column of the figure we show the R-B profile asymmetry in those spectral windows at a range of velocities
Finally, we consider the long-exposure SOT observations of the plage region. In Figure 5 we show the raw Ca iiH intensity (top left) and background subtracted (bottom left) images at the start of the sequence. The background subtracted images (formed by subtracting a 10 pixel -1.3″ square- boxcar smoothed image from the original) allow us to isolate the dynamic type-II spicules without time-differencing the data and losing some of their temporal behavior. The background subtraction leaves only structures that are small or thin, such as spicules
The right panel of Fig. 5 shows a sample x-t plot from the background-subtracted SOT Ca iiH timeseries (the red dashed lines in the left column panels, offset to approximately the same location as that used in Fig. 3 accounting for differential rotation). The panel is speckled with short (2″) bright streaks that are slightly inclined to the horizontal, revealing the typical high velocities of type II spicules, which originate near the small (but strong) magnetic flux concentrations that comprise the plage region. For reference we show a red solid line which indicates a projected velocity of 100km/s. We can easily see that there are many streaks that are faster and not many that are slower, but it appears that there is a range of Type-II spicule velocities consistent with those determined in De Pontieu et al. (2007b). These velocities are, in turn, consistent with those observed with EIS (De Pontieu et al. 2009), XRT (Sakao et al. 2007) and TRACE (Schrijver et al. 1999). We also note the rapid recurrence of Type-II all over the plot, but in particular towards the end of the sequence just to the left of the dashed vertical line.
In the previous section we established the intensity and velocity similarities of blobs observed in image timeseries of the corona from XRT (see also Sakao et al. 2007, who call these “continuous outflows”) and TRACE with low-amplitude, high-velocity, upflows revealed by novel spectroscopic diagnosis of emission lines in the upper transition region and corona. We finally connected these blobs to dynamic spicules visible in the upper chromosphere. In each temperature regime we have noted the correspondence in location (footpoints or extending away from them), amplitude (of order a few percent of background intensity) and the range of velocities (of order 100km/s) of these events. Our data is compatible with the following scenario: collections of fine jets initialized in the chromosphere (Type-II spicules) are heated and observed as faint high-velocity upflows at the footpoints of coronal loops (with EIS spectra), which are seen as discrete plasma blobs projecting upward and outward from the footpoints along coronal loops when the angle between the line-of-sight and magnetic field direction increases (XRT/TRACE images). In other words, we suggest that the flows observed by Sakao et al. (2007) are episodic in nature (blobs), that they are similar to the blobs observed by TRACE described in Schrijver et al. (1999), and provide a direct window into the process by which hot plasma is propelled from chromospheric to coronal heights in association with dynamic spicules driven from below (De Pontieu et al. 2009).
The ubiquitous presence of faint, but strongly blueshifted plasma at the footpoints of the loop structures that are seen to carry the outflowing blobs provides strong evidence for our interpretation of the data. First, the faint blue component is observed with EIS to be at the level of 1-5% of the background intensity of the line core. Given the strong upflows of order 100 km/s, and the inevitably changing viewing angle, our scenario thus naturally predicts that coronal imaging will show faint perturbations that propagate away from the footpoints at velocities of order 100 km/s, with amplitudes of order 1-5% of the background intensity of the loops - exactly what TRACE and XRT observe. The similarity in velocities, location and amplitudes provide strong evidence for a connection between all of these phenomena.
The visibility in spectral line profiles of the spicule-associated coronal outflows will strongly depend on the background emission of the loops that carry the flows. For example, along dark, long, coronal loops that have weak background emission, the spicular outflows at high speeds may contribute more to the overall spectral line profile, so that line-centroiding will sample more of the outflowing velocity distribution of the blobs, whereas the line widths will be enhanced. This is exactly what we observe for the blob-carrying loops in the maps of Doppler velocity (line centroid) and line width of the Fe xii 195 Å and Fe xiv 264 Å lines (bottom panels of Fig 1). The line centroid outflows of 20-30 km/s in the blob-carrying loops (reported previously by Harra et al. (2008) for the Fe xii line) imply that the physical velocities associated with the RB maps (of order 100 km/s) could be higher by 20-30 km/s. Our analysis suggests that the single Gaussian approach to profile fitting of a line dominated by multiple, unresolved components, will yield an ensemble velocity that, while strong (e.g., Harra et al. 2008; Sakao et al. 2007; Doschek et al. 2008) does not reflect the full range of velocities revealed by the RB analysis.
The presence of strong spectroscopically determined upflows associated with the blobs sheds light on another issue. Schrijver et al. (1999) interpreted these blobs as episodic flows and suggested that they were driven upward from low heights by Lorentz forces acting in the strong field regions. This scenario fits in well with that proposed by De Pontieu et al. (2009). However, some of these upward propagating disturbances in TRACE time sequences have been interpreted as propagating slow-mode magneto-acoustic waves (e.g., de Moortel et al. 2002a, b; McIntosh et al. 2008; Wang et al. 2009). Because of the similarities in periods, the waves above plage regions were thought to result from p-mode leakage from the photosphere (De Pontieu et al. 2005). The wave interpretation was compatible with the lack of strong Doppler shifts (which would have been expected from significant flows). However, our discovery of high speed (if faint) quasi-periodic upflows at the footpoints suggests that at least a fraction of these “waves” may, in fact, be better interpreted as flows. This highlights some of the poorly understood issues regarding p-mode leakage. For example, given the ubiquity of the p-mode leakage process (e.g., De Pontieu et al. 2004; Hansteen et al. 2006; Jefferies et al. 2006; De Pontieu et al. 2007a; Rouppe van der Voort et al. 2007), it is surprising that 5 min “waves” are seen (with TRACE) in only a small subset of loops, and usually only on the longest coronal loops (that are faint and cool as a result, e.g., Rosner et al. 1978). In addition, many of the “waves” show, at best, quasi-periodicity and some have periods much longer than the typical p-mode spectrum (10-25 minutes). This begs the question: what portion of the wave-like phenomena observed by TRACE are the result of p-mode leakage and which are just the result of observations with long exposure times (causing apparent lengthening of the blob) and high local image contrast (because of reduced background emission) of heating events triggered quasi-periodically in the lower atmosphere? The higher signal-to-noise and wide temperature coverage offered by SDO/AIA may help us establish the relative role of waves and episodic flows.
The quasi-periodicity of these events may also provide an intriguing window into the driver of chromospheric spicules (coronal blobs) and ubiquitous Alfvénic perturbations of the corona (Tomczyk et al. 2007). Current suggestions for the formation of type-II spicules indicate a role for reconnection (De Pontieu et al. 2007b) and their rooting in mostly unipolar flux regions may favor some form of “component” reconnection mechanism (Langangen et al. 2008). It is tempting to speculate that, whatever the formation mechanism, the quasi-periodicity on timescales of 5-10 minutes may be related to granular dynamics that drive the build-up of magnetic field stresses that lead to the formation of the spicules (and blobs). Work remains to investigate how clumps of concurrent type-II spicules “merge” to form the blobs observed by TRACE and XRT and how the plasma from these apparently discrete heating and injection events merges with the pre-existing coronal plasma.
We have shown that high velocity upflow events are visible in the chromospheric and coronal footpoints of loops that carry plasma blobs observed by TRACE and XRT. These velocities correspond spatially and spectrally to weak, blue-wing asymmetries that are observed across a range of temperatures in many EIS lines. We suggest that the weak upflows are related to episodes of dynamic spicule activity. Further, we suggest that the observed coronal outflows seen in long-exposure TRACE and XRT image sequences of the corona are the result of the heating of chromospheric material to coronal temperatures in dynamic Type-II spicules and, as such, may be the signature of discrete episodic coronal heating and mass injection events. The nature of these Type-II spicules, in tracing out the roots of the coronal magnetic topology and mass loading, suggests a more appropriate (or descriptive) name for these spicules: “radices” (singular: “radix” - Latin for “root”).
- affiliation: High Altitude Observatory, National Center for Atmospheric Research, P.O. Box 3000, Boulder, CO 80307
- affiliation: Lockheed Martin Solar and Astrophysics Lab, 3251 Hanover St., Org. ADBS, Bldg. 252, Palo Alto, CA 94304
- The movies in the online edition of the journal show the complete set of R-B measures in each line.
- The resulting scale of the background subtracted image is 50 Data Numbers where the original intensity ranged from 300 to 1500 DN.
- Beckers, J. M., 1968, Sol. Phys., 3, 367
- Culhane, J. L., et al., 2007, Sol. Phys., 243, 19
- de Moortel, I., et al., 2002a, Solar Phys., 209, 61
- de Moortel, I., et al., 2002b, Solar Phys., 209, 89
- De Pontieu, B., Erdelyi, R., James, S.P., 2004, Nature, 430, 536
- De Pontieu, B., Erdelyi, R., De Moortel, I., 2005, ApJ, 624, 61
- De Pontieu, B., et al., 2007a, ApJ, 655, 624
- De Pontieu, B., et al., 2007b, Pub. Ast. Soc. Jap., 59, 655
- De Pontieu, B., et al. 2007c, Science, 318, 1574
- De Pontieu, B., McIntosh, S. W., et al. 2009, 702, 1
- Doschek, G. A., et al., 2008, ApJ, 686, 1362
- Edlén, B., 1943, Zeits. für Astrophys., 22, 30
- Golub, L., et al., 2007, Sol. Phys.., 243, 63
- Handy, B. N., et al. 1999, Sol. Phys., 187, 229
- Hansteen, V. H., et al. 2006, ApJ, 647, L73
- Harra, L. K., et al. 2008, ApJ, 676, L147
- Jefferies, S.M., et al., 2006, ApJ, 648, L151
- Kosugi, T., et al., 2007, Sol. Phys., 243, 3
- Langangen, Ø, et al., 2008, ApJ, 679, 167
- Mazzotta, P., et al., 1998, A&ASupp., 133, 403
- McIntosh, S. W., et al. 2007, ApJ, 654, 650
- McIntosh, S. W., De Pontieu, B. & Tomczyk, S., 2008, Sol. Phys., 252, 321
- McIntosh, S. W., De Pontieu, B., 2009, submitted ApJ[Available on request!]
- Roberts, W. O., 1945, ApJ, 101, 136
- Rosner, R., Tucker, W. H., & Vaiana, G. S. 1978, ApJ, 220, 643
- Rouppe van der Voort, L., et al., 2007, ApJ, 660, L169
- Sakao, T., et al., 2007, Science, 318, 1585
- Sakamoto, Y., Tsuneta, S., & Vekstein, G., 2008, ApJ, 689, 1421
- Schrijver, C. J., et al., 1999, Sol. Phys., 187, 261
- Schrijver, C. J. & Zwaan, C., 2000, “Solar and Stellar Magnetic Activity”, Cambridge
- Tomczyk, S., et al. 2007, Science, 317, 1192
- Tsuneta, S., et al., 2008, Sol. Phys., 249, 167
- Wang, T. J., Ofman, L., Davila J. M., 2009, ApJ, 696, 1448