We compare the coronal magnetic energy and helicity of two solar active regions (ARs), prolific in major eruptive (AR~11158) and confined (AR~12192) flaring, and analyze the potential of deduced proxies to forecast upcoming flares. Based on nonlinear force-free (NLFF) coronal magnetic field models with a high degree of solenoidality, and applying three different computational methods to investigate the coronal magnetic helicity, we are able to draw conclusions with a high level of confidence. Based on real observations of two solar ARs we checked trends regarding the potential eruptivity of the active-region corona, as suggested earlier in works that were based on numerical simulations, or solar observations. Our results support that the ratio of current-carrying to total helicity, $|H_\mathrm{J}|/|H_\mathrm{V}|$, shows a strong ability to indicate the eruptive potential of a solar AR. However, $|H_\mathrm{J}|/|H_\mathrm{V}|$ seems not to be indicative for the magnitude or type of an upcoming flare (confined or eruptive). Interpreted in context with earlier observational studies, our findings furthermore support that the total relative helicity normalized to the magnetic flux at the NLFF model's lower boundary, $H_\mathrm{V}/\phi^2$, represents no indicator for the eruptivity.

Context. In September 2017 the largest X-class flare of Solar Cycle 24 occurred from the most active region (AR) of this cycle, AR 12673. The AR attracted much interest because of its unique morphological and evolution characteristics. Among the parameters examined in the AR was magnetic helicity, but either only approximately, and/or intermittently. Aims. This work is interested in studying the evolution of the relative magnetic helicity and of the two components of its decomposition, the non-potential, and the volume-threading one, in the time interval around the highest activity of AR 12673. Special emphasis is given on the study of the ratio of the non-potential to total helicity, that was recently proposed as an indicator of ARs eruptivity. Methods. For these, we first approximate the coronal magnetic field of the AR with two different optimization-based extrapolation procedures, and choose the one that produces the most reliable helicity value at each instant. Moreover, in one of these methods, we weight the optimization by the uncertainty estimates derived from the Helioseismic and Magnetic Imager (HMI) instrument, for the first time. We then follow an accurate method to compute all quantities of interest. Results. The first observational determination of the evolution of the non-potential to total helicity ratio seems to confirm the quality it has in indicating eruptivity. This ratio increases before the major flares of AR 12673, and afterwards it relaxes to smaller values. Additionally, the evolution patterns of the various helicity, and energy budgets of AR 12673 are discussed and compared with other works.

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.

Context. Magnetic helicity is an important quantity in studies of magnetized plasmas as it provides a measure of the geometrical complexity of the magnetic field in a given volume. A more detailed description of the spatial distribution of magnetic helicity is given by the field line helicity that expresses the amount of helicity associated to individual field lines, rather than in the full analysed volume. Aims. Magnetic helicity is not a gauge-invariant quantity in general, unless it is computed with respect to a reference field, yielding the so called relative magnetic helicity. The field line helicity corresponding to the relative magnetic helicity has only been examined under specific conditions so far. This work aims to define the field line helicity corresponding to relative magnetic helicity in the most general way. In addition to its general form, we provide the expression for the relative magnetic field line helicity in a few commonly used gauges, and reproduce known results as a limit of our general formulation. Methods. By starting from the definition of relative magnetic helicity, we derive the corresponding field line helicity, and we note the assumptions it is based on. Results. We check that the developed quantity reproduces relative magnetic helicity by using three different numerical simulations. For these cases we also show the morphology of field line helicity in the volume, and on the photospheric plane. As an application to solar situations, we compare the morphology of field line helicity on the photosphere with that of the connectivity-based helicity flux density in two reconstructions of an active region's magnetic field. We discuss how the derived relative magnetic field line helicity has a wide range of applications, notably in solar physics and magnetic reconnection studies.

The relative magnetic helicity is a quantity that is often used to describe the level of entanglement of non-isolated magnetic fields, such as the magnetic field of solar active regions.The aim of this paper is to investigate how different kinds of photospheric boundary flows accumulate relative magnetic helicity in the corona and if and how-well magnetic helicity related quantities identify the onset of an eruption. We use a series of three-dimensional, parametric magnetohydrodynamic simulations of the formation and eruption of magnetic flux ropes. All the simulations are performed on the same grid, using the same parameters, but they are characterized by different driving photospheric flows, i.e., shearing, convergence, stretching, peripheral- and central- dispersion flows. For each of the simulations, the instant of the onset of the eruption is carefully identified by using a series of relaxation runs. We find that magnetic energy and total relative helicity are mostly injected when shearing flows are applied at the boundary, while the magnetic energy and helicity associated with the coronal electric currents increase regardless of the kind of photospheric flows. We also find that, at the onset of the eruptions, the ratio between the non-potential magnetic helicity and the total relative magnetic helicity has the same value for all the simulations, suggesting the existence of a threshold in this quantity. Such threshold is not observed for other quantities as, for example, those related to the magnetic energy.

In the solar corona, magnetic helicity slowly and continuously accumulates in response to plasma flows tangential to the photosphere and magnetic flux emergence through it. Analyzing this transfer of magnetic helicity is key for identifying its role in the dynamics of active regions (ARs). The connectivity-based helicity flux density method was recently developed for studying the 2D and 3D transfer of magnetic helicity in ARs. The method takes into account the 3D nature of magnetic helicity by explicitly using knowledge of the magnetic field connectivity, which allows it to faithfully track the photospheric flux of magnetic helicity. Because the magnetic field is not measured in the solar corona, modeled 3D solutions obtained from force-free magnetic field extrapolations must be used to derive the magnetic connectivity. Different extrapolation methods can lead to markedly different 3D magnetic field connectivities, thus questioning the reliability of the connectivity-based approach in observational applications. We address these concerns by applying this method to the isolated and internally complex AR 11158 with different magnetic field extrapolation models. We show that the connectivity-based calculations are robust to different extrapolation methods, in particular with regards to identifying regions of opposite magnetic helicity flux. We conclude that the connectivity-based approach can be reliably used in observational analyses and is a promising tool for studying the transfer of magnetic helicity in ARs and relate it to their flaring activity.

