© ESO, 2013

1. Introduction

Prominences consist of relatively extended, cool, and overdense plasma seen in the lower corona above the solar limb (Martin 1998; Tandberg-Hanssen 1995; Mackay et al. 2010). Their structure and composition are exceedingly complicated. The plasma mainly resides in highly tangled magnetic fields (van Ballegooijen & Cranmer 2010). Coronal cavities are often observed to have cooler prominence plasma at their bases (Hudson et al. 1999; Gibson et al. 2006; Régnier et al. 2011). The cavity is a region of relatively low-density high-temperature plasma (Gibson et al. 2010; Habbal et al. 2010). Berger et al. (2011) have proposed that the prominence and its cavity are a form of magneto-thermal convective structure, macroscopically stable, but internally in a constant state of ubiquitous motions (Low et al. 2012a,b).

Fig. 1 Region around the tornado on 25 September 2011: a) EUVI-A 195 Å; b) AIA 171 Å images. On the EUVI image, the black dashed line shows the position of the AIA limb, and the long diagonal arrow is the epipolar line for the prominence position marked with + on the AIA image. The orange dashed line outlines the position of the southern edge of the active region (AR) corona.

Viewed at the limb, quiescent prominences mostly appear as curtains of vertical thread-like structures (Berger et al. 2008). Occasionally they look like tornados with rotations along, or of, their magnetic structures (Pettit 1943; Liggett & Zirin 1984). At first glance it is not obvious why tornados do not erupt but continue to rotate for several hours before quietening down. Recently, images from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamic Observatory (SDO) have allowed the study of solar tornados (Li et al. 2012; Su et al. 2012; Liu et al. 2012a). The driving mechanism may be a coupling and expansion of a twisted flux rope into the coronal cavity (Li et al. 2012) and/or be related to photospheric vortices at the footpoint of the tornado (Attie et al. 2009; Wedemeyer-Böhm et al. 2012; Su et al. 2012).

We study the impressive solar tornado observed by the AIA on 25 September 2011. Li et al. (2012) described the formation and disappearance of the related prominence by analysing SDO data from 24−26 September. The evolution of the tornado was attributed to the expansion of helical structures into the cavity. But the reason behind the expansion remained an open question.

Here we analyse observations of the mechanism that led to the expansion of the prominence. We show with Solar TErrestrial RElations Observatory (STEREO) Extreme UltraViolet Imager (EUVI) observations that three strong flares in a neighbouring active region coincided with phases of the tornado activation. Each flare was associated with a coronal mass ejection (CME) and an EUV wave. Particularly after the third flare, a slow but significant expansion of the overlying prominence cavity was observed together with more rapid rotation at the top of the prominence. We describe the observations in Sect. 2. The evolution of the tornado, its relationship to the solar flares and the EUV waves are described in Sect. 3. In Sect. 4, we summarize our observations and speculate on the link between the flares and the cavity expansion.

2. Observations

The solar tornado on 25 September 2011 was observed on the southwest limb by SDO/AIA (Lemen et al. 2012). In images from EUVI on the Ahead spacecraft of STEREO (Howard et al. 2008), it appeared on the solar disk around 15° E. The separation angle between SDO and STEREO-A was 103°. We studied the tornado with combined observations from the two directions.

The AIA takes high spatial resolution (0.6′′ pixel^-1) full-disk images with a cadence of 12 s. For the analysis, we selected images from the 171 Å channel, which is centred on the Fe ix line formed around 0.63 MK. Images from this channel show the hot, bright, emitting and the cold, dense, absorbing parts of the prominence, as well as the dark, low-density cavity and the surrounding coronal loops. To enhance the visibility of the cavity and faint loops, we removed an average background from the 171 Å images by taking the median of two months of data – September and October, 2011.

During the period of 01:30−13:30 UT, SDO was in eclipse from 06:02-07:13 UT. This gap is covered by the SWAP instrument on PROBA-2 (Halain et al. 2010). SWAP provides a full Sun 174 Å image every 2−3 min with a spatial resolution of 3.16′′ pixel^-1. The images are not as detailed as the AIA images, but they are essential for checking the behaviour of the prominence during the data gap.

The EUVI-A 195 Å images have a time cadence of 5 min and a resolution of 1.6′′ pixel^-1. The 195 Å emission is mainly Fe xii formed at  ~1.2 MK. To confirm the prominence position in the EUVI images, we obtained the three-dimensional coordinates of the prominence with the routine SCC$\begin{array}{}\end{array}$MEASURE (Thompson 2009) available in the SolarSoft library.

2.1. Overview

