Thomas Wiegelmann

A Critical Assessment of Nonlinear Force-Free Field Modeling of the Solar Corona for Active Region 10953

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

Optimization code with weighting function for the reconstruction of coronal magnetic fields

We developed a code for the reconstruction of nonlinear force-free and non-force-free coronal magnetic fields. The 3D magnetic field is computed numerically with the help of an optimization principle. The force-free and non-force-free codes are compiled in one program. The force-free approach needs photospheric vector magnetograms as input. The non-force-free code additionally requires the line-of-sight integrated coronal density distribution in combination with a tomographic inversion code. Previously the optimization approach has been used to compute magnetic fields using all six boundaries of a computational box. Here we extend this method and show how the coronal magnetic field can be reconstructed only from the bottom boundary, where the boundary conditions are measured with vector magnetographs. The program is planed for use within the Stereo mission.

Nonlinear Force‐free Field Modeling of a Solar Active Region around the Time of a Major Flare and Coronal Mass Ejection

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and are 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 preexisting 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 AR 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with four 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 ~1032 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.

Evolution of Magnetic Field and Energy in a Major Eruptive Active Region Based on SDO/HMI Observation

We report the evolution of magnetic field and its energy in NOAA active region 11158 over 5 days based on a vector magnetogram series from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). Fast flux emergence and strong shearing motion led to a quadrupolar sunspot complex that produced several major eruptions, including the first X-class flare of Solar Cycle 24. Extrapolated non-linear force-free coronal fields show substantial electric current and free energy increase during early flux emergence near a low-lying sigmoidal filament with sheared kilogauss field in the filament channel. The computed magnetic free energy reaches a maximum of ∼2.6 × 10 32 erg, about 50% of which is stored below 6 Mm. It decreases by ∼0.3 × 10 32 erg within 1 hour of the X-class flare, which is likely an underestimation of the actual energy loss. During the flare, the photospheric field changed rapidly: horizontal field was enhanced by 28% in the core region, becoming more inclined and more parallel to the polarity inversion line. Such change is consistent with the conjectured coronal field “implosion”, and is supported by the coronal loop retraction observed by the Atmospheric Imaging Assembly (AIA). The extrapolated field becomes more “compact” after the flare, with shorter loops in the core region, probably because of reconnection. The coronal field becomes slightly more sheared in the lowest layer, relaxes faster with height, and is overall less energetic.

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.

Coexisting Flux Rope and Dipped Arcade Sections Along One Solar Filament

We compute the three-dimensional magnetic field of an active region in order to study the magnetic configuration of active region filaments. The nonlinear force-free field model is adopted to compute the magnetic field above the photosphere, where the vector magnetic field was observed by THEMIS/MTR on 2005 May 27. We propose a new method to remove the 180 ◦ ambiguity of the transverse field. Next, we analyze the implications of the preprocessing of the data by minimizing the total force and torque in the observed vector fields. This step provides a consistent bottom boundary condition for the nonlinear force-free field model. Then, using the optimization method to compute the coronal field, we find a magnetic flux rope along the polarity inversion line. The magnetic flux rope aligns well with part of an Hα filament, while the total distribution of the magnetic dips coincides with the whole Hα filament. This implies that the magnetic field structure in one section of the filament is a flux rope, while the other is a sheared arcade. The arcade induced a left-bearing filament in the magnetic field of negative helicity, which is opposite to the chirality of barbs that a flux rope would induce in a magnetic field of the same helicity sign. The field strength in the center of the flux rope is about 700 G, and the twist of the field lines is ∼1.4 turns.

Coronal plasma flows and magnetic fields in solar active regions: Combined observations from SOHO and NSO/Kitt Peak

During the early days of the SOHO mission, SUMER observed a few active regions (ARs) connected with sunspots on the Sun and took their images and spectra in various EUV emission lines. In addition to these spectroscopic data magne- tograms of the photospheric footpoint regions of the AR loops were available from the MDI on SOHO and the National Solar Observatory/Kitt Peak (NSO/KP), data which here are used to construct the coronal magnetic field of the ARs by force-free- field extrapolation. The combined data set is analysed with respect to the large-scale circulation of coronal matter, which means that the Dopplershifts of various lines used as tracers of the plasma flow are investigated in close connection with the ambient magnetic field, which is found to be either closed or open in the coronal volume considered. The Dopplershift pattern is found to be clearly linked with the field topology, and several regions of strong velocity shear are identified. We also estimate the coronal currents. We discuss the results of this mainly phenomenological correlative study with the perspective to understand coronal heating and mass supply to the extended corona, and with respect to the role played by the field in guiding and constraining plasma flows.

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.

Driving Mechanism and Onset Condition of a Confined Eruption

We study a confined eruption accompanied by an M1.1 flare in solar active region (AR) NOAA 10767 on 2005 May 27, where a pre-eruptive magnetic flux rope was reported in a nonlinear force-free field (NLFFF) extrapolation. The observations show a strong writhing motion of the erupting structure, suggesting that a flux rope was indeed present and converted some of its twist into writhe in the course of the eruption. Using the NLFFF extrapolation, we calculate the twist of the pre-eruptive flux rope and find that it is in very good agreement with thresholds of the helical kink instability found in numerical simulations. We conclude that the activation and rise of the flux rope were triggered and driven by the instability. Using a potential field extrapolation, we also estimate the height distribution of the decay index of the external magnetic field in the AR 1 hr prior to the eruption. We find that the decay index stays below the threshold for the torus instability for a significant height range above the erupting flux rope. This provides a possible explanation for the confinement of the eruption to the low corona.

Comparing magnetic field extrapolations with measurements of magnetic loops

We compare magnetic field extrapolations from a photospheric magnetogram with the observationally inferred struc- ture of magnetic loops in a newly developed active region. This is the first time that the reconstructed 3D-topology of the magnetic field is available to test the extrapolations. We compare the observations with potential fields, linear force-free fields and non-linear force-free fields. This comparison reveals that a potential field extrapolation is not suitable for a reconstruction of the magnetic field in this young, developing active region. The inclusion of field-line-parallel electric currents, the so called force-free approach, gives much better results. Furthermore, a non-linear force-free computation reproduces the observations better than the linear force-free approximation, although no free parameters are available in the former case.