Nonlinear Force-free Field Extrapolation Research Articles

STUDY OF TWO SUCCESSIVE THREE-RIBBON SOLAR FLARES ON 2012 JULY 6

This Letter reports two rarely observed three-ribbon flares (M1.9 and C9.2) on 2012 July 6 in NOAA AR 11515, which we found with Halpha observations of 0.1" resolution from the New Solar Telescope and CaII H images from Hinode. The flaring site is characterized with an intriguing "fish-bone-like" morphology evidenced by both Halpha images and a nonlinear force-free field (NLFFF) extrapolation, where two semi-parallel rows of low-lying, sheared loops connect an elongated, parasitic negative field with the sandwiching positive fields. The NLFFF model also shows that the two rows of loops are asymmetric in height and have opposite twists, and are enveloped by large-scale field lines including open fields. The two flares occurred in succession in half an hour and are located at the two ends of the flaring region. The three ribbons of each flare run parallel to the PIL, with the outer two lying in the positive field and the central one in the negative field. Both flares show surge-like flows in Halpha apparently toward the remote region, while the C9.2 flare is also accompanied by EUV jets possibly along the open field lines. Interestingly, the 12-25 keV hard X-ray sources of the C9.2 flare first line up with the central ribbon then shift to concentrate on the top of the higher branch of loops. These results are discussed in favor of reconnection along the coronal null-line producing the three flare ribbons and the associated ejections.

NONLINEAR FORCE-FREE EXTRAPOLATION OF THE CORONAL MAGNETIC FIELD BASED ON THE MAGNETOHYDRODYNAMIC RELAXATION METHOD

We develop a nonlinear force-free field (NLFFF) extrapolation code based on the magnetohydrodynamic (MHD) relaxation method. We extend the classical MHD relaxation method in two important ways. First, we introduce an algorithm initially proposed by cite{2002JCoPh.175..645D} to effectively clean the numerical errors associated with $nabla cdot vec{B}$. Second, the multi-grid type method is implemented in our NLFFF to perform direct analysis of the high-resolution magnetogram data. As a result of these two implementations, we successfully extrapolated the high resolution force-free field introduced by cite{1990ApJ...352..343L} with better accuracy in a drastically shorter time. We also applied our extrapolation method to the MHD solution obtained from the flux-emergence simulation by cite{2012ApJ...748...53M}. We found that NLFFF extrapolation may be less effective for reproducing areas higher than a half-domain, where some magnetic loops are found in a state of continuous upward expansion. However, an inverse S shaped structure consisting of the sheared and twisted loops formed in the lower region can be captured well through our NLFFF extrapolation method. We further discuss how well these sheared and twisted fields are reconstructed by estimating the magnetic topology and twist quantitatively.

EVIDENCE FOR SOLAR TETHER-CUTTING MAGNETIC RECONNECTION FROM CORONAL FIELD EXTRAPOLATIONS

Magnetic reconnection is one of the primary mechanisms for triggering solar eruptive events, but direct observation of its rapid process has been of challenge. In this Letter we present, using a nonlinear force-free field (NLFFF) extrapolation technique, a visualization of field line connectivity changes resulting from tether-cutting reconnection over about 30 minutes during the 2011 February 13 M6.6 flare in NOAA AR 11158. Evidence for the tether-cutting reconnection was first collected through multiwavelength observations and then by the analysis of the field lines traced from positions of four conspicuous flare 1700 A footpoints observed at the event onset. Right before the flare, the four footpoints are located very close to the regions of local maxima of magnetic twist index. Especially, the field lines from the inner two footpoints form two strongly twisted flux bundles (up to ~1.2 turns), which shear past each other and reach out close to the outer two footpoints, respectively. Immediately after the flare, the twist index of regions around the footpoints greatly diminish and the above field lines become low lying and less twisted (~0.6 turns), overarched by loops linking the later formed two flare ribbons. About 10% of the flux (~3x10^19 Mx) from the inner footpoints has undergone a footpoint exchange. This portion of flux originates from the edge regions of the inner footpoints that are brightened first. These rapid changes of magnetic field connectivity inferred from the NLFFF extrapolation are consistent with the tether-cutting magnetic reconnection model.

