Free Access
Issue
A&A
Volume 583, November 2015
Article Number A47
Number of page(s) 9
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201323284
Published online 26 October 2015

© ESO, 2015

1. Introduction

On 2013, September 29 a large solar filament rose smoothly, without any previous flare, vigorous dynamics, or significant brightening arcade. This filament lifted and afterwards created a two-ribbon flare, and its X-ray flux reached the category C1.2 by GOES 1–8 Å1. The eruption led to an Earth-directed halo coronal mass ejection (CME). The interplanetary plasma interacted with the Earth’s magnetic field creating a moderate geomagnetic storm with effects on the ground (Dst geomagnetic index reached –75 nT, and storm classification follows Gonzalez et al. 1994) on 2013 October 22. However, since the interplanetary medium and terrestrial effects are beyond the scope of this paper, they will be treated elsewhere. This eruption went unnoticed for flare early alerts and predictions and it produced a geomagnetic storm. This case is a very pertinent example, which can be used to obtain input for simulations and for solar and space weather predictions and warnings.

thumbnail Fig. 1

Context figures of the filament. The grid spaces mark 15°. Left: image in AIA 193 Å, on 2013 Sept. 29, at 00:00 UT. Image is clipped from 50 to 5000 DN. The filament is located on 20\hbox{$\deg$}N 15\hbox{$\deg$}W, with two coronal holes surrounding it on coordinates 40\hbox{$\deg$}N 20\hbox{$\deg$}W (CH1) and 0\hbox{$\deg$}N 7\hbox{$\deg$}W (CH2). Right: context image showing a derotated AIA 304 Å image displaying the filament, and the positive polarity (in green) and negative (blue). Different colour shades are explained in caption of Fig. 2. The bipolar flux emergence occurrence is located at [250′′, 200′′], enclosed in a white square. A movie available online shows the temporal evolution of the whole visible hemisphere from 21:00 to 22:09 UT in the AIA 304 Å channel. In the movie, the FOV of the right panel of this figure is marked by a black rectangle, and the emergence by a white square. The green and blue contours have the same meaning as in the right panel.

In this paper, we focus on the solar source and trigger of this filament eruption, which seems to be a supergranular-sized magnetic flux emergence occurrence underneath the filament. We investigate the reasons for this kind of an undisturbed rise using Solar Dynamics Observatory (SDO) data by means of the Atmospheric Imaging Assembly (AIA) at various EUV wavelengths (AIA/SDO 304 Å, 193 Å, 211Å and 94 Å). Also, to explain the event we obtained a close-up using Helioseismic and Magnetic Imager (HMI) magnetograms, continuum, and linear polarization images. We also took advantage of the fully inverted HMI data in this work.

In the review of Linker et al. (2003) it is stated that only 20% of the CMEs are associated with a large flare. Flareless CMEs have been reported by Song et al. (2013). Gosain et al. (2012) state that most of the filament eruptions related to CMEs pass undetected because of their low emission and, consequently, are not associated with flares. Therefore, trusting only in flare occurrence may be misleading and the region can actually be CME-productive.

Some very recent and comprehensive summaries of the various models that can trigger the eruption of filaments and CMEs are Aulanier et al. (2010), Schmieder et al. (2013), Aulanier (2014), Parenti (2014) and references therein; the “magnetic breakout model” (Antiochos 1998; Antiochos et al. 1999); a highly twisted flux rope, as in the “torus instability” (Török & Kliem 2003, 2005; Kliem & Török 2006), different instabilities, e.g. MHD (Klimchuk & Sturrock 1989; Sturrock 1989), catastrophic loss of equilibrium (Forbes & Isenberg 1991); approaching polarities (Forbes & Priest 1995), moving features and magnetic cancellation (van Ballegooijen & Martens 1989; Forbes & Isenberg 1991); and the “tether cutting” model, implying reconnection in the arcades below the flux rope (Moore & Roumeliotis 1992; Moore et al. 2001). Considering the magnetic flux emergence, filament eruptions due to flux emergence in the surroundings have been analysed in simulations by Chen & Shibata (2000). The magnetic configuration of active and quiescent filaments and their CME productivity have been studied in Feynman & Martin (1995). Slow and fast eruptive phases have been studied by Sterling & Moore (2005), Sterling et al. (2007). We consider the observational features of these models to choose the most plausible model for this phenomenon.

More particularly, different works are revised in the review by Chen (2011): Zhang et al. (2008) studied events of flux emergence and cancellation 12 h prior to a CME initiation, finding that 91% are related. Wang & Sheeley (1999) also found filament eruptions close to flux emergence regions. Jing et al. (2004) detected that more than half of the eruptive filaments were related to flux emergence. Zhang et al. (2001) related relatively small flux cancellation areas with the filament eruption of the Bastille Day event.

However, many of these works were qualitative, only with the presence and closeness of the flux emergence to relate with other magnitudes. Here we aim to describe the filament quantitatively and to assess its instabilities, and we describe how a flux emergence can influence a filament lift-off.

This paper is organized as follows: first, the observations are described in Sect. 2, describing the flux emergence and filament dynamics, along with the surrounding coronal holes. All calculations are included in Sect. 2. The discussion and conclusions are presented in Sects. 3 and 4, respectively.

thumbnail Fig. 2

Top row: frame sequence of AIA 304 Å images of a filament’s barb (darker areas) with LOS-magnetic field superimposed. Different shades of green (from lighter to darker) contours mark different levels, as 90, 200, 500, 600, and 700 G, while a gamut of blue contours mark –90, –200, –500, –700, –900 G for the negative polarities. The frame size is 100′′× 100′′. Sequence timing aims to show the magnetic flux emergence (panel 1), pre-flare (panel 2), flare (panel 3), and decay (panel 4), from left to right. Bottom row: same temporal sequence in AIA 211 and 94 Å. An associated movie is available online.

2. Observations of the emerging region and filament lift-off

