Coronal Flux Rope Equilibria in Closed Magnetic Fields

Using a 2.5-dimensional ideal MHD model in Cartesian coordinates,we investigate the equilibrium properties of coronal magnetic flux ropes in background magnetic fields that are completely closed.The background fields are produced by a dipole,a quadrupole,and an octapole,respectively,located below the photosphere at the same depth.A magnetic flux rope is then launched from below the photo-sphere,and its magnetic properties,i.e,the annular magnetic fluxφp and the axial magnetic fluxφz,are controlled by a single emergence parameter.The whole sys-tem eventually evolves into equilibrium,and the resultant flux rope is characterized by three geometrical parameters:the height of the rope axis,the half-width of the rope,and the length of the vertical current sheet below the rope.It is found that the geometrical parameters increase monotonically and continuously with increasing φp and φz:no catastrophe occurs.Moreover,there exists a steep segment in the profiles of the geometrical parameters versus either φp or φz,and the faster the background field decays with height,the larger both the gradient and the growth amplitude within the steep segment will be.     

Magnetic Energy of Force-Free Fields with Detached Field Lines

Using an axisymmetrical ideal MHD model in spherical coordinates, we present a numerical study of magnetic configurations characterized by a levitating flux rope embedded in a bipolar background field whose normal field at the solar surface is the same or very close to that of a central dipole. The characteristic plasma β (the ratio between gas pressure and magnetic pressure) is taken to be sosmall (β= 10^-4) that the magnetic field is close to being force-free. The system as a whole is then let evolve quasi-statically with a slow increase of either the annular magnetic flux or the axial magnetic flux of the rope, and the total magneticenergy of the system grows accordingly. It is found that there exists an energy threshold: the flux rope sticks to the solar surface in equilibrium if the magneticenergy of the system is below the threshold, whereas it loses equilibrium if the threshold is exceeded. The energy threshold is found to be larger than that of thecorresponding fully-open magnetic field by a factor of nearly 1.08 irrespective as towhether the background field is completely closed or partly open, or whether the magnetic energy is enhanced by an increase of annular or axial flux of the rope.This gives an example showing that a force-free magnetic field may have an energy larger than the corresponding open field energy if part of the field lines is allowed tobe detached from the solar surface. The implication of such a conclusion in coronal mass ejections is briefly discussed and some comments are made on the maximum energy of force-free magnetic fields.     

Catastrophe of coronal magnetic flux ropes in partially open magnetic fields

Using a 2.5-D, time-dependent ideal MHD model in Cartesian coordinates, a numerical study is carried out to find equilibrium solutions associated with a magnetic flux rope in the corona. The ambient magnetic field is partially open, consisting of a closed arcade in the center and an open field at the flank. The coronal magnetic flux rope is characterized by its magnetic properties, including the axial and annular magnetic fluxes and the magnetic helicity, and its geometrical features, including the height of the rope axis, the halfwidth of the rope and the length of the vertical current sheet below the rope. It is shown that for a given partially open ambient magnetic field, the dependence of the geometrical features on the magnetic properties displays a catastrophic behavior, namely, there exists a certain critical point, across which an infinitesimal enhancement of the magnetic parameters causes a finite jump of the geometrical parameters for the rope. The amplitude of the jump depends on the extent to which the ambient magnetic field in open, and approaches to zero when the ambient magnetic field becomes completely closed. The implication of such a catastrophe in solar active phenomena is briefly discussed.     

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…     

Transient Coronal Sigmoids and Rotating Erupting Flux Ropes

To determine the relationship between transient coronal (soft X-ray or EUV) sigmoids and erupting flux ropes, we analyse four events in which a transient sigmoid could be associated with a filament whose apex rotates upon eruption and two further events in which the two phenomena were spatially but not temporally coincident. We find the helicity sign of the erupting field and the direction of filament rotation to be consistent with the conversion of twist into writhe under the ideal MHD constraint of helicity conservation, thus supporting our assumption of flux rope topology for the rising filament. For positive (negative) helicity the filament apex rotates clockwise (counterclockwise), consistent with the flux rope taking on a reverse (forward) S shape, which is opposite to that observed for the sigmoid. This result is incompatible with two models for sigmoid formation: one identifying sigmoids with upward arching kink-unstable flux ropes and one identifying sigmoids with a current layer between two oppositely sheared arcades. We find instead that the observations agree well with the model by Titov and Démoulin (Astron. Astrophys. 351, 707, 1999), which identifies transient sigmoids with steepened current layers below rising flux ropes.     