In Fig. 1 we show the EUVI-A 195 Å and AIA 171 Å images of the filament/prominence and the active region NOAA 11303 (marked with an arrow) as seen from the two angles. The position of the cross on the stem of the prominence lies along the epipolar line drawn as a long black arrow on the EUVI image. The dark lane crossed by the epipolar line is the filament channel and the tornado is the bright region that coincides with the epipolar line and the filament channel, indicated with a black arrow. The core of the active region is separated by 300′′ from the filament channel. However, the southern edge of the active-region corona (outlined by the orange dashed line) reaches to within 50′′ of the filament channel.

Fig. 2 EUVI-A 195 Å running-ratio images showing the onset of the three flares (a), b), and c)) and their associated EUV waves (d), e), and f)). In a) the epipolar line shown in Fig. 1a is drawn as a long white arrow, and the solar limb is drawn as a dashed white line. The small white arrow in all images points to the tornado site. The dot-dashed/dashed lines in b) mark the positions of time series shown in Fig. 4, and the orange dashed line along the EUV wave front in e) is the position of the time series in Fig. 3.

2.2. The three flares

There were flares at approximately 02:45, 07:00, 09:40 UT from the nearby active region. They were all associated with CMEs and EUV waves. The third was GOES class M1.4 according to the NOAA flare catalogue¹. The GOES class of the first two are not given in the catalogue. We therefore checked the hard X-ray quicklook images from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). Unfortunately, both had their peaks during RHESSI data gaps. The first RHESSI image in the flares’ decay phase showed that these flares were the brightest RHESSI sources. At the time of the first and second flares the GOES 1−8 Å fluxes reached the M4.4 and M1.0 level, so it is possible that these were M-class flares as well.

The flares’ onsets and related EUV waves can be seen in the 195 Å running ratio images shown in Fig. 2. The running ratio images are the log of the intensity at the time shown in the images divided by the intensity 5 min earlier. Because this reflects relative changes, the faint EUV waves show up against bright inactive active regions. The faint EUV wavefronts sweep over the filament channel about half an hour after the flare (bottom row). The small white arrows point to the tornado site, which showed increased activity along the filament channel after the flares.

When we take a time series along the orange dashed line drawn in Fig. 2e, we are able to see the relation between the EUV waves and the filament activity (Fig. 3). In this time series image, the waves associated with the three flares are marked by F1, F2, and F3 and a red arrow points to the filament activity. The filament activity associated with F2 increased before the arrival of the EUV wave. Further more, each EUV wave front was followed by a series of oscillations with a period of 20 min at the edge of the active region and along the filament channel. Similar post-EUV cavity oscillations, seen at the limb in the AIA images, have been reported by Liu et al. (2012b). EUVI images have a much lower cadence than those from the AIA, and may not be resolving the wave trains.

As mentioned above, the EUV wave reached the filament channel after the filament activation. To investigate whether the activation could have been triggered by the flare, we took time series through the flare sites (vertical line Fig. 2b) and along the filament channel across the activated filament (horizontal line in Fig. 2b). We observe that there are perturbations in the filament plasma immediately after each flare (yellow arrows in Fig. 4b). Large-scale changes along the filament channel started after F2 and increased further after the F3.

The reaction of the prominence to the flares is illustrated in Fig. 7. The positions of the three time series are drawn in Fig. 6, an AIA image of the active region and the tornado. The oscillations of the prominence at the time of F1 (Fig. 7a) are consistent with a flare trigger from the north since it moved first towards the pole. The wobbling in the prominence after the flare is also clearly visible in movie6a (Fig. 6).

Fig. 3 EUVI 195 Å running-ratio space-time image along the dashed line in Fig. 2e. The three flares are labelled F1, F2, and F3. The red arrow points to the filament.

The second flare, F2, was at  ~06:45 UT when SDO was in eclipse. Since this phase is important for understanding the tornado, we checked the SWAP data to see if the prominence activation started before F2. The SWAP data are shown in the movie attached to Fig. 5. Although the SWAP images are not as sharp as the AIA images, they show quite clearly that increased prominence activity occurred after F2 and that the first AIA images after eclipse caught the expansion of the flare loops and the rapid growth in prominence activity. The middle image in Fig. 7b, taken along the dashed line in Fig. 6, shows an arm of prominence plasma that reaches out in the direction of the active region. There is also a very slight expansion of the cavity boundary towards the active region.

Fig. 4 EUVI 195 Å running-difference time-series along a) the vertical and b) the horizontal line in Fig. 2b. The red dashed lines are drawn at the times of flare F1, F2, and F3. In b) the yellow arrows point to the changes in filament structure after the flares.