We used data provided by the space solar facility Solar Dynamics Observatory (SDO, Pesnell et al. 2012). Data were acquired with the Atmospheric Imaging Assembly (AIA, Lemen et al. 2012) and the Helioseismic and Magnetic Imager (HMI, Scherrer et al. 2012). The AIA multi-wavelength data comprise observations in 304 and 193 Å, which sample the high chromosphere, transition region, and corona, respectively. We also included AIA 211 Å data, which are useful as an intermediate wavelength measurement since the filament opacity, compared to 304 Å, is somewhat reduced. The AIA 94 Å data also probe the inner corona. The HMI data sample the Fe i line at 6173 Å to generate line-of-sight (hereafter, LOS) photospheric magnetograms. This line samples a photospheric height around ~100 km (Fleck et al. 2011). The data cadence for the general study is chosen to be 12 min from September 29 to October 1 on both AIA and HMI instruments. In other particular cases of the study, the cadence is reduced to 1 min or increased up to 180 min (3 h). Images are co-spatial and quasi-synchronized (within seconds). Pixel sampling for both instruments is 0′′6.

All images are reduced through the SolarSoftWare (SSW) read_sdo.pro and aia_prep.pro routines. Next, frames are corrected from rotation and aligned.

Figure 1 shows context images. The left panel indicates the whole solar hemisphere, showing the filament and two surrounding coronal holes in AIA 193 Å. The right panel shows the filament in AIA 304 Å as a dark, inverse-S-like, elongated structure located between two opposite polarities of NOAA AR11850 and the following active region. The filament was visible close to the trailing boundary of AR 11850, with its southern part close to its positive facular polarity, and the spine located at the polarity inversion line (PIL), with a facular negative polarity on the left and another positive on the right. As seen in AIA 304 and 193 Å, in close to five days prior to the eruption, the filament developed from a short structure (Sept. 25, 21:00 UT) to a inverse-S-shaped elongated in just one day. This fully developed filament evolved in a dark area in 193 Å, which is an open field region. The length was about ~0.40 R on Sept. 25 to ~0.63 R two days later. This filament showed very active mass motion along the spine and barb. The southern part is more elevated than the northern part at this stage. The heights look similar at Sept. 28, 10:00 UT. This situation seems to reverse on Sept. 28 19:00 UT. The magnetic flux emergence started around Sept. 29 at 00:00 UT and the appearance of pushing the filament from below is conspicuous. This emergence displayed small pores in continuum. There was a positive polarity where one of the barbs was located, and this region seems to be attaching the filament to the surface. The apex of the filament rose, placed somewhat northern from the barb, around 21:00 UT. Eventually the filament developed a two-ribbon flare on Sept. 29 at 21:43 UT, which led to a maximum flux of X-ray of C1.2 in GOES. The southern leg was detached early in the flare, and the northern leg produced a shower of coronal rain on Sept. 29, at 23:29 UT, lasting around one hour. Thread-like features of the filament feet appeared to fall down.

The subsequent CME appeared in LASCO at 22:24 UT, with a plane-of-the-sky speed (POS) of ~600 km s-1. The halo CME is reported in the CACTUS catalog3 (Robbrecht & Berghmans 2004; Robbrecht et al. 2009). However, the reported speed in the LASCO catalog4 is 1179 km s-1. This discrepancy is due to the external CME shell, which is faster than the filament.

thumbnail Fig. 3

Left: computed magnetic flux on positive patches (solid line), on negative patches (dashed line), signed magnetic flux (dotted-dashed line), unsigned magnetic flux (dotted line). A vertical line marks the time of the flare at 21:43 UT. Right: mean values of the magnetic flux density of the positive (solid) and negative (dashed) regions.

2.1. Magnetic flux emergence

We selected AIA 304 Å images and HMI magnetograms with a cadence of 12 min, from Sept. 29 00:00 UT to Oct 1 10:00 UT. We analysed AIA 304 Å close-up images with overlaid HMI magnetograms of the aforementioned filament barb (located on ~18°N, 15°W in the left panel of Fig. 1), which appears to be the location where the photospheric magnetic flux emergence happens. The panels in Fig. 2 show the evolution, displaying selected frames5. This barb seems to be close to a positive polarity, on Sept. 29 at 00:00 UT. A large magnetic supergranular extent appeared, as seen in the HMI magnetograms. This large supergranular emergence exhibited mixed polarities at first, and later on the positive and negative polarities drifted, rearranging as ordinary network elements. Three hours later, the main negative polarity started developing a pore. The positive polarities also started developing small pores in the following hours. Six hours later, this emergence extent, consisting of positive polarities (green) on the left and negative polarities (blue) on the right, moved apart up to ~25′′ (panel 1 in Fig. 2), and reached a maximum size of ~33′′ around 13:00 UT. At this moment, the negative polarity peaked –1000 G, while the positive was about 800 G. Therefore, the magnetic field of both polarities increased from around +(–)200 to +(–)900 G in only six hours. A negative polarity on Sept. 30, 11:24 UT, emerged and drifted to the right side with other negative polarity elements. These magnetic rearrangements are visible for the whole period and the loops in AIA 304 Å mark the polarity correspondence; however, no pronounced change in polarity is noticeable in HMI magnetograms, neither shear in AIA 304 nor 193 Å data. At the end of the sequence, the polarities diminished their flux and approached each other. The magnetic emergence seems to push the barb of the filament from below during the whole event.

In the second row of panels in Fig. 2, we plot AIA 211 Å images in mainly two lines: Fe xiv, corresponding to log T ~ 6.30; and Fe xiii, corresponding to log T ~ 6.25 (O’Dwyer et al. 2010). This wavelength happened to become very adequate to observe flux emergence patches under filaments, since the signal-to-noise ratio is good, samples high temperatures, and the opacity of the filament becomes reduced but not totally suppressed.

We also present evidence of an EUV counterpart of the flux emergence in 94 Å in the third row of panels in Fig. 2. These images sample a temperature range in the corona (O’Dwyer et al. 2010) with two lines: the inner quiet corona, sampled by Fe x, corresponding to log T ~ 6.05, and the flaring range, sampled by Fe xviii, corresponding to log T ~ 6.85. The region started to increase its intensity from 6 DN to 30 DN. In O’Dwyer et al. (2010), the quiet sun intensity is synthetically defined as 6 DN. Therefore, the intensity of this emerging region is five times that of the quiet sun, and 10% of the maximum intensity in the flare of the full disk in this period.