Solar flares and coronal mass ejections (CMEs) are among the most energetic events in the solar system, impacting the near-Earth environment. Flare productivity is empirically known to be correlated with the size and complexity of active regions. Several indicators, based on magnetic-field data from active regions, have been tested for flare forecasting in recent years. None of these indicators, or combinations thereof, have yet demonstrated an unambiguous eruption or flare criterion. Furthermore, numerical simulations have been only barely used to test the predictability of these parameters. In this context, we used the 3D parametric MHD numerical simulations of the self-consistent formation of the flux emergence of a twisted flux tube, inducing the formation of stable and unstable magnetic flux ropes of Leake (2013, 2014). We use these numerical simulations to investigate the eruptive signatures observable in various magnetic scalar parameters and provide highlights on data analysis processing. Time series of 2D photospheric-like magnetograms are used from parametric simulations of stable and unstable flux emergence, to compute a list of about 100 different indicators. This list includes parameters previously used for operational forecasting, physical parameters used for the first time, as well as new quantities specifically developed for this purpose. Our results indicate that only parameters measuring the total non-potentiality of active regions associated with magnetic inversion line properties, such as the Falconer parameters $L_{ss}$, $WL_{ss}$, $L_{sg}$ and $WL_{sg}$, as well as the new current integral $WL_{sc}$ and length $L_{sc}$ parameters, present a significant ability to distinguish the eruptive cases of the model from the non-eruptive cases, possibly indicating that they are promising flare and eruption predictors.

We study the writhe, twist and magnetic helicity of different magnetic flux ropes, based on models of the solar coronal magnetic field structure. These include an analytical force-free Titov--D\'emoulin equilibrium solution, non force-free magnetohydrodynamic simulations, and nonlinear force-free magnetic field models. The geometrical boundary of the magnetic flux rope is determined by the quasi-separatrix layer and the bottom surface, and the axis curve of the flux rope is determined by its overall orientation. The twist is computed by the Berger--Prior formula that is suitable for arbitrary geometry and both force-free and non-force-free models. The magnetic helicity is estimated by the twist multiplied by the square of the axial magnetic flux. We compare the obtained values with those derived by a finite volume helicity estimation method. We find that the magnetic helicity obtained with the twist method agrees with the helicity carried by the purely current-carrying part of the field within uncertainties for most test cases. It is also found that the current-carrying part of the model field is relatively significant at the very location of the magnetic flux rope. This qualitatively explains the agreement between the magnetic helicity computed by the twist method and the helicity contributed purely by the current-carrying magnetic field.

We present a detailed study of a confined circular flare dynamics associated with 3 UV late phases in order to understand more precisely which topological elements are present and how they constrain the dynamics of the flare. We perform a non-linear force free field extrapolation of the confined flare observed with the HMI and AIA instruments onboard SDO. From the 3D magnetic field we compute the squashing factor and we analyse its distribution. Conjointly, we analyse the AIA EUV light curves and images in order to identify the post-flare loops, their temporal and thermal evolution. By combining both analysis we are able to propose a detailed scenario that explains the dynamics of the flare. Our topological analysis shows that in addition to a null-point topology with the fan separatrix, the spine lines and its surrounding Quasi-Separatix Layers halo (typical for a circular flare), a flux rope and its hyperbolic flux tube (HFT) are enclosed below the null. By comparing the magnetic field topology and the EUV post-flare loops we obtain an almost perfect match 1) between the footpoints of the separatrices and the EUV 1600~\AA{} ribbons and 2) between the HFT's field line footpoints and bright spots observed inside the circular ribbons. We showed, for the first time in a confined flare, that magnetic reconnection occured initially at the HFT, below the flux rope. Reconnection at the null point between the flux rope and the overlying field is only initiated in a second phase. In addition, we showed that the EUV late phase observed after the main flare episode are caused by the cooling loops of different length which have all reconnected at the null point during the impulsive phase.

The discovery of clear criteria that can deterministically describe the eruptive state of a solar active region would lead to major improvements on space weather predictions. Using series of numerical simulations of the emergence of a magnetic flux rope in a magnetized coronal, leading either to eruptions or to stable configurations, we test several global scalar quantities for the ability to discriminate between the eruptive and the non-eruptive simulations. From the magnetic field generated by the three-dimensional magnetohydrodynamical simulations, we compute and analyse the evolution of the magnetic flux, of the magnetic energy and its decomposition into potential and free energies, and of the relative magnetic helicity and its decomposition. Unlike the magnetic flux and magnetic energies, magnetic helicities are able to markedly distinguish the eruptive from the non-eruptive simulations. We find that the ratio of the magnetic helicity of the current-carrying magnetic field to the total relative helicity presents the highest values for the eruptive simulations, in the pre-eruptive phase only. We observe that the eruptive simulations do not possess the highest value of total magnetic helicity. In the framework of our numerical study, the magnetic energies and the total relative helicity do not correspond to good eruptivity proxies. Our study highlights that the ratio of magnetic helicities diagnoses very clearly the eruptive potential of our parametric simulations. Our study shows that magnetic-helicity-based quantities may be very efficient for the prediction of solar eruptions.

