Nonlinear Force-free Coronal Magnetic Field Research Articles

Deducing the reliability of relative helicities from nonlinear force-free coronal models

Aims. We study the relative helicity of active region (AR) NOAA 12673 during a ten-hour time interval centered around a preceding X2.2 flare (SOL2017-09-06T08:57) and also including an eruptive X9.3 flare that occurred three hours later (SOL2017-09-06T11:53). In particular, we aim for a reliable estimate of the normalized self-helicity of the current-carrying magnetic field, the so-called helicity ratio, |HJ|/|H𝒱|, a promising candidate to quantity the eruptive potential of solar ARs. Methods. Using Solar Dynamics Observatory Helioseismic and Magnetic Imager vector magnetic field data as an input, we employ nonlinear force-free (NLFF) coronal magnetic field models using an optimization approach. The corresponding relative helicity, and related quantities, are computed using a finite-volume method. From multiple time series of NLFF models based on different choices of free model parameters, we are able to assess the spread of |HJ|/|H𝒱|, and to estimate its uncertainty. Results. In comparison to earlier works, which identified the non-solenoidal contribution to the total magnetic energy, Ediv/E, as selection criterion regarding the required solenoidal quality of magnetic field models for subsequent relative helicity analysis, we propose to use in addition the non-solenoidal contribution to the free magnetic energy, |Emix|/EJ, s. As a recipe for a reliable estimate of the relative magnetic helicity (and related quantities), we recommend to employ multiple NLFF models based on different combinations of free model parameters, to retain only those that exhibit smallest values of both Ediv/E and |Emix|/EJ, s at a certain time instant, to subsequently compute mean estimates, and to use the spread of the individually contributing values as an indication for the uncertainty.

Erratum: “On the Reliability of Magnetic Energy and Helicity Computations Based on Nonlinear Force-free Coronal Magnetic Field Models” (2019, ApJL, 880, L6)

Erratum: “On the Reliability of Magnetic Energy and Helicity Computations Based on Nonlinear Force-free Coronal Magnetic Field Models” (2019, ApJL, 880, L6)

Magnetic Helicity Budget of Solar Active Regions Prolific of Eruptive and Confined Flares

Abstract We compare the coronal magnetic energy and helicity of two solar active regions (ARs), prolific in major eruptive (AR 11158) and confined (AR 12192) flaring, and analyze the potential of deduced proxies to forecast upcoming flares. Based on nonlinear force-free (NLFF) coronal magnetic field models with a high degree of solenoidality, and applying three different computational methods to investigate the coronal magnetic helicity, we are able to draw conclusions with a high level of confidence. Based on real observations of two solar ARs we checked trends regarding the potential eruptivity of the active-region corona, as suggested earlier in works that were based on numerical simulations, or solar observations. Our results support that the ratio of current-carrying to total helicity, , shows a strong ability to indicate the eruptive potential of a solar AR. However, does not seem to be indicative for the magnitude or type of an upcoming flare (confined or eruptive). Interpreted in the context of earlier observational studies, our findings furthermore support that the total relative helicity normalized to the magnetic flux at the NLFF model’s lower boundary, , represents no indicator for the eruptivity.

On the Reliability of Magnetic Energy and Helicity Computations Based on Nonlinear Force-free Coronal Magnetic Field Models

Abstract We demonstrate the sensitivity of magnetic energy and helicity computations regarding the quality of the underlying coronal magnetic field model. We apply the method of Wiegelmann & Inhester to a series of Solar Dynamics Observatory/Helioseismic and Magnetic Imager vector magnetograms, and discuss nonlinear force-free (NLFF) solutions based on two different sets of the free model parameters. The two time series differ from each other concerning their force-free and solenoidal quality. Both force- and divergence-freeness are required for a consistent NLFF solution. Full satisfaction of the solenoidal property is inherent in the definition of relative magnetic helicity in order to ensure gauge independence. We apply two different magnetic helicity computation methods to both NLFF time series and find that the output is highly dependent on the level to which the NLFF magnetic fields satisfy the divergence-free condition, with the computed magnetic energy being less sensitive than the relative helicity. Proxies for the nonpotentiality and eruptivity derived from both quantities are also shown to depend strongly on the solenoidal property of the NLFF fields. As a reference for future applications, we provide quantitative thresholds for the force- and divergence-freeness, for the assurance of reliable computation of magnetic energy and helicity, and of their related eruptivity proxies.

