Preprocessing the Photospheric Vector Magnetograms for an NLFFF Extrapolation Using a Potential-Field Model and an Optimization Method

Numerical reconstruction/extrapolation of the coronal nonlinear force-free magnetic field (NLFFF) usually takes the photospheric vector magnetogram as input at the bottom boundary. The magnetic field observed at the photosphere, however, contains a force that is in conflict with the fundamental assumption of the force-free model. It also contains measurement noise, which hinders the practical computation. Wiegelmann, Inhester, and Sakurai (Solar Phys. 233, 215, 2006) have proposed to preprocess the raw magnetogram to remove the force and noise to provide better input for NLFFF modeling. In this paper we develop a new code of magnetogram preprocessing that is consistent with our extrapolation method CESE–MHD–NLFFF (Jiang, Feng, and Xiang in Astrophys. J. 755, 62, 2012; Jiang and Feng in Astrophys. J. 749, 135, 2012a). Based on the magnetic-splitting rule that a magnetic field can be split into a potential-field part and a non-potential part, we split the magnetogram and dealt with the two parts separately. The preprocessing of the magnetogram’s potential part is based on a numerical potential-field model, and the non-potential part is preprocessed using the similar optimization method of Wiegelmann, Inhester, and Sakurai (2006). The code was applied to the SDO/HMI data, and results show that the method can remove the force and noise efficiently and improve the extrapolation quality.     

How Should One Optimize Nonlinear Force-Free Coronal Magnetic Field Extrapolations from SDO/HMI Vector Magnetograms?

The Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) provides photospheric vector magnetograms with a high spatial and temporal resolution. Our intention is to model the coronal magnetic field above active regions with the help of a nonlinear force-free extrapolation code. Our code is based on an optimization principle and has been tested extensively with semianalytic and numeric equilibria and applied to vector magnetograms from Hinode and ground-based observations. Recently we implemented a new version which takes into account measurement errors in photospheric vector magnetograms. Photospheric field measurements are often affected by measurement errors and finite nonmagnetic forces inconsistent for use as a boundary for a force-free field in the corona. To deal with these uncertainties, we developed two improvements: i) preprocessing of the surface measurements to make them compatible with a force-free field, and ii) new code which keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which must be optimized for use with data from SDO/HMI. In this work we describe the corresponding analysis method and evaluate the force-free equilibria by how well force-freeness and solenoidal conditions are fulfilled, by the angle between magnetic field and electric current, and by comparing projections of magnetic field lines with coronal images from the Atmospheric Imaging Assembly (SDO/AIA). We also compute the available free magnetic energy and discuss the potential influence of control parameters.     

A Comparison Between Nonlinear Force-Free Field and Potential Field Models Using Full-Disk SDO/HMI Magnetogram

Measurements of magnetic fields and electric currents in the pre-eruptive corona are crucial to the study of solar eruptive phenomena, like flares and coronal mass ejections (CMEs). However, spectro-polarimetric measurements of certain photospheric lines permit a determination of the vector magnetic field only at the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field above multiple active regions with the help of a potential field and a nonlinear force-free field (NLFFF) extrapolation code over the full solar disk using Helioseismic and Magnetic Imager (SDO/HMI) data as boundary conditions. We compare projections of the resulting magnetic field lines with full-disk coronal images from the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study has found that the NLFFF model reconstructs the magnetic configuration closer to observation than the potential field model for full-disk magnetic field extrapolation. We conclude that many of the trans-equatorial loops connecting the two solar hemispheres are current-free.     

Preprocessing of Vector Magnetograph Data for a Nonlinear Force-Free Magnetic Field Reconstruction

Knowledge regarding the coronal magnetic field is important for the understanding of many phenomena, like flares and coronal mass ejections. Because of the low plasma beta in the solar corona, the coronal magnetic field is often assumed to be force-free and we use photospheric vector magnetograph data to extrapolate the magnetic field into the corona with the help of a nonlinear force-free optimization code. Unfortunately, the measurements of the photospheric magnetic field contain inconsistencies and noise. In particular, the transversal components (say B _x and B _y) of current vector magnetographs have their uncertainties. Furthermore, the magnetic field in the photosphere is not necessarily force free and often not consistent with the assumption of a force-free field above the magnetogram. We develop a preprocessing procedure to drive the observed non–force-free data towards suitable boundary conditions for a force-free extrapolation. As a result, we get a data set which is as close as possible to the measured data and consistent with the force-free assumption.     