Determination of the Topology Skeleton of Magnetic Fields in a Solar Active Region

Magnetic topology has been a key to the understanding of magnetic energy re-lease mechanism. Based on observed vector magnetograms, we have determined the three-dimensional (3D) topology skeleton of the magnetic fields in the active region NOAA 10720.The skeleton consists of six 3D magnetic nulls and a network of corresponding spines, fans,and null-null lines. For the first time, we have identified a spiral magnetic null in Sun's corona.The magnetic lines of force twisted around the spine of the null, forming a 'magnetic wreath'with excess of free magnetic energy and resembling observed brightening structures at extra-ultraviolet (EUV) wavebands. We found clear evidence of topology eruptions which are re-ferred to as catastrophic changes of topology skeleton associated with a coronal mass ejection(CME) and an explosive X-ray flare. These results shed new lights on the structural complex-ity and its role in explosive magnetic activity. The concept of flux rope has been widely used in modelling explosive magnetic activity, although their observational identity is rather ob-scure or, at least, lacking of necessary details up to date. We suggest that the magnetic wreath associated with the 3D spiral null is likely an important class of the physical entity of flux ropes.     

Magnetic Helicity Injection in Solar Active Regions

We present the evolution of magnetic field and its relationship with mag- netic(current)helicity in solar active regions from a series of photospheric vector magnetograms obtained by Huairou Solar Observing Station,longitudinal magne- tograms by MDI of SOHO and white light images of TRACE.The photospheric current helicity density is a quantity reflecting the local twisted magnetic field and is related to the remaining magnetic helicity in the photosphere,even if the mean current helicity density brings the general chiral property in a layer of solar active regions.As new magnetic flux emerges in active regions,changes of photospheric cur- rent helicity density with the injection of magnetic helicity into the corona from the subatmosphere can be detected,including changes in sign caused by the injection of magnetic helicity of opposite sign.Because the injection rate of magnetic helicity and photospheric current helicity density have different means in the solar atmosphere, the injected magnetic helicity is probably not proportional to the current helicity den- sity remaining in the photosphere.The evidence is that rotation of sunspots does not synchronize exactly with the twist of photospheric transverse magnetic field in some active regions(such as,delta active regions).They represent different aspects of mag- netic chirality.A combined analysis of the observational magnetic helicity parameters actually provides a relative complete picture of magnetic helicity and its transfer in the solar atmosphere.     

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.     

Distribution of Helical Properties of Solar Magnetic Fields

We summarize studies of helical properties of solar magnetic fields such as current helicity and twist of magnetic fields in solar active regions (ARs), that are observational tracers of the alpha-effect in the solar convective zone (SCZ). Information on their spatial distribution is obtained by analysis of systematic mag-netographic observations of active regions taken at Huairou Solar Observing Station of National Astronomical Observatories of Chinese Academy of Sciences. The main property is that the tracers of the alpha-effect are antisymmetric about the solar equator. Identifying longitudinal migration of active regions with their individual rotation rates and taking into account the internal differential rotation law within the SCZ known from helioseismology, we deduce the distribution of the effect over depth. We have found evidence that the alpha-effect changes its value and sign near the bottom of the SCZ, and this is in accord with the theoretical studies and numerical simulations. We discuss     

中尺度磁绳结构的行星际观测

由磁绳结构主导、平均尺度约二、三十个小时的行星际磁云是日冕物质抛射在行星际膨胀、传播的体现。最近，Moldwin等人报道在太阳风中还观测到一些尺度在几十分钟的小尺度磁绳结构，并认为太阳风中的磁绳结构在尺度分布上可能具有双峰特征，在全面检视了WIND卫星(1995年-2000年)和ACE卫星(1998年-2000年)的观测资料后，发现了在行星际太阳风中一些尺度为几个小时的中尺度磁绳结构，利用初步整理的其中28个中尺度磁绳结构事件，认为太阳风中的磁绳结构在尺度分布上可能是连续的，这对行星际太阳风中磁绳结构物理起源的研究可能提出重要的物理限制。     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

Formation of Solar Delta Active Regions： Twist and Writhe of Magnetic Ropes