A Comparison Between Nonlinear Force-Free Field and Potential Field Models Using Full-Disk SDO/HMI Magnetogram

Measurements of magnetic fields and electric currents in the pre-eruptive corona are crucial to the study of solar eruptive phenomena, like flares and coronal mass ejections (CMEs). However, spectro-polarimetric measurements of certain photospheric lines permit a determination of the vector magnetic field only at the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field above multiple active regions with the help of a potential field and a nonlinear force-free field (NLFFF) extrapolation code over the full solar disk using Helioseismic and Magnetic Imager (SDO/HMI) data as boundary conditions. We compare projections of the resulting magnetic field lines with full-disk coronal images from the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study has found that the NLFFF model reconstructs the magnetic configuration closer to observation than the potential field model for full-disk magnetic field extrapolation. We conclude that many of the trans-equatorial loops connecting the two solar hemispheres are current-free.

First observational application of a connectivity-based helicity flux density

Measuring the magnetic helicity distribution in the solar corona can help in understanding the trigger of solar eruptive events because magnetic helicity is believed to play a key role in solar activity due to its conservation property. A new method for computing the photospheric distribution of the helicity flux was recently developed. This method takes into account the magnetic field connectivity whereas previous methods were based on photospheric signatures only. This novel method maps the true injection of magnetic helicity in active regions. We applied this method for the first time to an observed active region, NOAA 11158, which was the source of intense flaring activity. We used high-resolution vector magnetograms from the SDO/HMI instrument to compute the photospheric flux transport velocities and to perform a nonlinear force-free magnetic field extrapolation. We determined and compared the magnetic helicity flux distribution using a purely photospheric as well as a connectivity-based method. While the new connectivity-based method confirms the mixed pattern of the helicity flux in NOAA 11158, it also reveals a different, and more correct, distribution of the helicity injection. This distribution can be important for explaining the likelihood of an eruption from the active region. The connectivity-based approach is a robust method for computing the magnetic helicity flux, which can be used to study the link between magnetic helicity and eruptivity of observed active regions.

MAGNETOHYDRODYNAMIC SIMULATION OF A SIGMOID ERUPTION OF ACTIVE REGION 11283

Current magnetohydrodynamic (MHD) simulations of the initiation of solar eruptions are still commonly carried out with idealized magnetic field models, whereas the realistic coronal field prior to eruptions can possibly be reconstructed from the observable photospheric field. Using a nonlinear force-free field extrapolation prior to a sigmoid eruption in AR 11283 as the initial condition in an MHD model, we successfully simulate the realistic initiation process of the eruption event, as is confirmed by a remarkable resemblance to the SDO/AIA observations. Analysis of the pre-eruption field reveals that the envelope flux of the sigmoidal core contains a coronal null and furthermore the flux rope is prone to a torus instability. Observations suggest that reconnection at the null cuts overlying tethers and likely triggers the torus instability of the flux rope, which results in the eruption. This kind of simulation demonstrates the capability of modeling the realistic solar eruptions to provide the initiation process.

Recurrent coronal jets induced by repetitively accumulated electric currents

