Issue |
A&A
Volume 521, October 2010
|
|
---|---|---|
Article Number | A56 | |
Number of page(s) | 12 | |
Section | The Sun | |
DOI | https://doi.org/10.1051/0004-6361/200913275 | |
Published online | 20 October 2010 |
Trend of photospheric magnetic helicity flux in active regions generating halo coronal mass ejections![[*]](/icons/foot_motif.png)
A. Smyrli1,2 - F. Zuccarello1 - P. Romano3 - F. P. Zuccarello4 - S. L. Guglielmino1 - D. Spadaro3 - A. W. Hood2 - D. Mackay2
1 - Dipartimento di Fisica e Astronomia - Sezione Astrofisica, Università di Catania, via S. Sofia 78, 95123 Catania, Italy
2 -
Institute of Mathematics, University of St. Andrews, The North Haugh, St. Andrews, Fife KY 169SS,
Scotland, UK
3 -
INAF - Osservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy
4 -
Centre for Plasma Astrophysics, K. U. Leuven, Celestijnenlaan 200B, 3001 Leuven, Belgium
Received 10 September 2009 / Accepted 22 April 2010
Abstract
Context. Coronal mass ejections (CMEs) are very energetic events (1032 erg)
initiated in the solar atmosphere, resulting in the expulsion of
magnetized plasma clouds that propagate into interplanetary space. It
has been proposed that CMEs can play an important role in shedding
magnetic helicity, avoiding its endless accumulation in the corona.
Aims. The aim of this work is to investigate the behavior of
magnetic helicity accumulation in sites where the initiation of CMEs
occurred to determine whether and how changes in magnetic helicity
accumulation are temporally correlated with CME occurrence.
Methods. We used MDI/SOHO line-of-sight magnetograms to
calculate magnetic flux evolution and magnetic helicity injection in 10
active regions that gave rise to halo CMEs observed during the period
2000 February to 2003 June.
Results. The magnetic helicity injection does not have a unique
trend in the events analyzed: in 40% of the cases it shows a large
sudden and abrupt change that is temporally correlated with a CME
occurrence, while in the other cases it shows a steady monotonic trend,
with a slight change in magnetic helicity at CME occurrence.
Conclusions. The results obtained from the sample of events that
we have analyzed indicate that major changes in magnetic helicity flux
are observed in active regions characterized by emergence of new
magnetic flux and/or generating halo CMEs associated with X-class
flares or filament eruptions. In some of the analyzed cases the changes
in magnetic helicity flux follow the CME events and can be attributed
to a process of restoring a torque balance between the subphotospheric
and the coronal domain of the flux tubes.
Key words: Sun: activity - Sun: coronal mass ejections (CMEs) - magnetic fields
1 Introduction
Coronal mass ejections (CMEs) are sudden expulsions of magnetized plasma clouds from the solar atmosphere into the interplanetary medium. They are usually observed by coronagraphs above the occulting disk in projection on the plane of the sky, and therefore their connection to the magnetic region where they were initiated is not straightforward. In many cases, a twisted flux tube ejected in a CME appears in interplanetary space as a magnetic cloud in which the twisted nature remains well observable (Qiu et al. 2007; Longcope et al. 2007).
Depending on the velocity and acceleration values, it is possible to distinguish between two different classes of CMEs: a) gradual CMEs, characterized by velocities
-600 km s-1 and gradual acceleration (
-40 m s-2 within a distance from the Sun less than
30
), generally associated with eruptive quiescent filaments; b) impulsive CMEs with higher initial velocities, in the range
750-1000 km s-1, and decelerating at distances
2
(Sheeley et al., 1999), generally associated with flares and eruptive active region filaments.
MacQueen & Fisher (1983) suggested that the reason for this different dynamical behavior might be due to a different driver mechanism (or at least one with a different strength): gradual CMEs, which are events associated to eruptive prominence, are supposed to undergo a significant net propelling force over extended periods, while impulsive CMEs, which are flare-associated events, are supposed to arise from an impulsive input into the low corona (see also Low & Zhang 2002; Chen & Krall 2003).
Although several models for the initiation of CMEs have been developed (see e.g. Forbes 2000), there is a lack of observational evidence to provide the necessary constraints. In the breakout model (Antiochos et al. 1999) a quadrupolar flux distribution is considered consisting of an overlying arcade whose flux is oppositely directed with respect to the underlying central arcade and two lateral flux systems. In this scenario the emergence of new magnetic flux (Zuccarello et al. 2008) and/or the shearing motions along the magnetic inversion line eventually result in a CME. Priest & Forbes (1990) and Forbes & Priest (1995) showed that within a bipolar topology, converging motions can result in the loss of equilibrium of a magnetic flux-rope leading to a catastrophic eruption. Chen & Shibata (2000) showed that the injection of new magnetic flux from below the photosphere is also responsible for the eruption of a pre-existing flux-rope. More global MHD simulation of the emergence of a twisted flux tube have been carried out by Archontis & Török (2008). The authors found that a flux rope is formed within the expanding field due to shearing and reconnection of field lines at low atmospheric heights. Depending on if the corona is magnetized or not, the expanding flux rope experiences a full eruption or remains confined.
A useful quantity to distinguish between these cases is the magnetic helicity. It provides a measure of the connection of the magnetic field lines within a volume and also quantifies how much a flux tube is twisted or a magnetic arcade is sheared (Berger 1984; Heyvaerts & Priest 1984; Pariat et al. 2006; Démoulin & Pariat 2009). It is one of the few global quantities which is conserved also in resistive magnetohydrodynamics (MHD), on a timescale much smaller than the global diffusion timescale. However, when CMEs release magnetic energy stored over a large volume, they lead to the partial opening of the magnetic field lines and, as a result, the helicity does not remain conserved in a volume on the active region scale. In particular, the amount of magnetic helicity that can be expelled during a CME can be estimated by the variation of coronal helicity in the active region that generates a CME.
In the recent simulations of Yeates & Mackay (2009), the authors carried out global simulations of the solar corona based on photospheric observations for a period of six months in 1999. The model included the effects of flux emergence, shearing by surface motions and cancellation along with the formation and ejection of flux ropes due to loss of equilibrium. By varying the rate at which magnetic helicity was injected, the model produced 50% of the observed CME rate over the simulated period (see also Mackay & van Ballegooijen 2006). In contrast to this model, Cook et al. (2009) studied the number and locations of coronal null points, a key element in the magnetic breakout model, over two solar cycles. They found that the number of nulls followed a cycle variation, peaking at maximum of 15-17 per rotation. The vast majority of null points were located above active latitudes and lay at low heights in the corona.
Coronal mass ejections are prime candidates for carrying away the accumulated helicity (Heyvaerts & Priest 1984; Low 1996). They also provide a mechanism for shedding the required amount of helicity to avoid the saturation of the mean field
to a value much lower than the equipartion field strength, estimated as 1046 Mx2 per cycle from simulations (Brandenburg & Sandin, 2004).
Thus, it is crucial to estimate the amount of magnetic helicity that
can be expelled during a CME, for instance by measuring the variation
of coronal helicity in the active region that generates the CME.
Nindos & Andrews (2004) performed a statistical study of 133 big flares of M and X class and found that in the pre-flare phase the coronal magnetic helicity of active regions producing big flares not correlated with CMEs is smaller than the coronal magnetic helicity of active regions producing CME-associated big flares. Furthermore, Moon et al. (2002) found that an amount of magnetic helicity on the order of 1041 Mx2 was impulsively injected into the corona around the peak time of four eruptive flares.
On the basis of these results it is clear that the study of the magnetic helicity evolution in active regions can provide a useful tool to obtain important information about the mechanisms able to produce instabilities in the magnetic configuration. Moreover, the investigation of the magnetic helicity trend, before and after a CME occurrence, might contribute to our understanding of the role effectively played by these events in any eventual change in the magnetic helicity accumulation process (see, e.g., Nindos et al. 2003).
In this regard it is important to make a distinction between two cases. The former is related with the change in magnetic helicity prior to the CME occurrence, which should be related with phenomena able to destabilize the coronal configuration, like new flux emergence or shearing/twisting of the magnetic field lines. The latter is related with changes in magnetic helicity occurring after the CME that, as suggested by Longcope & Welsh (2000) and Chae et al. (2003), might be due to the unbalance of torque between the sub-photospheric part of the flux rope and the coronal field which just lost some stress via the CME launch.
Therefore we determined here the evolution of both magnetic flux and magnetic helicity flux in active regions where a flare or a filament eruption occurred close in time with a halo CME to investigate the above mentioned scenarios.
The layout of this paper is as follows: in the next section we recall the definition of magnetic helicity and related parameters, in Sect. 3 we describe the data, the methods of analysis and the results, in Sect. 4 we present the discussion and our conclusions.
2 Magnetic helicity
Magnetic helicity is given by a volume integral
![]() |
(1) |
where