By using a new method of forced-field extrapolation, we study the emerging flux region AR 11850 observed by the Interface Region Imaging Spectrograph (IRIS) and Solar Dynamical Observatory (SDO). Our results suggest that the bright points (BPs) in this emerging region have responses in lines formed from the upper photosphere to the transition region, with a relatively similar morphology. They have an oscillation of several minutes according to the Atmospheric Imaging Assembly (AIA) data at 1600 and 1700 A . The ratio between the BP intensities measured in 1600 A and 1700 A filtergrams reveals that these BPs are heated differently. Our analysis of the Helioseismic and Magnetic Imager (HMI) vector magnetic field and the corresponding topology in AR11850 indicates that the BPs are located at the polarity inversion line (PIL) and most of them related with magnetic reconnection or cancelation. The heating of the BPs might be different due to different magnetic topology. We find that the heating due to the magnetic cancelation would be stronger than the case of bald patch reconnection. The plasma density rather than the magnetic field strength could play a dominant role in this process. Based on physical conditions in the lower atmosphere, our forced-field extrapolation shows consistent results between the bright arcades visible in slit-jaw image (SJI) 1400 A and the extrapolated field lines that pass through the bald patches. It provides a reliable observational evidence for testing the mechanism of magnetic reconnection for the BPs and arcades in emerging flux region, as proposed in simulation works.

Three homologous C-class flares and one last M-class flare were observed by both the Solar Dynamics Observatory (SDO) and the Hinode EUV Imaging Spectrometer (EIS) in the AR 11429 on March 9, 2012. All the recurrent flares occurred within a short interval of time (less than 4 hours), showed very similar plasma morphology and were all confined, until the last one when a large-scale eruption occurred. The C-class flares are characterized by the appearance, at approximatively the same locations, of two bright and compact footpoint sources of $\approx$~3--10~MK evaporating plasma, and a semi-circular ribbon. During all the flares, the continuous brightening of a spine-like hot plasma ($\approx$~10~MK) structure is also observed. Spectroscopic observations with Hinode/EIS are used to measure and compare the blueshift velocities in the \fexxiii\ emission line and the electron number density at the flare footpoints for each flare. Similar velocities, of the order of 150--200~km~s$^{-1}$, are observed during the C2.0 and C4.7 confined flares, in agreement with the values reported by other authors in the study of the last M1.8 class flare. On the other hand, lower electron number densities and temperatures tend to be observed in flares with lower peak soft X-ray flux.In order to investigate the homologous nature of the flares, we performed a Non-Linear Force-Free Field (NLFFF) extrapolation of the 3D magnetic field configuration in the corona. The NLFFF extrapolation and the Quasi-Separatrix Layers (QSLs) provide the magnetic field context which explains the location of the kernels, spine-like and semi-circular brightenings observed in the (non-eruptive) flares. Given the absence of a coronal null point, we argue that the homologous flares were all generated by the continuous recurrence of bald patch reconnection.

Magnetic helicity is a conserved quantity of ideal magneto-hydrodynamics characterized by an inverse turbulent cascade. Accordingly, it is often invoked as one of the basic physical quantities driving the generation and structuring of magnetic fields in a variety of astrophysical and laboratory plasmas. We provide here the first systematic comparison of six existing methods for the estimation of the helicity of magnetic fields known in a finite volume. All such methods are reviewed, benchmarked, and compared with each other, and specifically tested for accuracy and sensitivity to errors. To that purpose, we consider four groups of numerical tests, ranging from solutions of the three-dimensional, force-free equilibrium, to magneto-hydrodynamical numerical simulations. Almost all methods are found to produce the same value of magnetic helicity within few percent in all tests. In the more solar-relevant and realistic of the tests employed here, the simulation of an eruptive flux rope, the spread in the computed values obtained by all but one method is only 3%, indicating the reliability and mutual consistency of such methods in appropriate parameter ranges. However, methods show differences in the sensitivity to numerical resolution and to errors in the solenoidal property of the input fields. In addition to finite volume methods, we also briefly discuss a method that estimates helicity from the field lines' twist, and one that exploits the field's value at one boundary and a coronal minimal connectivity instead of a pre-defined three-dimensional magnetic-field solution.

Context: A subclass of broad EUV and X-ray jets, called blowout jets, have become a topic of research since they could be the link between standard collimated jets and CMEs.} Aim: Our aim is to understand the origin of a series of broad jets, some accompanied by flares and associated with narrow and jet-like CMEs. Methods: We analyze observations of a series of recurrent broad jets observed in AR 10484 on 21-24 October 2003. In particular, one of them occurred simultaneously with an M2.4 flare on 23 October at 02:41 UT (SOLA2003-10-23). Both events were observed by ARIES H-alpha Solar Tower-Telescope, TRACE, SOHO, and RHESSI instruments. The flare was very impulsive and followed by a narrow CME. A local force-free model of AR 10484 is the basis to compute its topology. We find bald patches (BPs) at the flare site. This BP topology is present for at least two days before. Large-scale field lines, associated with the BPs, represent open loops. This is confirmed by a global PFSS model. Following the brightest leading edge of the H-alpha and EUV jet emission, we can temporarily associate it with a narrow CME. Results: Considering their characteristics, the observed broad jets appear to be of the blowout class. As the most plausible scenario, we propose that magnetic reconnection could occur at the BP separatrices forced by the destabilization of a continuously reformed flux rope underlying them. The reconnection process could bring the cool flux-rope material into the reconnected open field lines driving the series of recurrent blowout jets and accompanying CMEs. Conclusions: Based on a model of the coronal field, we compute the AR 10484 topology at the location where flaring and blowout jets occurred from 21 to 24 October 2003. This topology can consistently explain the origin of these events.