We analyze the process of formation of delta configuration in some well-known super active regions based on photospheric vector magnetogram observations. It is found that the magnetic field in the initial developing stage of some delta active regions shows a potential-like configuration in the solar atmosphere, the magnetic shear develops mainly near the magnetic neutral line with magnetic islands of opposite polarities, and the large-scale photospheric twisted field forming gradually later. Some results are obtained: (1) The analysis of magnetic writhe of whole active regions cannot be limited in the strong field of sunspots, because the contribution of the fraction of decayed magnetic field is non-negligible. (2) The magnetic model of kink magnetic ropes, supposed to be generated in the subatmosphere, is not consistent with the evolution of large-scale twisted photospheric transverse magnetic field and not entirely consistent with the relationship with magnetic shear in some delta active regions. (3) T     

Catastrophe of Coronal Magnetic Flux Ropes Caused by Photospheric Motions

Using a 2.5-D, time-dependent ideal MHD model in Cartesian coordinates, we carried out numerical simulations to investigate the equilibrium and evolution properties of a magnetic configuration that consists of a coronal magnetic flux rope and a partly open photospheric background field, which is equivalent to that produced by a two-patch magnetic source on the photospheric surface. The axial and annular magnetic fluxes of the flux rope are given and fixed. The global magnetic configuration evolves in response to three types of changes of the background field: decreasing of the distance between the two sources, shrinking of the size of each source, and increasing of the shear in the closed component of the background field. As a result, the geometrical parameters of the flux rope, i.e. the height of the rope axis, the half-width of the rope and the length of the vertical current sheet below the rope, change due to the variation of the background field. It is shown that for a given coronal magnetic flux rope in a partly open background field, the variation of the geometrical parameters of the flux rope displays a catastrophic behavior, namely, there exists a critical point for each case, at which an infinitesimal change of the background field leads to a loss of equilibrium, and thus a jump of the flux rope. The implication of such a catastrophe in solar active phenomena is briefly discussed.     

Numerical Experiments Base on the Catastrophe Model of Solar Eruptionstwo

On the basis of the catastrophe model developed by Isenberg et al., we have used the NIRVANA code to perform the magnetohydrodynamics (MHD) numerical experiments to look into the various behaviors of the coronal magnetic configuration that includes a current-carrying flux rope for modelling the prominence levitation in the corona. These behaviors include the evolution of the equilibrium height of magnetic flux rope with the background magnetic field, the corresponding interior equilibrium of magnetic flux rope, the dynamic properties of magnetic flux rope after the system loses equilibrium, as well as the impact of the reference radius on the equilibrium height of magnetic flux rope. In our calculations, an empirical model of the coronal density distribution given by Sittler & Guhathakurta is used, and the physical dissipation is included. Our experiments show that a deviation between the simulated equilibrium height of magnetic flux rope and the theoretical result of Isenberg et al. exists, but it is not apparent, and the evolutionary features of the two results are similar. If the magnetic flux rope is initially located at the stable branch of the theoretical equilibrium curve, the magnetic flux rope will quickly reach the equilibrium position after several rounds of oscillations as a result of the self-adjustment of the system; when the system is located at the critical point it will quickly lose equilibrium and evolve to the eruptive state; the impact of the variation of reference radius on the equilibrium height of magnetic flux rope is consistent with the prediction of the theory; in the eruptive state, the kinetic properties of magnetic flux rope are consistent with the results given by the Lin-Forbes model and observation, and the fast-mode shock in front of the magnetic flux rope is observed in our experiments; furthermore, because that the dissipation is included in our numerical experiments, the energy conversion from the magnetic energy to other forms of energy is very apparent in the eruptive process.     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Initiation of cmes by magnetic flux emergence

The initiation of solar Coronal Mass Ejections (CMEs) is studied in the framework of numerical magnetohydrodynamics (MHD). The initial CME model includes a magnetic flux rope in spherical, axi-symmetric geometry. The initial configuration consists of a magnetic flux rope embedded in a gravitationally stratified solar atmosphere with a background dipole magnetic field. The flux rope is in equilibrium due to an image current below the photosphere. An emerging flux triggering mechanism is used to make this equilibrium system unstable. When the magnetic flux emerges within the filament below the flux rope, this results in a catastrophic behavior similar to previous models. As a result, the flux rope rises and a current sheet forms below it. It is shown that the magnetic reconnection in the current sheet below the flux rope in combination with the outward curvature forces results in a fast ejection of the flux rope as observed for solar CMEs. We have done a parametric study of the emerging flux rate.     

