Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

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.

Between 2008/3/24 and 2008/4/2, the three active regions NOAA active regions 10987, 10988 and 10989 were observed daily by the Synoptic Optical Long-term Investigations of the Sun (SOLIS) Vector Spectro-Magnetograph (VSM) while they traversed the solar disk. We use these measurements and the nonlinear force-free magnetic field code XTRAPOL to reconstruct the coronal magnetic field for each active region and compare model field lines with images from the Solar Terrestrial RElations Observatory (STEREO) and Hinode X-ray Telescope (XRT) telescopes. Synoptic maps made from continuous, round-the-clock Global Oscillations Network Group (GONG) magnetograms provide information on the global photospheric field and potential-field source-surface models based on these maps describe the global coronal field during the Whole Heliospheric Interval (WHI) and its neighboring rotations. Features of the modeled global field, such as the coronal holes and streamer belt locations, are discussed in comparison with extreme ultra-violet and coronagraph observations from STEREO. The global field is found to be far from a minimum, dipolar state. From the nonlinear models we compute physical quantities for the active regions such as the photospheric magnetic and electric current fluxes, the free magnetic energy and the relative helicity for each region each day where observations permit. The interconnectivity of the three regions is addressed in the context of the potential-field source-surface model. Using local and global quantities derived from the models, we briefly discuss the different observed activity levels of the regions.

We consider a three-dimensional bipolar force-free magnetic field with non zero magnetic helicity, occupying a half-space, and study the problem of its evolution driven by an imposed photospheric flux decrease. For this specific setting of the Flux Cancellation Model describing coronal mass ejections occuring in active regions, we address the issues of the physical meaning of flux decrease, of the influence on field evolution of the size of the domain over which this decrease is imposed, and of the existence of an energetic criterion characterizing the possible onset of disruption of the configuration. We show that: (1) The imposed flux disappearance can be interpreted in terms of transport of positive and negative fluxes towards the inversion line, where they get annihilated. (2) For the particular case actually computed, in which the initial state is quite sheared, the formation of a twisted flux rope and the subsequent global disruption of the configuration are obtained when the flux has decreased by only a modest amount over a limited part of the whole active region. (3) The disruption is produced when the magnetic energy becomes of the order of the decreasing energy of a semi-open field, and then before reaching the energy of the associated fully open field. This suggests that the mechanism leading to the disruption is nonequilibrium as in the case where flux is imposed to decrease over the whole region.

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and powered by the partial dissipation of electrical currents in the solar atmosphere. A critical unanswered question is whether the currents involved are induced by the motion of pre-existing atmospheric magnetic flux subject to surface plasma flows, or whether these currents are associated with the emergence of flux from within the solar convective zone. We address this problem by applying state-of-the-art nonlinear force-free field (NLFFF) modeling to the highest resolution and quality vector-magnetographic data observed by the recently launched Hinode satellite on NOAA Active Region 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with 4 different codes and a variety of boundary conditions. We find that the model fields differ markedly in geometry, energy content, and force-freeness. We discuss the relative merits of these models in a general critique of present abilities to model the coronal magnetic field based on surface vector field measurements. For our application in particular, we find a fair agreement of the best-fit model field with the observed coronal configuration, and argue (1) that strong electrical currents emerge together with magnetic flux preceding the flare, (2) that these currents are carried in an ensemble of thin strands, (3) that the global pattern of these currents and of field lines are compatible with a large-scale twisted flux rope topology, and (4) that the ~10^32 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.