Three extreme-ultraviolet (EUV) jets recurred in about one hour on 2010 September 17 in the following magnetic polarity of active region 11106. The EUV jets were observed by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). The Helioseismic and Magnetic Imager (HMI) on board SDO measured the vector magnetic field, from which we derive the magnetic flux evolution, the photospheric velocity field, and the vertical electric current evolution. The magnetic configuration before the jets is derived by the nonlinear force-free field (NLFFF) extrapolation. We derive that the jets are above a pair of parasitic magnetic bipoles which are continuously driven by photospheric diverging flows. The interaction drove the build up of electric currents that we indeed observed as elongated patterns at the photospheric level. For the first time, the high temporal cadence of HMI allows to follow the evolution of such small currents. In the jet region, we found that the integrated absolute current peaks repetitively in phase with the 171 A flux evolution. The current build up and its decay are both fast, about 10 minutes each, and the current maximum precedes the 171 A by also about 10 minutes. Then, HMI temporal cadence is marginally fast enough to detect such changes. The photospheric current pattern of the jets is found associated to the quasi-separatrix layers deduced from the magnetic extrapolation. From previous theoretical results, the observed diverging flows are expected to build continuously such currents. We conclude that magnetic reconnection occurs periodically, in the current layer created between the emerging bipoles and the large scale active region field. It induced the observed recurrent coronal jets and the decrease of the vertical electric current magnitude.

The state of nonlinear force-free magnetic field extrapolation

Magnetic field extrapolation is the construction of a model solution for the coronal magnetic field in active regions from magnetic boundary data originating close to the Sun's surface. The nonlinear force-free model (in which the electric current density is parallel to the magnetic field) is often adopted to describe the coronal field. The solution of the nonlinear force-free equations is a challenging computational task, and the application of codes to available boundary data has revealed a number of significant problems with nonlinear force-free extrapolation. This paper summarises the present status of coronal field extrapolation, and describes some recent developments.

Modeling coronal magnetic field using spherical geometry: cases with several active regions

The magnetic fields in the solar atmosphere structure the plasma, store free magnetic energy and produce a wide variety of active solar phenomena, like flare and coronal mass ejections(CMEs). The distribution and strength of magnetic fields are routinely measured in the solar surface(photosphere). Therefore, there is considerable interest in accurately modeling the 3D structure of the coronal magnetic field using photospheric vector magnetograms. Knowledge of the 3D structure of magnetic field lines also help us to interpret other coronal observations, e.g., EUV images of the radiating coronal plasma. Nonlinear force-free field (NLFFF) models are thought to be viable tools for those task. Usually those models use Cartesian geometry. However, the spherical nature of the solar surface cannot be neglected when the field of view is large. In this work, we model the coronal magnetic field above multiple active regions using NLFFF extrapolation code using vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/ Vector Spectromagnetograph (VSM) as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with their respective coronal EUV-images from the Atmospheric Imaging Assembly (SDO/AIA) observed on October 11, 2011 and November 13, 2012. This study has found that the NLFFF model in spherical geometry reconstructs the magnetic configurations for several active regions which agrees with observations. During October 11, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.

Topological study of active region 11158

AbstractWith the cylindrical equal area (CEA) projection data from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), we reconstructed the three-dimensional (3D) magnetic fields in the corona, using a non-linear force-free field (NLFFF) extrapolation method every 12 minutes during five days, to calculate the squashing degree factor Q in the volume. The results show that this AR has an hyperbolic flux tube (HFT) configuration, a typical topology of quadrupole, which is stable even during the two large flares (M6.6 and X2.2 class flares).

MAGNETIC ENERGY PARTITION BETWEEN THE CORONAL MASS EJECTION AND FLARE FROM AR 11283