Fig. 5 SWAP 174 Å intensity image from 25 September 2011. The arrows show the positions of the cavity and the active region (AR). This frame is taken from the movie “MOVIE5.mp4”.

Fig. 6 AIA 171 Å intensity image. The diagonal long arrows show the positions of the time-series images through the prominence/tornado, cavity, active region hot loop (AR loop) and active-region (AR) shown in Fig. 7. This is the frame from the movie “MOVIE6b.mp4”. “MOVIE6a.mp4” shows the evolution of this region at earlier times.

Fig. 7 AIA 171 Å intensity time-series along the diagonal arrows in Fig. 6: a) oscillations of the prominence stem triggered by F1 taken along the bottom white arrow; b) activity triggered by F2 and F3 taken along the middle white arrow; c) tornado activity and cavity changes triggered by F2 and F3 taken along the top arrow. The arrow at the time of the flare (~09.40 UT) points to the movement of the cavity boundary.

The biggest change in the cavity occurred after F3 when the bright and complex top part of the prominence started to rotate faster. Movie6b, attached to Fig. 6, shows the evolution of the tornado and flare from F2 to a couple of hours after F3. The time series in Fig. 7c, taken along the top white arrow in Fig. 6, shows the main features. The flare erupted at 09:20 UT. The active-region loops started to expand at 09:40:48 UT. At 09:50:13 UT there was a sharp contraction, approximately 10 Mm, of the observed cavity. The cavity started to expand back towards the active region after  ~10:00 UT. It grew from 140 Mm to 167 Mm in 3.5 h. At the same time, the prominence grew in height and the bright and complex head of the prominence started to rotate faster, forming the tornado. The rotations at the top of the prominence/tornado lasted  ~3.5 h with a speed from 55 km s^-1 to 95 km s^-1 (Li et al. 2012).

3. Discussion

A solar tornado was observed by the AIA and EUVI-A on 25 September 2011. The tornado was reported by Li et al. (2012), who attributed it to the growth of a helical prominence system. The reason for the prominence growth remained unanswered. After careful analysis of the AIA full-disk images, we noticed that AR 11303, to the north of the prominence cavity, was flaring and suspected it of having influenced the prominence in some way. We observed three flares, which were all associated with EUV waves that swept over the prominence cavity. Each time the cavity was buffeted by the wave, the prominence plasma became more active. After the third flare a tornado had developed at the top of the prominence.

After the first flare, there were clear oscillations along the prominence stem. The second flare caused plasma to move along arm-like extensions projecting out from the stem in the direction of the active region. These extensions appear to have subsequently rotated back towards the prominence before the third flare. After the third flare the tornado became very strong and was stable for about 3.5 h. The third flare seemed to have the biggest impact on the surrounding plasma, or at least it produced the most visible (at 195 Å) EUV wave. It also caused large active-region loops to move towards the edge of the prominence cavity. Afterwards, the loops swayed back and the cavity started to slowly (~2.5 km s^-1) but visibly expand over the next 3.5 h. During this time the prominence rose into the cavity and created the tornado.

To understand the relationship between the flare, the prominence, the cavity, and the tornado, several effects need to be considered. The expansion of the cavity is clearly associated with the rise of the prominence plasma. One possible explanation for the expansion is that the loss of free magnetic energy from the active region by flares and CMEs (Zhang & Low 2005) resulted in a contraction of the active region field, the Hudson effect (Hudson 2000; Zhang & Low 2003; Janse & Low 2007). Subsequently, the surrounding fields, including the cavity field, expanded to fill the vacated space.

An idealized illustration of the Hudson effect, in terms of the force-free equilibrium between two bipolar fields before and after one of the fields has lost its free energy, is given in Fig. 8. Initially, the active-region field on the right has excess free energy, and the cavity field on the left is taken to be potential to simplify the calculation. During the eruption, the active-region field loses free energy, decreasing its pressure since the magnetic energy density B² of a field B is also the magnetic pressure. This leads to an expansion of the cavity field to restore the pressure balance in the system.

