Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDeducing the reliability of relative helicities from nonlinear force-free coronal models
74
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study the relative helicity of active region (AR) NOAA~12673 during a ten-hour time interval centered around a preceding X2.2 flare (SOL2017-09-06T08:57) and also including an eruptive X9.3 flare that occurred three hours later (SOL2017-09-06T11:53). In particular, we aim for a reliable estimate of the normalized self-helicity of the current-carrying magnetic field, the so-called helicity ratio $|H_{mathrm{J}}|/|H_{mathcal{V}}|$, a promising candidate to quantity the eruptive potential of solar ARs. Using SDO/HMI vector magnetic field data as an input, we employ nonlinear force-free (NLFF) coronal magnetic field models using an optimization approach. The corresponding relative helicity, and related quantities, are computed using a finite-volume method. From multiple time series of NLFF models based on different choices of free model parameters, we are able to assess the spread of $|H_{mathrm{J}}|/|H_{mathcal{V}}|$, and to estimate its uncertainty. In comparison to earlier works, which identified the non-solenoidal contribution to the total magnetic energy, $E_{rm div}/E$, as selection criterion regarding the required solenoidal quality of magnetic field models for subsequent relative helicity analysis, we propose to use in addition the non-solenoidal contribution to the free magnetic energy, $|E_{rm mix}|/E_{mathrm{J,s}}$. As a recipe for a reliable estimate of the relative magnetic helicity (and related quantities), we recommend to employ multiple NLFF models based on different combinations of free model parameters, to retain only those that exhibit smallest values of both $E_{rm div}/E$ and $|E_{rm mix}|/E_{mathrm{J,s}}$ at a certain time instant, to subsequently compute mean estimates, and to use the spread of the individually contributing values as an indication for the uncertainty.
rate research
Read More
We demonstrate the sensitivity of magnetic energy and helicity computations regarding the quality of the underlying coronal magnetic field model. We apply the method of Wiegelmann & Inhester (2010) to a series of SDO/HMI vector magnetograms, and discuss nonlinear force-free (NLFF) solutions based on two different sets of the free model parameters. The two time series differ from each other concerning their force-free and solenoidal quality. Both force- and divergence-freeness are required for a consistent NLFF solution. Full satisfaction of the solenoidal property is inherent in the definition of relative magnetic helicity in order to insure gauge-independence. We apply two different magnetic helicity computation methods (Thalmann et al. 2011; Valori et al. 2012) to both NLFF time series and find that the output is highly dependent on the level to which the NLFF magnetic fields satisfy the divergence-free condition, with the computed magnetic energy being less sensitive than the relative helicity. Proxies for the non-potentiality and eruptivity derived from both quantities are also shown to depend strongly on the solenoidal property of the NLFF fields. As a reference for future applications, we provide quantitative thresholds for the force- and divergence-freeness, for the assurance of reliable computation of magnetic energy and helicity, and of their related eruptivity proxies.
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.
Presently, 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 (NLFFF) 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 at present is not used due to the mathematical complications of incorporating such input into numerical models. Projection effects and the limited number of usable loops further complicate the use of coronal information. We develop and test an algorithm to use images showing coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with field lines of constant-als 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 to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying method presented in citet{Wheatland2007}. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques citep{Schrijver2006, Schrijver2008} and compare the results with those obtained by presently available methods that rely 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 photospheric magnetic field with information from coronal images. Vector magnetograms over the full or partial photospheric boundary of the numerical domain could optionally be used.
Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the self-consistency procedure of Wheatland & Regnier (2009). In that procedure, the inconsistency is resolved by an iterative process of constructing two solutions based on the values of the force-free parameter alpha on the two polarities of the field in the boundary (the P and N polarities), and taking uncertainty-weighted averages of the boundary alpha values in the P and N solutions. When the alpha values in the P and N regions are very different, the self-consistent solution may lose high alpha values from the boundary conditions. We show how, by altering the weighting of the uncertainties in the P or N boundary conditions, we can preserve high alpha values in the self-consistent solution. The weighted self-consistent extrapolation method is demonstrated on an analytic bipole field and applied to vector magnetogram data taken by the Helioseismic and Magnetic Imager (HMI) instrument for NOAA active region AR 12017 on 2014 March 29.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