The Causes of Quasi-homologous CMEs

Abstract In this paper, we identified the magnetic source locations of 142 quasi-homologous (QH) coronal mass ejections (CMEs), of which 121 are from solar cycle (SC) 23 and 21 from SC 24. Among those CMEs, 63% originated from the same source location as their predecessor (defined as S-type), while 37% originated from a different location within the same active region as their predecessor (defined as D-type). Their distinctly different waiting time distributions, peaking around 7.5 and 1.5 hr for S- and D-type CMEs, suggest that they might involve different physical mechanisms with different characteristic timescales. Through detailed analysis based on nonlinear force-free coronal magnetic field modeling of two exemplary cases, we propose that the S-type QH CMES might involve a recurring energy release process from the same source location (by magnetic free energy replenishment), whereas the D-type QH CMEs can happen when a flux tube system is disturbed by a nearby CME.

COMPARISON OF FORCE-FREE CORONAL MAGNETIC FIELD MODELING USING VECTOR FIELDS FROMHINODEANDSOLAR DYNAMICS OBSERVATORY

Photospheric magnetic vector maps from two different instruments are used to model the nonlinear force-free coronal magnetic field above an active region. We use vector maps inferred from polarization measurements of the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (HMI) and the Solar Optical Telescope Spectropolarimeter (SP) aboard Hinode. Besides basing our model calculations on HMI data, we use both, SP data of original resolution and scaled down to the resolution of HMI. This allows us to compare the model results based on data from different instruments and to investigate how a binning of high-resolution data effects the model outcome. The resulting 3D magnetic fields are compared in terms of magnetic energy content and magnetic topology. We find stronger magnetic fields in the SP data, translating into a higher total magnetic energy of the SP models. The net Lorentz forces of the HMI and SP lower boundaries verify their force-free compatibility. We find substantial differences in the absolute estimates of the magnetic field energy but similar relative estimates, e.g., the fraction of excess energy and of the flux shared by distinct areas. The location and extension of neighboring connectivity domains differs and the SP model fields tend to be higher and more vertical. Hence, conclusions about the magnetic connectivity based on force-free field models are to be drawn with caution. We find that the deviations of the model solution when based on the lower-resolution SP data are small compared to the differences of the solutions based on data from different instruments.

Coronal Magnetic Field Structure and Evolution for Flaring AR 11117 and Its Surroundings

In this study, photospheric vector magnetograms obtained with the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS), are used as boundary conditions to model the three-dimensional nonlinear force-free (NLFF) coronal magnetic fields as a sequence of nonlinear force-free equilibria in spherical geometry. We study the coronal magnetic field structure inside active regions and its temporal evolution. We compare the magnetic field configuration obtained from NLFF extrapolation before and after flaring event in active region (AR) 11117 and its surroundings observed on 27 October 2010. We compare the magnetic field topologies and the magnetic energy densities and study the connectivities between AR 11117 and its surroundings. During the investigated time period, we estimate the change in free magnetic energy from before to after the flare to be 1.74x10^{32}erg which represents about 13.5% of nonlinear force-free magnetic energy before the flare. In this study, we find that electric currents from AR 11117 to its surroundings were disrupted after the flare.

Solar Force-free Magnetic Fields

The structure and dynamics of the solar corona is dominated by the magnetic field. In most areas in the corona magnetic forces are so dominant that all non-magnetic forces like plasma pressure gradient and gravity can be neglected in the lowest order. This model assumption is called the force-free field assumption, as the Lorentz force vanishes. This can be obtained by either vanishing electric currents (leading to potential fields) or the currents are co-aligned with the magnetic field lines. First we discuss a mathematically simpler approach that the magnetic field and currents are proportional with one global constant, the so-called linear force-free field approximation. In the generic case, however, the relation between magnetic fields and electric currents is nonlinear and analytic solutions have been only found for special cases, like 1D or 2D configurations. For constructing realistic nonlinear force-free coronal magnetic field models in 3D, sophisticated numerical computations are required and boundary conditions must be obtained from measurements of the magnetic field vector in the solar photosphere. This approach is currently of large interests, as accurate measurements of the photospheric field become available from ground-based (for example SOLIS) and space-born (for example Hinode and SDO) instruments. If we can obtain accurate force-free coronal magnetic field models we can calculate the free magnetic energy in the corona, a quantity which is important for the prediction of flares and coronal mass ejections. 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.