![]() |
(2) |
where





![]() |
(3) |
This relative magnetic helicity may vary due to the passage of magnetic field lines through the photospheric boundary (when the flux emerges), or through the sides (when the magnetic field expands to be in contact with neighboring active regions), or through the outer boundary (when a coronal mass ejection occurs). Unlike the magnetic helicity of fully enclosed magnetic fields, the helicity of open magnetic fields may change with time as it is transported through the boundaries. In this regard, it is worthwhile mentioning that most of the helicity injection in an active region occurs during magnetic flux emergence (Tian & Alexander 2008).
Therefore, the magnetic helicity in the coronal volume of an active
region is a function of time. In order to determine this time
variation, it is possible to use the following equation for the
helicity flux (or injection rate) dH/dt through a planar surface
(indicating the photospheric portion of the surface S bounding the volume V) (Berger and Field, 1984)
![]() |
(4) |
where the magnetic field





The magnetic helicity content of an active region at a specific time can be determined by integrating the helicity injection rate over time since its birth to the specified time, provided that no CME occurred which would have carried away helicity. According to Eq. (4), the magnetic helicity in an open volume can change either by the shuffling horizontal motions of field lines on the surface (first term, called shear term) or by the passage of twisted or sheared fields through the solar surface (second term, called emergence term).
![]() |
Figure 1: a) Difference image showing the halo CME observed by LASCO-C2 on 2002 july 15 at 20:30:05 UT, superposed on the EIT difference image relevant to 20:24 UT; b) SOHO/EIT full disk image acquired on 2002 july 15 at 20:24:10 UT, showing the location on the disk of the CME source region. |
Open with DEXTER |
![]() |
Figure 2: a) Zoom of the active region NOAA 10030 where the CME was initiated; the image was acquired by SOHO/EIT during the main phase of the flare; b) MDI/SOHO magnetogram showing the magnetic configuration of NOAA 10030 on 15 July 2002. |
Open with DEXTER |
Table 1: Data on the halo CMEs and analyzed relevant eruptive events.
Chae (2001) showed that when magnetic fields are normal to the surface, the flux distribution on the surface
can
be exploited to determine the non-steady transverse motion of the
field-line footpoints, and proposed to use the displacement velocity
determined by the technique of the local cross-correlation tracking (LCT) as a proxy for the transverse motion
.
He also introduced a Fourier transform method for the fast computation of
from the
.
Moreover, following Démoulin & Berger (2003), it is possible to write
![]() |
(5) |
and the photospheric helicity flux (or injection rate) simplifies to
![]() |
(6) |
These equations put together the shear and emergence terms, as Démoulin & Berger (2003) realized that the emergence of a twisted flux system would cause what may appear as shear flows.
However, Pariat et al. (2005) showed
that Eq. (6) introduces artificial signals because helicity flux
densities per unit surface are not physical quantities. Therefore they
presented a new estimation of helicity flux per unit of elementary
magnetic flux. In this case the helicity flux is given by
![]() |
(7) |
where S = S' is the integration surface,