On 2011 September 6, an X-class flare and a halo coronal mass ejection (CME) were observed from Earth erupting from the same active region AR 11283. The magnetic energy partition between them has been investigated. SDO/HMI vector magnetograms were used to obtain the coronal magnetic field using the nonlinear force-free field (NLFFF) extrapolation method. The free magnetic energies before and after the flare were calculated to estimate the released energy available to power the flare and the CME. For the flare energetics, thermal and nonthermal energies were derived using the RHESSI and GOES data. To obtain the radiative output, SDO/EVE data in the 0.1-37 nm waveband were utilized. We have reconstructed the three-dimensional (3D) periphery of the CME from the coronagraph images observed by STEREO-A, B, and SOHO. The mass calculations were then based on a more precise Thomson-scattering geometry. The subsequent estimate of the kinetic and potential energies of the CME took advantage of the more accurate mass, and the height and speed in a 3D frame. The released free magnetic energy resulting from the NLFFF model is about 6.4x10(31) erg, which has a possible upper limit of 1.8x10(32) erg. The thermal and nonthermal energies are lower than the radiative output of 2.2x10(31) erg from SDO/EVE for this event. The total radiation covering the whole solar spectrum is probably a few times larger. The sum of the kinetic and potential energy of the CME could go up to 6.5x10(31) erg. Therefore, the free energy is able to power the flare and the CME in AR 11283. Within the uncertainty, the flare and the CME may consume a similar amount of free energy.

Evidence for partial Taylor relaxation from changes in magnetic geometry and energy during a solar flare

Solar flares are powered by energy stored in the coronal magnetic field, a portion of which is released when the field reconfigures into a lower energy state. Investigation of sunspot magnetic field topology during flare activity is useful to improve our understanding of flaring processes. Here we investigate the deviation of the non-linear field configuration from that of the linear and potential configurations, and study the free energy available leading up to and after a flare. The evolution of the magnetic field in NOAA region 10953 was examined using data from Hinode/SOT-SP, over a period of 12 hours leading up to and after a GOES B1.0 flare. Previous work on this region found pre- and post-flare changes in photospheric vector magnetic field parameters of flux elements outside the primary sunspot. 3D geometry was thus investigated using potential, linear force-free, and non-linear force-free field extrapolations in order to fully understand the evolution of the field lines. Traced field line geometrical and footpoint orientation differences show that the field does not completely relax to a fully potential or linear force-free state after the flare. Magnetic and free magnetic energies increase significantly ~ 6.5-2.5 hours before the flare by ~ 10^31 erg. After the flare, the non-linear force-free magnetic energy and free magnetic energies decrease but do not return to pre-flare 'quiet' values. The post-flare non-linear force-free field configuration is closer (but not equal) to that of the linear force-free field configuration than a potential one. However, the small degree of similarity suggests that partial Taylor relaxation has occurred over a time scale of ~ 3-4 hours.

Full-disk nonlinear force-free field extrapolation of SDO/HMI and SOLIS/VSM magnetograms

Extrapolation codes in Cartesian geometry for modelling the magnetic field in the corona do not take the curvature of the Sun's surface into account and can only be applied to relatively small areas, e.g., a single active region. We compare the analysis of the photospheric magnetic field and subsequent force-free modeling based on full-disk vector maps from Helioseismic and Magnetic Imager (HMI) on board solar dynamics observatory (SDO) and Vector Spectromagnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun (SOLIS). We use Helioseismic and Magnetic Imager and Vector Spectromagnetograph photospheric magnetic field measurements to model the force-free coronal field above multiple solar active regions, assuming magnetic forces to dominate. We solve the nonlinear force-free field equations by minimizing a functional in spherical coordinates over a full disk excluding the poles. After searching for the optimum modeling parameters for the particular data sets, we compare the resulting nonlinear force-free model fields. We compare quantities like the total magnetic energy content and free magnetic energy, the longitudinal distribution of the magnetic pressure and surface electric current density using our spherical geometry extrapolation code. The magnetic field lines obtained from nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager and Vector Spectromagnetograph data have good agreement. However, the nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager data have more contents of total magnetic energy, free magnetic energy, the longitudinal distribution of the magnetic pressure and surface electric current density compared to the one from Vector Spectromagnetograph data.

MODELING MAGNETIC FIELD STRUCTURE OF A SOLAR ACTIVE REGION CORONA USING NONLINEAR FORCE-FREE FIELDS IN SPHERICAL GEOMETRY

