Solar Flare Activity and Variability of Electric Current Helicity

1 INTRODUCTIONRecently Bao, Zhang, Ai, and Zhang (1999), using Huairou vector magnetograph data,have shown that the average current helicity (h.) or the curreflt helicity imbalance ph of activeregions change rapidly after so1ar flares. Up'an the onset of flares it tends to decrease for a fewhours and then to increase again, whereas ifQ some cases the flare promotes an increase in thecurrent helicity The observations led to tbe fol1owing conclusions: (1) raPid and substantialchanges of c…     

Basic Properties of Mutual Magnetic Helicity

We derive the magnetic helicity for configurations formed by flux tubes contained fully or only partially in the spatial domain considered (called closed and open configurations, respectively). In both cases, magnetic helicity is computed as the sum of mutual helicity over all possible pairs of magnetic flux tubes weighted by their magnetic fluxes. We emphasize that these mutual helicities have properties which are not those of mutual inductances in classical circuit theory. For closed configurations, the mutual helicity of two closed flux tubes is their relative winding around each other (known as the Gauss linkage number). For open configurations, the magnetic helicity is derived directly from the geometry of the interlaced flux tubes so it can be computed without reference to a ground state (such as a potential field). We derive the explicit expression in the case of a planar and spherical boundary. The magnetic helicity has two parts. The first one is given only by the relative positions of the flux tubes on the boundary. It is the only part if all flux tubes are arch-shaped. The second part counts the integer number of turns each pair of flux tubes wind about each other. This provides a general method to compute the magnetic helicity with discrete or continuous distributions of magnetic field. The method sets closed and open configurations on an equal level within the same theoretical framework.     

Multiwavelength Study of a Solar Eruption from AR NOAA 11112 I. Flux Emergence, Sunspot Rotation and Triggering of a Solar Flare

We analyze the multiwavelength observations of an M2.9/1N flare that occurred in the active region (AR) NOAA 11112 in the vicinity of a huge filament system on 16 October 2010. SDO/HMI magnetograms reveal the emergence of a bipole (within the existing AR) 50 hours prior to the flare event. During the emergence, both the positive and negative sunspots in the bipole show translational as well as rotational motion. The positive-polarity sunspot shows significant motion/rotation in the south-westward/clockwise direction, and we see continuously pushing/sliding of the surrounding opposite-polarity field region. On the other hand, the negative-polarity sunspot moves/rotates in the westward/anticlockwise direction. The positive-polarity sunspot rotates ≈?70^° within 30 hours, whereas the one with negative polarity rotates ≈?20^° within 10 hours. SDO/AIA 94 Å EUV images show the emergence of a flux tube in the corona, consistent with the emergence of the bipole in HMI. The footpoints of the flux tube were anchored in the emerging bipole. The initial brightening starts at one of the footpoints (western) of the emerging loop system, where the positive-polarity sunspot pushes/slides towards a nearby negative-polarity field region. A high speed plasmoid ejection (speed ≈?1197 km?s^?1) was observed during the impulsive phase of the flare, which suggests magnetic reconnection of the emerging positive-polarity sunspot with the surrounding opposite-polarity field region. The entire AR shows positive-helicity injection before the flare event. Moreover, the newly emerging bipole reveals the signature of a negative (left-handed) helicity. These observations provide unique evidence of the emergence of twisted flux tubes from below the photosphere to coronal heights, triggering a flare mainly due to the interaction between the emerging positive-polarity sunspot and a nearby negative-polarity sunspot by the shearing motion of the emerging positive sunspot towards the negative one. Our observations also strongly support the idea that the rotation can most likely be attributed to the emergence of twisted magnetic fields, as proposed by recent models.     

Role of Sunspot and Sunspot-Group Rotation in Driving Sigmoidal Active Region Eruptions