Jets are dynamic, impulsive, well-collimated plasma events that develop at many different scales and in different layers of the solar atmosphere. Jets are believed to be induced by magnetic reconnection, a process central to many astrophysical phenomena. Within the solar atmosphere, jet-like events develop in many different environments, e.g., in the vicinity of active regions as well as in coronal holes, and at various scales, from small photospheric spicules to large coronal jets. In all these events, signatures of helical structure and/or twisting/rotating motions are regularly observed. The present study aims to establish that a single model can generally reproduce the observed properties of these jet-like events. In this study, using our state-of-the-art numerical solver ARMS, we present a parametric study of a numerical tridimensional magnetohydrodynamic (MHD) model of solar jet-like events. Within the MHD paradigm, we study the impact of varying the atmospheric plasma $\beta$ on the generation and properties of solar-like jets. The parametric study validates our model of jets for plasma $\beta$ ranging from $10^{-3}$ to $1$, typical of the different layers and magnetic environments of the solar atmosphere. Our model of jets can robustly explain the generation of helical solar jet-like events at various $\beta \le 1$. This study introduces the new result that the plasma $\beta$ modifies the morphology of the helical jet, explaining the different observed shapes of jets at different scales and in different layers of the solar atmosphere. Our results allow us to understand the energisation, triggering, and driving processes of jet-like events. Our model allows us to make predictions of the impulsiveness and energetics of jets as determined by the surrounding environment, as well as the morphological properties of the resulting jets.

Coronal jets represent important manifestations of ubiquitous solar transients, which may be the source of significant mass and energy input to the upper solar atmosphere and the solar wind. While the energy involved in a jet-like event is smaller than that of "nominal" solar flares and coronal mass ejections (CMEs), jets share many common properties with these phenomena, in particular, the explosive magnetically driven dynamics. Studies of jets could, therefore, provide critical insight for understanding the larger, more complex drivers of the solar activity. On the other side of the size-spectrum, the study of jets could also supply important clues on the physics of transients close or at the limit of the current spatial resolution such as spicules. Furthermore, jet phenomena may hint to basic process for heating the corona and accelerating the solar wind; consequently their study gives us the opportunity to attack a broad range of solar-heliospheric problems.

Coronal-hole jets occur ubiquitously in solar coronal holes, at EUV and X-ray bright points associated with intrusions of minority magnetic polarity. The embedded-bipole model for these jets posits that they are driven by explosive, fast reconnection between the stressed closed field of the embedded bipole and the open field of the surrounding coronal hole. Previous numerical studies in Cartesian geometry, assuming uniform ambient magnetic field and plasma while neglecting gravity and solar wind, demonstrated that the model is robust and can produce jet-like events in simple configurations. We have extended these investigations by including spherical geometry, gravity, and solar wind in a nonuniform, coronal hole-like ambient atmosphere. Our simulations confirm that the jet is initiated by the onset of a kink-like instability of the internal closed field, which induces a burst of reconnection between the closed and external open field, launching a helical jet. Our new results demonstrate that the jet propagation is sustained through the outer corona, in the form of a traveling nonlinear Alfven wave front trailed by slower-moving plasma density enhancements that are compressed and accelerated by the wave. This finding agrees well with observations of white-light coronal-hole jets, and can explain microstreams and torsional Alfven waves detected in situ in the solar wind. We also use our numerical results to deduce scaling relationships between properties of the coronal source region and the characteristics of the resulting jet, which can be tested against observations.

The standard model for eruptive flares has in the past few years been extended to 3D. It predicts typical J-shaped photospheric footprints of the coronal current layer, forming at similar locations as the Quasi-Separatrix Layers (QSLs). Such a morphology is also found for flare ribbons observed in the EUV band, as well as in non-linear force-free field (NLFFF) magnetic field extrapolations and models. We study the evolution of the photospheric traces of the current density and flare ribbons, both obtained with the SDO instruments. We investigate the photospheric current evolution during the 6 September 2011 X-class flare (SOL2011-09-06T22:20) from observational data of the magnetic field obtained with HMI. This evolution is compared with that of the flare ribbons observed in the EUV filters of the AIA. We also compare the observed electric current density and the flare ribbon morphology with that of the QSLs computed from the flux rope insertion method/NLFFF model. The NLFFF model shows the presence of a fan-spine configuration of overlying field lines, due to the presence of a parasitic polarity, embedding an elongated flux rope that appears in the observations as two parts of a filament. The QSLs, evolved via a magnetofrictional method, also show similar morphology and evolution as both the current ribbons and the EUV flare ribbons obtained at several times during the flare. For the first time, we propose a combined analysis of the photospheric traces of an eruptive flare, in a complex topology, with direct measurements of electric currents and QSLs from observational data and a magnetic field model. The results, obtained by two different and independent approaches, 1) confirm previous results of current increase during the impulsive phase of the flare, 2) show how NLFFF models can capture the essential physical signatures of flares even in a complex magnetic field topology.