We consider the 2D Cartesian domain 0 < y < π, 0 < z < ∞, where x is an ignorable coordinate, taking the boundary to be a rigid perfect electrical conductor for simplicity. The initial state in Fig. 8a is a continuous global solution for the magnetic field, B = B_initial of the force-free equations $\mathrm{\nabla }\mathrm{\times }{B}\mathrm{=}\mathit{\alpha }{B}\mathrm{=}\mathrm{0}\mathrm{;} \mathrm{\nabla }\mathrm{·}{B}\mathrm{=}\mathrm{0}\mathit{,}$(1)which describe a bipolar (α = 0) potential field that represents the cavity that occupies the partial domain 0 < y < 2π/3,0 < z < ∞. This field continues into a constant-α force-free field that represents the active region that occupies the complementary partial domain 2π/3 < y < π,0 < z < ∞. Their interface y = 2π/3, is in force balance because the fields on its two sides exert equal magnetic pressure. Suppose the electric current of the constant-α force-free field is removed from the partial domain 2π/3 < y < π,0 < z < ∞ to represent a flare-like loss of free energy. The distribution of the normal component of the field along the rigid boundary and the magnetic foot-points on the base z = 0 cannot change. This boundary condition then determines the unique end-state, B_end, with α = 0 everywhere in the domain (Fig. 8b). The former bipolar field on the right has contracted downward and has withdrawn into a partial domain of a finite height, while the other bipolar field now occupies the whole infinite space above that height. The fluid-interface between the two fields is now an arc of finite length. These mathematical solutions, given in the appendix, serve to conceptually make specific a basic effect in the complex processes of a real energy release.

In addition to the active-region contraction, the CME/flare eruptions generated hydromagnetic disturbances, seen as EUV waves, which perturbed the surrounding low corona. They evidently caused swaying in the prominence and oscillations along the cavity boundary. These perturbations may have been intense enough to partially destabilize the cavity-prominence system, thereby contributing to, or possibly even initializing, the tornado-like activity.

Fig. 8 Idealized 2D Cartesian magnetic fields of the cavity and active region a) before and b) after the flare, in a domain 0 < y < π,0 < z < ∞, where z denotes the coronal height.

Our observational study supports the interpretation that the CME/flare activities in the active region had a causal relationship with the nearby tornado. Clearly, the tornado was dynamically associated with the expansion of the prominence’s helical field and its cavity. It is an open question whether the Hudson effect drove the cavity expansion and the helical field was responding dynamically to the expansion, or if the helical field lost its meta-stability as the result of hydromagnetic disturbances from the CME/flares from the active region. This physical question is worthy of an MHD time-dependent simulation. Evidence of the Hudson effect has been reported in recent solar observational studies (e.g. Sun et al. 2012; Liu et al. 2009). Our study shows the importance of this effect in relation to CMEs, flares, and prominences, and provides observational motivation for theoretical MHD modeling to address the specific questions posed by this relationship.

Online material

Movie of Fig.5 Access here

Movie a of Fig.6 Access here

Movie b of Fig.6 Access here

Acknowledgments

We are obliged to the SDO/AIA, STEREO/EUVI and SWAP teams. N.K.P. acknowledges the facilities provided by the MPS. We are thankful to Joan Burkepile, Thomas Berger, and Wei Liu for their useful inputs. The National Center for Atmospheric Research is sponsored by the US National Science Foundation.

References

Appendix A: Mathematical illustration of the Hudson effect