Nonlinear Force-Free Modeling of Coronal Magnetic Fields. II. Modeling a Filament Arcade and Simulated Chromospheric and Photospheric Vector Fields

We compare a variety of nonlinear force-free field (NLFFF) extrapolation algorithms, including optimization, magneto-frictional, and Grad – Rubin-like codes, applied to a solar-like reference model. The model used to test the algorithms includes realistic photospheric Lorentz forces and a complex field including a weakly twisted, right helical flux bundle. The codes were applied to both forced “photospheric” and more force-free “chromospheric” vector magnetic field boundary data derived from the model. When applied to the chromospheric boundary data, the codes are able to recover the presence of the flux bundle and the field’s free energy, though some details of the field connectivity are lost. When the codes are applied to the forced photospheric boundary data, the reference model field is not well recovered, indicating that the combination of Lorentz forces and small spatial scale structure at the photosphere severely impact the extrapolation of the field. Preprocessing of the forced photospheric boundary does improve the extrapolations considerably for the layers above the chromosphere, but the extrapolations are sensitive to the details of the numerical codes and neither the field connectivity nor the free magnetic energy in the full volume are well recovered. The magnetic virial theorem gives a rapid measure of the total magnetic energy without extrapolation though, like the NLFFF codes, it is sensitive to the Lorentz forces in the coronal volume. Both the magnetic virial theorem and the Wiegelmann extrapolation, when applied to the preprocessed photospheric boundary, give a magnetic energy which is nearly equivalent to the value derived from the chromospheric boundary, but both underestimate the free energy above the photosphere by at least a factor of two. We discuss the interpretation of the preprocessed field in this context. When applying the NLFFF codes to solar data, the problems associated with Lorentz forces present in the low solar atmosphere must be recognized: the various codes will not necessarily converge to the correct, or even the same, solution. On 07/07/2007, the NLFFF team was saddened by the news that Tom Metcalf had died as the result of an accident. We remain grateful for having had the opportunity to benefit from his unwavering dedication to the problems encountered in attempting to understand the Sun’s magnetic field; Tom had completed this paper several months before his death, leading the team through the many steps described above.     

Modeling coronal magnetic field using spherical geometry: cases with several active regions

The magnetic fields in the solar atmosphere structure the plasma, store free magnetic energy and produce a wide variety of active solar phenomena, like flare and coronal mass ejections (CMEs). The distribution and strength of magnetic fields are routinely measured in the solar surface (photosphere). Therefore, there is considerable interest in accurately modeling the 3D structure of the coronal magnetic field using photospheric vector magnetograms. 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. Nonlinear force-free field (NLFFF) models are thought to be viable tools for those task. Usually those models use Cartesian geometry. However, the spherical nature of the solar surface cannot be neglected when the field of view is large. In this work, we model the coronal magnetic field above multiple active regions using NLFFF extrapolation code using vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/Vector Spectromagnetograph (VSM) as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with their respective coronal EUV-images from the Atmospheric Imaging Assembly (SDO/AIA) observed on October 15, 2011 and November 13, 2012. This study has found that the NLFFF model in spherical geometry reconstructs the magnetic configurations for several active regions which agrees to some extent with observations. During October 15, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.     

Global Solar Free Magnetic Energy and Electric Current Density Distribution of Carrington Rotation 2124