The LCT method has some important practical advantages. First of
all, it makes use of the line-of-sight magnetograms that are commonly
available. Secondly, it is simple to implement, robust and flexible.
However, it is possible to introduce large systematic errors into the
results by not choosing interpolation algorithms with sufficient care (Potts et al., 2003)
or by not choosing the correct size for the two input parameters of
this method: FWHM (full width at half maximum) that defines the width
of the apodizing window, and the time interval
between two images of comparison.
The accumulated magnetic helicity
(or total magnetic helicity injection) due to changes across the
photospheric boundary can be determined by integrating the values of
the helicity injection rate over the entire time interval
analyzed:
![]() |
(8) |
Finally, indicating with



3 Data analysis
The method used to carry out our analysis can be described by the
following steps: 1) selection (using SOHO/LASCO data) of halo CMEs; 2)
search for the flare or eruptive filament closest in time for each
selected CME; 3) attempt to identify the active region (AR) where the
CME was initiated by means of a double cross-check based on the
flaring-eruptive activity and the use of SOHO/EIT difference images; 4)
determination of the magnetic flux trend, magnetic helicity injection
rate dH/dt, accumulated magnetic helicity
and ratio
for the selected active regions.
More precisely, we selected a sample of halo CMEs from the SOHO/LASCO on-line catalog (http://cdaw.gsfc.nasa.gov/CME_list/). The choice to select halo CME is related to the fact that these events are generally initiated in areas located in the central part of the solar disk: this aspect is quite important when an analysis of line-of-sight magnetograms is carried out, because in this situation possible problems related to projection effects can be avoided. There were a total of 38 selected events in the time interval 2000 February 1st-2003 June 30; we show in Fig. 1a an example of a halo CME and the associated flare observed by SOHO/EIT in Fig. 1b.
We successively compared the time inferred by a back-extrapolation in time of the CME height-time trajectory with the time of the flare or filament eruption recorded in a time interval between one hour before and one hour after the CME, using the NOAA flare catalog (http://www.solarmonitor.org/index.php), in an attempt to obtain a first association between each CME and the corresponding flare or filament eruption.
In a next step, we used SOHO/EIT 171 Å images acquired during a time interval of two hours before and two hours after the CME occurrence and performed difference images between successive images to determine a possible increase or decrease of brightness. The double-cross check between the eruptive event registered in the NOAA report and the brightness variation inferred by the difference maps provided us a further indication about the active region that was the site of the initiation of the halo CME.
As an example, we show in Fig. 2a the EIT/SOHO image of the active region NOAA 10030 which was identified by the previously mentioned double check as the source region of the halo CME shown in Fig. 1a. Figure 2b shows the MDI/SOHO magnetogram relevant to NOAA 10030, acquired at 20:30 UT.
For some active regions that were very close to other active regions, it was not possible to draw the respective boundaries (e.g., NOAA 9114, NOAA 9393, NOAA 10069), therefore to avoid to neglect potentially significant flux that could influence the helicity flux measurements, we merged the above mentioned active regions with their neighboring ones (e.g., NOAA 9114+9115, NOAA 9393+9394, NOAA 10069 + 10077) and calculated the total helicity flux .
In order to minimize problems related to projection effects, we focused our attention on active regions located between
of longitude from the central meridian. With this choice, we limited
our analysis to 10 active regions and, correspondingly, to 12 CMEs,
because both NOAA 10030 and 10365 were the source regions of two halo
CMEs.
The halo CME events selected and the relevant active regions found by this method are reported in Table 1, where we can distinguish two classes of active regions: the first containing ARs giving rise to gradual halo CMEs, i.e. with positive acceleration (hereinafter called Class I) and the second relevant to ARs that were sources of impulsive halo CMEs, i.e. characterized by negative acceleration (Class II).
We determined the magnetic flux evolution and the above mentioned
parameters concerning the magnetic helicity injection for these two
sets of active regions by means of MDI/SOHO low-resolution magnetograms
taken with a 96 min cadence. We estimated the vertical field
strength
by multiplying the MDI/SOHO line-of-sight component by 1/cos
,
where
is the heliocentric angle.
We studied the magnetic evolution of each active region with respect to
the time of the first observation of the CME in the LASCO/C2
field-of-view, which is considered below as the reference-time point.
The collection of the line-of-sight magnetograms for each active
region was carried out over a period of time that was chosen on the
basis of the following two characteristics: 1) the time period expands
over at least 24 h before and after the reference time; 2) during
this time interval the longitude of the active region must be within 35
.
We spatially aligned the sequence of magnetograms relevant to each active region by successive non-linear mapping, which takes into account the solar differential rotation effect.
Then we calculated the horizontal velocities of magnetic field
lines at photospheric footpoints by measuring the local displacements
of magnetic flux concentrations between two successive magnetograms
with the LCT technique. We selected
and
min.
![]() |
Figure 3:
From top to bottom: evolution of the magnetic flux (continuous
line indicates the average total flux, marked line the positive flux
and dotted line the negative flux), magnetic helicity injection rate,
total helicity injection, ratio
|
Open with DEXTER |
![]() |
Figure 4: Same as in Fig. 3 for ARs 10030 and 10069+10077, belonging to Class II (impulsive CMEs). |
Open with DEXTER |
We used the method developed by Pariat et al. (2005) to estimate the magnetic helicity flux dH/dt and to reduce spurious signals (see also Pariat et al. 2007, and references therein; Démoulin & Pariat 2009).
Then the accumulated amount of helicity
could be estimated by integrating the measured dH/dt
from the start of the observing run to a specific time. Finally, the
total, positive and negative magnetic flux crossing the photosphere and
the normalized helicity
were computed.
3.1 Class I: gradual CMEs
In Fig. 3 we report the results obtained for two active regions related to gradual CMEs. Similar plots, concerning the other ARs belonging to the same Class, are reported in Fig. A.1 of Appendix A.
More specifically, the plots in Fig. 3 report the evolution of total, positive, and negative magnetic flux expressed in 1022 Mx (first row), helicity injection rate dH/dt in 1042 Mx2 h-1 (second row), accumulated helicity injection
in 1042 Mx2 (third row), ratio
(fourth row), positive (H+) and negative (H-) accumulated helicity injection (fifth row) for NOAA 8858 and NOAA 9415.
NOAA 8858 (Fig. 3,
left panel): the total magnetic flux shows a very small decrease during
the analyzed period (more precisely, the positive flux remains almost
constant, while the negative flux slightly decreases). The helicity
injection rate is generally negative, with the exception of few points
with positive values. The accumulated helicity injection
shows a continuous increase (in absolute value), with a significant change in the curve slope after the CME event (third row). This change is also evident in the
curve (fourth row), where it seems possible to distinguish two
different regimes. The inspection of the graph shown in the fifth row,
which gives the trend of the accumulated negative (asterisks) and
positive (crosses) helicity injection, reveals that H- increases and changes its trend after the CME occurrence, while H+
is very small and almost constant during the analyzed period. Note that
the gradual CME was characterized by a velocity of 944 km s-1 and an acceleration of 11.4 m s-2, and was associated with a C7.3 class flare.
NOAA 9415 (Fig. 3,
right panel): the total magnetic flux is constant during the analyzed
period, whereas the helicity injection rate shows alternating negative
and positive values, with some more abrupt changes just before the CME.
The accumulated helicity injection shows predominantly negative values
prior to the CME occurrence and a successive rise to positive values
after it. The ratio
has a very irregular behavior, with a sudden change immediately before the CME and a clear increasing trend after it.
The values of H- and H+
show a very peculiar behavior, characterized by alternate phases where
one component remains constant, while the other increases.
Interestingly, one of the changes is temporarily correlated with the
CME occurrence. The velocity and acceleration of the CME were
2411 km s-1 and 211.6 m s-2, respectively, and the associated flare was of X2.3 class.
Table 2: Characteristics of the active regions, halo CMEs and analyzed relevant eruptive events.
3.2 Class II: impulsive CMEs
Figure 4 shows the same parameters as in Fig. 3 for NOAA 10030 and 10069+10077, where impulsive CMEs occurred, while results concerning NOAA 10162, 10229 and 10365, belonging to the same Class, are reported in the Appendix (Fig. A.2).
NOAA 10030 (Fig. 4,
left panel): during the analyzed period (five days), the total magnetic
flux (as well as the positive and negative fluxes) shows a continuous
increase, the helicity injection rate varies initially from positive to
negative values and shows sudden changes in coincidence with the
occurrence of two halo CMEs; in particular, after the second event, it
again assumes positive values.
shows a continuous decrease after an initial increase to a maximum value of
Mx2 h-1, and an abrupt variation immediately after the second CME. However, when we look at the curve showing the ratio
(fourth row), we notice that there is a very abrupt change immediately before the first CME and immediately after
the second CME. The first change is due to the change in the sign of
the accumulated helicity injection, while the latter might be explained
by the changes in the trends of H- and H+. In particular, concerning the trend of H- and H+, reported in the fifth row, we can see that H-
is initially constant, then (some hours before the first CME) it starts
to increase and this increase abruptly stops just before the second
CME. The behavior of H+ is almost specular: it is increasing at the beginning of the period examined, then it stops (approximately when H- starts to rise) and again suddenly starts to increase in coincidence with the CME (when H-
becomes constant). Note that the first CME, associated with an X3.0
class flare, was characterized by a velocity, as measured by LASCO C2,
of 1151 km s-1 and a deceleration of 25.6 m s-2, whereas the second CME had a velocity of 1636 km s-1, a deceleration of 41.0 m s-2, and was associated to an erupting filament.
NOAA 10069+10077 (Fig. 4,
right panel): the magnetic flux shows an increase during the analyzed
period, while the helicity injection rate (second row) is initially
characterized by negative values and, after the CME occurrence, by
positive values. This trend is reflected in the curve representing the
accumulated helicity injection ,
where the abrupt variation immediately after the CME is evident. The same behavior is shown by the
values (fourth row). The graph shown in the fifth row shows that H- and H+ have a similar trend to that observed for NOAA 10030: H- initially increases while H+ remains constant and after the CME they invert their trends: H- remains constant while H+ starts to increase. In this case, the CME velocity was 1585 km s-1, the deceleration was 67.1 m s-2, and the associated M5.2 flare was related to an erupting filament.
We note that the
values obtained in our analysis for both Classes of CMEs are lower than those obtained by Démoulin & Pariat (2009):
this is probably due to the limited time interval we used for the
measurements of helicity accumulation, which in most of the cases
excluded the emergence phase.
4 Discussion and conclusions
The role of magnetic helicity evolution in active regions producing eruptive events is considered to be very important and it is believed that the accurate measurement and modeling of the magnetic helicity is essential to understand the coupling between the photosphere, the corona and the heliosphere.
The study presented here was aimed at determining the temporal variations of magnetic flux and magnetic helicity injection in a sample of active regions generating halo CMEs, and how these variations could possibly be correlated with the CME occurrence.
We investigated the evolution of magnetic flux, magnetic helicity injection rate, and total magnetic helicity injection for ten active regions leading to twelve halo CMEs.
As a general result, we can deduce from the graphs shown in Figs. 3, 4, Figs. A.1, A.2
and from Table 2 that there is no unique behavior in magnetic
helicity injection accompanying halo CME occurrence. In fact, in some
cases there is an abrupt change in
correlated in time with the CME event (see, e.g., NOAA 8858, 9415,
10030 [second CME], 10069+10077, 10365 (second CME)), in others there
is a slight change in the curve slope (NOAA 9114+9115, 9393+9394,
10162), and in some no significant variation is recorded (NOAA 8948,
10229, 10365 (first CME)). Therefore we cannot infer a general rule
from this analysis.
We can, however, try to determine what conditions could give rise to these different behaviors and we discuss some possible relationships between the parameters characterizing the sample of ARs we analyzed.
A first conclusion that is possible to draw from the graphs shown in Figs. 3, 4 and Figs. A.1, A.2 is that the presence (or absence) of an abrupt change in cannot
be ascribed to the fact that the CME is gradual or impulsive: we
actually found these changes in both classes (even if the most
conspicuous
variations are associated with impulsive CMEs).
Nevertheless, an interesting result that can be deduced from the analysis of the total magnetic flux concerning the two classes is that the latter remains almost constant or shows a very slight increase for all the ARs that give rise to gradual CMEs, while it is significantly increasing or decreasing (with the exception of NOAA 10162) for impulsive CMEs. Moreover, on average the variation in helicity for the impulsive CMEs is more significant than for gradual CMEs. So the contribution in magnetic helicity injection due to flux emergence is quite significant also for more mature active regions, analogous to the results obtained by Tian & Alexander (2008) for new-born active regions.
In this regard, we recall that as already stated by other authors (see, e.g. Démoulin & Pariat 2009, and references therein), an increase of the magnetic flux, indicating the emergence of new flux tubes from the sub-photospheric layers, is the most important effect in causing the increase in the magnetic helicity flux. It is also possible that this emergence destabilizes the coronal magnetic configuration and leads to an eruptive event. On the basis of our results, we can conclude that this process was probably at work in some of the analyzed active regions, that is NOAA 9114+9115, NOAA 9393+9394 and NOAA 10030 (first CME), taking into account that in these active regions the change in magnetic helicity flux was observed before the CMEs. Note that these CME events are associated with X-class flares or erupting filaments.
Another effect, which causes a change in magnetic helicity flux after the CME, is related to a possible lack of torque balance between the subphotospheric and the coronal domain of a flux rope, which can occur once the CME has removed the accumulated stress in the magnetic field (Longcope & Welsh 2000; Chae et al. 2003; Petsov et al. 2003; Petsov et al. 2008).
The imbalance of magnetic torque can induce a transfer of helicity from the convective zone below the active region by torsional Alfvén waves, which is present until the coronal field has a torque large enough to balance the sub-photospheric torque.
Our analysis indicates that the observed change in the magnetic helicity flux occurred after the CMEs could be ascribed to this process for NOAA 8858, NOAA 9415, NOAA 10030 (second CME), NOAA 10069+10070 and probably also for NOAA 10365 (second CME).
However, in some of these cases, a flux emergence process also occurs, so it is not possible to state definitely what the main mechanism at work is. What we can notice from the trend of H+ and H- is though that for some of the active regions just mentioned these parameters show an alternate trend. In these active regions the cancellation of magnetic helicity should be taking place. Therefore it would be worthwhile to extend this study to other active regions that show a change in the magnetic helicity flux after a CME to obtain further information on this process.
Concerning NOAA 8948 and NOAA 10229, both characterized by decreasing magnetic flux and no change in the magnetic helicity trend both before and after the CME, we recall that also the photospheric shearing motions can be responsible for increasing the stress in the magnetic field. However, the new, less energetic configuration reached by these active regions after the eruptive events was probably not characterized by a torque imbalance as in the other cases. Interestingly, both events are associated with M-class flares, while the majority of the events quoted in the second category (i.e., those associated with the torque unbalance) are associated with X-class flares (characterized by a higher energy output than M-class flares) or with erupting filaments, whose expulsion could greatly change the coronal magnetic configuration of the active region.
Another interesting result concerns the hemispheric helicity
rule, which implies that H should be negative for ARs in the Northern
hemisphere and positive for those in the Southern hemisphere (Seehafer, 1990). Figures 3, 4, Figs. A.1, A.2 and Table 2
show that the analyzed ARs have a helicity accumulation in agreement
with this rule, with the exception of NOAA 9415 and NOAA 10030 (showing
both positive and negative H
values during the analyzed time interval), 10069+10077 (showing
negative helicity) and 10162 (characterized by positive helicity).
Therefore 4/5 of the active regions of Class I have a helicity
accumulation in agreement with the hemispheric helicity rule, while 3/5
of the active regions associated with impulsive CMEs seem to be related
with an exception to this rule.
In conclusion, this study showed that magnetic helicity accumulation in our sample of active regions generating halo CMEs can exhibit significant changes when there are indications of new magnetic flux emergence. This emergence can change the helicity injection rate and also destabilize the coronal magnetic configuration, leading to an eruptive event.
Significant changes of the helicity injection rate can apparently also occur without the involvement of magnetic flux emergence. In most of these cases the observed changes follow the CME events and can be attributed to a process of restoring a torque balance between the subphotospheric and the coronal portion of the flux tubes.
Finally, we emphasize the result concerning the magnetic helicity injection with alternating sign found in those active regions where a considerable change of magnetic helicity was found, that was correlated in time with the CME occurrence. As in these active regions a cancellation of magnetic helicity should take place, we plan to further investigate this process both in these ARs as well in others, where a similar trend could be singled out.
AcknowledgementsThe authors thank the anonymous referee for helpful comments and suggestions. This work was supported by the European Commission through the SOLAIRE Network (MTRN-CT-2006-035484), by the Istituto Nazionale di Astrofisica (INAF), by the Agenzia Spaziale Italiana (contract I/015/07/0) and by the Università degli Studi di Catania. The CME catalog is generated and maintained at the CDAW Data Center by NASA and The CatholicUniversity of America in cooperation with the Naval Research Laboratory. SOHO is a project of international cooperation between ESA and NASA.
References
- Alissandrakis, C. E. 1981, A & A, 100, 197 [Google Scholar]
- Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485 [Google Scholar]
- Archontis, V., & Török, T. 2008, A&A, 492, L38 [Google Scholar]
- Berger, T. E., & Lites, B. W. 2003, Sol. Phys., 213, 213 [Google Scholar]
- Berger, M. A. 1984, Geophys. Astrophys. Fluid. Dyn., 30, 79 [Google Scholar]
- Berger, M. A., & Field, G. B. 1984, J. Fluid Mech., 147, 133 [NASA ADS] [CrossRef] [Google Scholar]
- Brandenburg, A. 2007, Highlights Astron., 14, 291 [NASA ADS] [Google Scholar]
- Brandenburg, A., & Sandin, C. 2004, A&A, 427, 13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chae, J. 2001, ApJ, 560, L95 [NASA ADS] [CrossRef] [Google Scholar]
- Chae, J., Moon, Y.-J., Rust, D. M., Wang, H., & Goode, P. R. 2003, JKoAS, 36, 33 [Google Scholar]
- Chen, J. 1996, J. Geophys. Res., 101, 27499 [NASA ADS] [CrossRef] [Google Scholar]
- Chen, J., & Krall, J. 2003, J. Geophys. Res., 108, SSH 2-1 [Google Scholar]
- Chen, P. F., & Shibata, K. 2000, ApJ, 545, 524 [Google Scholar]
- Cook, G. R., Mackay, D. H., & Nandy, D., 2009, ApJ, 704, 1021 [NASA ADS] [CrossRef] [Google Scholar]
- Démoulin, P., & Berger, M. A. 2003, Sol. Phys., 215, 203 [NASA ADS] [CrossRef] [Google Scholar]
- Démoulin, P., & Pariat, E. 2009, Adv. Space Res., 43, 1013 [NASA ADS] [CrossRef] [Google Scholar]
- Forbes, T. G. 2000, J. Geophys. Res., 105, 23153 [NASA ADS] [CrossRef] [Google Scholar]
- Forbes, T. G., & Priest, E. R. 1995, ApJ, 446, 377 [NASA ADS] [CrossRef] [Google Scholar]
- Heyvaerts, J., & Priest, E. R. 1984, A&A, 137, 63 [NASA ADS] [Google Scholar]
- Longcope, D. W., & Welsch, B. T. 2000, ApJ, 545, 1089 [NASA ADS] [CrossRef] [Google Scholar]
- Longcope, D., Beveridge, C., Qiu, J., et al. 2007, Sol. Phys., 244, 45 [NASA ADS] [CrossRef] [Google Scholar]
- Low, B. C. 1996, Solar Phys., 167, 217 [Google Scholar]
- Low, B. C., & Zhang, M. 2002, ApJ, 564, L53 [NASA ADS] [CrossRef] [Google Scholar]
- Mackay, D. H., & van Ballegooijen, A. A. 2006, ApJ, 641, 577 [NASA ADS] [CrossRef] [Google Scholar]
- MacQueen, R. M., & Fisher, R. R. 1983, Sol. Phys., 89, 89 [NASA ADS] [CrossRef] [Google Scholar]
- Moon, Y. J., Chae, J., Wang, H., Choe, G. S., & Park, Y. D. 2002, ApJ, 580, 528 [NASA ADS] [CrossRef] [Google Scholar]
- Nindos, A., & Andrews, M. D. 2004, ApJ, 616, L175 [NASA ADS] [CrossRef] [Google Scholar]
- Nindos, A., Zhang, J., & Zhang, H. 2003, ApJ, 594, 1033 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pariat, E., Démoulin, P., & Berger, M. A. 2005, A&A, 439, 1191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pariat, E., Nindos, A., Demoulin, P., & Berger, M. A. 2006, A&A, 452, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pariat, E., Démoulin, P., & Nindos, A. 2007, Adv. Space Res., 39, 1706 [NASA ADS] [CrossRef] [Google Scholar]
- Petsov, A. A. 2008, J. Astrophys. Astron., 29, 49 [NASA ADS] [CrossRef] [Google Scholar]
- Petsov, A. A., Maleev, V. M., & Longcope, D. W. 2003, ApJ, 593, 1217 [NASA ADS] [CrossRef] [Google Scholar]
- Potts, H. E., Barrett, R. K., & Diver, D. A. 2003, Sol. Phys., 217, 69 [NASA ADS] [CrossRef] [Google Scholar]
- Priest, E. R., & Forbes, T. G. 1990, Sol. Phys., 126, 319 [NASA ADS] [CrossRef] [Google Scholar]
- Qiu, J., Hu, Q., Howard, T. A., & Yurchyshyn, V. B. 2007, ApJ, 659, 758 [NASA ADS] [CrossRef] [Google Scholar]
- Seehafer, N. 1990, Sol. Phys., 125, 219 [NASA ADS] [CrossRef] [Google Scholar]
- Sheeley, N. R., Jr., Walters, J. H., Wang, Y.-M., & Howard, R. A. 1999, J. Geophys. Res., 104, 24739 [Google Scholar]
- Tian, L., & Alexander, D. 2008, ApJ, 673, 532 [NASA ADS] [CrossRef] [Google Scholar]
- Yeates, A. R., & Mackay, D. H. 2009, ApJ, 699, 1024 [Google Scholar]
- Zuccarello, F. P., Soenen, A., Poedts, S., Zuccarello, F., & Jacobs, C. 2008, ApJ, 689, L157 [NASA ADS] [CrossRef] [Google Scholar]
Online Material
Appendix A: Magnetic flux and magnetic helicity evolution
The plots in Figs. A.1 and A.2 report the evolution of total, positive, and negative magnetic flux expressed in 1022 Mx (first row), helicity injection rate dH/dt in 1042 Mx2 h-1 (second row), accumulated helicity injection
in 1042 Mx2 (third row), ratio
(fourth row), positive (H+) and negative (H-)
accumulated helicity injection (fifth row) for NOAA 8948, 9114+9115 and
9393+9394, belonging to the Class I and for NOAA 10162, 10229 and
10365, belonging to Class II.
NOAA 8948 (Fig. A.1,
leftt panel): the total, positive, and negative magnetic flux are
slightly decreasing during the period analyzed. However, if we limit
our attention to the time interval preceding the CME, we note a slight
increase of both positive and negative magnetic fluxes. The helicity
injection rate dH/dt (second row) is prevalently
positive, with the exception of the last part of the analyzed period.
The accumulated helicity injection (third row) is initially
characterized by increasing positive values and in the last part of the
analyzed period by a decreasing trend. The
ratio (fourth row) shows an increasing trend till the time of the corresponding increase in H-
(see the graph in the fifth row), but not correlated with the CME
occurrence. The graph in the fifth row shows that the main contribution
to the helicity injection is due to the increase of H+, while there is an increase of H- only during the last part of the analyzed period. The CME had a velocity of 409 km s-1, an acceleration of 2.5 m s-2, and was associated to an M3.1 flare.
NOAA 9114+9115 (Fig. A.1,
middle panel): the total magnetic flux is slightly increasing during
the analyzed period and in particular just before the CME occurrence;
the helicity injection rate (second row), which is always characterized
by negative values, shows an evident increasing trend before the CME
occurrence (till the value of
Mx2) and a phase characterized by lower values (till
Mx2) after
it. The accumulated helicity injection (third row) shows a different
trend before and after the CME occurrence. This behavior is also
evident in the curve showing the ratio
.
The graph reported in the fifth row indicates that the contribution to
is only due to H-. Note that the CME had a velocity of 720 km s-1, an acceleration of 2.8 m s-2 and it was associated with a C2.3 flare and an erupting filament.
NOAA 9393+9394 (Fig. A.1,
right panel): the total magnetic flux shows an increase during the
analyzed period (in this case, the positive flux increases, whereas the
negative flux decreases); the helicity injection rate (second row) is
predominantly characterized by negative values. The curve showing the
accumulated helicity injection
indicates that there is a continuous increase of helicity, always
characterized by negative values, also after the CME. Interestingly,
the trend shows a slight variation a few hours before the CME. A
similar trend is shown by the
values (fourth row). The graph reported in the fifth row indicates that the main contribution to
is due to H-, and that H+
starts to increase some hours before the CME occurrence. In this case,
the CME, associated to an X1.7 flare, had a velocity of
942 km s-1 and an acceleration of 3.5 m s-2.
![]() |
Figure A.1: Same as in Fig. 3 for ARs 8948, 9114+9115 and 9393+9394, belonging to Class I (gradual CMEs). |
Open with DEXTER |
NOAA 10162 (Fig. A.2,
left panel): the total magnetic flux is almost constant during the
analyzed period (this is due to a smooth increase of the positive flux
and to a contemporary decrease of the negative flux); the helicity
injection rate shows a very spiky behavior, generally characterized by
positive values. The accumulated helicity shows an almost monotonic
increase, with a small plateau (reflecting the dH/dt value) at the CME
occurrence. The ratio
shows an increasing trend with a similar plateau in coincidence with the CME occurrence. H+ has a trend similar to H
and
,
while H- remains constant (assuming values very close to zero). In this case, the CME had a velocity of 870 km s-1, a deceleration of 23.4 m s-2 and was associated with an erupting filament.
NOAA 10229 (Fig. A.2,
middle panel): in this active region the total magnetic flux decreases
during the analyzed period (the negative flux shows a stronger decrease
than the positive flux), whereas the helicity injection rate generally
shows negative values, with some abrupt changes to positive values. The
accumulated helicity, and correspondingly the ratio
are monotonically increasing. In this case H- (see the graph in the fifth row) shows a continuous increase, while H+
has only a slight increase just before the CME occurrence and again
some hours later. The velocity and deceleration of the CME were
1092 km s-1 and 36.2 m s-2, respectively, and the associated flare was of M2.7 class.
NOAA 10365 (Fig. A.2,
right panel): in the first row we can see that the total magnetic flux
shows an evident increase during the analyzed period. However, a more
careful inspection of this graph indicates that there is a clear
increase before the first CME and almost no increase before the second
CME. The helicity injection rate dH/dt is initially
characterized by positive values and does not show significant
variation after the first CME, while it changes, assuming also negative
values, after the second CME. This behavior is clearly shown in the
accumulated helicity (third row) and even more evidently in the
trend, especially as far as the second event is concerned. Similarly, H+
(see the graph in the fifth row), which increases during the entire
analyzed period, shows a small plateau after the second CME, and later
starts to increase again, while H-, which was almost zero until this event, starts to increase. The first CME had a velocity of 509 km s-1, a deceleration of 14.6 m s-2, and was associated with an M1.6 flare, while the second CME presented a velocity of 964 km s-1, a deceleration of 9.6 m s-2, and was associated with an X1.3 flare.
![]() |
Figure A.2: Same as in Fig. 3 for ARs 10162, 10229 and 10365, belonging to Class II (impulsive CMEs). |
Open with DEXTER |
Footnotes
- ... ejections
- Appendix is only available in electronic form at http://www.aanda.org
All Tables
Table 1: Data on the halo CMEs and analyzed relevant eruptive events.
Table 2: Characteristics of the active regions, halo CMEs and analyzed relevant eruptive events.
All Figures
![]() |
Figure 1: a) Difference image showing the halo CME observed by LASCO-C2 on 2002 july 15 at 20:30:05 UT, superposed on the EIT difference image relevant to 20:24 UT; b) SOHO/EIT full disk image acquired on 2002 july 15 at 20:24:10 UT, showing the location on the disk of the CME source region. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: a) Zoom of the active region NOAA 10030 where the CME was initiated; the image was acquired by SOHO/EIT during the main phase of the flare; b) MDI/SOHO magnetogram showing the magnetic configuration of NOAA 10030 on 15 July 2002. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
From top to bottom: evolution of the magnetic flux (continuous
line indicates the average total flux, marked line the positive flux
and dotted line the negative flux), magnetic helicity injection rate,
total helicity injection, ratio
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: Same as in Fig. 3 for ARs 10030 and 10069+10077, belonging to Class II (impulsive CMEs). |
Open with DEXTER | |
In the text |
![]() |
Figure A.1: Same as in Fig. 3 for ARs 8948, 9114+9115 and 9393+9394, belonging to Class I (gradual CMEs). |
Open with DEXTER | |
In the text |
![]() |
Figure A.2: Same as in Fig. 3 for ARs 10162, 10229 and 10365, belonging to Class II (impulsive CMEs). |
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.