Magnetic Connectivity Between Active Regions 10987, 10988, and 10989 by Means of Nonlinear Force-Free Field Extrapolation

Extrapolation codes for modelling the magnetic field in the corona in cartesian geometry do not take the curvature of the Sun's surface into account and can only be applied to relatively small areas, \textit{e.g.}, a single active region. We apply a method for nonlinear force-free coronal magnetic field modelling of photospheric vector magnetograms in spherical geometry which allows us to study the connectivity between multi-active regions. We use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/Vector Spectromagnetograph(VSM) to model the coronal magnetic field, where we study three neighbouring magnetically connected active regions (ARs: 10987, 10988, 10989) observed on 28, 29, and 30 March 2008, respectively. We compare the magnetic field topologies and the magnetic energy densities and study the connectivities between the active regions(ARs). We have studied the time evolution of magnetic field over the period of three days and found no major changes in topologies as there was no major eruption event. From this study we have concluded that active regions are much more connected magnetically than the electric current.

Nonlinear force-free field extrapolation in spherical geometry: improved boundary data treatment applied to a SOLIS/VSM vector magnetogram

Routine measurements of the solar magnetic field are mainly carried out in the photosphere. Therefore, one has to infer the field strength in the higher layers of the solar atmosphere from the measured photospheric field based on the assumption that the corona is force-free. Meanwhile, those measured data are inconsistent with the above force-free assumption. Therefore, one has to apply some transformations to these data before nonlinear force-free extrapolation codes can be applied. 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. Here we apply a method for nonlinear force-free coronal magnetic field modelling and preprocessing of photospheric vector magnetograms in spherical geometry using the optimization procedure.We solve the nonlinear force-free field equations by minimizing a functional in spherical coordinates over a restricted area of the Sun. We extend the functional by an additional term, which allows to incorporate measurement error and treat regions with lacking observational data. We use vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS) to model the coronal magnetic field. We study two neighbouring magnetically connected active regions observed on May 15 2009. For vector magnetograms with variable measurement precision and randomly scattered data gaps (e.g., SOLIS/VSM) the new code yields field models which satisfy the solenoidal and force-free condition significantly better as it allows deviations between the extrapolated boundary field and observed boundary data within measurement errors. Data gaps are assigned to an infinite error. We extend this new scheme to spherical geometry and apply it for the first time to real data.

How to deal with measurement errors and lacking data in nonlinear force-free coronal magnetic field modelling?

Context. The measured solar photospheric magnetic field vector is extrapolated into the solar corona under the assumption of a force-free plasma. In the generic case this problem is nonlinear. Aims. We aim to improve an algorithm for computing the nonlinear force-free coronal magnetic field. We are in particular interested to incorporate measurement errors and to handle lacking data in the boundary conditions.Methods. We solve the nonlinear force-free field equations by minimizing a functional. Within this work we extend the functional by an additional term, which allows us to incorporate measurement errors and treat regions with lacking observational data. We test the new code with the help of a well known semi-analytic test case. We compare coronal magnetic field extrapolations from ideal boundary conditions and boundary conditions where the transversal magnetic field information is lacking or has a poor signal-to-noise ratio in weak field regions. Results. For ideal boundary conditions the new code gives the same result as the old code. The advantage of the new approach, which includes an error matrix, is visible only for non-ideal boundary conditions. The force-free and solenoidal conditions are fulfilled significantly better and the solutions agrees somewhat better with the exact solution. The new approach also relaxes the boundary and allows a deviation from the boundary data in poor signal-to-noise ratio areas. Conclusions. The incorporation of measurement errors in the updated extrapolation code significantly improves the quality of nonlinear force-free field extrapolation from imperfect boundary conditions.