This observational fact may suggest rapid connectivity between the photospheric longitudinal and transverse fields and the corona, since the signal starts to increase from 01:00 UT onwards, as shown in Fig. 2 and in the movie added as online material.

After the flare, we can see that opposite magnetic polarities rearrange and sometimes approach and disappear, displaying a brightening between them as a cancellation signature. In the post-flare panels, the last one is a very clear case.

Since we want to measure the magnetic flux emergence and check whether magnetic cancellation takes place, we compute the magnetic flux and magnetic flux density to this small area in HMI LOS-magnetograms. We set a threshold of +(–)90 G, that is, 10σ for 45-s magnetograms (Liu et al. 2012). This threshold is high enough to account for the flux computations in the inter-network and network, avoiding noise effects. Under all considerations, the filling factor is equal to 1.

thumbnail Fig. 4

Left and centre: AIA 94 Å image at different instants with LOS-magnetic field superimposed, with the same contours as Fig. 2. Orange and red contours mark levels of 180, 200, 250, 400, 500 G in transverse field. Right: linear polarization map at 03:48 UT, plot in logarithmic colour scale.

We calculated the signed magnetic (positive and negative) flux separately, as in the left panel of Fig. 3. The magnetic flux has been corrected from heliocentric angle (μ). The vertical bars mark the flare onset. The starting time is Sept. 29 at 00:00 UT. The positive flux reaches a maximum value at Δt = 600 min, but there is a plateau from Δt = 600 min to Δt = 960 min (from 10–16h UT). Both positive and negative fluxes start to decrease around Sept. 29 at 12:00 UT. The maximum unsigned flux is ~3.4 × 1020 Mx, and the unsigned flux rate, 4.29 ± 0.01 × 1017 Mx min-1. Calculating the emergence flux rate (up to the middle of the plateau) yields +2.23 ± 0.07 × 1017 Mx min-1. Considering the first maximum, the rate is +2.60 ± 0.08 × 1017 Mx min-1. We performed the same operation for the negative polarity, obtaining –2.52 ± 0.07 × 1017 Mx min-1. For the negative polarity, we can see a clear decay in both magnetic flux and flux density, as shown in Fig. 3. Calculating the decay rate of this negative area yields 5.2 ± 1.2 × 1016 Mx min-1, which is much slower than the emergence rate.

The signed magnetic flux densities, shown in the right panel of Fig. 3, are both around ~+(–)250 G and follow the trend as the signed magnetic flux. There are two peaks on the positive magnetic flux density some time after the flare, probably caused by recurrent magnetic intensification in the area (see Sect. 3).