We test a nonlinear force-free field (NLFFF) optimization code in spherical geometry using an analytical solution from Low and Lou. Several tests are run, ranging from idealized cases where exact vector field data are provided on all boundaries, to cases where noisy vector data are provided on only the lower boundary (approximating the solar problem). Analytical tests also show that the NLFFF code in the spherical geometry performs better than that in the Cartesian one when the field of view of the bottom boundary is large, say, $20^\circ \times 20^\circ$. Additionally, We apply the NLFFF model to an active region observed by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) both before and after an M8.7 flare. For each observation time, we initialize the models using potential field source surface (PFSS) extrapolations based on either a synoptic chart or a flux-dispersal model, and compare the resulting NLFFF models. The results show that NLFFF extrapolations using the flux-dispersal model as the boundary condition have slightly lower, therefore better, force-free and divergence-free metrics, and contain larger free magnetic energy. By comparing the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the Atmospheric Imaging Assembly (AIA) on board SDO, we find that the NLFFF performs better than the PFSS not only for the core field of the flare productive region, but also for large EUV loops higher than 50 Mm.

A NON-RADIAL ERUPTION IN A QUADRUPOLAR MAGNETIC CONFIGURATION WITH A CORONAL NULL

We report one of several homologous non-radial eruptions from NOAA active region (AR) 11158 that are strongly modulated by the local magnetic field as observed with the Solar Dynamic Observatory (SDO). A small bipole emerged in the sunspot complex and subsequently created a quadrupolar flux system. Non-linear force-free field (NLFFF) extrapolation from vector magnetograms reveals its energetic nature: the fast-shearing bipole accumulated ~2e31 erg free energy (10% of AR total) over just one day despite its relatively small magnetic flux (5% of AR total). During the eruption, the ejected plasma followed a highly inclined trajectory, over 60 degrees with respect to the radial direction, forming a jet-like, inverted-Y shaped structure in its wake. Field extrapolation suggests complicated magnetic connectivity with a coronal null point, which is favorable of reconnection between different flux components in the quadrupolar system. Indeed, multiple pairs of flare ribbons brightened simultaneously, and coronal reconnection signatures appeared near the inferred null. Part of the magnetic setting resembles that of a blowout-type jet; the observed inverted-Y structure likely outlines the open field lines along the separatrix surface. Owing to the asymmetrical photospheric flux distribution, the confining magnetic pressure decreases much faster horizontally than upward. This special field geometry likely guided the non-radial eruption during its initial stage.

GUIDING NONLINEAR FORCE-FREE MODELING USING CORONAL OBSERVATIONS: FIRST RESULTS USING A QUASI-GRAD-RUBIN SCHEME

Presently, many models of the coronal magnetic field rely on photospheric vector magnetograms but these data have been shown to be problematic as the sole boundary information for nonlinear force-free field (NLFFF) extrapolations. Magnetic fields in the corona manifest themselves in high-energy images (X-rays and EUV) in the shapes of coronal loops, providing an additional constraint that at present is not used due to the mathematical complications of incorporating such input into numerical models. Projection effects and the limited number of usable loops further complicate the use of coronal information. We develop and test an algorithm to use images showing coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with field lines of constant-\als force-free fields to approximate the three-dimensional distribution of currents in the corona along a sparse set of trajectories. We then apply a Grad-Rubin-like iterative technique to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying method presented in \citet{Wheatland2007}. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques \citep{Schrijver2006, Schrijver2008} and compare the results with those obtained by presently available methods that rely only on the photospheric data. We conclude that we have developed a functioning method of modeling the coronal magnetic field by combining the line-of-sight component of photospheric magnetic field with information from coronal images. Vector magnetograms over the full or partial photospheric boundary of the numerical domain could optionally be used.

Understanding the white-light flare on 2012 March 9: evidence of a two-step magnetic reconnection