Solar eruptive phenomena, like flares and coronal mass ejections (CMEs), are governed by magnetic fields. To describe the structure of these phenomena one needs information on the magnetic flux density and the electric current density vector components in three dimensions throughout the atmosphere. However, current spectro-polarimetric measurements typically limit the determination of the vector magnetic field to only the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field for global solar atmosphere using nonlinear force-free field (NLFFF) extrapolation codes implemented to a synoptic maps of photospheric vector magnetic field synthesized from the Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) as boundary condition. Using the resulting three-dimensional magnetic field, we calculate the three-dimensional electric current density and magnetic energy throughout the solar atmosphere for Carrington rotation 2124 using our global extrapolation code. We found that spatially, the low-lying, current-carrying core field demonstrates a strong concentration of free energy in the active-region core, from the photosphere to the lower corona (about 70 Mm). The free energy density appears largely co-spatial with the electric current distribution.     

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests

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 (Astrophys. J. 352, 343, 1990) model. The forward-fitting tests demonstrate: i) a satisfactory convergence behavior (with typical misalignment angles of μ≈1^°?–?10^°), ii) relatively fast computation times (from seconds to a few minutes), and iii) the high fidelity of retrieved force-free α-parameters (α _fit/α _model≈0.9?–?1.0 for simulations and α _fit/α _model≈0.7±0.3 for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-force-free magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing nonlinear force-free field (NLFFF) codes based on photospheric magnetograms without coronal constraints.     

A Magnetostatic Grad�CRubin Code for Coronal Magnetic Field Extrapolations

The coronal magnetic field cannot be directly observed, but, in principle, it can be reconstructed from the comparatively well observed photospheric magnetic field. A?popular approach uses a nonlinear force-free model. Non-magnetic forces at the photosphere are significant, meaning the photospheric data are inconsistent with the force-free model, and this causes problems with the modeling (De Rosa et?al., Astrophys. J. 696, 1780, 2009). In this paper we present a numerical implementation of the Grad?CRubin method for reconstructing the coronal magnetic field using a magnetostatic model. This model includes a pressure force and a non-zero magnetic Lorentz force. We demonstrate our implementation on a simple analytic test case and obtain the speed and numerical error scaling as a function of the grid size.     

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.     

Including stereoscopic information in the reconstruction of coronal magnetic fields

We present a method to include stereoscopic information about the three-dimensional structure of flux tubes into the reconstruction of the coronal magnetic field. Due to the low plasma beta in the corona we can assume a force-free magnetic field, with the current density parallel to the magnetic field lines. Here we use linear force-free fields for simplicity. The method uses the line-of-sight magnetic field on the photosphere as observational input. The value of is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force-free solution constrained by both the photospheric magnetogram and the observed plasma structures. As an example we apply our method to SOHO MDI/EIT data of an active region. In the future it is planned to apply the method to analyse data from the SECCHI instrument aboard the STEREO mission.     

The Force-Free Electrodynamics Method for the Extrapolation of Coronal Magnetic Fields from Vector Magnetograms

We present a new improved version of our force-free electrodynamics (FFE) numerical code in spherical coordinates that extrapolates the magnetic field in the inner solar corona from a photospheric vector magnetogram. The code satisfies the photospheric boundary condition and the condition ??B=0 to machine accuracy. The performance of our method is evaluated with standard convergence parameters, and is found to be comparable to that of other nonlinear force-free extrapolations.     

Computation of inner coronal magnetic fields from longitudinal field components on a spherical photosphere

A method is developed which for a certain day permits the approximate calculation of closed small and large scale magnetic field lines. From the photospheric longitudinal components of the magnetic field measured at this day normal components are derived taking into account the curvature of the solar surface. The magnetic fields are assumed to be potential or force-free fields.The method is applied to observations of September 5 and September 7, 1973. The projected magnetic field lines are compared with the loop structures which are visible in XUV pictures taken on these days. In the cases where no good agreement could be obtained for potential fields, force-free fields are calculated and fitted to the observed structures.     

The effects of viewing angle on the inference of magnetic shear in preflare active regions

The magnetic shear at a point within an active region field configuration can be defined (Hagyard et al., 1984b) as the difference in angle between the observed photospheric transverse field and that of a reference potential field calculated using the observed line-of-sight field as a boundary condition. Using analytic models for non-potential (but force-free) fields representative of preflaring active regions, we calculate the degree of magnetic shear along the magnetic neutral line that such fields would exhibit, as a function of the location and orientation of the active region on the solar disk. We find that, except for regions close to disk center, the position of the inferred neutral line (zero line-of-sight field) is significantly different from the actual neutral line (zero radial field), and that the calculated reference potential field also varies significantly with the position of the region. Thus the inferred degree of shear can vary significantly with the position and orientation of the region, due to (a) straightforward geometric projection effects, (b) the shift of the inferred neutral line relative to its true position, and (c) variations in the reference potential field. The significance of these results for flare prediction is considered.Presidential Young Investigator.     

Computation of solar magnetic fields from photospheric observations

The observational difficulties of obtaining the magnetic field distribution in the chromosphere and corona of the Sun has led to methods of extending photospheric magnetic measurements into the solar atmosphere by mathematical procedures. A new approach to this problem presented here is that a constant alpha force-free field can be uniquely determined from the tangential components of the measured photospheric flux alone. The vector magnetographs now provide measurements of both the solar photospheric tangential and the longitudinal magnetic field. This paper presents derivations for the computation of the solar magnetic field from these type of measurements. The fields considered are assumed to be a constant alpha force-free fields or equivalent, producing vanishing Lorentz forces. Consequently, magnetic field lines and currents are related by a constant and hence show an identical distribution. The magnetic field above simple solar regions are described from the solution of the field equations.     

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. I. Theory

We derive an analytical approximation of nonlinear force-free magnetic field solutions (NLFFF) that can efficiently be used for fast forward-fitting to solar magnetic data, constrained either by observed line-of-sight magnetograms and stereoscopically triangulated coronal loops, or by 3D vector-magnetograph data. The derived NLFFF solutions provide the magnetic field components B _x (x), B _y (x), B _z (x), the force-free parameter α(x), the electric current density j(x), and are accurate to second-order (of the nonlinear force-free α-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.     

Solar Cycle Variation of Magnetic Flux Ropes in a Quasi-Static Coronal Evolution Model

The structure of electric current and magnetic helicity in the solar corona is closely linked to solar activity over the 11-year cycle, yet is poorly understood. As an alternative to traditional current-free “potential-field” extrapolations, we investigate a model for the global coronal magnetic field which is non-potential and time-dependent, following the build-up and transport of magnetic helicity due to flux emergence and large-scale photospheric motions. This helicity concentrates into twisted magnetic flux ropes, which may lose equilibrium and be ejected. Here, we consider how the magnetic structure predicted by this model – in particular the flux ropes – varies over the solar activity cycle, based on photospheric input data from six periods of cycle 23. The number of flux ropes doubles from minimum to maximum, following the total length of photospheric polarity inversion lines. However, the number of flux rope ejections increases by a factor of eight, following the emergence rate of active regions. This is broadly consistent with the observed cycle modulation of coronal mass ejections, although the actual rate of ejections in the simulation is about a fifth of the rate of observed events. The model predicts that, even at minimum, differential rotation will produce sheared, non-potential, magnetic structure at all latitudes.     

Probing Twisted Magnetic Field Using Microwave Observations in an M Class Solar Flare on 11 February, 2014

This work demonstrates the possibility of magnetic-field topology investigations using microwave polarimetric observations. We study a solar flare of GOES M1.7 class that occurred on 11 February, 2014. This flare revealed a clear signature of spatial inversion of the radio-emission polarization sign. We show that the observed polarization pattern can be explained by nonthermal gyrosynchrotron emission from the twisted magnetic structure. Using observations of the Reuven Ramaty High Energy Solar Spectroscopic Imager, Nobeyama Radio Observatory, Radio Solar Telescope Network, and Solar Dynamics Observatory, we have determined the parameters of nonthermal electrons and thermal plasma and identified the magnetic structure where the flare energy release occurred. To reconstruct the coronal magnetic field, we use nonlinear force-free field (NLFFF) and potential magnetic-field approaches. Radio emission of nonthermal electrons is simulated by the GX Simulator code using the extrapolated magnetic field and the parameters of nonthermal electrons and thermal plasma inferred from the observations; the model radio maps and spectra are compared with observations. We have found that the potential-magnetic-field approach fails to explain the observed circular polarization pattern; on the other hand, the Stokes-\(V\) map is successfully explained by assuming nonthermal electrons to be distributed along the twisted magnetic structure determined by the NLFFF extrapolation approach. Thus, we show that the radio-polarization maps can be used for diagnosing the topology of the flare magnetic structures where nonthermal electrons are injected.     

RECONSTRUCTION OF THE THREE-DIMENSIONAL CORONAL MAGNETIC FIELD

Studies of solar flares indicate that the mechanism of flares is magnetic in character and that the coronal magnetic field is a key to understanding solar high-energy phenomena. In our ongoing research we are conducting a systematic study of a large database of observations which includes both coronal structure (from the Soft X-ray Telescope on the Yohkoh spacecraft) and photospheric vector magnetic fields (from the Haleakala Stokes Polarimeter at Mees Solar Observatory). We compare the three-dimensional nonlinear force-free coronal magnetic field, computed from photospheric boundary data, to images of coronal structure. In this paper we outline our techniques and present results for active region AR 7220/7222. We show that the computed force-free coronal magnetic field agrees well with Yohkoh X-ray coronal loops, and we discuss the properties of the coronal magnetic field and the soft X-ray loops.