Refine
Has Fulltext
- no (14) (remove)
Document Type
- Article (14)
Language
- English (14) (remove)
Is part of the Bibliography
- yes (14)
Keywords
- Sun: corona (14) (remove)
Institute
The onset of a solar eruption is formulated here as either a magnetic catastrophe or as an instability. Both start with the same equation of force balance governing the underlying equilibria. Using a toroidal flux rope in an external bipolar or quadrupolar field as a model for the current-carrying flux, we demonstrate the occurrence of a fold catastrophe by loss of equilibrium for several representative evolutionary sequences in the stable domain of parameter space. We verify that this catastrophe and the torus instability occur at the same point; they are thus equivalent descriptions for the onset condition of solar eruptions.
In this paper, we address the formation of a magnetic flux rope (MFR) that erupted on 2012 July 12 and caused a strong geomagnetic storm event on July 15. Through analyzing the long-term evolution of the associated active region observed by the Atmospheric Imaging Assembly and the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, it is found that the twisted field of an MFR, indicated by a continuous S-shaped sigmoid, is built up from two groups of sheared arcades near the main polarity inversion line a half day before the eruption. The temperature within the twisted field and sheared arcades is higher than that of the ambient volume, suggesting that magnetic reconnection most likely works there. The driver behind the reconnection is attributed to shearing and converging motions at magnetic footpoints with velocities in the range of 0.1-0.6 km s(-1). The rotation of the preceding sunspot also contributes to the MFR buildup. Extrapolated three-dimensional non-linear force-free field structures further reveal the locations of the reconnection to be in a bald-patch region and in a hyperbolic flux tube. About 2 hr before the eruption, indications of a second MFR in the form of an S-shaped hot channel are seen. It lies above the original MFR that continuously exists and includes a filament. The whole structure thus makes up a stable double-decker MFR system for hours prior to the eruption. Eventually, after entering the domain of instability, the high-lying MFR impulsively erupts to generate a fast coronal mass ejection and X-class flare; while the low-lying MFR remains behind and continuously maintains the sigmoidicity of the active region.
Context. For magnetically driven events, the magnetic energy of the system is the prime energy reservoir that fuels the dynamical evolution. In the solar context, the free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. A trustworthy estimation of the magnetic energy is therefore needed in three-dimensional (3D) models of the solar atmosphere, e. g., in coronal fields reconstructions or numerical simulations.
Aims. The expression of the energy of a system as the sum of its potential energy and its free energy (Thomson's theorem) is strictly valid when the magnetic field is exactly solenoidal. For numerical realizations on a discrete grid, this property may be only approximately fulfilled. We show that the imperfect solenoidality induces terms in the energy that can lead to misinterpreting the amount of free energy present in a magnetic configuration.
Methods. We consider a decomposition of the energy in solenoidal and nonsolenoidal parts which allows the unambiguous estimation of the nonsolenoidal contribution to the energy. We apply this decomposition to six typical cases broadly used in solar physics. We quantify to what extent the Thomson theorem is not satisfied when approximately solenoidal fields are used.
Results. The quantified errors on energy vary from negligible to significant errors, depending on the extent of the nonsolenoidal component of the field. We identify the main source of errors and analyze the implications of adding a variable amount of divergence to various solenoidal fields. Finally, we present pathological unphysical situations where the estimated free energy would appear to be negative, as found in some previous works, and we identify the source of this error to be the presence of a finite divergence.
Conclusions. We provide a method of quantifying the effect of a finite divergence in numerical fields, together with detailed diagnostics of its sources. We also compare the efficiency of two divergence-cleaning techniques. These results are applicable to a broad range of numerical realizations of magnetic fields.
The structure of the coronal magnetic field prior to eruptive processes and the conditions for the onset of eruption are important issues that can be addressed through studying the magnetohydrodynamic (MHD) stability and evolution of nonlinear force-free field (NLFFF) models. This paper uses data-constrained NLFFF models of a solar active region (AR) that erupted on 2010 April 8 as initial conditions in MHD simulations. These models, constructed with the techniques of flux rope insertion and magnetofrictional relaxation (MFR), include a stable, an approximately marginally stable, and an unstable configuration. The simulations confirm previous related results of MFR runs, particularly that stable flux rope equilibria represent key features of the observed pre-eruption coronal structure very well, and that there is a limiting value of the axial flux in the rope for the existence of stable NLFFF equilibria. The specific limiting value is located within a tighter range, due to the sharper discrimination between stability and instability by the MHD description. The MHD treatment of the eruptive configuration yields a very good agreement with a number of observed features, like the strongly inclined initial rise path and the close temporal association between the coronal mass ejection and the onset of flare reconnection. Minor differences occur in the velocity of flare ribbon expansion and in the further evolution of the inclination; these can be eliminated through refined simulations. We suggest that the slingshot effect of horizontally bent flux in the source region of eruptions can contribute significantly to the inclination of the rise direction. Finally, we demonstrate that the onset criterion, formulated in terms of a threshold value for the axial flux in the rope, corresponds very well to the threshold of the torus instability in the considered AR.