A&A 382, 650-665 (2002)
P. Démoulin1 - C. H. Mandrini2 - L. van Driel-Gesztelyi 1,3,4,5 - B. J. Thompson6 - S. Plunkett7 - Zs. Kovári4 - G. Aulanier1 - A. Young6
1 - Observatoire de Paris, section Meudon, DASOP, URA 2080 (CNRS), 92195 Meudon Principal Cedex, France
2 - Instituto de Astronomía y Física del Espacio, IAFE, CC. 67 Suc. 28, 1428 Buenos Aires, Argentina
3 - Mullard Space Science Laboratory, Univ. College London, UK
4 - Konkoly Observatory, Hungary
5 - Centre for Plasma Astrophysics, K.U. Leuven, Celestijnenlaan 200B, 3001 Heverlee, Belgium
6 - Goddard Space Flight Center - NASA, USA
7 - USRA, Naval Research Laboratory, Washington, DC 20375, USA
Received 3 October 2001 / Accepted 15 November 2001
An isolated active region (AR) was observed on the Sun during seven rotations, starting from its birth in July 1996 to its full dispersion in December 1996. We analyse the long-term budget of the AR relative magnetic helicity. Firstly, we calculate the helicity injected by differential rotation at the photospheric level using MDI/SoHO magnetograms. Secondly, we compute the coronal magnetic field and its helicity selecting the model which best fits the soft X-ray loops observed with SXT/Yohkoh. Finally, we identify all the coronal mass ejections (CMEs) that originated from the AR during its lifetime using LASCO and EIT/SoHO. Assuming a one to one correspondence between CMEs and magnetic clouds, we estimate the magnetic helicity which could be shed via CMEs. We find that differential rotation can neither provide the required magnetic helicity to the coronal field (at least a factor 2.5 to 4 larger), nor to the field ejected to the interplanetary space (a factor 4 to 20 larger), even in the case of this AR for which the total helicity injected by differential rotation is close to the maximum possible value. However, the total helicity ejected is equivalent to that of a twisted flux tube having the same magnetic flux as the studied AR and a number of turns in the interval [0.5,2.0]. We suggest that the main source of helicity is the inherent twist of the magnetic flux tube forming the active region. This magnetic helicity is transferred to the corona either by the continuous emergence of the flux tube for several solar rotations (i.e. on a time scale much longer than the classical emergence phase), or by torsional Alfvén waves.
Key words: Sun: corona - Sun: coronal mass ejections (CMEs) - Sun: magnetic fields - solar-terrestrial relations
CMEs are thought to be the last physical mechanism in the chain of processes which transfer magnetic flux and helicity from the base of the convective zone (from the dynamo region) into the interplanetary medium. Both magnetic flux and helicity are injected from the convective zone into the corona in a systematic, though bursty, way. Magnetic flux reverses its sign with the solar cycle while magnetic helicity does not (in a given solar hemisphere). Because magnetic helicity dissipation is a very slow process (Berger 1984), the Sun can avoid endless accumulation of helicity in the solar corona only by ejecting the accumulated helicity via the launch of CMEs (Low 1997).
Theoretical MHD models of coronal arcades have shown that magnetic shear exceeding a critical threshold can lead to their destabilization (Mikic et al. 1994; Mikic & Linker 1997), indicating that they can be considered candidates for the initiation of CMEs. Reconnection occurring in a current sheet within the arcade could then produce not only newly closed arcades, but also a twisted flux rope which is expelled from the corona. Several other CME models invoke the presence of twisted flux ropes in the corona before the ejection. These models show that a slow photospheric evolution (Forbes & Priest 1995; Lin et al. 1998), or a fast injection of new magnetic flux (Chen 1996), can lead to the ejection of the magnetic flux rope. All these models, thus, suggest that the observed interplanetary magnetic clouds are due to the ejection of magnetic flux ropes from the corona (seen as CMEs).
Interplanetary magnetic clouds are formed by twisted magnetic flux tubes, as revealed by "in situ'' measurements of their magnetic fields (Burlaga 1981; Burlaga et al. 1988; Lepping et al. 1990; Hidalgo et al. 2000; Shimazu & Marubashi 2000). These clouds have often been associated with disappearing filaments on the Sun (Marubashi 1986; Burlaga et al. 1998; Srivastava et al. 1998). The magnetic flux rope, together with the disappearing filament, is thought to be ejected from the Sun in the form of a CME that traverses the interplanetary medium carrying magnetic flux and helicity with it. The computation of the helicity in magnetic clouds can give us clues concerning the amount of helicity shed by CMEs, a quantity which is useful to compare to the helicity injected and stored in the corona.
The life of active regions (ARs) reflects the evolution of the magnetic field that emerges through the photosphere, interacts with the pre-existing magnetic environment, becomes dispersed by large-scale flows and supergranular or granular convection and, furthermore, is sheared by the differential rotation. Once an AR appears at the photospheric level, it reaches its maximum development in a few days. Compared to this short emergence period, its decay time can be much longer, of the order of months (see the review by van Driel-Gesztelyi 1998). During this decay phase flaring, which might have been prolific and violent during the AR growth, diminishes considerably, while CMEs may well occur accompanied by the eruption of long filaments. The total number of CMEs produced by one single active region during its entire lifetime, from the cradle to the grave, can give us information about the total amount of magnetic helicity the AR contains. However, during solar maximum, both because there are several active regions on the Sun and because a significant part of the CMEs are launched outside ARs, it is difficult to make a full CME count for a source region, especially, when it is behind the limb of the Sun. Periods of low solar activity provide a better opportunity for such study.
Solar activity during the second half of 1996 was governed by a single AR, which was located on the southern hemisphere. This region was identified as NOAA AR 7978 when it first appeared on July 6, 1996. From this month and for the next five solar rotations (until November 1996), only a few small, short-lived ARs were observed on the solar disk besides AR 7978. Thus, solar activity at that period came mainly from AR 7978, and a direct comparison between the evolution of the AR magnetic field and its flare/CME production can be made. This long-term (six months) study gives us the unique opportunity to investigate the origin of the helicity in the field, as generated by differential rotation or/and brought up from the convective zone by twisted flux emergence. Furthermore, the simple bipolar configuration of AR 7978, which is maintained throughout its life, makes the results clearer, since the studied configuration is relatively close to theoretical configurations: in the convective zone (emergence of a single flux tube), in the corona (no interaction with neighboring active regions) and for the injection of helicity by differential rotation (a bipolar field is the simplest configuration, still twist and writhe magnetic helicity injected by differential rotation nearly cancel as shown by Démoulin et al. 2001).
This paper is organized as follows: in Sect. 2
we describe the data used in our long-term analysis; then we present a
global description of the evolution of AR 7978 at different
atmospheric levels, together with a study of the evolution of the
magnetic shear based on models of the coronal magnetic field
(preliminary results can be found in van Driel-Gesztelyi et al.
1999; Mandrini et al. 2000). Magnetic
helicity is defined in Sect. 3. In
Sect. 4, we briefly describe the injection
of magnetic helicity in the corona by differential rotation. After
that, we present the computation of magnetic helicity in the coronal
field (Sect. 5) and in interplanetary magnetic
clouds (Sect. 6). We apply the results to AR 7978
and Sect. 7 describes its magnetic helicity
budget. Finally, in Sect. 8, we discuss our
results. We find that differential rotation is neither sufficient to
generate the helicity found in the coronal field, nor the helicity
ejected from the Sun.
|Figure 1: Evolution of the longitudinal magnetic field during the emergence of AR 7978 (white: positive polarity, black: negative polarity). The spatial coordinates have been corrected from projection effects. The negative polarity (black) which grow in front of the main positive polarity is the classical evolution of the longitudinal magnetic field as an AR rotates towards the west limb.|
|Open with DEXTER|
|Figure 2: Evolution of the longitudinal magnetic field during the second rotation of AR 7978 over the solar disk (same convention as in Fig. 1). The remnant field of the second emergence is present as a tilted bipole (the first bipole being oriented nearly east-west).|
|Open with DEXTER|
There was a second episode of flux emergence in the AR, which started while it was on the invisible side of the Sun and continued for a few days after the AR rotated onto the disk (Fig. 2). The new flux emergence increased the magnetic flux measured during the first rotation (total unsigned flux of Mx) by about 50% ( Mx), see Fig. 4. The peak flux was reached during the third rotation ( Mx), this was due to an increase of the negative flux creating a roughly 10% flux imbalance in the AR. The imbalance was plainly due to the fact that the AR emerged in a dominantly negative polarity magnetic environment and, as its flux was getting diffused and its magnetic area grew, more and more of this pre-existing negative flux was included in the measurements (van Driel-Gesztelyi et al. 1999). Later the total flux decreased slowly, and even increased during the last two rotations due to an increasing level of small-scale new flux emergence and the ever-increasing area.
After the first two rotations the global deformation of the flux
concentrations were mainly due to the effect of diffusion and
differential rotation. Figure 3, right column,
illustrates the long-term evolution of the magnetic field of AR 7978
at successive central meridian passages (CMPs). In the left column of
the same figure, we show magnetic maps rotated to the time of the
following CMP map, applying differential rotation. Comparing maps in
the left and right columns, we observe the same tilt of the magnetic
inversion line implying that the global deformation of the AR, as seen
by the longitudinal magnetograms of MDI, was just the result of
|Figure 3: On the right: SOHO/MDI magnetic maps showing the long-term evolution of AR 7978. On the left: magnetograms of the right column (first date) transformed by differential rotation after one solar rotation (second date). The difference between magnetic maps of the same line represents the evolution of the AR without the distortion of the differential rotation.|
|Open with DEXTER|
|Figure 4: Evolution of the total unsigned (continuous line), positive and negative (dashed lines) photospheric magnetic flux in AR 7978. The longitudinal magnetic field of MDI at the central meridian passages is used.|
|Open with DEXTER|
|No.||first day||Flares (GOES)||CMEs||CMEs|
|of rot.||of rot.||X||M||C||B||observed||corrected|
|1st||7 July 1996||1||2||14||16||8||11|
|2nd||3 August 1996||-||-||1||17||5||5|
|3rd||30 August 1996||-||-||-||2||2||3|
|4th||25 September 1996||-||-||-||-||5||5|
|5th||23 October 1996||-||-||-||1||3||4|
|6th||19 November 1996||-||-||-||-||3||3|
|Figure 5: Yohkoh/SXT observations and linear force-free computations of the coronal field using SOHO/MDI magnetic maps. The top rows correspond to Aug. 3 and 30, while bottom rows to Sep. 25 and Oct. 23, respectively. The value is selected to have the best global match between the SXT loops and the computed field lines (see Sect. 2.5). We have drawn with thick (thin) continuous lines the field lines computed taking the highest (lowest) value of , and have represented them on the same figure (see values in Table 3). The abscissas and ordinates in each map are expressed in Mm. Isocontours (100, 500, 1000 G) of the line of sight magnetic field have been drawn with continuous/dashed lines for positive/negative values.|
|Open with DEXTER|
|Date||angular size||Importance||Remarks||Rotation phase|
|1996/07/15 04:35 (10:46)||14||43||2||0.30|
|1996/07/16 05:12 (08:27)||29||47||4||0.33|
|1996/08/19 14:09 (11:18)||63||82||9||1.61|
|1996/10/05 09:45 (09:37)||253||161||9||T||3.35|
|1996/10/19 15:00 (17:39)||192||206||10||3.88|
|1996/10/31 07:00 (06:25)||13||2||4.30|
|1996/11/07 09:00 (07:25)||22||36||3||4.56|
|1996/11/07 19:05 (23:20)||360||10||H||4.58|
|1996/11/24 17:40 (17:52)||57||77||4||5.22|
|1996/11/26 16:57 (15:35)||39||41||10||T||5.29|
|1996/11/28 17:00 (16:50)||79||78||9||T||5.37|
CME activity, which was at first related in several cases to flare events, continued at a high level until the seventh rotation of the AR (CMP on Dec. 16, 1996). CMEs originating from AR 7978 were carefully identified in full-disc EIT images taking as proxies the presence of dimmings and/or loops in expansion. Furthermore, LASCO C2 high resolution movies were constructed in order to confirm the existence of a CME associated with the AR, and to determine its characteristics. When AR 7978 was on the far side of the Sun, we identified the related CMEs in a more subjective way, based on the position angle and the characteristics of the CME in the coronagraph observations compared to the location of the AR. This task resulted easier because of the low solar activity and the fact that, in the period from August to October 1996, AR 7978 was the only major AR on the Sun. Still, a large fraction (2/3) of CMEs observed during this period did not originate from AR 7978; so, though this AR was alone on the Sun, we still needed a careful identification of the CME source region for each and every CME.
Table 2 lists the observed CMEs linked to AR 7978 from its appearance and until the seventh rotation. The angular size of a CME as viewed with LASCO gives an indication of the CME spatial size, but it is biased by projection effects. These include an increase of the angular size, as CMEs are launched closer to the centre disk (on both sides of the Sun), and the orientation of magnetic configuration relative to the line of sight (see Gibson & Low 1998 for a model with a twisted flux tube viewed with different observing directions). Then, we also consider the brightness of the CME to estimate its importance. The subjective combination of the angular size and of the brightness is noted from 1 to 10 in Table 2. It can be seen that the AR produced 15 medium to large (>5) CMEs and 11 minor (<5) CMEs. This classification is not critical for the helicity budget, but it is rather a detail whether to take into account, or not, the "smaller'' CMEs in the budget.
For future space weather studies, we add the following complementary information. Most of the CMEs which originated from AR 7978 were not geoeffective. However, comparing our CME list with data corresponding to the geomagnetic index Dst (see ftp://ftp.ngdc.noaa.gov/STP/GEOMAGNETIC.DATA/INDICES/DST/), whose perturbation can be indicative of a geomagnetic storm (González et al. 1994, 1999), we have found that a very intense ( ) geomagnetic storm occurred four days after the CME on 19 October, 1996. This CME was associated with the disparition brusque of a filament as shown in H images obtained at Meudon Observatory. Figure 6 shows the evolution in LASCO C2 of this CME. A similar CME eruption linked to this AR has been studied by van Driel-Gesztelyi et al. (1998).
|Figure 6: An example of a CME which originated from AR 7978 as observed with LASCO C2. Running difference images show the evolution of the CME on October 19, 1996. Yohkoh/SXT images, taken with time difference of a few (10) min. to the times of the C2 observations, show the lower coronal evolution during the CME. Note the formation of a cusp structure in AR7978.|
|Open with DEXTER|
|Date||- South||- North|
|7 July, 1996||-||1.00|
|3 August, 1996||0.30||0.75|
|30 August, 1996||0.90||1.00|
|25 September, 1996||1.00||1.40|
|23 October, 1996||0.90||1.40|
|18 November, 1996||0.90||0.90|
We have found that not all SXT loops observed at a given time can be represented using the same value of . In general, a shear gradient exists in the North-South direction. Figure 5 illustrates the temporal evolution of the coronal shear from the second to the fifth rotation and, at the same time, the just mentioned shear gradient. In Table 3 we have listed the maximum and minimum values of determined for each rotation; these values correspond to the North and South portion of the AR, respectively. It can be seen that increases from the first to the fourth rotation, after that it starts decreasing.
The magnetic helicity of a field
within a volume V is
Berger & Fields (1984) have shown that for cases where
one can define a relative magnetic helicity ()
helicity of a reference field ,
having the same
distribution of Bn on S:
Because of the very different physical conditions
present in these three regions, a different procedure is required to
estimate the gauge-independent magnetic helicity
in each of
them: for twisted flux tubes in the convection zone see
Appendix A, for the corona see
Sect. 5 and for magnetic clouds see
Sect. 6. Magnetic helicity is also transferred by
plasma flows from one region to a neighbor one. We analyze in the next
section the transfer from the convective zone to the corona.
|No. of rot.||Date|
|1st||0.2 (2)||[16, 64]||(7)|
|3 Aug.||[5, 11]|
|30 Aug.||[17, 23]|
|25 Sep.||[9, 12]|
|23 Oct.||[4, 6]|
is well preserved under solar conditions, the
only way helicity can be modified inside V is because of helicity
flux crossing the boundary S. The change of relative helicity is
written as (Berger & Fields 1984):
For simplicity, let's assume that the photosphere is locally planar
at the scale size of an AR (see Berger & Ruzmaikin 2000
for a spherical
as a function of the Bndistribution, Berger (1984, 1988) derived an
that depends only on observable photospheric quantities
(Bn and ):
Considering only horizontal motions, the elementary fluxes
are time-independent (though they can
still change in shape or field strength) and Eq. (10)
can be integrated:
The total helicity injected into AR 7978 by differential rotation is Mx2. This value is in close agreement with the estimate of the maximum total helicity injected in a dipole that is initially parallel to the equator (1043 Mx2) having a magnetic flux of 1022 Mx (see Fig. 5 in Démoulin et al. 2001). In fact, AR 7978 keeps its magnetic coherence long enough so that differential rotation has nearly the time to reach the maximum injection of helicity; is one order of magnitude smaller during the 6th (Nov.) rotation compared to its value during the 2nd and 3th rotations. These results are in agreement with DeVore (2000).
), Eq. (12) shows that
grows proportionally with
for a given magnetogram.
becomes very close to ,
the factor l at the
denominator decreases (this effect being strongest for harmonics
(1,0) and (0,1)); then, the helicity increases rapidly (see Fig. 1 in Berger 1985, though this figure corresponds to a spherical
geometry the effect is the same in Cartesian geometry). In the limit
the helicity is infinite. Therefore,
when the value of
used in the model of the coronal field is
close to ,
the computed helicity can be an overestimation
of the real helicity value. A way to avoid this is to use a
linearized form for :
Apart from computing helicity values within an interval (see Sect. 5.2), it is difficult to estimate an error bar for the computed helicity because the observed coronal field is not fully relaxed to a linear force-free field. The general property of a linear force-free field is to make the shorter (resp. longer) field lines less (resp. more) sheared than the observed coronal loops. This implies that the linear force-free field has very low (resp. too large) magnetic helicity in the small scales (resp. large scales) compared to the coronal field (see e.g. Schmieder et al. 1996). This property is rooted in the inverse cascade (towards the large spatial scales) of magnetic helicity which occurs when a magnetic field relaxes to a linear force-free field (minimizing the magnetic energy under the constraint of a fixed total relative helicity). Presently, one can only suppose that the over/under estimate of helicity at large/small scales compensates well enough to keep the right order of magnitude because the value of is selected so that the linear force-free field best represent the observed coronal loops.
We have computed the coronal of AR 7978 for the magnetic field models obtained from MDI magnetograms and shown in Fig. 5. The calculation has been done considering a computational box centered in the AR and having the same extension in all cases. As shown in Fig. 5 and Table 3, the value of varies along the AR (from North to South), and also from rotation to rotation. In the third column of Table 4 we list an interval for per rotation, where the minimum (resp. maximum) values have been obtained from the minimum (resp. maximum) values of in each rotation in the linearized approach (Eq. (13)). Therefore, these values are likely to be a lower bound in the estimation of the coronal . However, we want to stress that in those cases where was close to , the computed from Eq. (12) was at most a factor 1.5 larger than the linearized value computed from Eq. (13). The value between parenthesis for Nov. 18 rotation is only indicative. It was obtained for the coronal field model of the southern portion of AR 7978, the SXT loops in the northern portion were so faint that it was not possible to model them.
In order to estimate the degree of confidence one can have in the given by the linear force-free field, we have done the following test. On August 30 rotation, AR 7978 has such a simple magnetic configuration that applying a hypothetic shearing motion parallel to its inversion line, and opposite to the observed shear, we can bring it to an approximately potential configuration. Doing so, we can estimate the amount of helicity coming from this shearing motion; we have found that this last value is similar to the computed from the linear force-free field.
Can the helicity injected in AR 7978 by differential rotation be the origin of the estimated coronal helicity? On Aug. 3 rotation the value of the minimum coronal helicity, Mx2(Table 4), is still higher than the maximum value ( Mx2) given by differential rotation acting on AR 7978 and its neighbor decaying ARs (from July 7 to Aug. 3). When we restrict the computation of the injected helicity only to AR 7978 the obtained result is one order of magnitude smaller. Later on, the minimum coronal helicity is Mx2 on Aug. 30, increasing in Mx2 compared with Aug. 3, while the injection by differential rotation during that period is only of Mx2, so a factor 4 smaller. If we rather consider for the variation of the coronal helicity the average between the minimum and maximum values (see Table 4), we find comparing the same two cases (July 7 to Aug. 3 and Aug. 3 to Aug. 30) that the helicity injected by differential rotation is always a factor 4 lower. After that, the coronal helicity decreases, but we have also to consider the helicity ejected from the Sun (see Sect. 6) when analyzing the helicity budget of the AR.
All these results indicate that differential rotation does not inject enough magnetic helicity into the coronal field to account for the helicity found there, even if we take the minimum of the estimates for the coronal helicity.
A series of studies, related to the magnetic field in
interplanetary clouds, has shown that they are basically twisted flux
tubes expanding and moving outwards in the interplanetary medium.
Moreover, "in situ'' measurements of the magnetic field are globally
well-described by a linear-force free field, suggesting that
Taylor's relaxation (1974) was efficient after their launch from the
corona (Lepping et al. 1990, 1997; Marubashi
1997; Osherovich & Burlaga 1997;
Shimazu & Marubashi 2000). Indeed, only one harmonic of the
linear force-free field
in cylindrical geometry is enough to describe these measurements of
! This clearly contrasts with the large number of harmonics
needed to compute the coronal field. The magnetic field of
interplanetary clouds, in general, is well-described in cylindrical
First, let us consider the evolution from the second to the sixth rotation when AR 7978 was dipole-like and nearly alone on the Sun. AR 7978 ejected at least 18 CMEs during that period (see Table 2). Multiplying that number of CMEs by the above estimate of Mx2 per CME, we find that a magnetic helicity Mx2 was ejected from AR 7978 from the second to the sixth rotation while differential rotation only produced Mx2 (see Table 4). Then, the magnetic helicity ejected is a factor 4 larger than the helicity injected by differential rotation. We have also computed the injection of helicity for AR 7978 and its surrounding during the 1st (July) rotation (Table 4). Including this rotation, we see that the deficit of helicity injected by differential rotation is still higher compared to the estimated ejected helicity (at most Mx2 compared to Mx2ejected for 26 CMEs). So we get at least a factor 5 when we consider their total budgets over the full AR lifetime.
We now refine this estimate taking into account the time-dependence of the CME productivity. Both histograms of Fig. 7 show a clear peak around the West limb (longitudes 45-135 degrees) while a minimum is observed on the East limb (longitudes 225-315 degrees). The physics of the light scattering is the same in both cases so this difference should be an intrinsic property of AR7978. The first flux emergence episode (in the first half of July) occurred dominantly when AR 7978 was on the western hemisphere and behind the West limb (see Sect. 2.2). Then the AR produced many CMEs, with a mean frequency of 0.6 CMEday-1. During the second flux emergence, which started while the AR was behind the limb but probably already on the eastern hemisphere, and lasted for about 10 days, we did not observe CMEs at all. However, there were regular data gaps during this period lasting 9 hours per day on average and adding up to 4.5 days. Supposing a similar CME frequency as during the first flux emergence period, we probably missed 2-3 CMEs around the East limb. During the decay phase of the AR the CME frequency decreased to 0.2 CMEday-1. Since the total length of the data gaps during this period was 11.5 days, we may have missed 2 CMEs due to the gaps. All together, we find that of all CMEs associated with AR7978 are likely to have been missed.
While we clearly find this intrinsic temporal dependence in the
CME production, the number of CMEs launched both towards (longitudes
315-360 and 0- 45 degrees) and away (longitudes 135-225 degrees) from
Earth is approximately the mean of the number of CMEs launched in the
eastern and western quadrants! We conclude that we have no evidence
that CMEs are under-observed for some phase angle. Additionally, CMEs
from behind the solar limb, or CMEs which were not able to be
definitely identified with the AR were not included. Therefore,
considering both the data gaps and uncertainty in identification, the
total number of CMEs identified may be underestimated by 1/3.
|Figure 7: Top figure. Thick line: histogram of the observed CMEs coming from AR 7978 summed in the period July-December 1996. In the abscissa the longitude of AR is plotted with 12 bins of 30 degrees (bin 1: 15 degrees around CMP, bin 4: 15 degrees around the West limb). The histogram with thin lines has been corrected because of the data gaps assuming that the frequency of the CMEs is the same during data gaps than during observing times (see Sect. 6.3). Bottom figure. Histogram for the fraction of time when no LASCO data were available (defined as periods with no data during more than 3 hours), so when CMEs could have been missed.|
|Open with DEXTER|
The largest uncertainty in the helicity ejected indeed comes from the average helicity that one CME can carry away. Presently, this estimate can only be done from the magnetic measurements in magnetic clouds. It is certainly an important hypothesis to suppose that the mean helicity contained in the CMEs of AR 7978 corresponds to the mean helicity contained in a set of 18 magnetic clouds observed at a different period of time. For example, if only the medium to large CMEs (15 observed associated to AR 7978, see Sect. 2.4) are large enough to correspond to the magnetic clouds, the helicity ejected by the CMEs would be reduced by a factor 0.6. However, we have presently no quantitative means to relate the properties of the observed CMEs to those of the associated magnetic clouds. Then, consistent with DeVore (2000), we assume that the CMEs launched from AR 7978 carry in average a magnetic helicity corresponding to the 18 magnetic clouds analysed by Lepping et al. (1990). Additionally, none of the magnetic clouds observed by WIND/MFI correspond to the CMEs in our study, so explicit helicity comparisons are not possible.
Besides the intrinsic limitation of the model used to derive the global cloud parameters (see e.g. Vandas & Geranios 2001, for a comparative study of a set of models fitted to the data of a magnetic cloud), the length of the magnetic cloud is an unknown. Above we have used the conventional value of 0.5 AU following DeVore's suggestion. But indeed, bidirectional flows observed in interplanetary clouds indicate that they are probably still rooted in the Sun when observed at 1 AU (e.g. Richardson et al. 1991; Richardson 1997; Shodhan et al. 2000). This implies that the length of clouds should be at least 2 AU, a factor 4 larger than DeVore's assumption. At this point it is worth remembering that the measurements of the magnetic field in clouds are local, and such local measurements are extrapolated to the full magnetic cloud. We may expect a uniform distribution of the twist along the cloud as a consequence of both magnetic flux and current conservation along the flux tube. Therefore, considering that magnetic clouds are still rooted in the Sun we can estimate a total ejected helicity Mx2, a factor 20 larger than what differential rotation provides at most!
In Sects. 5.2, 6.2 and 6.3 we have shown that differential rotation
cannot generate enough magnetic helicity to provide either the coronal
helicity or the helicity ejected into the interplanetary medium.
Another source of helicity is then required. From the AR evolution
(Figs. 1-3) there are no obvious
horizontal photospheric motions (like shearing ones) which can provide
the needed helicity. We conclude that the helicity comes from the
last term on the right hand side of Eq. (7). The budget of
magnetic helicity per rotation writes:
The variation of the coronal helicity, , was first estimated from the minimum and maximum values derived with Eq. (13). Both values are coherent in magnitude, then we average them and subtract the contribution of the differential rotation. The last column of Table 4 gives this estimation of . Then, the sum of the last two columns of Table 4 gives an estimate for . It indicates the need of a major injection of magnetic helicity during the first two rotations. This is coherent with the two episodes of flux emergence observed in July and August and tells us that most of the magnetic helicity indeed emerged from the convective zone.
The budget of helicity for the period September to November is less certain. Remembering the large errors that we may have on and , it is plausible that, during that period, the launch of CMEs simply depletes the coronal helicity without any need of injecting helicity from below ( ). This would require that the helicity ejected on average in CMEs be low (lower than Mx2). If indeed magnetic clouds are more extended than 0.5 AU at 1 AU (because they are still linked to the Sun), their magnetic helicity could be much higher. In such a case the input of helicity from the convective zone should continue well after the "classical'' emergence period (characterized by an increase of photospheric magnetic flux). There is no observational evidence in AR 7978 to support or reject this hypothesis.
An alternative mechanism to this extended emergence, which does not require that the flux tube remains buoyant in the convective zone, has been considered by Longcope & Welsch (2000). The twisted flux tube in the convective zone creates a net torque on the photospheric plasma which cannot be compensated by the plasma forces. After a transient period, during which torsional Alfvén waves are transferring twist to the corona, a steady state is reached, where the torque of the sub-photospheric flux tube is balanced by the torque of the coronal magnetic field. In short, for the steady state, the coronal helicity is determined by the amount of twist present in the sub-photospheric part of the flux tube. Let us extend this model qualitatively to include CMEs. If the coronal helicity of the steady state is above the threshold of the global instability for the coronal field, a CME will occur, removing part of the helicity. Next, the imbalance of torque will again charge the coronal field with helicity typically in a day (see the time estimate of Longcope & Welsch 2000) and, the process will start again. In this approach, the coronal field is permanently, but in a time dependent manner, filled by magnetic helicity from the convective zone until the flux tube twist is exhausted or the flux tube is destroyed by convective motions.
At present, observations do not let us decide whether the injection of magnetic helicity is limited or not to the "classical'' emergence period. In order to probe this, a much more precise determination of the magnetic helicity is needed, in particular, of that contained in magnetic clouds. This could clearly constrain the mechanisms of flux emergence, but also the physics of CMEs, e.g. does the long-term persistence of the CME launch occur only because of the dispersion the AR photospheric magnetic field (leading to a coronal expansion of the magnetic configuration) or does it need a nearly continuous feeding from the convective zone?
The long-term evolution of a single active region on the Sun, with a large amount of data coverage during its six month evolution, was a unique opportunity to evaluate its magnetic helicity budget. The isolation of AR7978, together with its simple bipolar nature, allows us to analyse an AR with a configuration closest to theoretical models. This clearly facilitates a quantitative analysis of the helicity budget (e.g. minimizing the exchange of helicity with a neighbor AR, and avoiding the possible caveat of multipolar configurations where several magnetic bipoles interact). We put emphasis on the magnetic helicity because, even in resistive MHD, it is dissipated on time scales much longer than six months. Magnetic helicity can then be used as a trace to understand the transport of the magnetic field from sub-photospheric layers to the corona, and then its ejection into the interplanetary medium.
The main question addressed in this study is whether magnetic helicity mainly is injected into the corona by horizontal photospheric motions or not. Considering the magnetic field evolution, the main candidate is the differential rotation. It provides a helicity injection which has the right sign for AR 7978 (compared to the inferred coronal helicity). We have estimated the magnetic helicity present in the coronal magnetic field of the AR. Even with the minimum values, the injection of helicity by differential rotation is too small, by at least a factor 2.5, to explain the coronal helicity (Sect. 5.2). Then, we have estimated the magnetic helicity carried away by CMEs assuming a one to one correspondence to magnetic clouds and using the mean values measured for the physical parameters of the clouds. We find that differential rotation provides a magnetic helicity too small by at least a factor of 4 (Sect. 6.2). This factor can be as high as 20 if the magnetic clouds are still attached to the Sun when their apex reach 1 AU (Sect. 6.3). These two factors take into account only the CMEs observed with LASCO; their number might be slightly higher, up to 1/3, due the data gaps and identification problems. We conclude that, contrary to DeVore (2000), differential rotation is not able to provide the magnetic helicity ejected from AR 7978, but that this helicity has to come dominantly from the convective zone.
We agree with DeVore (2000) on the amount of magnetic helicity injected by differential rotation and our conclusion, opposite to DeVore, comes only from the number of CMEs associated to AR 7978. This AR is indeed a typical AR, it has a magnetic flux close to the average flux value for an active region (1022 Mx, as used by Devore), but AR 7978 has a factor 5 more CMEs than the average value used by DeVore. This difference is mainly due to the much higher sensitivity of LASCO compared to previous instruments (e.g. aboard Skylab, P78-1 and Solar Maximum Mission), to the nearly continuous observations of LASCO (only of data gaps) and also that at the minimum of the cycle different CME events can be more clearly distinguished.
The above estimates of the magnetic helicity show that the differential rotation can neither provide the required magnetic helicity to the coronal field (factor 2.5 to 4), nor to the field ejected to the interplanetary space (factor 4 to 20). It is instructive to rewrite the helicity in the "natural'' units of , where Mx is the magnetic flux of the active region. In this units, the maximum coronal helicity is in the interval [0.2,0.3], while the total helicity ejected is [0.5,2.0] and the total helicity input from differential rotation in AR 7978 is 0.1 (which is nearly the maximum possible value). In the same units, the magnetic helicity of a twisted flux tube making N turns is simply N (see Eq. (A.5)). Then, the ejected magnetic helicity could originate from the twist of the flux tube which formed AR7978.
Recent studies (e.g. Emonet & Moreno-Insertis 1998) have shown that a minimum critical twist is needed so that a buoyant 2.5-D flux-tube is not destroyed during its rise by the hydrodynamic wakes that develop behind it. The pinch angle of the field should be typically larger than when it starts leaving the base of the convective zone. More recent studies do not significantly modify this result (see Abbett et al. 2000, and references therein). From this, one can estimate the length L of the flux tube which has become unstable at the base of the convective zone and which later on has formed AR 7978, , where N is the number of turns and a is the flux tube radius. The radius a can be simply estimated by using a typical field strength of 105 G for flux tubes becoming unstable at the tachocline (recalling that the magnetic flux Mx). We find for the unstable flux tube a length Mm, showing that an N of the order of a few units is indeed plausible with the mechanism of the ondulary (Parker) instability.
We conclude from the budget of magnetic helicity that the magnetic flux tube forming AR 7978 continues to emerge, or at least to provide an upward propagation of helicity to the corona (see the end of Sect. 7), for a much longer time than the observed photospheric increase of magnetic flux. One way to further constrain the duration of the emergence of flux tubes is to improve magnetic helicity estimations.
We thank C. R. DeVore and J. A. Klimchuk for useful discussions. P.D. and C.H.M. acknowledge financial support from ECOS (France) and SETCIP (Argentina) through their cooperative science program (A97U01). L.v.D.G. and Zs.K. were supported by the Hungarian Government grants OTKA T-026165, T-032846. P.D. L.v.D.G. and Zs.K. acknowledge the Hungarian-French S&T cooperative program. The authors thank the SOHO/MDI, SOHO/EIT and SOHO/LASCO consortia for the SOHO data and the MSSL SURF for YOHKOH/SXT data. SOHO is a joint project by ESA and NASA. The CME catalog is generated and maintained by the Center for Solar Physics and Space Weather, The Catholic University of America in cooperation with the Naval Research Laboratory and NASA.
We can now compute the magnetic helicity in the volume V,
which is a section of the cylinder (
). Because the magnetic field crosses the boundary
of the volume (at
), one should use
Eq. (5) that defines the relative helicity. It is worth
noting that using rather the initial definition of Eq. (3)
(only suited for fields fully contained in a volume) gives the
obviously wrong result of a null magnetic helicity (independently of
the amount of twist)! However, Eq. (5) requires the
reference helicity of the potential field. It is much easier to use a
different, but equivalent, formulation for the relative helicity
following Berger (1999). Let's write
where, in Berger's
notation, "open'' means a field crossing the boundary S of V and
"closed'' means a field confined inside the volume V.
Berger (1999) shows that the relative magnetic helicity should be