We attempt to understand the white-light flare (WLF) that was observed on 2012 March 9 with a newly constructed multi-wavelength solar telescope called the Optical and Near-infrared Solar Eruption Tracer (ONSET). We analyzed WLF observations in radio, H-alpha, white-light, ultraviolet, and X-ray bands. We also studied the magnetic configuration of the flare via the nonlinear force-free field (NLFFF) extrapolation and the vector magnetic field observed by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). Continuum emission enhancement clearly appeared at the 3600 angstrom and 4250 angstrom bands, with peak contrasts of 25% and 12%, respectively. The continuum emission enhancement closely coincided with the impulsive increase in the hard X-ray emission and a microwave type III burst at 03:40 UT. We find that the WLF appeared at one end of either the sheared or twisted field lines or both. There was also a long-lasting phase in the H-alpha and soft X-ray bands after the white-light emission peak. In particular, a second, yet stronger, peak appeared at 03:56 UT in the microwave band. This event shows clear evidence that the white-light emission was caused by energetic particles bombarding the lower solar atmosphere. A two-step magnetic reconnection scenario is proposed to explain the entire process of flare evolution, i.e., the first-step magnetic reconnection between the field lines that are highly sheared or twisted or both, and the second-step one in the current sheet, which is stretched by the erupting flux rope. The WLF is supposed to be triggered in the first-step magnetic reconnection at a relatively low altitude.

Nonlinear force-free magnetic field extrapolation in AR 11158 by GPU-DBIE

AbstractA GPU-based acceleration for the direct boundary integral equation method (GPU-DBIE) to extrapolate solar coronal magnetic fields is developed, which is about 1000 times faster than the original DBIE. The 3-d coronal magnetic field is reconstructed for NOAA 11158 on 14-Feb-2011 with the SDO/HMI vector magnetogram as bottom boundary condition. The extrapolated results agree well with the projected SDO/AIA, EUV loops and the STEREO EUV sideviews, which verifies the correctness of our GPU-DBIE method. It is also found that the group of bright EUV loops along magnetic neutral lines agree well with current lines, which may have played an important role in the flare process of the active region.

Tests and applications of nonlinear force-free field extrapolations in spherical geometry

AbstractWe test a nonlinear force-free field (NLFFF) optimization code in spherical geometry with an analytical solution from Low and Lou. The potential field source surface (PFSS) model is served as the initial and boundary conditions where observed data are not available. The analytical solution can be well recovered if the boundary and initial conditions are properly handled. Next, we discuss the preprocessing procedure for the noisy bottom boundary data, and find that preprocessing is necessary for NLFFF extrapolations when we use the observed photospheric magnetic field as bottom boundaries. Finally, we apply the NLFFF model to a solar area where four active regions interacting with each other. An M8.7 flare occurred in one active region. NLFFF modeling in spherical geometry simultaneously constructs the small and large scale magnetic field configurations better than the PFSS model does.

A NEW CODE FOR NONLINEAR FORCE-FREE FIELD EXTRAPOLATION OF THE GLOBAL CORONA

Reliable measurements of the solar magnetic field are still restricted to the photosphere, and our present knowledge of the three-dimensional coronal magnetic field is largely based on extrapolation from photospheric magnetogram using physical models, e.g., the nonlinear force-free field (NLFFF) model as usually adopted. Most of the currently available NLFFF codes have been developed with computational volume like Cartesian box or spherical wedge while a global full-sphere extrapolation is still under developing. A high-performance global extrapolation code is in particular urgently needed considering that Solar Dynamics Observatory (SDO) can provide full-disk magnetogram with resolution up to $4096\times 4096$. In this work, we present a new parallelized code for global NLFFF extrapolation with the photosphere magnetogram as input. The method is based on magnetohydrodynamics relaxation approach, the CESE-MHD numerical scheme and a Yin-Yang spherical grid that is used to overcome the polar problems of the standard spherical grid. The code is validated by two full-sphere force-free solutions from Low & Lou's semi-analytic force-free field model. The code shows high accuracy and fast convergence, and can be ready for future practical application if combined with an adaptive mesh refinement technique.