To rule out any possible contribution from other areas to the flux emergence in the region, we have also analysed three different areas where magnetic emergence is detected (close to the main sunspot of the active region [700′′, 90′′]. This area was actually flare-productive on Sept. 29 at 05:11 UT, with the onset of a C1.6 flare. In the active region, the positive polarity decreased from 2 × 1021 Mx about 25%, and the negative polarity increased from –1.5 × 1021 Mx by 16%, becoming both fluxes equal 2.5 h before the flare. Furthermore, we studied the areas around coordinates [400′′, 150′′] and [400′′, 30′′], corresponding to supergranular voids, where magnetic flux emergence is also observed. For these regions and the threshold set, the flux emergence rate is negligible. We computed the total magnetic flux of the whole AR as well, from September 24 00:00 UT to Sept. 29 00:00 UT. The positive flux started at 0.9 × 1022 Mx, peaking on Sept. 26 with 1.3 × 1022 Mx. The negative flux started at –1.1 × 1022 Mx, peaking at the same time. At the time of the flare, the fluxes were similar to Sept. 24, keeping the decreasing trend. The rate of the total flux emergence respect to the AR is around 1%.

We also measured the expansion velocity of the emergence, with two main opposite polarities, from Sept. 29 01:00 UT to 13:00 UT with different methods. The first method is the magnetic centroids, as detailed in Balmaceda et al. (2010), but the distance was computed between the positive and negative polarities above the threshold of +(–)90 G. With this method, the expansion velocity is 0.165 ± 0.006 km s-1. The second method we used is visually measuring the distance between the two main polarities (as a diameter and obtaining the radius), getting 0.100 ± 0.003 km s-1. The third procedure is the one described in Palacios et al. (2012), measuring the area and assuming it as circular, subsequently calculating the radius from it. This yields a velocity of 0.157 ± 0.004 km s-1.

thumbnail Fig. 5

Left: height versus time of the filament apex. Two plots are shown: triangles represent the uncorrected height data; and the solid line indicates the height after correcting of heliocentric angle projection. Error bars are computed considering the visual error on height and angle correction. Time starts from 21:00 UT and covers 1 h. In the last 38 min, the height starts increasing. A vertical bar marks the flare onset, at 21:43 UT. Right: height-time diagram, including LASCO C2 (diamonds) and C3 (asterisks) data.

2.1.1. Linear polarization and transverse field of the flux emergence

SDO/HMI retrieves a full-disk imaging of Stokes parameters, in the form of HMI full-disk Stokes spectrograms (HMI.S _720s data series). We used Stokes U and Q images to create linear polarization images. These images are created with six points (line sampling of about 69 mÅ from the centre line; Schou et al. 2012; Centeno et al. 2011) in the Fe i line at 6173 Å. Using the standard SSW routine for HMI/SP reduction, images are accumulated and processed. Profiles are extracted from the images and normalized by the red-most wavelength point in Stokes I (I6). For each frame, we have computed the linear polarization using the following formula (e.g. Solanki et al. 2010; Martínez Pillet et al. 2011): Lpol=16i=16Ui2+Qi2I6·\begin{eqnarray} L_{\rm pol}=\frac{1}{6}\sum^{6}_{i=1}\frac{ \sqrt{U_{i}^{2}+Q_{i}^{2}}}{ I_{6}}\cdot \end{eqnarray}(1)In this stage, we also used also HMI inverted data by the Milne-Eddington inversion code VFISV (Borrero et al. 2011; Centeno et al. 2014). These data consist of the magnetic field, azimuth, and inclination to compute and plot the transverse field through B·sin(γ) with γ the inclination with respect to the LOS, as described in Hoeksema et al. (2014). For a general comparison between the LOS-magnetograms and the LOS-magnetic field obtained by weak field approximation through the Stokes profiles in active regions, see Palacios et al. (2014).

We studied the transverse field magnetic flux density. This includes a threshold of 180 G, which is double the longitudinal flux. The magnetic flux density is included in Fig. 4. In the left and central panels, we show both longitudinal (green and blue) and transverse field (shades of orange), with a 94 Å background, in two different instants: 03:48 and 12:00 UT. The rightmost panel shows the linear polarization at 03:48 UT plotted in logarithmic colour scale. The maximum linear polarization is 1.5%. The structure of the transverse field in the blob is not uniform, but creates bridges between the main positive and negative polarities. In the first four hours, this situation is conspicuous. After that, the transverse field patches remain in the main polarities because some inclination of the mainly vertical magnetic field. Some showers of cosmic rays appear at some moments, altering the most sensitive part of the polarization images, i.e. the linear polarization, translated into transverse magnetic field and making it look noisy.

2.2. Filament rising

We calculated the rising speed on the apex of the filament on de-rotated AIA 304 Å images, also provided in the online material, and whose detailed reference frame is shown in the right panel of Fig. 1. The cadence of this series is 1 min from 21:00 to 22:00 UT. We examined the previous frames from 20:00–21:00 UT, but the rising is negligible at that period. As shown in Fig. 5, the filament kept a constant height during 22 min, and the lift-off took around 38 min. The lift-off of the filament started much before the flare onset at 21:43 UT. The height was estimated on de-rotated images, considering a fiducial point [340′′, 340′′] that lies in the line of the axis defined by the two-ribbon flare middle point and the footpoints. The distance between footpoints is around 0.8 R.The apex height was measured, albeit the southern part of the filament was more active on rising the first 20 min. These measurements are indicated with triangles in Fig. 5. The visual error for the uncorrected height is considered as 3 pixels of a rebinned image of 1024 × 1024 pixels, translated as 5220 km (not shown). We corrected from projection of the heliocentric angle, i.e. dividing by the heliocentric angle sine, considering the eruption follows the local surface normal direction, as shown in the plots of Fig. 5. The mean velocity from the flare to the end at 22:00 UT is 188 ± 5 km s-1, and the subsequent acceleration in this period is 0.065 ± 0.007 km s-2. The maximum velocity is 312 km s-1. (Considering the 38-min period that the rising lasts, the velocity is 115 ± 5 km s-1 and acceleration is 0.049 ± 0.001 km s-2 ).

For this event, we fit the ascent height to different functions, since this can shed light on the physical mechanism that triggered the eruption. We tried exponential, log −log , and parabolic function. Then we evaluated them with the merit function χ2 (unreduced). The exponential fit for the time-height (uncorrected) yields χ2 = 1.7, while the log −log  fitting yields χ2 = 5.9. In the case of projection-corrected height, the fitting is even better: χ2 = 0.9, while the log −log  fitting yields χ2 = 3.4. Parabolic fittings get very large χ2.

We also completed the study with SOHO (Domingo et al. 1995) LASCO C2 and C3 level 1 data (Brueckner et al. 1995). After centring and rotating, we performed image base differences and then, applied a Sobel operator for enhancement. We measured the apex in the filament images. The plate scale for C2 is 11.9′′ pix-1 and for C3 is 56′′ pix-1. We show the whole time-height scale for AIA304, C2 and C3 in Fig 5 (bottom). The χ2 for the whole sequence is 3.80.

2.3. Potential field modelling and critical decay index for the filament region and emergence

We used the package Potential Field Source Surface (PFSS) available in SSW (Schrijver & De Rosa 2003) to assess the critical height of the filament through the decay index of the external field n, since we aim to know whether the filament was already unstable before the emergence or the emergence made the filament unstable. The underlying model for the decay index n is the torus instability. The input of the PFSS model is a synoptic magnetogram. We used the highest resolution input data, the available daily quasi-synoptic HMI magnetograms, obtained at 12:00 UT, on Sept. 25, 26, 28, and 29 (27 not available). These synoptic maps were not derotated, so care had to be taken to identify the areas. We set the PFSS computation grid as 800 points in solar longitude, 400 in latitude, and 66 points in height, covering from 1 R to 2 R. Therefore, the spatial sampling of the grid is 0.45° per grid point in longitude and latitude, and around 10 600 km per grid point on height. We included a context image of the PFSS loop set for this region on Sept. 29, 18:00 UT (top panel of Fig. 6). The height rendered is in agreement with the estimated filament height.

Then, with the decay index of the external field n, used recently in Kliem & Török (2006), Török & Kliem (2007), Zuccarello et al. (2014,and references therein), we estimated this magnitude as n=Rln(Bex)∂R·\begin{eqnarray} n=-R \frac{\partial{\ln(B_{\rm ex})}}{\partial{R}}\cdot \end{eqnarray}(2)We generated datacubes of n for Sept. 25, 26, 28, and 29. We selected a box of 10 × 65 × 66 grid points, equivalent to 4°5 × 29°25 × 1 R (solar longitude, latitude, radial height). The Z axis actually goes from ~20 Mm to 700 Mm, since we removed the points on the very bottom. This final volume of ~55 Mm × 357 Mm × 680 Mm includes the filament and flux emergence. The Y axis covers approximately the filament’s length, but considering its shape, it cannot be considered as wholly lying in the solar longitude selected plane. We clipped the n values from 1.2 to 2.0. We cut this box along and across the exact location where the flux emergence happens, as shown in Fig. 6.

The decay index n is represented in shades of purple-blue on the volume where it is torus-unstable (n becomes a critical factor for instability, ncrit, when values ranges from 1.3 to 1.5, as in Zuccarello et al. 2014). Considering the low impulsivity of the ejection, this threshold is also in agreement with Török & Kliem (2007).

During the previous days, some network points are rooted in the area, but the most important variation are from Sept. 28 12:00 UT to Sept. 29 12:00 UT. The height of the filament may vary and is only marked for Sept. 29, with a hollow circle the projection-uncorrected height, and the corrected height with a cross. On Sept. 28, n increased on the emergence area. Black areas are those where n< 1.2, and these areas decreased from 25 onwards. On Sept. 26, a large network area also made n increase, but it was not located in same coordinates of the emergence.

thumbnail Fig. 6

Top left: PFSS magnetic rendering corresponding to the filament region rendered by PFSS for Sept. 29 data. A low-resolution magnetogram was used, with the same FOV as Fig. 1 (right). Top right: transverse slice of the decay index n across the volume that encloses the filament and emergence. The estimated height of the filament for Sept. 29 is marked with a black cross (corrected), or small circumference (uncorrected). Bottom: longitudinal slice of the decay index n along the volume that encloses the filament and emergence, and same symbols for the filament height. This figure belongs to Sect. 2.3.

3. Discussion

The emergence of this relatively small-scale region follows Hale-Nicholson’s law of polarities (the same leading and trailing polarity in the same hemisphere, Hale et al. 1919); in this case, negative is leading, and positive is trailing, as the leading AR11850. However, in the right panel of Fig. 1, we notice a different magnetic configuration. The filament is located between the negative polarity of the next active region (to the left) and the positive polarity from AR11850 between where we can locate the polarity inversion line (PIL). Therefore, the emergence exhibits a normal configuration, but inverse to the filament configuration, which can help the magnetic instability. The situation is similar to the simulations that originate eruptions in Chen & Shibata (2000), as an opposite polarity configuration emerges beneath or at the side of a filament. Observing in different coronal wavelengths, as 193 Å and 94 Å, there is no evidence of a large magnetic canopy over the filament in the form of bright loops. LASCO images also reveal a leading edge and the filament inside this edge. These images seem to be in agreement with the standard model of a twisted flux rope surrounded by a less twisted longitudinal magnetic field, but without including a large-scale magnetic arcade. The inverse-S shape and barbs would define it as dextral chirality, and the longitudinal field of the filament would then point to the solar north direction. The twisted field lines would exhibit negative helicity.

We have to note that this filament configuration (dextral) follows the hemispheric pattern; however, exceptions to this rule are reported (e.g. see Mackay et al. 2014, and references therein). Furthermore, under a full polarimetry analysis, the chirality definition may change, from a topological orientation of features to a real estimation of the magnetic field direction and sense (Hanaoka & Sakurai 2014).

The barb is very close to the positive polarity from before the supergranular flux emergence, and seems fundamental to the mass dynamics of the filament. Since the barb is very close to the emergence, it may be the point where the instability is transmitted.

The rate of flux emergence, which is about 2 × 1017 Mx min-1, is at least one order of magnitude larger than those reported by Palacios et al. (2012). However, the expansion velocity of the flux emergence is much slower, at 0.16 km s-1, than in Palacios et al., which is approximately 0.6 km s-1. In this case, the higher magnetic fields may slow down the velocity of the upflowing plasma. The expansion velocity calculated with centroids is similar.

Expecting dramatic changes in photospheric magnetic flux as large reconnection features, we only found the bright blob in 94 Å, being conspicuously bright in 211 Å. The magnetic flux emergence lasts around 12 h and smoothly decays, probably due to the reconnection that started a bit after the flux emergence detection. We included spectropolarimetric data to clarify the magnetic structure by adding the transverse component during the whole emergence. We may expect a larger transverse field structure, but that may be smeared out because of the averaging over 12 min, or some circular to linear cross-talk. Combined images between AIA 94 Å images and longitudinal magnetograms suggest the magnetic link between the emergence in the photosphere and the coronal counterpart, in the form of a blob, brighter than the surroundings. This coronal counterpart is also visible in 304 and 193 Å, but somewhat occulted by the barb due to its opacity. The blob is best seen in 211 Å since the emergence patch is better seen through the filament barb, as the latter is less opaque.

Compared to simulations of Kusano et al. (2012) and Toriumi et al. (2013), this emergence region resembles an OP-type configuration with azimuth angle φe ~ 160°, and shear θ0 ~ 60° (φe is the angle between the PIL and the orientation of the bipolar region. With this shear and according to these simulations, the injected kinetic energy might be large enough to provoke an eruption. Bamba et al. (2013) states that OP-type may influence the destabilization of a flux rope via reconnection.

We followed the evolution of the filament for about six days, and the phenomenon known as the “sliding door effect” (Okamoto et al. 2008, 2009; Kuckein et al. 2012) is not evident, since the distance between the positive and negative facular boundaries remains the same, ~220′′. The filament was already there at the beginning of the study, and the part that grew from there does not exhibit any distinctive pattern of the sliding door effect.

Development of pores have been reported by Kuckein et al. (2012), and this is considered tracers of the filament emergence. However, Vargas Domínguez et al. (2012) had studied the event of Okamoto et al. (2008, 2009) and considered these photospheric signatures as necessary but not sufficient pieces of evidence for supporting the hypothesis of emergence. We also find very small appearance of pores. However, it is probably more related to the intensification of the magnetic field (Grossmann-Doerth et al. 1998; Bellot Rubio et al. 2001; Danilovic et al. 2010) rather than characteristics of another flux rope emerging, since the magnetic field density increases from –200 to –1000 G. In the right panel of Fig. 3, other peaks appear in the positive magnetic field density, probably caused by recurrent magnetic field intensification in the region.

Another point is that the CME speed was reported as 600 km s-1. In the case of CME velocities, we always should have in mind the solar escape velocity, that is 617 km s-1. Any difference from this velocity can account for the projection angle or the kinematic energy injection to the CME.

Furthermore, we checked that CH1 becomes darker visually, which might be an indicator of more open field lines there. Some of these features will be treated in another paper.

Regarding the filament eruption, the current instability models might explain or rule out the mechanism that causes the eruption. In the review of, e.g. Aulanier (2014), different models are presented. The breakout model (Antiochos 1998; Antiochos et al. 1999) depends on significant shearing, and involving reconnection above the flux rope (Schmieder et al. 2013). Reconnection is strongly enclosed in some of these models (Zuccarello et al. 2013). Loss of equilibrium models (Forbes & Isenberg 1991; Forbes & Priest 1995) may explain the slow filament rising (Schmieder et al. 2013). These models rely on photospheric motions, polarity approaching, or vortical motion (Amari et al. 2003; Aulanier et al. 2010). Another model to consider is “tether-weakening”, with observational examples presented by Sterling & Moore (2005), Sterling et al. (2007). This weakening is usually due to emergence off the filament channel (Moore & Roumeliotis 1992).

In addition to the observational signatures of the models, the height-time fitting may provide hints about the physical mechanism. In Schrijver et al. (2008), the expansion is fit with different functional forms (polynomial, exponential, power-law...) and evaluated with the merit function χ2. Time-height exponential fitting is also used in very rapid CMEs, like in the case of Gallagher et al. (2003). Gosain et al. (2012) also uses fitting to exponential. The time-height fittings presented here gets better results for exponential functions than any other, and even better to projection-corrected height than uncorrected. Actually, an exponential ascent may be a characteristic of different ideal instabilities, e.g. “torus instability” at the first stages (Schrijver et al. 2008; Gosain et al. 2012). Observing the filament expansion and its shape at several solar radii, it may reveal a kink instability, as in Török & Kliem (2005).

We checked the decay index n of the area from Sept. 25 to 29, similar to the study by Zuccarello et al. (2014). The parameter ncrit is always between 1.3 and 1.5 on the filament and emergence volume, which may indicate a torus instability. However, the lower part of the volume exhibits an important change of n in the area corresponding to the flux emergence, proving useful to identify ongoing instabilities. Importantly, ncrit ranging the mentioned values may also imply an exponential-to-linear expansion (Kliem & Török 2006).

After exploring the main models of eruption in the above paragraphs, ideal instabilities (torus instability and also kink) can be the most plausible mechanism. We may establish a possible scenario as follows: the filament is destabilized, and a large magnetic flux emergence below the filament spine can help the destabilization, elevating against gravity and magnetic tension, increasing the magnetic pressure.

The filament may be located close by a neutral point, which might be metastable (the paler shade of blue in Fig. 6, which is the minimum corresponding n in the Z direction). This point may allow the topology of the magnetic field to be maintained without 3D reconnection with the lower part of the filament, and may help the torus instability to develop due to the flux emergence beneath this point. These kind of eruptions may be slow, leading to CMEs, and with one/several flux emergence patches underneath.

4. Summary and conclusions

This filament, located in the boundaries of the active region AR11850, evolves during almost five days of observations and eventually erupts. On Sept. 29, 00:00 UT, a supergranular flux area started emerging and about 20 h later, the filament started rising. The most conspicuous rising phase lasted about 38 min and the C1.2 flare onset started at 21:43 UT, while the two-ribbon arcade started developing around 15 min later.

The height-time profile corresponds to an exponential profile, and the theoretical model that fits this smooth ascent most is the torus instability, plausibly related to a large magnetic flux emergence, and supported by the height-time profile and the succession of events (first the filament eruption and the flare afterwards).

Online material

Movie of Fig. 1 Access here

Movie of Fig. 4 Access here


5

The temporal evolution shown in Fig. 2 is available in the online edition.

Acknowledgments

The authors want to acknowledge SDO/AIA and SDO/HMI Data Science Centers and Teams, and the GOES and LASCO teams. The LASCO catalog is a CME catalogue that is generated and maintained at the CDAW Data Center by NASA and The Catholic University of America in cooperation with the Naval Research Laboratory. SOHO is a project of international cooperation between ESA and NASA. We would like to thank to the Virtual Solar Observatory (VSO) and Helioviewer for data acquisition and visualisation, respectively. Also, we acknowledge data use from WDC from Geomagnetism, Kyoto, and LMSAL Solarsoft. This research has made use of NASA’s Astrophysics Data System. We would like to thank funding from the Spanish project PPII10-0183-7802 from “Junta de Comunidades de Castilla – La Mancha” and MINECO project code AYA2013-47735-P. We want to thank the anonymous referee for useful comments.

References

  1. Amari, T., Luciani, J. F., Aly, J. J., Mikic, Z., & Linker, J. 2003, ApJ, 595, 1231 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Antiochos, S. K. 1998, ApJ, 502, L181 [NASA ADS] [CrossRef] [Google Scholar]
  3. Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485 [Google Scholar]
  4. Aulanier, G. 2014, in IAU Symp. 300, 184 [Google Scholar]
  5. Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2010, ApJ, 708, 314 [Google Scholar]
  6. Balmaceda, L., Vargas Domínguez, S., Palacios, J., Cabello, I., & Domingo, V. 2010, A&A, 513, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bamba, Y., Kusano, K., Yamamoto, T. T., & Okamoto, T. J. 2013, ApJ, 778, 48 [Google Scholar]
  8. Bellot Rubio, L. R., Rodríguez Hidalgo, I., Collados, M., Khomenko, E., & RuizCobo, B. 2001, ApJ, 560, 1010 [NASA ADS] [CrossRef] [Google Scholar]
  9. Borrero, J. M., Tomczyk, S., Kubo, M., et al. 2011, Sol. Phys., 273, 267 [Google Scholar]
  10. Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, Sol. Phys., 162, 357 [NASA ADS] [CrossRef] [Google Scholar]
  11. Centeno, R., Tomczyk, S., Borrero, J. M., et al. 2011, in Solar Polarization 6, eds. J. R. Kuhn, D. M. Harrington, H. Lin, ASP Conf. Ser., 437, 147 [Google Scholar]
  12. Centeno, R., Schou, J., Hayashi, K., et al. 2014, Sol. Phys., 289, 3531 [Google Scholar]
  13. Chen, P. F. 2011, Living Rev. Sol. Phys., 8, 1 [Google Scholar]
  14. Chen, P. F., & Shibata, K. 2000, ApJ, 545, 524 [Google Scholar]
  15. Danilovic, S., Schüssler, M., & Solanki, S. K. 2010, A&A, 509, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Domingo, V., Fleck, B., & Poland, A. I. 1995, Space Sci. Rev., 72, 81 [Google Scholar]
  17. Feynman, J., & Martin, S. F. 1995, J. Geophys. Res., 100, 3355 [NASA ADS] [CrossRef] [Google Scholar]
  18. Fleck, B., Couvidat, S., & Straus, T. 2011, Sol. Phys., 271, 27 [NASA ADS] [CrossRef] [Google Scholar]
  19. Forbes, T. G., & Isenberg, P. A. 1991, ApJ, 373, 294 [Google Scholar]
  20. Forbes, T. G., & Priest, E. R. 1995, ApJ, 446, 377 [NASA ADS] [CrossRef] [Google Scholar]
  21. Gallagher, P. T., Lawrence, G. R., & Dennis, B. R. 2003, ApJ, 588, L53 [Google Scholar]
  22. Gonzalez, W. D., Joselyn, J. A., Kamide, Y., et al. 1994, J. Geophys. Res., 99, 5771 [Google Scholar]
  23. Gosain, S., Schmieder, B., Artzner, G., Bogachev, S., & Török, T. 2012, ApJ, 761, 25 [NASA ADS] [CrossRef] [Google Scholar]
  24. Grossmann-Doerth, U., Schuessler, M., & Steiner, O. 1998, A&A, 337, 928 [NASA ADS] [Google Scholar]
  25. Hale, G. E., Ellerman, F., Nicholson, S. B., & Joy, A. H. 1919, ApJ, 49, 153 [NASA ADS] [CrossRef] [Google Scholar]
  26. Hanaoka, Y., & Sakurai, T. 2014, in IAU Symp. 300, 515 [Google Scholar]
  27. Hoeksema, J. T., Liu, Y., Hayashi, K., et al. 2014, Sol. Phys., 289, 3483 [Google Scholar]
  28. Jing, J., Yurchyshyn, V. B., Yang, G., Xu, Y., & Wang, H. 2004, ApJ, 614, 1054 [NASA ADS] [CrossRef] [Google Scholar]
  29. Kliem, B., & Török, T. 2006, Phys. Rev. Lett., 96, 255002 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  30. Klimchuk, J. A., & Sturrock, P. A. 1989, ApJ, 345, 1034 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kuckein, C., Martínez Pillet, V., & Centeno, R. 2012, A&A, 539, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Kusano, K., Bamba, Y., Yamamoto, T. T., et al. 2012, ApJ, 760, 31 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
  34. Linker, J. A., Mikić, Z., Riley, P., Lionello, R., & Odstrcil, D. 2003, in Solar Wind Ten, eds. M. Velli, R. Bruno, F. Malara, & B. Bucci, AIP Conf. Ser., 679, 703 [Google Scholar]
  35. Liu, Y., Hoeksema, J. T., Scherrer, P. H., et al. 2012, Sol. Phys., 279, 295 [Google Scholar]
  36. Mackay, D. H., Gaizauskas, V., & Yeates, A. R. 2014, in IAU Symp. 300, 445 [Google Scholar]
  37. Martínez Pillet, V., Del Toro Iniesta, J. C., Álvarez-Herrero, A., et al. 2011, Sol. Phys., 268, 57 [NASA ADS] [CrossRef] [Google Scholar]
  38. Moore, R. L., & Roumeliotis, G. 1992, in Eruptive Solar Flares, eds. Z. Svestka, B. V. Jackson, & M. E. Machado, (Berlin: Springer Verlag), IAU Colloq. 133, Lect. Notes Phys., 399, 69 [Google Scholar]
  39. Moore, R. L., Sterling, A. C., Hudson, H. S., & Lemen, J. R. 2001, ApJ, 552, 833 [NASA ADS] [CrossRef] [Google Scholar]
  40. O’Dwyer, B., Del Zanna, G., Mason, H. E., Weber, M. A., & Tripathi, D. 2010, A&A, 521, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Okamoto, T. J., Tsuneta, S., Lites, B. W., et al. 2008, ApJ, 673, L215 [NASA ADS] [CrossRef] [Google Scholar]
  42. Okamoto, T. J., Tsuneta, S., Lites, B. W., et al. 2009, ApJ, 697, 913 [NASA ADS] [CrossRef] [Google Scholar]
  43. Palacios, J., Blanco Rodríguez, J., Vargas Domínguez, S., et al. 2012, A&A, 537, A21 [Google Scholar]
  44. Palacios, J., Balmaceda, L. A., & Vieira, L. E. 2014, in Polarimetry: From the Sun to Stars and Stellar Environments, Proc. IAU, 10, 97 [Google Scholar]
  45. Parenti, S. 2014, Liv. Rev. Sol. Phys., 11, 1 [Google Scholar]
  46. Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 [NASA ADS] [CrossRef] [Google Scholar]
  47. Robbrecht, E., & Berghmans, D. 2004, A&A, 425, 1097 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Robbrecht, E., Berghmans, D., & Van der Linden, R. A. M. 2009, ApJ, 691, 1222 [NASA ADS] [CrossRef] [Google Scholar]
  49. Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207 [Google Scholar]
  50. Schmieder, B., Démoulin, P., & Aulanier, G. 2013, Adv. Space Res., 51, 1967 [NASA ADS] [CrossRef] [Google Scholar]
  51. Schou, J., Scherrer, P. H., Bush, R. I., et al. 2012, Sol. Phys., 275, 229 [Google Scholar]
  52. Schrijver, C. J., & De Rosa, M. L. 2003, Sol. Phys., 212, 165 [Google Scholar]
  53. Schrijver, C. J., Elmore, C., Kliem, B., Török, T., & Title, A. M. 2008, ApJ, 674, 586 [NASA ADS] [CrossRef] [Google Scholar]
  54. Solanki, S. K., Barthol, P., Danilovic, S., et al. 2010, ApJ, 723, L127 [NASA ADS] [CrossRef] [Google Scholar]
  55. Song, H. Q., Chen, Y., Ye, D. D., et al. 2013, ApJ, 773, 129 [NASA ADS] [CrossRef] [Google Scholar]
  56. Sterling, A. C., & Moore, R. L. 2005, ApJ, 630, 1148 [Google Scholar]
  57. Sterling, A. C., Harra, L. K., & Moore, R. L. 2007, ApJ, 669, 1359 [NASA ADS] [CrossRef] [Google Scholar]
  58. Sturrock, P. A. 1989, Sol. Phys., 121, 387 [NASA ADS] [CrossRef] [Google Scholar]
  59. Toriumi, S., Iida, Y., Bamba, Y., et al. 2013, ApJ, 773, 128 [NASA ADS] [CrossRef] [Google Scholar]
  60. Török, T., & Kliem, B. 2003, A&A, 406, 1043 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Török, T., & Kliem, B. 2005, ApJ, 630, L97 [Google Scholar]
  62. Török, T., & Kliem, B. 2007, Astron. Nachr., 328, 743 [NASA ADS] [CrossRef] [Google Scholar]
  63. van Ballegooijen, A. A., & Martens, P. C. H. 1989, ApJ, 343, 971 [NASA ADS] [CrossRef] [Google Scholar]
  64. Vargas Domínguez, S., MacTaggart, D., Green, L., van Driel-Gesztelyi, L., & Hood, A. W. 2012, Sol. Phys., 278, 33 [NASA ADS] [CrossRef] [Google Scholar]
  65. Wang, Y.-M., & Sheeley, Jr., N. R. 1999, ApJ, 510, L157 [NASA ADS] [CrossRef] [Google Scholar]
  66. Zhang, J., Wang, J., Deng, Y., & Wu, D. 2001, ApJ, 548, L99 [NASA ADS] [CrossRef] [Google Scholar]
  67. Zhang, Y., Zhang, M., & Zhang, H. 2008, Sol. Phys., 250, 75 [NASA ADS] [CrossRef] [Google Scholar]
  68. Zuccarello, F., Balmaceda, L., Cessateur, G., et al. 2013, J. Space Weather and Space Climate, 3, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Zuccarello, F. P., Seaton, D. B., Mierla, M., et al. 2014, ApJ, 785, 88 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Context figures of the filament. The grid spaces mark 15°. Left: image in AIA 193 Å, on 2013 Sept. 29, at 00:00 UT. Image is clipped from 50 to 5000 DN. The filament is located on 20\hbox{$\deg$}N 15\hbox{$\deg$}W, with two coronal holes surrounding it on coordinates 40\hbox{$\deg$}N 20\hbox{$\deg$}W (CH1) and 0\hbox{$\deg$}N 7\hbox{$\deg$}W (CH2). Right: context image showing a derotated AIA 304 Å image displaying the filament, and the positive polarity (in green) and negative (blue). Different colour shades are explained in caption of Fig. 2. The bipolar flux emergence occurrence is located at [250′′, 200′′], enclosed in a white square. A movie available online shows the temporal evolution of the whole visible hemisphere from 21:00 to 22:09 UT in the AIA 304 Å channel. In the movie, the FOV of the right panel of this figure is marked by a black rectangle, and the emergence by a white square. The green and blue contours have the same meaning as in the right panel.

In the text
thumbnail Fig. 2

Top row: frame sequence of AIA 304 Å images of a filament’s barb (darker areas) with LOS-magnetic field superimposed. Different shades of green (from lighter to darker) contours mark different levels, as 90, 200, 500, 600, and 700 G, while a gamut of blue contours mark –90, –200, –500, –700, –900 G for the negative polarities. The frame size is 100′′× 100′′. Sequence timing aims to show the magnetic flux emergence (panel 1), pre-flare (panel 2), flare (panel 3), and decay (panel 4), from left to right. Bottom row: same temporal sequence in AIA 211 and 94 Å. An associated movie is available online.

In the text
thumbnail Fig. 3

Left: computed magnetic flux on positive patches (solid line), on negative patches (dashed line), signed magnetic flux (dotted-dashed line), unsigned magnetic flux (dotted line). A vertical line marks the time of the flare at 21:43 UT. Right: mean values of the magnetic flux density of the positive (solid) and negative (dashed) regions.

In the text
thumbnail Fig. 4

Left and centre: AIA 94 Å image at different instants with LOS-magnetic field superimposed, with the same contours as Fig. 2. Orange and red contours mark levels of 180, 200, 250, 400, 500 G in transverse field. Right: linear polarization map at 03:48 UT, plot in logarithmic colour scale.

In the text
thumbnail Fig. 5

Left: height versus time of the filament apex. Two plots are shown: triangles represent the uncorrected height data; and the solid line indicates the height after correcting of heliocentric angle projection. Error bars are computed considering the visual error on height and angle correction. Time starts from 21:00 UT and covers 1 h. In the last 38 min, the height starts increasing. A vertical bar marks the flare onset, at 21:43 UT. Right: height-time diagram, including LASCO C2 (diamonds) and C3 (asterisks) data.

In the text
thumbnail Fig. 6

Top left: PFSS magnetic rendering corresponding to the filament region rendered by PFSS for Sept. 29 data. A low-resolution magnetogram was used, with the same FOV as Fig. 1 (right). Top right: transverse slice of the decay index n across the volume that encloses the filament and emergence. The estimated height of the filament for Sept. 29 is marked with a black cross (corrected), or small circumference (uncorrected). Bottom: longitudinal slice of the decay index n along the volume that encloses the filament and emergence, and same symbols for the filament height. This figure belongs to Sect. 2.3.

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.