We study active region NOAA 9684 (N06L285) which produced an X1.0/3B flare on November 4, 2001 associated with a fast CME (1810 km s⁻¹) and the largest proton event (31 700 pfu) in cycle 23. SOHO/MDI continuum image data show that a large leading sunspot rotated counter-clockwise around its umbral center for at least 4 days prior to the flare. Moreover, it is found from SOHO/MDI 96 m line-of-sight magnetograms that the systematic tilt angle of the bipolar active region, a proxy for writhe of magnetic fluxtubes, changed from a positive value to a negative one. This signifies a counter-clockwise rotation of the spot-group as a whole. Using vector magnetograms from Huairou Solar Observing Station (HSOS), we find that the twist of the active region magnetic fields is dominantly left handed (α_best = −0.03), and that the vertical current and current helicity are predominantly negative, and mostly distributed within the positive rotating sunspot. The active region exhibits a narrow inverse S-shaped H_α filament and soft X-ray sigmoid distributed along the magnetic neutral line. The portion of the filament which is most closely associated with the rotating sunspot disappeared on November 4, and the corresponding portion of the sigmoid was observed to erupt, producing the flare and initiating the fast CME and proton event. These results imply that the sunspot rotation is a primary driver of helicity production and injection into the corona. We suggest that the observed active region dynamics and subsequent filament and sigmoid eruption are driven by a kink instability which occurred due to a large amount of the helicity injection.     

The rotation of sunspots in the solar active region NOAA 10930

The rotation of sunspots in the solar active region NOAA 10930 was investigated on the basis of the data on the longitudinal magnetic field and the Doppler velocities using magnetograms and dopplergrams taken with the Solar Optical Telescope installed aboard the HINODE mission. Under the assumption of axial symmetry, areally-mean vertical, radial, and azimuthal components of the magnetic field and velocity vectors were calculated in both sunspots. The plasma in the sunspots rotated in opposite directions: in the leading sunspot, clockwise, and in the following sunspot, counterclockwise. The magnetic flux tubes that formed sunspots of the active region on the solar surface were twisted in one direction, clockwise. Electric currents generated as a result of the rotation and twisting of magnetic flux tubes were also flowing in one direction. Azimuthal components of magnetic and velocity fields of both sunspot umbrae reached their maximum on December 11, 2006. By the start of the X3.4 flare (December 13, 2006), their values became practically equal to zero.     

On the possibility of investigation of the magnetic field structure in subphotospheric layers of the Sun according to observations of sunspot torsional oscillations

A method of investigation of the magnetic field structure in subphotospheric layers of the Sun has been developed. The method is based on observations of the torisonal oscillations of single sunspots. Characteristics of the torsional oscillations have been obtained from observations of the longitudinal magnetic field and radial velocities of seven single sunspots in the photospheric line Fe I λ5253 Å. The parameters of the torsional oscillations and magnetic tubes in the deep layers have been determined. The radius of the cross section of a magnetic flux tube forming a sunspot is greatest near the Sun’s surface and is approximately equal to the radius of a sunspot umbra. Down to the deeper layers, it decreases quite quickly. The longitudinal electric current appearing in the magnetic tube changes direction. The typical time of the current changes is determined by the period of the torsional oscillations. The intensity of the longitudinal magnetic field in the tube increases with depth. The Alfven wave velocity averaged over the length of a magnetic tube is tens or hundreds of times less than this velocity in a sunspot umbra. It decreases with an increase in the period of oscillations. A decrease in the Alfven wave velocity leads to an increase in the twisting angle of magnetic field lines.     

Magnetic Helicity Budget of Solar-Active Regions from the Photosphere to Magnetic Clouds

We have analyzed the long-term evolution of two active regions (ARs) from their emergence through their decay using observations from several instruments on board SoHO (MDI, EIT and LASCO) and Yohkoh/SXT. We have computed the evolution of the relative coronal magnetic helicity combining data from MDI and SXT with a linear force-free model of the coronal magnetic field. Next, we have computed the injection of helicity by surface differential rotation using MDI magnetic maps. To estimate the depletion of helicity we have counted all the CMEs of which these ARs have been the source, and we have evaluated their magnetic helicity assuming a one to one correspondence with magnetic clouds with an average helicity contain. When these three values (variation of coronal magnetic helicity, injection by differential rotation and ejection via CMEs) are compared, we find that surface differential rotation is a minor contributor to the helicity budget since CMEs carry away at least 10 times more helicity than the one differential rotation can provide. Therefore, the magnetic helicity flux needed in the global balance should come from localized photospheric motions that, at least partially, reflect the emergence of twisted flux tubes. We estimate that the total helicity carried away in CMEs can be provided by the end-to-end helicity of the flux tubes forming these ARs. Therefore, we conclude that most of the helicity ejected in CMEs is generated below the photosphere and emerges with the magnetic flux.     

