Variations of the magnetic field within solar coronal mass ejections (CMEs) in the heliosphere depend on the CME`s magnetic structure as it leaves the solar corona and its subsequent evolution through interplanetary space. To account for this evolution, we developed a new numerical model of the inner heliosphere that simulates the propagation of a CME through a realistic background solar wind and allows various CME magnetic topologies. To this end, we incorporate the Gibson-Low CME model within our global MHD model of the inner heliosphere, GAMERA-Helio. We apply the model to study the propagation of the geoeffective CME that erupted on 3 April, 2010 with the aim to reproduce the temporal variations of the magnetic field vector during the CME passage by Earth. Parameters of the Gibson-Low CME are informed by STEREO white-light observations near the Sun. The magnetic topology for this CME - the tethered flux rope - is informed by in-situ magnetic field observations near Earth. We performed two simulations testing different CME propagation directions. For an in-ecliptic direction, the simulation shows a rotation of all three magnetic field components within the CME, as seen at Earth, similar to that observed. With a southward propagation direction, suggested by coronal imaging observations, the modeled By and Bz components are consistent with the ACE data, but the Bx component lacks the observed change from negative to positive. In both cases, the model favors the East-West orientation of the CME flux rope, consistent with the orientation previously inferred from the STEREO/HI heliospheric images.

Coronal loops, seen in solar coronal images, are believed to represent emission from magnetic flux tubes with compact cross sections. We examine the 3D structure of plasma above an active region in a radiative magnetohydrodynamic simulation to locate volume counterparts for coronal loops. In many cases, a loop cannot be linked to an individual thin strand in the volume. While many thin loops are present in the synthetic images, the bright structures in the volume are fewer and of complex shape. We demonstrate that this complexity can form impressions of thin bright loops, even in the absence of thin bright plasma strands. We demonstrate the difficulty of discerning from observations whether a particular loop corresponds to a strand in the volume, or a projection artifact. We demonstrate how apparently isolated loops could deceive observers, even when observations from multiple viewing angles are available. While we base our analysis on a simulation, the main findings are independent from a particular simulation setup and illustrate the intrinsic complexity involved in interpreting observations resulting from line-of-sight integration in an optically thin plasma. We propose alternative interpretation for strands seen in Extreme Ultraviolet images of the corona. The "coronal veil" hypothesis is mathematically more generic, and naturally explains properties of loops that are difficult to address otherwise-such as their constant cross section and anomalously high density scale height. We challenge the paradigm of coronal loops as thin magnetic flux tubes, offering new understanding of solar corona, and by extension, of other magnetically confined bright hot plasmas.

Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientation) similar to weather prediction on Earth.

Magnetism defines the complex and dynamic solar corona. Coronal mass ejections (CMEs) are thought to be caused by stresses, twists, and tangles in coronal magnetic fields that build up energy and ultimately erupt, hurling plasma into interplanetary space. Even the ever-present solar wind possesses a three-dimensional morphology shaped by the global coronal magnetic field, forming geoeffective corotating interaction regions. CME evolution and the structure of the solar wind depend intimately on the coronal magnetic field, so comprehensive observations of the global magnetothermal atmosphere are crucial both for scientific progress and space weather predictions. Although some advances have been made in measuring coronal magnetic fields locally, synoptic measurements of the global coronal magnetic field are not yet available. We conclude that a key goal for 2050 should be comprehensive, ongoing 3D synoptic maps of the global coronal magnetic field. This will require the construction of new telescopes, ground and space-based, to obtain complementary, multiwavelength observations sensitive to the coronal magnetic field. It will also require development of inversion frameworks capable of incorporating multi-wavelength data, and forward analysis tools and simulation testbeds to prioritize and establish observational requirements on the proposed telescopes.

Solar and stellar flares are the most intense emitters of X-rays and extreme ultraviolet radiation in planetary systems(1,2). On the Sun, strong flares are usually found in newly emerging sunspot regions(3). The emergence of these magnetic sunspot groups leads to the accumulation of magnetic energy in the corona. When the magnetic field undergoes abrupt relaxation, the energy released powers coronal mass ejections as well as heating plasma to temperatures beyond tens of millions of kelvins. While recent work has shed light on how magnetic energy and twist accumulate in the corona4 and on how three-dimensional magnetic reconnection allows for rapid energy release(5,6), a self-consistent model capturing how such magnetic changes translate into observable diagnostics has remained elusive. Here, we present a comprehensive radiative magnetohydrodynamics simulation of a solar flare capturing the process from emergence to eruption. The simulation has sufficient realism for the synthesis of remote sensing measurements to compare with observations at visible, ultraviolet and X-ray wavelengths. This unifying model allows us to explain a number of well-known features of solar flares(7), including the time profile of the X-ray flux during flares, origin and temporal evolution of chromospheric evaporation and condensation, and sweeping of flare ribbons in the lower atmosphere. Furthermore, the model reproduces the apparent non-thermal shape of coronal X-ray spectra, which is the result of the superposition of multi-component super-hot plasmas(8) up to and beyond 100 million K.

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 describe the variation of the accumulated coronal helicity derived from the magnetic helicity flux through the photosphere in active region (AR) NOAA 10365, where several large flares and coronal mass ejections (CMEs) occurred. We used SOHO/MDI full-disk line-of-sight magnetograms to measure the helicity flux, and the integral of GOES X-ray flux as a proxy of the coronal energy variations due to flares or CMEs. Using the linear force-free field model, we transformed the accumulated helicity flux into a time sequence of the force-free parameter α accounting for flares or CMEs via the proxy derived from GOES observations. This method can be used to derive the value of α at different times during the AR evolution, and is a partial alternative to the commonly used match of field lines with EUV loops. By combining the accumulated helicity obtained from the observations with the linear force-free theory, we describe the main phases of the emergence process of the AR, and relate them temporally with the occurrence of flares or CMEs. Additionally, a comparison with the loop-matching method of fixing alpha at each time independently shows that the proposed method may be helpful in avoiding unrealistic or undetermined values of alpha that may originate from an insufficient quality of the image used to identify coronal loops at a given time. For the relative intensity of the considered events, the linear force-free field theory implies that there is a direct correlation between the released energy on the one hand and the product of the coronal helicity with the variation of α due to the event on the other. Therefore, the higher the value of the accumulated coronal helicity, the smaller the force-free parameter variation required to produce the same decrease in the free energy during the CMEs.

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. 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 nonlinear force-free field, and (3) thermodynamic models from a fully compressible, 3D MHD simulation of active region (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 Hinode X-ray Telescope 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.

The shapes of solar coronal loops are sensitive to the presence of electrical currents that are the carriers of the non-potential energy available for impulsive activity. We use this information in a new method for modeling the coronal magnetic field of active region (AR) 11158 as a nonlinear force-free field (NLFFF). The observations used are coronal images around the time of major flare activity on 2011 February 15, together with the surface line-of-sight magnetic field measurements. The data are from the Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. The model fields are constrained to approximate the coronal loop configurations as closely as possible, while also being 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 x 10(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 approximate to 0.9 full turns about its axis.

Most one-dimensional hydrodynamic models of plasma confined to magnetic flux tubes assume circular tube cross sections. 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 numerical experiments. We demonstrate that lifting this assumption in favor of realistic, non-circular loops makes the 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 result would impose a bias toward 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 heights even if they are resolved by current instruments.

At present, 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 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 is not at present used as constraints in the computational domain, directly influencing the evolution of the model. This is in part due to the mathematical complications of incorporating such input into numerical models. Projection effects, confusion due to overlapping loops (the coronal plasma is optically thin), and the limited number of usable loops further complicate the use of information from coronal images. We develop and test a new algorithm to use images of coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with those of constant-alpha 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, which uses these trajectories as volume constraints on the values of a, to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying a code and method presented by Wheatland. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques and compare the results with those obtained by currently available methods relying 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 the photospheric magnetic field with information from coronal images. Whereas we focus on the use of coronal loop information in combination with line-of-sight magnetograms, the method is readily extended to incorporate vector-magnetic data over any part of the photospheric boundary.

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 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.

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.

Nonlinear 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 nonlinear force-free fields, described in Low & Lou. 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 α), distribution of twist in the model photosphere, and strength of the magnetic field. We find that (1) for these trial fields, the method reconstructs twist with a mean absolute deviation of at most 15% of the range of photospheric twist, (2) heights of the loops are reconstructed with a mean absolute deviation of at most 5% of the range of trial heights, and (3) the magnitude of non-potential contribution to a photospheric field is reconstructed with a mean absolute deviation of at most 10% of the maximal value.

We introduce two different generalizations of relative helicity which may be applied to a portion of the coronal volume. Such a quantity is generally referred to as the self-helicity of the field occupying the subvolume. Each definition is a natural application of the traditional relative helicity but relative to a different reference field. One of the generalizations, which we term additive self-helicity, can be considered a generalization of twist helicity to volumes which are neither closed nor thin. It shares with twist the property of being identically zero for any portion of a potential magnetic field. The other helicity, unconfined self-helicity, is independent of the shape of the volume occupied by the field portion and is therefore akin to the sum of twist and writhe helicity. We demonstrate how each kind of self-helicity may be evaluated in practice. The set of additive self-helicities may be used as a constraint in the minimization of magnetic energy to produce a piecewise constant-α equilibrium. This class of fields falls into a hierarchy, along with the flux-constrained equilibria and potential fields, of fields with monotonically decreasing magnetic energies. Piecewise constant-α fields generally have fewer unphysical properties than genuinely constant-α fields, whose twist α is uniform throughout the entire corona.