Search

Solar jets are collimated plasma flows moving along magnetic field lines and accelerated at low altitude following magnetic reconnection. Several of them originate from anemone-shaped low-lying arcades and the most impulsive ones tend to be relatively wider and display untwisting motions. We aim to establish typical behaviours and observational signatures in the low atmosphere that can occur in response to the coronal development of such impulsive jets. We analysed an observed solar jet associated with a circular flare ribbon, using high-resolution observations from SST coordinated with IRIS and SDO. We related specifically identified features with those developing in a generic 3D line-tied numerical simulation of reconnection driven jets, performed with the ARMS code. We identified three features in the SST observations: the formation of a hook along the circular ribbon, the gradual widening of the jet through the apparent displacement of its kinked edge towards, and not away from the presumed reconnection site, and the falling back of some of the jet plasma towards a footpoint offset from that of the jet itself. The 3D numerical simulation naturally accounts for these features which were not imposed a priori. Our analyses allow to interpret them in the context of the 3D geometry of the asymmetric swirled anemone loops and their sequences of reconnection with ambient coronal loops. Given the relatively-simple conditions in which the observed jet occurred, together with the generic nature of the simulation that comprised minimum assumptions, we predict that the specific features that we identified and interpreted are probably typical of every impulsive jet, Comment: 13 pages, 8 figures, Accepted for publication in Astronomy and Astrophysics

Magnetic free energy powers solar flares and coronal mass ejections (CMEs), and the buildup of magnetic helicity might play a role in the development of unstable structures that subsequently erupt. To better understand the roles of energy and helicity in large flares and eruptions, we have characterized the evolution of magnetic energy and helicity associated with 21 X-class flares from 2010 to 2017. Our sample includes both confined and eruptive events, with 6 and 15 in each category, respectively. Using HMI vector magnetic field observations from several hours before to several hours after each event, we employ (a) the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) to determine the photospheric fluxes of energy and helicity, and (b) non-linear force-free field (NLFFF) extrapolations to estimate the coronal content of energy and helicity in source-region fields. Using Superposed Epoch analysis (SPE), we find, on average: (1) decreases in both magnetic energy and helicity, in both photospheric fluxes and coronal content, that persist for a few hours after eruptions, but no clear changes, notably in relative helicity, for confined events; (2) significant increases in the twist of photospheric fields in eruptive events, with twist uncertainties too large in confined events to constrain twist changes (and lower overall twist in confined events); and (3) on longer time scales (event time +12 hours), replenishment of free magnetic energy and helicity content to near pre-event levels for eruptive events. For eruptive events, magnetic helicity and free energy in coronal models clearly decrease after flares, with the amounts of decrease proportional to each region's pre-flare content., Comment: 25 pages, 15 figures, 1 table, accepted for publication in the Astrophysical Journal

Spectropolarimetric reconstructions of the photospheric vector magnetic field are intrinsically limited by the so-called 180$^\circ$ ambiguity in the orientation of the transverse component. The successful launch and operation of Solar Orbiter has made the removal of the 180$^\circ$ ambiguity possible using solely observations obtained from two different vantage points. While the exploitation of such a possibility is straightforward in principle, it is less so in practice and it is therefore important to assess the accuracy and limitations, as a function of both the satellites orbits and measurement principles. In this work we present a stereoscopic disambiguation method (SDM) and discuss a thorough testing of its accuracy in applications to modeled active regions and quiet Sun observations. In a first series of tests, we employ magnetograms extracted from three different numerical simulations as test fields, and model observations of the magnetograms from different angles and distances. In these more idealized tests, the SDM is proven to to reach a 100% disambiguation accuracy when applied to moderately-to-well resolved fields. Even in the case of disambiguation of quiet Sun magnetograms with significant under-resolved scale, the SDM provides an accuracy between 82% and 98% depending on the field strength. The accuracy of the SDM is found to be mostly sensitive to the variable resolution of Solar Orbiter on its highly elliptic orbit, as well as to the intrinsic scale of the observed field. Finally, as a more realistic test, we consider magnetograms that are obtained using a radiative transfer inversion code and the SOPHISM instrument simulator applied to a 3D simulation of a pore, and present a preliminary discussion of the effect of the viewing angle on the observed field., Comment: 41 pages; 18 Figures, accepted for publication in Solar Physics