Vector magnetogram and dopplergram observation of magnetic flux emergence and its explanation

During 23–28 August 1988, at the Huairou Solar Observation Station of Beijing Observatory, the full development process of the region HR 88059 was observed. It emerged near the center of the solar disk and formed a medium active region. A complete series of vector magnetograms and photospheric and chromospheric Dopplergrams was obtained. From an analysis of these data, combined with some numerical simulations, the following conclusions can be drawn. (1) The emergence of new magnetic flux from enhanced networks followed by sunspot formation is an interesting physical process which can be simply described by MHD numerical simulation. The phenomena accompanying it occur according to a definite law summarized by Zwaan (1985). The condition for gas cooling and sunspot formation seems to be transverse field strength > 50 G together with longitudinal field strength > 700 G. For a period of 4 to 5 hours, the orientation of the transverse field shows little change. The configuration of field lines may be derived from vector magnetograms. The arch filament system can be recognized as an MHD shock. (2) New opposite bipolar features emerge within the former bipolar field with an identical strength which will develop a sunspot group complex. Also, arch filament systems appear there located in the position of flux emergence. The neutral line is often pushed aside and curved, leading to faculae heating and the formation of a current sheet. In spite of complicated Dopplergrams, the same phenomena occur at the site of flux emergence as usual: upward flow appears at the location of the emerging and rapidly varying flux near the magnetic neutral line, and downdraft occurs over large parts of the legs of the emerging flux tubes. The age of magnetic emerging flux (or a sunspot) can be estimated in terms of transverse field strengths: when 50 G < transverse field < 200 G, the longitudinal magnetogram and Dopplergram change rapidly, which indicates a rigourously emerging magnetic flux. When the transverse field is between 200 and 400 G, the area concerned is in middle age, and some of the new flux is still emerging there. When the transverse field > 400 G, the variation of the longitudinal magnetogram slows down and the emerging arch becomes relatively stable and a photospheric Evershed flow forms at the penumbra of the sunspot.     

The X10 Flare on 29 October 2003: Was It Triggered by Magnetic Reconnection between Counter-Helical Fluxes?

Vector magnetograms taken at Huairou Solar Observing Station (HSOS) and Mees Solar Observatory (MSO) reveal that the super active region (AR) NOAA 10486 was a complex region containing current helicity flux of opposite signs. The main positive sunspots were dominated by negative helicity fields, while positive helicity patches persisted both inside and around the main positive sunspots. Based on a comparison of two days of deduced current helicity density, pronounced changes associated with the occurrence of an X10 flare that peaked at 20:49 UT on 29 October 2003 were noticed. The average current helicity density (negative) of the main sunspots decreased significantly by about 50%. Accordingly, the helicity densities of counter-helical patches (positive) were also found to decay by the same proportion or more. In addition, two hard X-ray (HXR) “footpoints” were observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) during the flare in the 50 – 100 keV energy range. The cores of these two HXR footpoints were adjacent to the positions of two patches with positive current helicity that disappeared after the flare. This strongly suggested that the X10 flare on 29 October 2003 resulted from reconnection between magnetic flux tubes having opposite current helicity. Finally, the global decrease of current helicity in AR 10486 by ∼50% can be understood as the helicity launched away by the halo coronal mass ejection (CME) associated with the X10 flare. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.     

On Formation of the Sigmoidal Structure in Solar Active Region NOAA 8100

