The shapes of solar coronal loops are sensitive to the presence of electrical currents that are the carriers of the nonpotential energy available for impulsive activity. We use this information in a new method for modeling the coronal magnetic field
of AR 11158 as a nonlinear force-free field (NLFFF). The observations used are coronal images around time of major flare activity on 2011/02/15, together with the surface line-of-sight magnetic field measurements. The data are from the Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly (HMI and AIA, respectively) onboard the Solar Dynamics Observatory (SDO). The model fields are constrained to approximate the coronal loop configurations as closely as possible, while also subject to the force-free constraints. The method does not use transverse photospheric magnetic field components as input, and is thereby distinct from methods for modeling NLFFFs based on photospheric vector magnetograms. We validate the method using observations of AR 11158 at a time well before major flaring, and subsequently review the field evolution just prior to and following an X2.2 flare and associated eruption. The models indicate that the energy released during the instability is about $1times10^{32}$ erg, consistent with what is needed to power such a large eruptive flare. Immediately prior to the eruption the model field contains a compact sigmoid bundle of twisted flux that is not present in the post-eruption models, which is consistent with the observations. The core of that model structure is twisted by $approx0.9$ full turns about its axis.
Most 1d hydrodynamic models of plasma confined to magnetic flux tubes assume circular cross-section of these tubes. We use potential field models to show that flux tubes in circumstances relevant to the solar corona do not in general maintain the sam
e cross-sectional shape through their length and therefore the assumption of a circular cross-section is rarely true. We support our hypothesis with mathematical reasoning and numeric experiments. We demonstrate that lifting this assumption in realistic non-circular loops make apparent expansion of magnetic flux tubes consistent with that of observed coronal loops. We propose that in a bundle of ribbon-like loops those that are viewed along the wide direction would stand out against those that are viewed across the wide direction, due to the difference in their column depths. That would impose a bias towards selecting loops that appear not to be expanding seen projected in the plane of sky. An implication of this selection bias is that the preferentially selected non-circular loops would appear to have increased pressure scale height even if they are resolved by current instruments.
In the present work we study evolution of magnetic helicity in the solar corona. We compare the rate of change of a quantity related to the magnetic helicity in the corona to the flux of magnetic helicity through the photosphere and find that the two
rates are similar. This gives observational evidence that helicity flux across the photosphere is indeed what drives helicity changes in solar corona during emergence. For the purposes of estimating coronal helicity we neither assume a strictly linear force-free field, nor attempt to construct a non-linear force-free field. For each coronal loop evident in Extreme Ultraviolet (EUV) we find a best-matching line of a linear force-free field and allow the twist parameter alpha to be different for each line. This method was introduced and its applicability was discussed in Malanushenko et. al. (2009). The object of the study is emerging and rapidly rotating AR 9004 over about 80 hours. As a proxy for coronal helicity we use the quantity <alpha_i*L_i/2> averaged over many reconstructed lines of magnetic field. We argue that it is approximately proportional to flux-normalized helicity H/Phi^2, where H is helicity and Phi is total enclosed magnetic flux of the active region. The time rate of change of such quantity in the corona is found to be about 0.021 rad/hr, which is compatible with the estimates for the same region obtained using other methods Longcope et. al. (2007), who estimated the flux of normalized helicity of about 0.016 rad/hr.
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.
Non-linear force-free fields are the most general case of force-free fields, but the hardest to model as well. There are numerous methods of computing such fields by extrapolating vector magnetograms from the photosphere, but very few attempts have s
o far made quantitative use of coronal morphology. We present a method to make such quantitative use of X-Ray and EUV images of coronal loops. Each individual loop is fit to a field line of a linear force-free field, allowing the estimation of the field lines twist, three-dimensional geometry and the field strength along it. We assess the validity of such a reconstruction since the actual corona is probably not a linear force-free field and that the superposition of linear force-free fields is generally not itself a force-free field. To do so, we perform a series of tests on non-linear force-free fields, described in Low & Lou (1990). For model loops we project field lines onto the photosphere. We compare several results of the method with the original field, in particular the three-dimensional loop shapes, local twist (coronal alpha), distribution of twist in the model photosphere and strength of the magnetic field. We find that, (i) for these trial fields, the method reconstructs twist with mean absolute deviation of at most 15% of the range of photospheric twist, (ii) that heights of the loops are reconstructed with mean absolute deviation of at most 5% of the range of trial heights and (iii) that the magnitude of non-potential contribution to photospheric field is reconstructed with mean absolute deviation of at most 10% of the maximal value.
In this paper we propose that additive self helicity, introduced by Longcope and Malanushenko (2008), plays a role in the kink instability for complex equilibria, similar to twist helicity for thin flux tubes (Hood and Priest (1979), Berger and Field
(1984)). We support this hypothesis by a calculation of additive self helicity of a twisted flux tube from the simulation of Fan and Gibson (2003). As more twist gets introduced, the additive self helicity increases, and the kink instability of the tube coincides with the drop of additive self helicity, after the latter reaches the value of $H_A/Phi^2approx 1.5$ (where $Phi$ is the flux of the tube and $H_A$ is additive self helicity). We compare additive self helicity to twist for a thin sub-portion of the tube to illustrate that $H_A/Phi^2$ is equal to the twist number, studied by Berger and Field (1984), when the thin flux tube approximation is applicable. We suggest, that the quantity $H_A/Phi^2$ could be treated as a generalization of a twist number, when thin flux tube approximation is not applicable. A threshold on a generalized twist number might prove extremely useful studying complex equilibria, just as twist number itself has proven useful studying idealized thin flux tubes. We explicitly describe a numerical method for calculating additive self helicity, which includes an algorithm for identifying a domain occupied by a flux bundle and a method of calculating potential magnetic field confined to this domain. We also describe a numerical method to calculate twist of a thin flux tube, using a frame parallelly transported along the axis of the tube.
