No Arabic abstract
The SDO/HMI instruments provide 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 semi-analytic and numeric equilibria and been applied before to vector magnetograms from Hinode and ground based observations. Recently we implemented a new version which takes measurement errors in photospheric vector magnetograms into account. Photospheric field measurements are often due to measurement errors and finite nonmagnetic forces inconsistent as a boundary for a force-free field in the corona. In order to deal with these uncertainties, we developed two improvements: 1.) Preprocessing of the surface measurements in order to make them compatible with a force-free field 2.) The new code keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which have to be optimized for use with data from SDO/HMI. Within this work we describe the corresponding analysis method and evaluate the force-free equilibria by means of how well force-freeness and solenoidal conditions are fulfilled, the angle between magnetic field and electric current and by comparing projections of magnetic field lines with coronal images from SDO/AIA. We also compute the available free magnetic energy and discuss the potential influence of control parameters.
Context: Knowledge about the coronal magnetic field is important to the understanding the structure of the solar corona. We compute the field in the higher layers of the solar atmosphere from the measured photospheric field under the assumption that the corona is force-free. Aims: Here we develop a method for nonlinear force-free coronal magnetic field medelling and preprocessing of photospheric vector magnetograms in spherical geometry using the optimization procedure. Methods: We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates over a restricted area of the Sun. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free equations in the solar corona. We develop a preprocessing procedure in spherical geometry to drive the observed non-force-free data towards suitable boundary conditions for a force-free extrapolation. Results: We test the code with the help of a semi-analytic solution and assess the quality of our reconstruction qualitatively by magnetic field line plots and quantitatively with a number of comparison metrics for different boundary conditions. The reconstructed fields from the lower boundary data with the weighting function are in good agreement with the original reference fields. We added artificial noise to the boundary conditions and tested the code with and without preprocessing. The preprocessing recovered all main structures of the magnetogram and removed small-scale noise. The main test was to extrapolate from the noisy photospheric vector magnetogram with and without preprocessing. The preprocessing was found to significantly improve the agreement between the extrapolated and the exact field.
The solar magnetic field is key to understanding the physical processes in the solar atmosphere. Nonlinear force-free codes have been shown to be useful in extrapolating the coronal field from underlying vector boundary data [see Schrijver et al. 2006 for an overview]. However, we can only measure the magnetic field vector routinely with high accuracy in the photosphere with, e.g., Hinode/SOT, and unfortunately these data do not fulfill the force-free consistency condition as defined by Aly (1989). We must therefore apply some transformations to these data before nonlinear force-free extrapolation codes can be legitimately applied. To this end, we have developed a minimization procedure that uses the measured photospheric field vectors as input to approximate a more chromospheric like field The method was dubbed preprocessing. See Wiegelmann et al. 2006 for details]. The procedure includes force-free consistency integrals and spatial smoothing. The method has been intensively tested with model active regions [see Metcalf et al. 2008] and been applied to several ground based vector magnetogram data before. Here we apply the preprocessing program to photospheric magnetic field measurements with the Hinode/SOT instrument.
In this paper, we came to conclusion that there is a significant systematic error in the SDO/HMI vector magnetic data, which reveals itself in a significant deviation of the lines of the knot magnetic fields from the radial direction. The value of this deviation demonstrates a clear dependence on the distance to the disk center. This paper suggests a method for correction of the vector magnetograms that eliminates the detected systematic error.
In the quiet solar photosphere, the mixed polarity fields form a magnetic carpet, which continuously evolves due to dynamical interaction between the convective motions and magnetic field. This interplay is a viable source to heat the solar atmosphere. In this work, we used the line-of-sight (LOS) magnetograms obtained from the Helioseismic and Magnetic Imager (HMI) on the textit{Solar Dynamics Observatory} (textit{SDO}), and the Imaging Magnetograph eXperiment (IMaX) instrument on the textit{Sunrise} balloon-borne observatory, as time dependent lower boundary conditions, to study the evolution of the coronal magnetic field. We use a magneto-frictional relaxation method, including hyperdiffusion, to produce time series of three-dimensional (3D) nonlinear force-free fields from a sequence of photospheric LOS magnetograms. Vertical flows are added up to a height of 0.7 Mm in the modeling to simulate the non-force-freeness at the photosphere-chromosphere layers. Among the derived quantities, we study the spatial and temporal variations of the energy dissipation rate, and energy flux. Our results show that the energy deposited in the solar atmosphere is concentrated within 2 Mm of the photosphere and there is not sufficient energy flux at the base of the corona to cover radiative and conductive losses. Possible reasons and implications are discussed. Better observational constraints of the magnetic field in the chromosphere are crucial to understand the role of the magnetic carpet in coronal heating.
The minimum-energy configuration for the magnetic field above the solar photosphere is curl-free (hence, by Amperes law, also current-free), so can be represented as the gradient of a scalar potential. Since magnetic fields are divergence free, this scalar potential obeys Laplaces equation, given an appropriate boundary condition (BC). With measurements of the full magnetic vector at the photosphere, it is possible to employ either Neumann or Dirichlet BCs there. Historically, the Neumann BC was used with available line-of-sight magnetic field measurements, which approximate the radial field needed for the Neumann BC. Since each BC fully determines the 3D vector magnetic field, either choice will, in general, be inconsistent with some aspect of the observed field on the boundary, due to the presence of both currents and noise in the observed field. We present a method to combine solutions from both Dirichlet and Neumann BCs to determine a hybrid, least-squares potential field, which minimizes the integrated square of the residual between the potential and actual fields. This has advantages in both not overfitting the radial field used for the Neumann BC, and maximizing consistency with the observations. We demonstrate our methods with SDO/HMI vector magnetic field observations of AR 11158, and find that residual discrepancies between the observed and potential fields are significant, and are consistent with nonzero horizontal photospheric currents. We also analyze potential fields for two other active regions observed with two different vector magnetographs, and find that hybrid potential fields have significantly less energy than the Neumann fields in every case --- by more than 10^(32) erg in some cases. This has major implications for estimates of free magnetic energy in coronal field models, e.g., non-linear force-free field extrapolations.