We studied the magnetic topology of active region 12158 on 2014 September 10 and compared it with the observations before and early in the flare which begins at 17:21 UT (SOL2014-09-10T17:45:00). Our results show that the sigmoidal structure and flare ribbons of this active region observed by SDO/AIA can be well reproduced from a Grad-Rubin non linear force free field extrapolation method. Various inverse-S and -J shaped magnetic field lines, that surround a coronal flux rope, coincide with the sigmoid as observed in different extreme ultraviolet wavelengths, including its multi-threaded curved ends. Also, the observed distribution of surface currents in the magnetic polarity where it was not prescribed is well reproduced. This validates our numerical implementation and set-up of the Grad-Rubin method. The modeled double inverse-J shaped Quasi-Separatrix Layer (QSL) footprints match the observed flare ribbons during the rising phase of the flare, including their hooked parts. The spiral-like shape of the latter may be related to a complex pre-eruptive flux rope with more than one turn of twist, as obtained in the model. These ribbon-associated flux-rope QSL-footprints are consistent with the new standard flare model in 3D, with the presence of a hyperbolic flux tube located below an inverse tear drop shaped coronal QSL. This is a new step forward forecasting the locations of reconnection and ribbons in solar flares, and the geometrical properties of eruptive flux ropes.

The dynamic properties of flare ribbons and the often associated filament eruptions can provide crucial information on the flaring coronal magnetic field. This Letter analyzes the GOES-class X1.0 flare on 2014 March 29 (SOL2014-03-29T17:48), in which we found an asymmetric eruption of a sigmoidal filament and an ensuing circular flare ribbon. Initially both EUV images and a preflare nonlinear force-free field model show that the filament is embedded in magnetic fields with a fan-spine-like structure. In the first phase, which is defined by a weak but still increasing X-ray emission, the western portion of the sigmoidal filament arches upward and then remains quasi-static for about five minutes. The western fan-like and the outer spine-like fields display an ascending motion, and several associated ribbons begin to brighten. Also found is a bright EUV flow that streams down along the eastern fan-like field. In the second phase that includes the main peak of hard X-ray (HXR) emission, the filament erupts, leaving behind two major HXR sources formed around its central dip portion and a circular ribbon brightened sequentially. The expanding western fan-like field interacts intensively with the outer spine-like field, as clearly seen in running difference EUV images. We discuss these observations in favor of a scenario where the asymmetric eruption of the sigmoidal filament is initiated due to an MHD instability and further facilitated by reconnection at a quasi-null in corona; the latter is in turn enhanced by the filament eruption and subsequently produces the circular flare ribbon.

There is a recurring question in solar physics about whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Another source of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net vs. neutralized currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi-static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: (1) produce both {\it direct} and {\it return} currents, (2) induce very weak compression currents - not observed in 2.5D - in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the absence of magnetic shear at the photospheric polarity inversion line - a special condition rarely observed. We conclude that, as magnetic flux emergence, photospheric flows can build up net currents in the solar atmosphere, in agreement with recent observations. These results thus provide support for eruption models based on pre-eruption magnetic fields possessing a net coronal current.

Magnetic helicity has the remarkable property of being a conserved quantity of ideal magnetohydrodynamics (MHD). Therefore, it could be used as an effective tracer of the magnetic field evolution of magnetized plasmas. Theoretical estimations indicate that magnetic helicity is also essentially conserved with non-ideal MHD processes, e.g. magnetic reconnection. This conjecture has however been barely tested, either experimentally or numerically. Thanks to recent advances in magnetic helicity estimation methods, it is now possible to test numerically its dissipation level in general three-dimensional datasets. We first revisit the general formulation of the temporal variation of relative magnetic helicity on a fully bounded volume when no hypothesis on the gauge are made. We introduce a method to precisely estimate its dissipation independently of the type of non-ideal MHD processes occurring. In a solar-like eruptive event simulation, using different gauges, we compare its estimation in a finite volume with its time-integrated flux through the boundaries, hence testing the conservation and dissipation of helicity. We provide an upper bound of the real dissipation of magnetic helicity: It is quasi-null during the quasi-ideal MHD phase. Even when magnetic reconnection is acting the relative dissipation of magnetic helicity is also very small (<2.2%), in particular compared to the relative dissipation of magnetic energy (>30 times larger). We finally illustrate how the helicity-flux terms involving velocity components are gauge dependent, hence limiting their physical meaning.

In this paper we present a topological magnetic field investigation of seven two-ribbon flares in sigmoidal active regions observed with Hinode, STEREO, and SDO. We first derive the 3D coronal magnetic field structure of all regions using marginally unstable 3D coronal magnetic field models created with the flux rope insertion method. The unstable models have been shown to be a good model of the flaring magnetic field configurations. Regions are selected based on their pre-flare configurations along with the appearance and observational coverage of flare ribbons, and the model is constrained using pre-flare features observed in extreme ultraviolet and X-ray passbands. We perform a topology analysis of the models by computing the squashing factor, Q, in order to determine the locations of prominent quasi-separatrix layers (QSLs). QSLs from these maps are compared to flare ribbons at their full extents. We show that in all cases the straight segments of the two J-shaped ribbons are matched very well by the flux-rope-related QSLs, and the matches to the hooked segments are less consistent but still good for most cases. In addition, we show that these QSLs overlay ridges in the electric current density maps. This study is the largest sample of regions with QSLs derived from 3D coronal magnetic field models, and it shows that the magnetofrictional modeling technique that we employ gives a very good representation of flaring regions, with the power to predict flare ribbon locations in the event of a flare following the time of the model.