Using soft X-ray images taken by the Soft X-ray Telescope on board Yohkoh, line-of-sight magnetograms taken by SOHO/MDI and vector magnetograms taken at Beijing Astronomical Observatory, we have studied the formation of the sigmoidal structure in active region NOAA 8100 on 3–4 November 1997. The sigmoidal structure appeared after the occurrences of a series of flares accompanied by new magnetic flux emergence. This implies that reconnection may play a role in formation of this sigmoid structure. We calculated the self-helicity (twist) and mutual helicity of the active region before and after the formation of the sigmoidal structure and found that the mutual helicity decreased. The twist of the sigmoidal structure was higher than the twist of the emerging magnetic flux and exceeded the critical twist for kink in stability. This result suggests that the reconnection increased the twist of magnetic flux tubes by converting mutual helicity to self-helicity, supporting the previous studies by Berger (1998, 1999).     

Solar Cycle Variation of Magnetic Flux Ropes in a Quasi-Static Coronal Evolution Model

The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free “potential-field” extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model – in particular the flux ropes – varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.     

A New Model for Propagating Parts of EIT Waves: A Current Shell in a CME

EIT waves are observed in EUV as bright fronts. Some of these bright fronts propagate across the solar disk. EIT waves are all associated with a flare and a CME and are commonly interpreted as fast-mode magnetosonic waves. Propagating EIT waves could also be the direct signature of the gradual opening of magnetic field lines during a CME. We quantitatively addressed this alternative interpretation. Using two independent 3D MHD codes, we performed nondimensional numerical simulations of a slowly rotating magnetic bipole, which progressively result in the formation of a twisted magnetic flux tube and its fast expansion, as during a CME. We analyse the origins, the development, and the observability in EUV of the narrow electric currents sheets that appear in the simulations. Both codes give similar results, which we confront with two well-known SOHO/EIT observations of propagating EIT waves (7 April and 12 May 1997), by scaling the vertical magnetic field components of the simulated bipole to the line of sight magnetic field observed by SOHO/MDI and the sign of helicity to the orientation of the soft X-ray sigmoids observed by Yohkoh/SXT. A large-scale and narrow current shell appears around the twisted flux tube in the dynamic phase of its expansion. This current shell is formed by the return currents of the system, which separate the twisted flux tube from the surrounding fields. It intensifies as the flux tube accelerates and it is co-spatial with weak plasma compression. The current density integrated over the altitude has the shape of an ellipse, which expands and rotates when viewed from above, reproducing the generic properties of propagating EIT waves. The timing, orientation, and location of bright and faint patches observed in the two EIT waves are remarkably well reproduced. We conjecture that propagating EIT waves are the observational signature of Joule heating in electric current shells, which separate expanding flux tubes from their surrounding fields during CMEs or plasma compression inside this current shell. We also conjecture that the bright edges of halo CMEs show the plasma compression in these current shells.     

Soft X-ray line profiles in the impulsive phase of electron-heated solar flares

We have measured the motion of facular points and granules in the same region near a decaying sunspot. It is found that both features move away across the moat surrounding the sunspot. The mean speed of facular points is larger than that of granules: 0.65 km s^–1 and 0.4 km s^–1, respectively. These results are consistent with previous measurements of the speed of bright network features and moving magnetic fields, as well as of non-magnetic photospherical material. They support models in which a decaying sunspot is at the center of a supergranule, whose horizontal motions sweep out granules and magnetic flux tubes associated to the facular points. It is also found that granules are dragged by supergranular motions away of the moat.Contributions from the Kwasan and Hida Observatories, University of Kyoto.A part of this work was done while one of the authors (R.M.) was staying at the Kwasan and Hida Observatories, University of Kyoto, Japan, as a JSPS research fellow.     

On the Shape of Force-Free Field Lines in the Solar Corona