FREE MAGNETIC ENERGY AND FLARE PRODUCTIVITY OF ACTIVE REGIONS

In this study, the photospheric vector magnetograms, obtained with the Spectro-Polarimeter of the Solar Optical Telescope on board Hinode, are used as the boundary conditions to extrapolate the three-dimensional nonlinear force-free (NLFF) coronal magnetic fields. The observed non-force-free photospheric magnetic fields are preprocessed toward the nearly force-free chromospheric magnetic fields. The performance of the preprocessing procedure is evaluated by comparing with chromospheric magnetic fields obtained by the Vector SpectroMagnetograph instrument located on the Synoptic Optical Long-term Investigations of the Sun Tower. Then, the weighted optimization method is applied to the preprocessed boundary data to extrapolate the NLFF fields with which we are able to estimate the free magnetic energy stored in the active regions. The magnitude scaling correlation between the free magnetic energy and the soft X-ray flare index (FI) of active regions is then studied. The latter quantifies the impending flare production of active regions over the subsequent 1, 2, and 3 day time windows. Based on 75 samples, we find a positive correlation between the free energy and the FI. We also study the temporal variation of free magnetic energy for three active regions, of which two are flare-active and one is flare-quiet during the observation over a period of several days. While the magnitude of free magnetic energy unambiguously differentiates between the flare-active and the flare-quiet regions, the temporal variation of free magnetic energy does not exhibit a clear and consistent pre-flare pattern. This may indicate that the trigger mechanism of flares is as important as the energy storage in active regions.

Nonlinear force-free coronal magnetic field modelling and preprocessing of vector magnetograms in spherical geometry

Context: Knowledge about the coronal magnetic field is important to the understanding the structure of the solar corona. We compute the field in the higher layers of the solar atmosphere from the measured photospheric field under the assumption that the corona is force-free. Aims: Here we develop a method for nonlinear force-free coronal magnetic field medelling and preprocessing of photospheric vector magnetograms in spherical geometry using the optimization procedure. Methods: We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates over a restricted area of the Sun. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free equations in the solar corona. We develop a preprocessing procedure in spherical geometry to drive the observed non-force-free data towards suitable boundary conditions for a force-free extrapolation. Results: We test the code with the help of a semi-analytic solution and assess the quality of our reconstruction qualitatively by magnetic field line plots and quantitatively with a number of comparison metrics for different boundary conditions. The reconstructed fields from the lower boundary data with the weighting function are in good agreement with the original reference fields. We added artificial noise to the boundary conditions and tested the code with and without preprocessing. The preprocessing recovered all main structures of the magnetogram and removed small-scale noise. The main test was to extrapolate from the noisy photospheric vector magnetogram with and without preprocessing. The preprocessing was found to significantly improve the agreement between the extrapolated and the exact field.

Optimization approach for the computation of magnetohydrostatic coronal equilibria in spherical geometry

Context: This paper presents a method which can be used to calculate models of the global solar corona from observational data. Aims: We present an optimization method for computing nonlinear magnetohydrostatic equilibria in spherical geometry with the aim to obtain self-consistent solutions for the coronal magnetic field, the coronal plasma density and plasma pressure using observational data as input. Methods: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma solves the non-force-free magnetohydrostatic equilibria using an optimization method. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements in Cartesian and spherical geometry, and magnetostatic-equilibria in Cartesian geometry. We test our code with the help of a known analytic 3D equilibrium solution of the magnetohydrostatic equations. The detailed comparison between the numerical calculations and the exact equilibrium solutions is made by using magnetic field line plots, plots of density and pressure and some of the usual quantitative numerical comparison measures. Results: We find that the method reconstructs the equilibrium accurately, with residual forces of the order of the discretisation error of the analytic solution. The correlation with the reference solution is better than 99.9% and the magnetic energy is computed accurately with an error of <0.1%. Conclusions: We applied the method so far to an analytic test case. We are planning to use this method with real observational data as input as soon as possible.

Computing Nonlinear Force-Free Coronal Magnetic Fields in Spherical Geometry