Relative magnetic helicity is conserved by magneto-hydrodynamic evolution even in the presence of moderate resistivity. For that reason, it is often invoked as the most relevant constraint to the dynamical evolution of plasmas in complex systems, such as solar and stellar dynamos, photospheric flux emergence, solar eruptions, and relaxation processes in laboratory plasmas. However, such studies often indirectly imply that relative magnetic helicity in a given spatial domain can be algebraically split into the helicity contributions of the composing subvolumes, i.e., that it is an additive quantity. A limited number of very specific applications have shown that this is not the case. Progress in understanding the non-additivity of relative magnetic helicity requires removal of restrictive assumptions in favour of a general formalism that can be used both in theoretical investigations as well as in numerical applications. We derive the analytical gauge-invariant expression for the partition of relative magnetic helicity between contiguous finite-volumes, without any assumptions on either the shape of the volumes and interface, or the employed gauge. The non-additivity of relative magnetic helicity in finite volumes is proven in the most general, gauge-invariant formalism, and verified numerically. More restrictive assumptions are adopted to derive known specific approximations, yielding a unified view of the additivity issue. As an example, the case of a flux rope embedded in a potential field shows that the non-additivity term in the partition equation is, in general, non-negligible. The relative helicity partition formula can be applied to numerical simulations to precisely quantify the effect of non-additivity on global helicity budgets of complex physical processes., Comment: 14 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Based on a decomposition of the magnetic field into potential and nonpotential components, magnetic energy and relative helicity can both also be decomposed into two quantities: potential and free energies, and volume-threading and current-carrying helicities. In this study, we perform a coupled analysis of their behaviors in a set of parametric 3D magnetohydrodynamic (MHD) simulations of solar-like eruptions. We present the general formulations for the time-varying components of energy and helicity in resistive MHD. We calculated them numerically with a specific gauge, and compared their behaviors in the numerical simulations, which differ from one another by their imposed boundary-driving motions. Thus, we investigated the impact of different active regions surface flows on the development of the energy and helicity-related quantities. Despite general similarities in their overall behaviors, helicities and energies display different evolutions that cannot be explained in a unique framework. While the energy fluxes are similar in all simulations, the physical mechanisms that govern the evolution of the helicities are markedly distinct from one simulation to another: the evolution of volume-threading helicity can be governed by boundary fluxes or helicity transfer, depending on the simulation. The eruption takes place for the same value of the ratio of the current-carrying helicity to the total helicity in all simulations. However, our study highlights that this threshold can be reached in different ways, with different helicity-related processes dominating for different photospheric flows. This means that the details of the pre-eruptive dynamics do not influence the eruption-onset helicity-related threshold. Nevertheless, the helicity-flux dynamics may be more or less efficient in changing the time required to reach the onset of the eruption., Comment: 18 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

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., Comment: 8 pages, 4 Figures, 1 Table, accepted for publication in ApJ Letters

Relative magnetic helicity is a gauge invariant quantity suitable for the study of the magnetic helicity content of heliospheric plasmas. Relative magnetic helicity can be decomposed uniquely into two gauge invariant quantities, the magnetic helicity of the non-potential component of the field, and a complementary volume-threading helicity. Recent analysis of numerical experiments simulating the generation of solar eruptions have shown that the ratio of the non-potential helicity to the total relative helicity is a clear marker of the eruptivity of the magnetic system, and that the high value of that quantity could be a sufficient condition for the onset of the instability generating the eruptions. The present study introduces the first analytical examination of the time variations of these non-potential and volume-threading helicities. The validity of the analytical formulas derived are confirmed with analysis of three-dimensional (3D) magnetohydrodynamics (MHD) simulations of solar coronal dynamics. Both the analytical investigation, and the numerical application show that, unlike magnetic helicity, the non-potential and the volume-threading helicities are not conserved quantities, even in the ideal MHD regime. A term corresponding to the transformation between the non-potential and volume-threading helicities frequently dominates their dynamics. This finding has an important consequence for their estimation in the solar corona: unlike with relative helicity, their volume coronal evolution cannot be ascertained by the flux of these quantities through the volume's boundaries. Only techniques extrapolating the 3D coronal field will enable both the proper study of the non-potential and volume-threading helicities, and the observational analysis of helicity-based solar-eruptivity proxies., Comment: 24 pages, 14 figures, accepted for publication in the Astrophysical Journal

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., Comment: 11 pages, 6 figures, in press in the Astrophysical Journal

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., Comment: 46 pages, 16 figures, accepted for publication in Space Weather and Space Climate on June, 8th

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.

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., Comment: Accepted for publication in Space Science Reviews. First paper by the ISSI International Team on 'Magnetic Helicity estimations in models and observations of the solar magnetic field', 54 pages, 12 figure. Data available at http://www.issibern.ch/teams/magnetichelicity/

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., Comment: Accepted for publication in ApJ

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., Comment: 7 pages, 5 figures, accepted to ApJ Letters

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., Comment: Accepted by ApJ

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., Comment: 15 pages, 13 figures

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., Comment: 8 pages, 8 figures, published in Astronomy and Astrophysics