Coronal mass ejections (CMEs) erupt and expand in a magnetically structured solar corona. Various indirect observational pieces of evidence have shown that the magnetic field of CMEs reconnects with surrounding magnetic fields, forming, e.g., dimming regions distant from the CME source regions. Analyzing Solar Dynamics Observatory (SDO) observations of the eruption from AR 11226 on 2011 June 7, we present the first direct evidence of coronal magnetic reconnection between the fields of two adjacent ARs during a CME. The observations are presented jointly with a data-constrained numerical simulation, demonstrating the formation/intensification of current sheets along a hyperbolic flux tube (HFT) at the interface between the CME and the neighbouring AR 11227. Reconnection resulted in the formation of new magnetic connections between the erupting magnetic structure from AR 11226 and the neighboring active region AR 11227 about 200 Mm from the eruption site. The onset of reconnection first becomes apparent in the SDO/AIA images when filament plasma, originally contained within the erupting flux rope, is re-directed towards remote areas in AR 11227, tracing the change of large-scale magnetic connectivity. The location of the coronal reconnection region becomes bright and directly observable at SDO/AIA wavelengths, owing to the presence of down-flowing cool, dense (10^{10} cm^{-3}) filament plasma in its vicinity. The high-density plasma around the reconnection region is heated to coronal temperatures, presumably by slow-mode shocks and Coulomb collisions. These results provide the first direct observational evidence that CMEs reconnect with surrounding magnetic structures, leading to a large-scale re-configuration of the coronal magnetic field.

We present for the first time the evolution of the photospheric electric currents during an eruptive X-class flare, accurately predicted by the standard 3D flare model. We analyze this evolution for the February 15, 2011 flare using HMI/SDO magnetic observations and find that localized currents in \J-shaped ribbons increase to double their pre-flare intensity. Our 3D flare model, developed with the OHM code, suggests that these current ribbons, which develop at the location of EUV brightenings seen with AIA imagery, are driven by the collapse of the flare's coronal current layer. These findings of increased currents restricted in localized ribbons are consistent with the overall free energy decrease during a flare, and the shape of these ribbons also give an indication on how much twisted the erupting flux rope is. Finally, this study further enhances the close correspondence obtained between the theoretical predictions of the standard 3D model and flare observations indicating that the main key physical elements are incorporated in the model.

We studied the temporal evolution of the magnetic topology of the active region (AR) 11158 based on the reconstructed three-dimensional magnetic fields in the corona. The \nlfff\ extrapolation method was applied to the 12 minutes cadence data obtained with the \hmi\ (HMI) onboard the \sdo\ (SDO) during five days. By calculating the squashing degree factor Q in the volume, the derived quasi-separatrix layers (QSLs) show that this AR has an overall topology, resulting from a magnetic quadrupole, including an hyperbolic flux tube (HFT) configuration which is relatively stable at the time scale of the flare ($\sim 1-2$ hours). A strong QSL, which corresponds to some highly sheared arcades that might be related to the formation of a flux rope, is prominent just before the M6.6 and X2.2 flares, respectively. These facts indicate the close relationship between the strong QSL and the high flare productivity of AR 11158. In addition, with a close inspection of the topology, we found a small-scale HFT which has an inverse tear-drop structure above the aforementioned QSL before the X2.2 flare. It indicates the existence of magnetic flux rope at this place. Even though a global configuration (HFT) is recognized in this AR, it turns out that the large-scale HFT only plays a secondary role during the eruption. In final, we dismiss a trigger based on the breakout model and highlight the central role of the flux rope in the related eruption.

To study the build up of a magnetic flux rope before a major flare and coronal mass ejection (CME), we compute the magnetic helicity injection, twist accumulation, and the topology structure of the three dimensional magnetic field, which is derived by the nonlinear force-free field model. The Extreme-ultraviolet Imaging Telescope on board the Solar and Heliospheric Observatory observed a series of confined flares without any CME before a major flare with a CME at 23:02 UT on 2005 January 15 in active region NOAA 10720. We derive the vector velocity at eight time points from 18:27 UT to 22:20 UT with the differential affine velocity estimator for vector magnetic fields, which were observed by the Digital Vector Magnetograph at Big Bear Solar Observatory. The injected magnetic helicity is computed with the vector magnetic and velocity fields. The helicity injection rate was (-16.47 \pm 3.52) \times 10^{40} Mx^2/hr. We find that only about 1.8% of the injected magnetic helicity became finally the internal helicity of the magnetic flux rope, whose twist increasing rate was -0.18 \pm 0.08 Turns/hr. The quasi-separatrix layers (QSLs) of the three dimensional magnetic field are computed by evaluating the squashing degree, Q. We find that the flux rope was wrapped by QSLs with large Q values, where the magnetic reconnection induced by the continuously injected magnetic helicity further produced the confined flares. We suggest that the flux rope was built up and heated by the magnetic reconnection in the QSLs.

Measuring the magnetic helicity distribution in the solar corona can help in understanding the trigger of solar eruptive events because magnetic helicity is believed to play a key role in solar activity due to its conservation property. A new method for computing the photospheric distribution of the helicity flux was recently developed. This method takes into account the magnetic field connectivity whereas previous methods were based on photospheric signatures only. This novel method maps the true injection of magnetic helicity in active regions. We applied this method for the first time to an observed active region, NOAA 11158, which was the source of intense flaring activity. We used high-resolution vector magnetograms from the SDO/HMI instrument to compute the photospheric flux transport velocities and to perform a nonlinear force-free magnetic field extrapolation. We determined and compared the magnetic helicity flux distribution using a purely photospheric as well as a connectivity-based method. While the new connectivity-based method confirms the mixed pattern of the helicity flux in NOAA 11158, it also reveals a different, and more correct, distribution of the helicity injection. This distribution can be important for explaining the likelihood of an eruption from the active region. The connectivity-based approach is a robust method for computing the magnetic helicity flux, which can be used to study the link between magnetic helicity and eruptivity of observed active regions.

