Magnetic reconnection occurs when new flux emerges into the corona and becomes incorporated into the existing coronal field. A new active region (AR) emerging in the vicinity of an existing AR provides a convenient laboratory in which reconnection of this kind can be quantified. We use high time-cadence 171 $\AA$ data from SDO/AIA focused on new/old active region pair 11147/11149, to quantify reconnection. We identify new loops as brightenings within a strip of pixels between the regions. This strategy is premised on the assumption that the energy brightening a loop originates in magnetic reconnection. We catalog 301 loops observed in the 48-hour time period beginning with the emergence of AR 11149. The rate at which these loops appear between the two ARs is used to calculate the reconnection rate between them. We then fit these loops with magnetic field, solving for each loop's field strength, geometry, and twist (via its proxy, coronal $\alpha$). We find the rate of newly-brightened flux overestimates the flux which could be undergoing reconnection. This excess can be explained by our finding that the interconnecting region is not at its lowest energy (constant-$\alpha$) state; the extrapolations exhibit loop-to-loop variation in $\alpha$. This flux overestimate may result from the slow emergence of AR 11149, allowing time for Taylor relaxation internal to the domain of the reconnected flux to bring the $\alpha$ distribution towards a single value which provides another mechanism for brightening loops after they are first created.

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

We test the feasibility of 3D coronal-loop tracing in stereoscopic EUV image pairs, with the ultimate goal of enabling efficient 3D reconstruction of the coronal magnetic field that drives flares and coronal mass ejections (CMEs). We developed an automated code designed to perform triangulation of coronal loops in pairs (or triplets) of EUV images recorded from different perspectives. The automated (or blind) stereoscopy code includes three major tasks: (i) automated pattern recognition of coronal loops in EUV images, (ii) automated pairing of corresponding loop patterns from two different aspect angles, and (iii) stereoscopic triangulation of 3D loop coordinates. We perform tests with simulated stereoscopic EUV images and quantify the accuracy of all three procedures. In addition we test the performance of the blind stereoscopy code as a function of the spacecraft-separation angle and as a function of the spatial resolution. We also test the sensitivity to magnetic non-potentiality. The automated code developed here can be used for analysis of existing {\sl Solar TErrestrial RElationship Observatory (STEREO)} data, but primarily serves for a design study of a future mission with dedicated diagnostics of non-potential magnetic fields. For a pixel size of 0.6\arcsec (corresponding to the {\sl Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA)} spatial resolution of 1.4\arcsec), we find an optimum spacecraft-separation angle of $\alpha_s \approx 5^\circ$.

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a non-linear force-free field, and (3) thermodynamic models from a fully-compressible, 3D MHD simulation of AR corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and XRT data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.

Sunspots are a canonical marker of the Sun's internal magnetic field which flips polarity every ~22-years. The principal variation of sunspots, an ~11-year variation in number, modulates the amount of magnetic field that pierces the solar surface and drives significant variations in our Star's radiative, particulate and eruptive output over that period. This paper presents observations from the Solar and Heliospheric Observatory and Solar Dynamics Observatory indicating that the 11-year sunspot variation is intrinsically tied it to the spatio-temporal overlap of the activity bands belonging to the 22-year magnetic activity cycle. Using a systematic analysis of ubiquitous coronal brightpoints, and the magnetic scale on which they appear to form, we show that the landmarks of sunspot cycle 23 can be explained by considering the evolution and interaction of the overlapping activity bands of the longer scale variability.

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 $1\times10^{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 same 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.

Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-fitting to six potential field and six nonpotential field cases simulated with our analytical model, as well as by forward-fitting to an exactly force-free solution of the Low and Lou (1990) model. The forward-fitting tests demonstrate: (i) a satisfactory convergence behavior (with typical misalignment angles of $\mu \approx 1^\circ-10^\circ$), (ii) relatively fast computation times (from seconds to a few minutes), and (iii) the high fidelity of retrieved force-free $\alpha$-parameters ($\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.9-1.0$ for simulations and $\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.7\pm0.3$ for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-forcefree magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing {\sl nonlinear force-free field (NLFFF)} codes based on photospheric magnetograms without coronal constraints.

The three-dimensional (3D) coordinates of stereoscopically triangulated loops provide strong constraints for magnetic field models of active regions in the solar corona. Here we use STEREO/A and B data from some 500 stereoscopically triangulated loops observed in four active regions (2007 Apr 30, May 9, May 19, Dec 11), together with SOHO/MDI line-of-sight magnetograms. We measure the average misalignment angle between the stereoscopic loops and theoretical magnetic field models, finding a mismatch of $\mu=19^\circ-46^\circ$ for a potential field model, which is reduced to $\mu=14^\circ-19^\circ$ for a non-potential field model parameterized by twist parameters. The residual error is commensurable with stereoscopic measurement errors ($\mu_{SE} \approx 8^\circ-12^\circ$). We developed a potential field code that deconvolves a line-of-sight magnetogram into three magnetic field components $(B_x, B_y, B_z)$, as well as a non-potential field forward-fitting code that determines the full length of twisted loops ($L \approx 50-300$ Mm), the number of twist turns (median $N_{twist}=0.06$), the nonlinear force-free $\alpha$-parameter (median $\alpha \approx 4 \times 10^{-11}$ cm$^{-1}$), and the current density (median $j_z \approx 1500$ Mx cm$^{-2}$ s$^{-1}$). All twisted loops are found to be far below the critical value for kink instability, and Joule dissipation of their currents is found be be far below the coronal heating requirement. The algorithm developed here, based on an analytical solution of nonlinear force-free fields that is accurate to second order (in the force-free parameter $\alpha$), represents the first code that enables fast forward-fitting to photospheric magnetograms and stereoscopically triangulated loops in the solar corona.

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.

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.

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 so 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 line's 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^2\approx 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.