We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free equations in the solar corona. The method is based on an optimization principle and the heritage of the newly developed code is a corresponding method in Cartesian geometry. We test the newly developed code with the help of a semi-analytic solution and rate the quality of our reconstruction qualitatively by magnetic field line plots and quantitatively with a number of comparison metrics. We find that we can reconstruct the original test field with high accuracy. The method is fast if the computation is limited to low co-latitudes (say 30°≤θ≤150°), but it becomes significantly slower if the polar regions are included.

An optimization principle for the computation of MHD equilibria in the solar corona

AIMS: We develop an optimization principle for computing stationary MHD equilibria. METHODS: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma uses non-force-free MHD equilibria. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements. The program uses photospheric vector magnetograms and coronal EUV images as input. We tested our reconstruction code with the help of a semi-analytic MHD-equilibrium. The quality of the reconstruction was judged by comparing the exact and reconstructed solution qualitatively by magnetic field-line plots and EUV-images and quantitatively by several different numerical criteria. RESULTS: Our code is able to reconstruct the semi-analytic test equilibrium with high accuracy. The stationary MHD optimization code developed here has about the same accuracy as its predecessor, a non-linear force-free optimization code. The computing time for MHD-equilibria is, however, longer than for force-free magnetic fields. We also extended a well-known class of nonlinear force-free equilibria to the non-force-free regime for purposes of testing the code. CONCLUSIONS: We demonstrate that the code works in principle using tests with analytical equilibria, but it still needs to be applied to real data.

Extrapolation of a nonlinear force-free field containing a highly twisted magnetic loop

The stress-and-relax method for the extrapolation of nonlinear force-free coronal magnetic fields from photospheric vector magnetograms is formulated and implemented in a manner analogous to the evolutionary extrapolation method. The technique is applied to a numerically constructed force-free equilibrium that has a simple bipolar structure of the normal field component in the bottom (magnetogram) plane but contains a highly twisted loop and a shear (current) layer, with a smooth but strong variation of the force-free parameter α in the magnetogram. A standard linear force-free extrapolation of this magnetogram, using the so-called αbest value, is found to fail in reproducing the twisted loop (or flux rope) and the shear layer; it yields a loop pair instead and the shear is not concentrated in a layer. With the nonlinear extrapolation technique, the given equilibrium is readily reconstructed to a high degree of accuracy if the magnetogram is sufficiently resolved. A parametric study quantifies the requirements on the resolution for a successful nonlinear extrapolation. Permitting magnetic reconnection by a controlled use of resistivity improved the extrapolation at a resolution comparable to the smallest structures in the magnetogram.

Reconstruction of the Three-Dimensional Coronal Magnetic Field

Studies of solar flares indicate that the mechanism of flares is magnetic in character and that the coronal magnetic field is a key to understanding solar high-energy phenomena. In our ongoing research we are conducting a systematic study of a large database of observations which includes both coronal structure (from the Soft X-ray Telescope on the Yohkoh spacecraft) and photospheric vector magnetic fields (from the Haleakala Stokes Polarimeter at Mees Solar Observatory). We compare the three-dimensional nonlinear force-free coronal magnetic field, computed from photospheric boundary data, to images of coronal structure. In this paper we outline our techniques and present results for active region AR 7220/7222. We show that the computed force-free coronal magnetic field agrees well with Yohkoh X-ray coronal loops, and we discuss the properties of the coronal magnetic field and the soft X-ray loops.

On the reconstruction of the nonlinear force-free coronal magnetic field from boundary data

Using a simple model in which the corona is represented by the half-space domain Ω = {z > 0} and the photosphere by the boundary plane ∂Ω = {z = 0}, we discuss some important aspects of the general problem of the reconstruction of the magnetic field B in a small isolated coronal region from the values of the vector B¦ ∂Ω measured by a magnetograph over its whole basis. Assuming B to be force-free in Ω: (i) we derive a series of relations which must be necessarily satisfied by the boundary field B¦ ∂Ω , and then by the magnetograph data if the force-free assumption is actually correct; (ii) we show how to extract directly from the measured B¦ ∂Ω some useful informations about the energy of B in Ω and the topological structure of its field lines; (iii) we present a critical discussion of the two methods which have been proposed so far for computing effectively B in Ω from B¦ ∂Ω .