Magnetic helicity quantifies how globally sheared and/or twisted is the magnetic field in a volume. This quantity is believed to play a key role in solar activity due to its conservation property. Helicity is continuously injected into the corona during the evolution of active regions (ARs). To better understand and quantify the role of magnetic helicity in solar activity, the distribution of magnetic helicity flux in ARs needs to be studied. The helicity distribution can be computed from the temporal evolution of photospheric magnetograms of ARs such as the ones provided by SDO/HMI and Hinode/SOT. Most recent analyses of photospheric helicity flux derive an helicity flux density proxy based on the relative rotation rate of photospheric magnetic footpoints. Although this proxy allows a good estimate of the photospheric helicity flux, it is still not a true helicity flux density because it does not take into account the connectivity of the magnetic field lines. For the first time, we implement a helicity density which takes into account such connectivity. In order to use it for future observational studies, we test the method and its precision on several types of models involving different patterns of helicity injection. We also test it on more complex configurations - from magnetohydrodynamics (MHD) simulations - containing quasi-separatrix layers. We demonstrate that this connectivity-based helicity flux density proxy is the best to map the true distribution of photospheric helicity injection.

A standard model for eruptive flares aims at describing observational 3D features of the reconnecting coronal magnetic field. Extensions to the 2D model require the physical understanding of 3D reconnection processes at the origin of the magnetic configuration evolution. However, the properties of 3D reconnection without null point and separatrices still need to be analyzed. We focus on magnetic reconnection associated with the growth and evolution of a flux rope and associated flare loops during an eruptive flare. We aim at understanding the intrinsic characteristics of 3D reconnection in the presence of quasi-separatrix layers (QSLs), how QSL properties are related to the slip-running reconnection mode in general, and how this applies to eruptive flares in particular. We studied the slip-running reconnection of field lines in a magnetohydrodynamic simulation of an eruptive flare associated with a torus-unstable flux rope. Field lines associated with the flux rope and the flare loops undergo a continuous series of magnetic reconnection, which results in their super-Alfvenic slipping motion. The time profile of their slippage speed and the space distribution of the mapping norm are shown to be strongly correlated. We find that the motion speed is proportional to the mapping norm. Moreover, this slip-running motion becomes faster as the flux rope expands, since the 3D current layer evolves toward a current sheet, and QSLs to separatrices. The present analysis extends our understanding of the 3D slip-running reconnection regime. We identified a controlling parameter of the apparent velocity of field lines while they slip-reconnect, enabling the interpretation of the evolution of post flare loops. This work completes the standard model for flares and eruptions by giving its 3D properties.