We give the mathematical solutions for 2D magnetic fields in Fig. 8, which are straightforward to construct by standard techniques. Use the representation ${B}\mathrm{=}\left[\mathit{Q,}\frac{\mathit{\partial A}}{\mathit{\partial z}}\mathit{,}\mathrm{-}\frac{\mathit{\partial A}}{\mathit{\partial y}}\right]$(A.1)in terms of the scalar flux-function A(y,z) and the component B_x = Q(y,z) corresponding to the ignorable coordinate x. The initial field B_initial in Fig. 8a is given by A = A_initial and Q = Q_initial, where $\begin{array}{ccc}& & {\mathit{A}}_{\mathrm{initial}}\mathrm{=}{\mathit{A}}_{\mathrm{pot}}\mathrm{=}\mathrm{2}\sin \left(\frac{\mathrm{3}}{\mathrm{2}}\mathit{y}\right)\exp \left(\mathrm{-}\frac{\mathrm{3}}{\mathrm{2}}\mathit{z}\right)\mathrm{;}\\ & & {\mathit{Q}}_{\mathrm{initial}}\mathrm{\equiv }\mathrm{0}\mathit{,} \mathrm{in} \mathrm{0}\mathit{<}\mathit{y}\mathit{<}\mathrm{2}\mathit{\pi }\mathit{/}\mathrm{3}\mathit{,}\mathrm{0}\mathit{<}\mathit{z}\mathit{<}\mathrm{\infty }\\ & & {\mathit{A}}_{\mathrm{initial}}\mathrm{=}{\mathit{A}}_{\mathrm{fff}}\mathrm{=}\mathrm{-}\sin \left(\mathrm{3}\mathit{y}\right)\exp \left(\mathrm{-}\frac{\mathrm{3}}{\mathrm{2}}\mathit{z}\right)\mathrm{;}\\ & & {\mathit{Q}}_{\mathrm{initial}}\mathrm{=}\sqrt{\frac{\mathrm{27}}{\mathrm{3}}}\mathit{A,} \mathrm{in} \mathrm{2}\mathit{\pi }\mathit{/}\mathrm{3}\mathit{<}\mathit{y}\mathit{<}\mathit{\pi ,}\mathrm{0}\mathit{<}\mathit{z}\mathit{<}\mathrm{\infty }\mathit{.}\end{array}$\arraycolsep1.75ptThe end-state field B_end has potential field everywhere with A = A_end, where $\begin{array}{ccc}& & {\mathit{A}}_{\mathrm{end}}\mathrm{=}\sum _{\mathit{n}\mathrm{=}\mathrm{1}}^{\mathrm{\infty }}{\mathit{a}}_{\mathit{n}}\sin \mathit{ny}\exp \mathrm{(}\mathrm{-}\mathit{nz}\mathrm{)}\mathit{,}\\ & & {\mathit{a}}_{\mathrm{3}}\mathrm{=}\frac{\mathrm{1}}{\mathrm{3}} \mathrm{and} {\mathit{a}}_{\mathit{n}}\mathrm{=}\frac{\mathrm{2}\mathrm{\times }{\mathrm{3}}^{\mathrm{4}}}{\mathit{\pi }}\sin \frac{\mathrm{2}\mathit{\pi n}}{\mathrm{3}}\frac{\mathrm{1}}{\mathrm{(}\mathrm{4}{\mathit{n}}^{\mathrm{2}}\mathrm{-}\mathrm{9}\mathrm{)}\mathrm{(}{\mathit{n}}^{\mathrm{2}}\mathrm{-}\mathrm{9}\mathrm{)}}\mathit{,} \mathrm{for} \mathit{n}\ne \mathrm{3}\mathit{,} \end{array}$with Q ≡ 0. Both B_initial and B_end have the same distribution of normal field component along the boundary of the domain 0 < y < π, 0 < z < ∞.

All Figures

	Fig. 1 Region around the tornado on 25 September 2011: a) EUVI-A 195 Å; b) AIA 171 Å images. On the EUVI image, the black dashed line shows the position of the AIA limb, and the long diagonal arrow is the epipolar line for the prominence position marked with + on the AIA image. The orange dashed line outlines the position of the southern edge of the active region (AR) corona.
In the text

Fig. 2 EUVI-A 195 Å running-ratio images showing the onset of the three flares (a), b), and c)) and their associated EUV waves (d), e), and f)). In a) the epipolar line shown in Fig. 1a is drawn as a long white arrow, and the solar limb is drawn as a dashed white line. The small white arrow in all images points to the tornado site. The dot-dashed/dashed lines in b) mark the positions of time series shown in Fig. 4, and the orange dashed line along the EUV wave front in e) is the position of the time series in Fig. 3.

In the text

	Fig. 3 EUVI 195 Å running-ratio space-time image along the dashed line in Fig. 2e. The three flares are labelled F1, F2, and F3. The red arrow points to the filament.
In the text

	Fig. 4 EUVI 195 Å running-difference time-series along a) the vertical and b) the horizontal line in Fig. 2b. The red dashed lines are drawn at the times of flare F1, F2, and F3. In b) the yellow arrows point to the changes in filament structure after the flares.
In the text

	Fig. 5 SWAP 174 Å intensity image from 25 September 2011. The arrows show the positions of the cavity and the active region (AR). This frame is taken from the movie “MOVIE5.mp4”.
In the text

	Fig. 6 AIA 171 Å intensity image. The diagonal long arrows show the positions of the time-series images through the prominence/tornado, cavity, active region hot loop (AR loop) and active-region (AR) shown in Fig. 7. This is the frame from the movie “MOVIE6b.mp4”. “MOVIE6a.mp4” shows the evolution of this region at earlier times.
In the text

Fig. 7 AIA 171 Å intensity time-series along the diagonal arrows in Fig. 6: a) oscillations of the prominence stem triggered by F1 taken along the bottom white arrow; b) activity triggered by F2 and F3 taken along the middle white arrow; c) tornado activity and cavity changes triggered by F2 and F3 taken along the top arrow. The arrow at the time of the flare (~09.40 UT) points to the movement of the cavity boundary.

In the text

	Fig. 8 Idealized 2D Cartesian magnetic fields of the cavity and active region a) before and b) after the flare, in a domain 0 < y < π,0 < z < ∞, where z denotes the coronal height.
In the text