This paper studies the shape parameters of looped field lines in a linear force-free magnetic field. Loop structures with a sufficient amount of kinking are generally seen to form S or inverse S (Z) shapes in the corona (as viewed in projection). For a single field line, we can ask how much the field line is kinked (as measured by the writhe), and how much neighbouring flux twists about the line (as measured by the twist number). The magnetic helicity of a flux element surrounding the field line can be decomposed into these two quantities. We find that the twist helicity contribution dominates the writhe helicity contribution, for field lines of significant aspect ratio, even when their structure is highly kinked. These calculations shed light on some popular assumptions of the field. First, we show that the writhe of field lines of significant aspect ratio (the apex height divided by the footpoint width) can sometimes be of opposite sign to the helicity. Secondly, we demonstrate the possibility of field line structures which could be interpreted as Z-shaped, but which have a helicity value sign expected of an S-shaped structure. These results suggest that caution should be exercised in using two-dimensional images to draw conclusions on the helicity value of field lines and flux tubes.     

An essay on sunspots and solar flares

The presently prevailing theories of sunspots and solar flares rely on the hypothetical presence of magnetic flux tubes beneath the photosphere and the two subsequent hypotheses, their emergence above the photosphere and explosive magnetic reconnection, converting magnetic energy carried by the flux tubes for solar flare energy.In this paper, we pay attention to the fact that there are large-scale magnetic fields which divide the photosphere into positive and negative (line-of-sight) polarity regions and that they are likely to be more fundamental than sunspot fields, as emphasized most recently by McIntosh (1981). A new phenomenological model of the sunspot pair formation is then constructed by considering an amplification process of these largescale fields near their boundaries by shear flows, including localized vortex motions. The amplification results from a dynamo process associated with such vortex flows and the associated convergence flow in the largescale fields.This dynamo process generates also some of the familiar “force-free” fields or the “sheared” magnetic fields in which the magnetic field-aligned currents are essential. Upward field-aligned currents generated by the dynamo process are carried by downward streaming electrons which are expected to be accelerated by an electric potential structure; a similar structure is responsible for accelerating auroral electrons in the magnetosphere. Depending on the magnetic field configuration and the shear flows, the current-carrying electrons precipitate into different geometrical patterns, causing circular flares, umbral flares, two-ribbon flares, etc. Thus, it is suggested that “low temperature flares” are directly driven by the photospheric dynamo process.     

Magnetic field, helicity and the 2000 July 14 flare in solar active region 9077

In this paper we analyse the non-potential magnetic field and the relationship with current (helicity) in the active region NOAA 9077 in 2000 July, using photospheric vector magnetograms obtained at different solar observatories and also coronal extreme-ultraviolet 171-Å images from the TRACE satellite.
We note that the shear and squeeze of magnetic field are two important indices for some flare-producing regions and can be confirmed by a sequence of photospheric vector magnetograms and EUV 171-Å features in the solar active region NOAA 9077. Evidence on the release of magnetic field near the photospheric magnetic neutral line is provided by the change of magnetic shear, electric current and current helicity in the lower solar atmosphere. It is found that the 'Bastille Day' 3B/5.7X flare on 2000 July 14 was triggered by the interaction of the different magnetic loop systems, which is relevant to the ejection of helical magnetic field from the lower solar atmosphere. The eruption of the large-scale coronal magnetic field occurs later than the decay of the highly sheared photospheric magnetic field and also current in the active region.     

Evolution of magnetic flux distribution of solar bipolar active regions

In this paper the magnetic flux distribution of bipolar active regions at the sunspot development stage is analyzed. It is shown that the ratio of the total sunspot area in an active region to the maximum one can be used as a characteristic of the development phase. Such a procedure allows combining the data attributed to different active regions for studying evolutionary changes. The expressions describing the evolution of magnetic flux distribution of bipolar active region were obtained and their interpretation with rise and descent of loop like magnetic flux tube leading to active region formation was justified.     

The Magnetic Field, Horizontal Motion and Helicity in a Fast Emerging Flux Region which Eventually forms a Delta Spot