Projection Effects on Physical Parameters Obtained from Solar Vector Magnetograms

1 INTRODUCTION Magnetic field plays an important role in solar activity. The stressing and subsequent partialrelaxation of magnetic fields in the active regions are generally accepted to be the energy sourceof solar flares. To quantitatively study the extent of stressed magnetic field as distinct from itspotential field, Hagyard et al. (1984) defined a magnetic shear angle膖he azimuth differencebetween the observed transverse magnetic field vector and the computed potential field vectorth…     

Transequatorial Filament Eruption and Its Link to a Coronal Mass Ejection

We revisit the Bastille Day flare/CME Event of 2000 July 14, and demonstrate that this flare/CME event is not related to only one single active region (AR). Activation and eruption of a huge transequatorial filament are seen to precede the simultaneous filament eruption and flare in the source active region, NOAA AR 9077, and the full halo-CME in the high corona. Evidence of reconfiguration of large-scale magnetic structures related to the event is illustrated by SOHO EIT and Yohkoh SXT observations, as well as, the reconstructed 3D magnetic lines of force based on the force-free assumption. We suggest that the AR filament in AR9077 was connected to the transequatorial filament. The large-scale magnetic composition related to the transequatorial filament and its sheared magnetic arcade appears to be an essential part of the CME parent magnetic structure. Estimations show that the filament-arcade system has enough magnetic helicity to account for the helicity carried by the related CMEs. In addition, rather global magnetic connectivity, covering almost all the visible range in longitude and a huge span in latitude on the Sun, is implied by the Nancay Radioheliograph (NRH) observations. The analysis of the Bastille Day event suggests that although the triggering of a global CME might take place in an AR, a much larger scale magnetic composition seems to be the source of the ejected magnetic flux, helicity and plasma. The Bastille Day event is the first described example in the literature, in which a transequatorial filament activity appears to play a key role in a global CME. Many tens of halo-CME are found to be associated with transequatorial filaments and their magnetic environment.     

Numerical simulations of fast and slow coronal mass ejections

Solar coronal mass ejections (CMEs) show a large variety in their kinematic properties. CMEs originating in active regions and accompanied by strong flares are usually faster and accelerated more impulsively than CMEs associated with filament eruptions outside active regions and weak flares. It has been proposed more than two decades ago that there are two separate types of CMEs, fast (impulsive) CMEs and slow (gradual) CMEs. However, this concept may not be valid, since the large data sets acquired in recent years do not show two distinct peaks in the CME velocity distribution and reveal that both fast and slow CMEs can be accompanied by both weak and strong flares. We present numerical simulations which confirm our earlier analytical result that a flux‐rope CME model permits describing fast and slow CMEs in a unified manner. We consider a force‐free coronal magnetic flux rope embedded in the potential field of model bipolar and quadrupolar active regions. The eruption is driven by the torus instability which occurs if the field overlying the flux rope decreases sufficiently rapidly with height. The acceleration profile depends on the steepness of this field decrease, corresponding to fast CMEs for rapid decrease, as is typical of active regions, and to slow CMEs for gentle decrease, as is typical of the quiet Sun. Complex (quadrupolar) active regions lead to the fastest CMEs. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

How to Reach Superequipartition Field Strengths in Solar Magnetic Flux Tubes

A number of independent arguments indicate that the toroidal flux system responsible for the sunspot cycle is stored at the base of the convection zone in the form of flux tubes with field strength close to 10⁵ G. Although the evidence for such strong fields is quite compelling, how such field strength can be reached is still a topic of debate. Flux expulsion by convection should lead to about the equipartition field strength, but the magnetic energy density of a 10⁵-G field is two orders of magnitude larger than the mean kinetic energy density of convective motions. Line stretching by differential rotation (i.e., the “Ω effect” in the classical mean-field dynamo approach) probably plays an important role, but arguments based on energy considerations show that it does not seem feasible that a 10⁵-G field can be produced in this way. An alternative scenario for the intensification of the toroidal flux system in the overshoot layer is related to the explosion of rising, buoyantly unstable magnetic flux tubes, which opens a complementary mechanism for magnetic-field intensification. A parallelism is pointed out with the mechanism of “convective collapse” for the intensification of photospheric magnetic flux tubes up to field strengths well above equipartition; both mechanisms, which are fundamentally thermal processes, are reviewed.