For magnetically driven events, the magnetic energy of the system is the prime energy reservoir that fuels the dynamical evolution. In the solar context, the free energy is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. A trustworthy estimation of the magnetic energy is therefore needed in three-dimensional models of the solar atmosphere, eg in coronal fields reconstructions or numerical simulations. The expression of the energy of a system as the sum of its potential energy and its free energy (Thomson's theorem) is strictly valid when the magnetic field is exactly solenoidal. For numerical realizations on a discrete grid, this property may be only approximately fulfilled. We show that the imperfect solenoidality induces terms in the energy that can lead to misinterpreting the amount of free energy present in a magnetic configuration. We consider a decomposition of the energy in solenoidal and nonsolenoidal parts which allows the unambiguous estimation of the nonsolenoidal contribution to the energy. We apply this decomposition to six typical cases broadly used in solar physics. We quantify to what extent the Thomson theorem is not satisfied when approximately solenoidal fields are used. The quantified errors on energy vary from negligible to significant errors, depending on the extent of the nonsolenoidal component. We identify the main source of errors and analyze the implications of adding a variable amount of divergence to various solenoidal fields. Finally, we present pathological unphysical situations where the estimated free energy would appear to be negative, as found in some previous works, and we identify the source of this error to be the presence of a finite divergence. We provide a method of quantifying the effect of a finite divergence in numerical fields, together with detailed diagnostics of its sources.

Solar flares strongly affect the Sun's atmosphere as well as the Earth's environment. Quantifying the maximum possible energy of solar flares of the present-day Sun, if any, is thus a key question in heliophysics. The largest solar flares observed over the past few decades have reached energies of a few times 10^{32} ergs, possibly up to 10^{33} ergs. Flares in active Sun-like stars reach up to about 10^{36} ergs. In the absence of direct observations of solar flares within this range, complementary methods of investigation are needed. Using historical reports for solar active region, we scaled to observed solar values a realistic dimensionless 3D MHD simulation for eruptive flares, which originate from a highly sheared bipole. This enabled us to calculate the magnetic fluxes and flare energies in the model in a wide paramater space. Firstly, commonly observed solar conditions lead to modeled magnetic fluxes and flare energies that are comparable to those estimated from observations. Secondly, we evaluate from observations that 30% of the area of sunspot groups are typically involved in flares. This is related to the strong fragmentation of these groups, which naturally results from sub-photospheric convection. When the model is scaled to 30% of the area of the largest sunspot group ever reported, with its peak magnetic field being set to the strongest value ever measured in a sunspot, it produces a flare with a maximum energy of ~ 6x10^{33} ergs. The results of the model suggest that the Sun is able to produce flares up to about six times as energetic in total solar irradiance fluence as the strongest directly-observed flare from Nov 4, 2003. Sunspot groups larger than historically reported would yield superflares for spot pairs that would exceed tens of degrees in extent. We thus conjecture that superflare-productive Sun-like stars should have a much stronger dynamo than in the Sun.

We investigate the X-ray and UV emission detected by RHESSI and TRACE in the context of a solar flare on the 16th November 2002 with the goal of better understanding the evolution of the flare. We analysed the characteristics of the X-ray emission in the 12-25 and 25-50 keV energy range while we looked at the UV emission at 1600 {\AA}. The flare appears to have two distinct phases of emission separated by a 25-second time delay, with the first phase being energetically more important. We found good temporal and spatial agreement between the 25-50 keV X-rays and the most intense areas of the 1600 {\AA} UV emission. We also observed an extended 100-arcsecond < 25 keV source that appears coronal in nature and connects two separated UV ribbons later in the flare. Using the observational properties in X-ray and UV wavelengths, we propose two explanations for the flare evolution in relation to the spine/fan magnetic field topology and the accelerated electrons. We find that a combination of quasi separatrix layer reconnection and null-point reconnection is required to account for the observed properties of the X-ray and UV emission.

We present a new model to explain how particles (solar energetic particles; SEPs), accelerated at a reconnection site that is not magnetically connected to the Earth, could eventually propagate along the well-connected open flux tube. Our model is based on the results of a low-beta resistive magnetohydrodynamics simulation of a three-dimensional line-tied and initially current-free bipole, that is embedded in a non-uniform open potential field. The topology of this configuration is that of an asymmetric coronal null-point, with a closed fan surface and an open outer spine. When driven by slow photospheric shearing motions, field lines, initially fully anchored below the fan dome, reconnect at the null point, and jump to the open magnetic domain. This is the standard interchange mode as sketched and calculated in 2D. The key result in 3D is that, reconnected open field lines located in the vicinity of the outer spine, keep reconnecting continuously, across an open quasi-separatrix layer, as previously identified for non-open-null-point reconnection. The apparent slipping motion of these field lines leads to form an extended narrow magnetic flux tube at high altitude. Because of the slip-running reconnection, we conjecture that if energetic particles would be traveling through, or be accelerated inside, the diffusion region, they would be successively injected along continuously reconnecting field lines that are connected farther and farther from the spine. At the scale of the full Sun, owing to the super-radial expansion of field lines below 3 solar radii, such energetic particles could easily be injected in field lines slipping over significant distances, and could eventually reach the distant flux tube that is well-connected to the Earth.

We present and interpret observations of two morphologically homologous flares that occurred in active region (AR) NOAA 10501 on 20 November 2003. Both flares displayed four homologous H-alpha ribbons and were both accompanied by coronal mass ejections (CMEs). The central flare ribbons were located at the site of an emerging bipole in the center of the active region. The negative polarity of this bipole fragmented in two main pieces, one rotating around the positive polarity by ~ 110 deg within 32 hours. We model the coronal magnetic field and compute its topology, using as boundary condition the magnetogram closest in time to each flare. In particular, we calculate the location of quasiseparatrix layers (QSLs) in order to understand the connectivity between the flare ribbons. Though several polarities were present in AR 10501, the global magnetic field topology corresponds to a quadrupolar magnetic field distribution without magnetic null points. For both flares, the photospheric traces of QSLs are similar and match well the locations of the four H-alpha ribbons. This globally unchanged topology and the continuous shearing by the rotating bipole are two key factors responsible for the flare homology. However, our analyses also indicate that different magnetic connectivity domains of the quadrupolar configuration become unstable during each flare, so that magnetic reconnection proceeds differently in both events.

The geoeffective magnetic cloud (MC) of 20 November 2003, has been associated to the 18 November 2003, solar active events in previous studies. In some of these, it was estimated that the magnetic helicity carried by the MC had a positive sign, as well as its solar source, active region (AR) NOAA 10501. In this paper we show that the large-scale magnetic field of AR 10501 had a negative helicity sign. Since coronal mass ejections (CMEs) are one of the means by which the Sun ejects magnetic helicity excess into the interplanetary space, the signs of magnetic helicity in the AR and MC should agree. Therefore, this finding contradicts what is expected from magnetic helicity conservation. However, using for the first time correct helicity density maps to determine the spatial distribution of magnetic helicity injection, we show the existence of a localized flux of positive helicity in the southern part of AR 10501. We conclude that positive helicity was ejected from this portion of the AR leading to the observed positive helicity MC.

We report on the first stereoscopic observations of polar coronal jets made by the EUVI/SECCHI imagers on board the twin STEREO spacecraft. The significantly separated viewpoints ($\sim$ 11$^\circ$) allowed us to infer the 3D dynamics and morphology of a well-defined EUV coronal jet for the first time. Triangulations of the jet's location in simultaneous image pairs led to the true 3D position and thereby its kinematics. Initially the jet ascends slowly at $\approx$10-20 $\mathrm{{km} {s}^{-1}}$ and then, after an apparent 'jump' takes place, it accelerates impulsively to velocities exceeding 300 $\mathrm{{km} {s}^{-1}}$ with accelerations exceeding the solar gravity. Helical structure is the most important geometrical feature of the jet which shows evidence of untwisting. The jet structure appears strikingly different from each of the two STEREO viewpoints: face-on in the one viewpoint and edge-on in the other. This provides conclusive evidence that the observed helical structure is real and is not resulting from possible projection effects of single viewpoint observations. The clear demonstration of twisted structure in polar jets compares favorably with synthetic images from a recent MHD simulation of jets invoking magnetic untwisting as their driving mechanism. Therefore, the latter can be considered as a viable mechanism for the initiation of polar jets.