We studied the behavior of magnetic field, horizontal motion and helicity in a fast emerging flux region NOAA 10488 which eventually forms a δ spot. It is found that the rotation of photospheric footpoints forms in the earlier stage of magnetic flux emergence and the relative shear motion of different magnetic flux systems appears later in this active region (AR). Therefore the emerging process of the AR can be separated into two phases: rotation and shear. We have computed the magnetic helicity injected into the corona using the local correlation tracking (LCT) technique. Furthermore we determined the vertical component of current helicity density and the vertical component of induction electric fields E_z = (V× B)_z in the photosphere. Particularly we have presented the comparison of the injection rate of magnetic helicity and the variation of the current helicity density. The main results are as follows: (1) The strong shear motion (SSM) between the new emerging flux system and the old one brings more magnetic helicity into the corona than the twisting motions. (2) After the maturity of the main bipolar spots, their twist decreases and the SSM becomes dominant and the major contributor of magnetic non-potentiality in the solar atmosphere in this AR. (3) The positions of the maxima of E_z (about 0.1 ∼ 0.2 V cm⁻¹) shift from the twisting areas to the areas showing SSMs as the AR evolved from the rotation phase to the shear one, but no obvious correlation is found between the kernels of Hα flare and E_z for the M1.6 flare in this AR. (4) The coronal helicity inferred from the horizontal motion of this AR amounts to −6 × 10⁴³ Mx². It is comparable with the coronal helicity of ARs producing flares with coronal mass ejections (CMEs) or helicity carried away by magnetic clouds (MCs) reported in previous studies (Nindos, Zhang, and Zhang, 2003; Nindos and Andrews, 2004). In addition, the formation of the δ configuration in this AR belongs to the third formation type indicated by Zirin and Liggett (1987), i.e., collision of opposite polarities from different dipoles, and can be naturally explained by the SSM.     

Sunspot proper motions and behaviour of rising magnetic flux tubes

Greenwich data for proper motions of the components of sunspot groups are analysed. We have found that only 4% of the 408 examined sunspot groups do not show systematic proper motions of the components at the beginning of their life. Interpreting proper motions as the results of ascending sub-photospheric magnetic flux tubes, information on characteristics of the tubes have been deduced. The influence of proper motions on the evaluation of the sunspot rotation rates is discussed.     

Soho Observations of the Structure and Dynamics of Sunspot Region Atmospheres

We present results from a study of the spatial distributions of line emission and relative line-of-sight velocity in the atmosphere above 17 sunspot regions, from the chromosphere, through the transition region and into the corona, based on simultaneous observations of ten EUV emission lines with the Coronal Diagnostic Spectrometer – CDS on SOHO. We find that the spatial distributions are nonuniform over the sunspot region and introduce the notation 'sunspot loop' to describe an enhanced transition region emission feature that looks like a magnetic loop, extending from inside the sunspot to the surrounding regions. We find little evidence for the siphon flow. Attention is given to the time variations since we observe both a rapid variation with a characteristic time of a few to several minutes and a slow variation with a time constant of several hours to 1 day. The most prominent features in the transition region intensity maps are the sunspot plumes. We introduce an updated criterion for the presence of plumes and find that 15 out of 17 sunspots contain a plume in the temperature range logT5.2–5.6. The relative line-of-sight velocity in sunspot plumes is high and directed into the Sun in the transition region. Almost all the sunspot regions contain one or a few prominent, strongly redshifted velocity channels, several of the channels extend from the sunspot plume to considerable distances from the sunspot. The flow appears to be maintained by plasmas at transition region temperatures, moving from regions located at a greater height outside the sunspots and towards the sunspot. The spatial correlation is high to moderate between emission lines formed in the transition region lines, but low between the transition region lines and the coronal lines. From detailed comparisons of intensity and velocity maps we find transition region emission features without any sign of coronal emission in the vicinity. A possible explanation is that the emission originates in magnetic flux tubes that are too cold to emit coronal emission. The comparisons suggest that gas at transition region temperature occur in loops different from loops with coronal temperature. However, we cannot exclude the presence of transition region temperatures close to the footpoints of flux tubes emitting at coronal temperatures. Regions with enhanced transition region line emission tend to be redshifted, but the correlation between line emission and relative line-of-sight velocity is weak. We extend our conditional probability studies and confirm that there is a tendency for line profiles with large intensities and red shifts (blue shifts) above the average to constitute an increasing (decreasing) fraction of the profiles as the wavelength shift increases.