In order to improve our understanding on the pre-requisites of eruptive solar flares, we study and compare different measures that characterize the eruptive potential of solar active regions - the critical height for torus instability as a local measure and the helicity ratio as a global measure - with the structural properties of the underlying magnetic field, namely the altitude of the center of the current-carrying magnetic structure. Using time series of 3D optimization-based nonlinear force-free magnetic field models for 10 different active regions (ARs) around the time of large solar flares, we determine the altitudes of the current-weighted centers of the non-potential model structures. Based on the potential magnetic field, we inspect the decay index, n, in multiple vertical planes oriented along of or perpendicular to the flare-relevant polarity inversion line, and estimate the critical height (h_crit) for torus instability (TI) using different thresholds of n. The critical heights are interpreted with respect to the altitudes of the current-weighted centers of the associated non-potential structures, as well as the eruptive character of the associated flares, and the eruptive potential of the host AR, as characterized by the helicity ratio. Our most important findings are that (i) h_crit is more segregated in terms of flare type than the helicity ratio, and that (ii) coronal field configurations with a higher eruptive potential (in terms of the helicity ratio) also appear to be more prone to TI. Furthermore, we find no pronounced differences in the altitudes of the non-potential structures prior to confined and eruptive flares.
Context. Nonlinear force-free (NLFF) modeling is regularly used to indirectly infer the 3D geometry of the coronal magnetic field, which is not otherwise accessible on a regular basis by means of direct measurements. Aims. We study the effect of binning in time-series NLFF modeling of individual active regions (ARs) in order to quantify the effect of a different underlying spatial sampling on the quality of modeling as well as on the derived physical parameters. Methods. We apply an optimization method to sequences of Solar Dynamics Observatory (SDO) Helioseismic and Magnetic Imager (HMI) vector magnetogram data at three different plate scales for three solar active regions to obtain nine NLFF model time series. From the NLFF models, we deduce active-region magnetic fluxes, electric currents, magnetic energies, and relative helicities, and analyze those with respect to the underlying spatial sampling. We calculate various metrics to quantify the quality of the derived NLFF models and apply a Helmholtz decomposition to characterize solenoidal errors. Results. At a given spatial sampling, the quality of NLFF modeling is different for different ARs, and the quality varies along the individual model time series. For a given AR, modeling at a certain spatial sampling is not necessarily of superior quality compared to that performed with a different plate scale. Generally, the NLFF model quality tends to be higher for larger pixel sizes with the solenoidal quality being the ultimate cause for systematic variations in model-deduced physical quantities. Conclusions. Optimization-based modeling using SDO/HMI vector data binned to larger pixel sizes yields variations in magnetic energy and helicity estimates of ≲30% on overall, given that concise checks ensure the physical plausibility and high solenoidal quality of the tested model. Spatial-sampling-induced differences are relatively small compared to those arising from other sources of uncertainty, including the effects of applying different data calibration methods, those of using vector data from different instruments, or those arising from application of different NLFF methods to identical input data.
Context. In order to better understand the underlying processes and prerequisites for solar activity, it is essential to study the time evolution of the coronal magnetic field of solar active regions (ARs) associated with flare activity. Aims. We investigate the coronal magnetic energy and helicity budgets of ten solar ARs around the times of large flares. In particular, we are interested in a possible relation of the derived quantities to the particular type of the flares that the AR produces, namely, whether they are associated with a CME or whether they are confined (i.e., not accompanied by a CME). Methods. Using an optimization approach, we employed time series of 3D nonlinear force-free magnetic field models of ten ARs, covering a time span of several hours around the time of occurrence of large solar flares (GOES class M1.0 and larger). We subsequently computed the 3D magnetic vector potentials associated to the model 3D coronal magnetic field using a finite-volume method. This allows us to correspondingly compute the coronal magnetic energy and helicity budgets, as well as related (intensive) quantities such as the relative contribution of free magnetic energy, EF/E (energy ratio), the fraction of non-potential (current-carrying) helicity, |HJ|/|HV| (helicity ratio), and the normalized current-carrying helicity, |HJ|/ϕ′2. Results. The total energy and helicity budgets of flare-productive ARs (extensive parameters) cover a broad range of magnitudes, with no obvious relation to the eruptive potential of the individual ARs, that is, whether or not a CME is produced in association with the flare. The intensive eruptivity proxies, EF/E and |HJ|/|HV|, and |HJ|/ϕ′2, however, seem to be distinctly different for ARs that produce CME-associated large flares compared to those which produce confined flares. For the majority of ARs in our sample, we are able to identify characteristic pre-flare magnitudes of the intensive quantities that are clearly associated with subsequent CME-productivity. Conclusions. If the corona of an AR exhibits characteristic values of ⟨|HJ|/|HV|⟩ > 0.1, ⟨EF/E⟩ > 0.2, and ⟨|HJ|/ϕ′2⟩ > 0.005, then the AR is likely to produce large CME-associated flares. Conversely, confined large flares tend to originate from ARs that exhibit coronal values of ⟨|HJ|/|HV|⟩ ≲ 0.1, ⟨EF/E⟩ ≲ 0.1, and ⟨|HJ|/ϕ′2⟩ ≲ 0.002.
Aims. 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, |HJ|/|H𝒱|, a promising candidate to quantity the eruptive potential of solar ARs. Methods. Using Solar Dynamics Observatory Helioseismic and Magnetic Imager 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 |HJ|/|H𝒱|, and to estimate its uncertainty. Results. In comparison to earlier works, which identified the non-solenoidal contribution to the total magnetic energy, Ediv/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, |Emix|/EJ, 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 Ediv/E and |Emix|/EJ, 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.