Showing 1–50 of 50 results for author: Pariat, E

Generic low-atmosphere signatures of swirled-anemone jets

Authors: Reetika Joshi, Guillaume Aulanier, Alice Radcliffe, Luc Rouppe van der Voort, Etienne Pariat, Daniel Nóbrega-Siverio, Brigitte Schmieder

Abstract: 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 re… ▽ More 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 △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 13 pages, 8 figures, Accepted for publication in Astronomy and Astrophysics

On the use of relative field line helicity as an indicator for solar eruptivity

Authors: K. Moraitis, S. Patsourakos, A. Nindos, J. K. Thalmann, É. Pariat

Abstract: Context. Relative field line helicity (RFLH) is a recently developed quantity which can approximate the density of relative magnetic helicity. Aims. This paper aims to determine whether RFLH can be used as an indicator of solar eruptivity. Methods. Starting from magnetographic observations from the Helioseismic and Magnetic Imager instrument onboard the Solar Dynamic Observatory of a sample of sev… ▽ More Context. Relative field line helicity (RFLH) is a recently developed quantity which can approximate the density of relative magnetic helicity. Aims. This paper aims to determine whether RFLH can be used as an indicator of solar eruptivity. Methods. Starting from magnetographic observations from the Helioseismic and Magnetic Imager instrument onboard the Solar Dynamic Observatory of a sample of seven solar active regions (ARs), which comprises over 2000 individual snapshots, we reconstruct the AR's coronal magnetic field with a widely-used non-linear force-free method. This enables us to compute RFLH using two independent gauge conditions for the vector potentials. We focus our study around the times of strong flares in the ARs, above the M class, and in regions around the polarity inversion lines (PILs) of the magnetic field, and of RFLH. Results. We find that the temporal profiles of the relative helicity that is contained in the magnetic PIL follow those of the relative helicity that is computed by the accurate volume method for the whole AR. Additionally, the PIL relative helicity can be used to define a parameter similar to the well-known parameter R (Schrijver 2007), whose high values are related with increased flaring probability. This helicity-based R-parameter correlates well with the original one, showing in some cases even higher values, and additionally, it experiences more pronounced decreases during flares. This means that there exists at least one parameter deduced from RFLH, that has important value as a solar eruptivity indicator. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: accepted by Astronomy & Astrophysics

Stereoscopic disambiguation of vector magnetograms: first applications to SO/PHI-HRT data

Authors: G. Valori, D. Calchetti, A. Moreno Vacas, É. Pariat, S. K. Solanki, P. Löschl, J. Hirzberger, S. Parenti, K. Albert, N. Albelo Jorge, A. Álvarez-Herrero, T. Appourchaux, L. R. Bellot Rubio, J. Blanco Rodríguez, A. Campos-Jara, A. Feller, A. Gandorfer, P. García Parejo, D. Germerott, L. Gizon, J. M. Gómez Cama, L. Guerrero, P. Gutierrez-Marques, F. Kahil, M. Kolleck , et al. (12 additional authors not shown)

Abstract: Spectropolarimetric reconstructions of the photospheric vector magnetic field are intrinsically limited by the 180$^\circ$-ambiguity in the orientation of the transverse component. So far, the removal of such an ambiguity has required assumptions about the properties of the photospheric field, which makes disambiguation methods model-dependent. The basic idea is that the unambiguous line-of-sight… ▽ More Spectropolarimetric reconstructions of the photospheric vector magnetic field are intrinsically limited by the 180$^\circ$-ambiguity in the orientation of the transverse component. So far, the removal of such an ambiguity has required assumptions about the properties of the photospheric field, which makes disambiguation methods model-dependent. The basic idea is that the unambiguous line-of-sight component of the field measured from one vantage point will generally have a non-zero projection on the ambiguous transverse component measured by the second telescope, thereby determining the ``true'' orientation of the transverse field. Such an idea was developed and implemented in the Stereoscopic Disambiguation Method (SDM), which was recently tested using numerical simulations. In this work we present a first application of the SDM to data obtained by the High Resolution Telescope (HRT) onboard Solar Orbiter during the March 2022 campaign, when the angle with Earth was 27 degrees. The method is successfully applied to remove the ambiguity in the transverse component of the vector magnetogram solely using observations (from HRT and from the Helioseismic and Magnetic Imager), for the first time. The SDM is proven to provide observation-only disambiguated vector magnetograms that are spatially homogeneous and consistent. A discussion about the sources of error that may limit the accuracy of the method, and of the strategies to remove them in future applications, is also presented. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 32 pages, 12 figures, accepted in A&A on 09/07/2023

Journal ref: A&A 677, A25 (2023)

Comparison of magnetic energy and helicity in coronal jet simulations

Authors: E. Pariat, P. F. Wyper, L. Linan

Abstract: While free/non-potential magnetic energy is a necessary element of any active phenomenon in the solar corona, its role as a marker of the trigger of eruptive process remains elusive. Based on the unique decomposition of the magnetic field into potential and non-potential components, magnetic energy and helicity can also both be uniquely decomposed into two quantities. Using two 3D MHD parametric s… ▽ More While free/non-potential magnetic energy is a necessary element of any active phenomenon in the solar corona, its role as a marker of the trigger of eruptive process remains elusive. Based on the unique decomposition of the magnetic field into potential and non-potential components, magnetic energy and helicity can also both be uniquely decomposed into two quantities. Using two 3D MHD parametric simulations of a configuration that can produce coronal jets, we compare the dynamics of the magnetic energies and of the relative magnetic helicities. Both simulations share the same initial set-up and line-tied bottom-boundary driving profile. However, they differ by the duration of the forcing. In one simulation the system is driven sufficiently so that an helical jet is induced. The generation of the jet is however markedly delayed: a relatively long phase of lower-intensity reconnection takes place before the jet is eventually induced. In the other reference simulation, the system is driven during a shorter time, and no jet is produced. As expected, we observe that the Jet producing simulation contains a higher value of non-potential energy and non-potential helicity. Focusing on the phase between the end of the driving-phase and the jet generation, we note that magnetic energies remain relatively constant, while magnetic helicities have a noticeable evolution. During this post-driving phase, the ratio of the non-potential to total magnetic energy very slightly decreases while the helicity eruptivity index significantly increases. The jet is generated when the system is at the highest value of this helicity eruptivity index. This proxy critically decreases during the jet generation phase. The free energy also decreases but does not present any peak when the jet is being generated. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: 19 pages, 7 Figures, accepted in Astronomy and Astrophysics

Journal ref: A&A 669, A33 (2023)

Changes of Magnetic Energy and Helicity in Solar Active Regions from Major Flares

Authors: Yang Liu, Brian T. Welsch, Gherardo Valori, Manolis K. Georgoulis, Yang Guo, Etienne Pariat, Sung-Hong Park, Julia K. Thalmann

Abstract: 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 20… ▽ More 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. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 25 pages, 15 figures, 1 table, accepted for publication in the Astrophysical Journal

Disambiguation of Vector Magnetograms by Stereoscopic Observations from the Solar Orbiter/Polarimetric and Helioseismic Imager (PHI) and the Solar Dynamic Observatory (SDO)/Helioseismic and Magnetic Imager (HMI)

Authors: Gherardo Valori, Philipp Löschel, David Stansby, Etienne Pariat, Johann Hirzberger, Feng Chen

Abstract: 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 s… ▽ More 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. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: 41 pages; 18 Figures, accepted for publication in Solar Physics

Magnetic helicity estimations in models and observations of the solar magnetic field. Part IV: Application to solar observations

Authors: J. K. Thalmann, M. K. Georgoulis, Y. Liu, E. Pariat, G. Valori, S. Anfinogentov, F. Chen, Y. Guo, K. Moraitis, S. Yang, A. Mastrano

Abstract: In this ISSI-supported series of studies on magnetic helicity in the Sun, we systematically implement different magnetic helicity calculation methods on high-quality solar magnetogram observations. We apply finite-volume, discrete flux tube (in particular, connectivity-based) and flux-integration methods to data from Hinode's Solar Optical Telescope. The target is NOAA active region 10930 during a… ▽ More In this ISSI-supported series of studies on magnetic helicity in the Sun, we systematically implement different magnetic helicity calculation methods on high-quality solar magnetogram observations. We apply finite-volume, discrete flux tube (in particular, connectivity-based) and flux-integration methods to data from Hinode's Solar Optical Telescope. The target is NOAA active region 10930 during a ~1.5 day interval in December 2006 that included a major eruptive flare (SOL2006-12-13T02:14X3.4). Finite-volume and connectivity-based methods yield instantaneous budgets of the coronal magnetic helicity, while the flux-integration methods allow an estimate of the accumulated helicity injected through the photosphere. The objectives of our work are twofold: A cross-validation of methods, as well as an interpretation of the complex events leading to the eruption. To the first objective, we find (i) strong agreement among the finite-volume methods, (ii) a moderate agreement between the connectivity-based and finite-volume methods, (iii) an excellent agreement between the flux-integration methods, and (iv) an overall agreement between finite-volume and flux-integration based estimates regarding the predominant sign and magnitude of the helicity. To the second objective, we are confident that the photospheric helicity flux significantly contributed to the coronal helicity budget, and that a right-handed structure erupted from a predominantly left-handed corona during the X-class flare. Overall, we find that the use of different methods to estimate the (accumulated) coronal helicity may be necessary in order to draw a complete picture of an active-region corona, given the careful handling of identified data (preparation) issues, which otherwise would mislead the event analysis and interpretation. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: Accepted for publication in The Astrophysical Journal, 29 pages, 9 figures

The Flare Likelihood and Region Eruption Forecasting (FLARECAST) Project: Flare forecasting in the big data & machine learning era

Authors: M. K. Georgoulis, D. S. Bloomfield, M. Piana, A. M. Massone, M. Soldati, P. T. Gallagher, E. Pariat, N. Vilmer, E. Buchlin, F. Baudin, A. Csillaghy, H. Sathiapal, D. R. Jackson, P. Alingery, F. Benvenuto, C. Campi, K. Florios, C. Gontikakis, C. Guennou, J. A. Guerra, I. Kontogiannis, V. Latorre, S. A. Murray, S. -H. Park, S. von Stachelski , et al. (3 additional authors not shown)

Abstract: The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) differe… ▽ More The EU funded the FLARECAST project, that ran from Jan 2015 until Feb 2018. FLARECAST had a R2O focus, and introduced several innovations into the discipline of solar flare forecasting. FLARECAST innovations were: first, the treatment of hundreds of physical properties viewed as promising flare predictors on equal footing, extending multiple previous works; second, the use of fourteen (14) different ML techniques, also on equal footing, to optimize the immense Big Data parameter space created by these many predictors; third, the establishment of a robust, three-pronged communication effort oriented toward policy makers, space-weather stakeholders and the wider public. FLARECAST pledged to make all its data, codes and infrastructure openly available worldwide. The combined use of 170+ properties (a total of 209 predictors are now available) in multiple ML algorithms, some of which were designed exclusively for the project, gave rise to changing sets of best-performing predictors for the forecasting of different flaring levels. At the same time, FLARECAST reaffirmed the importance of rigorous training and testing practices to avoid overly optimistic pre-operational prediction performance. In addition, the project has (a) tested new and revisited physically intuitive flare predictors and (b) provided meaningful clues toward the transition from flares to eruptive flares, namely, events associated with coronal mass ejections (CMEs). These leads, along with the FLARECAST data, algorithms and infrastructure, could help facilitate integrated space-weather forecasting efforts that take steps to avoid effort duplication. In spite of being one of the most intensive and systematic flare forecasting efforts to-date, FLARECAST has not managed to convincingly lift the barrier of stochasticity in solar flare occurrence and forecasting: solar flare prediction thus remains inherently probabilistic. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: 67 pages, 14 figures; submitted

Journal ref: Journal of Space Weather and Space Climate, 2021

The Solar Orbiter Science Activity Plan: translating solar and heliospheric physics questions into action

Authors: I. Zouganelis, A. De Groof, A. P. Walsh, D. R. Williams, D. Mueller, O. C. St Cyr, F. Auchere, D. Berghmans, A. Fludra, T. S. Horbury, R. A. Howard, S. Krucker, M. Maksimovic, C. J. Owen, J. Rodriiguez-Pacheco, M. Romoli, S. K. Solanki, C. Watson, L. Sanchez, J. Lefort, P. Osuna, H. R. Gilbert, T. Nieves-Chinchilla, L. Abbo, O. Alexandrova , et al. (160 additional authors not shown)

Abstract: Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operat… ▽ More Solar Orbiter is the first space mission observing the solar plasma both in situ and remotely, from a close distance, in and out of the ecliptic. The ultimate goal is to understand how the Sun produces and controls the heliosphere, filling the Solar System and driving the planetary environments. With six remote-sensing and four in-situ instrument suites, the coordination and planning of the operations are essential to address the following four top-level science questions: (1) What drives the solar wind and where does the coronal magnetic field originate? (2) How do solar transients drive heliospheric variability? (3) How do solar eruptions produce energetic particle radiation that fills the heliosphere? (4) How does the solar dynamo work and drive connections between the Sun and the heliosphere? Maximising the mission's science return requires considering the characteristics of each orbit, including the relative position of the spacecraft to Earth (affecting downlink rates), trajectory events (such as gravitational assist manoeuvres), and the phase of the solar activity cycle. Furthermore, since each orbit's science telemetry will be downloaded over the course of the following orbit, science operations must be planned at mission level, rather than at the level of individual orbits. It is important to explore the way in which those science questions are translated into an actual plan of observations that fits into the mission, thus ensuring that no opportunities are missed. First, the overarching goals are broken down into specific, answerable questions along with the required observations and the so-called Science Activity Plan (SAP) is developed to achieve this. The SAP groups objectives that require similar observations into Solar Orbiter Observing Plans (SOOPs), resulting in a strategic, top-level view of the optimal opportunities for science observations during the mission lifetime. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 20 pages, 1 figure, accepted by Astronomy & Astrophysics

Journal ref: A&A 642, A3 (2020)

Additivity of relative magnetic helicity in finite volumes

Authors: Gherardo Valori, Pascal Démoulin, Etienne Pariat, Anthony Yeates, Kostas Moraitis, Luis Linan

Abstract: 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 of… ▽ More 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. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 14 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 643, A26 (2020)

Energy and helicity fluxes in line-tied eruptive simulations

Authors: Luis Linan, Étienne Pariat, Guillaume Aulanier, Kostas Moraitis, Gherardo Valori

Abstract: 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-li… ▽ More 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. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 18 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 636, A41 (2020)

Magnetic helicity budget of solar active regions prolific of eruptive and confined flares

Authors: J. K. Thalmann, K. Moraitis, L. Linan, E. Pariat, G. Valori, K. Dalmasse

Abstract: 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 investigat… ▽ More 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}/φ^2$, represents no indicator for the eruptivity. △ Less

Submitted 15 October, 2019; originally announced October 2019.

Comments: 11 pages, 3 figures, accepted for publication in ApJ

Magnetic helicity and eruptivity in active region 12673

Authors: K. Moraitis, X. Sun, E. Pariat, L. Linan

Abstract: 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 e… ▽ More 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. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: accepted by A&A

Journal ref: A&A 628, A50 (2019)

On the reliability of magnetic energy and helicity computations based on nonlinear force-free coronal magnetic field models

Authors: Julia K. Thalmann, Luis Linan, Etienne Pariat, Gherardo Valori

Abstract: 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 concern… ▽ More We demonstrate the sensitivity of magnetic energy and helicity computations regarding the quality of the underlying coronal magnetic field model. We apply the method of Wiegelmann & Inhester (2010) to a series of SDO/HMI vector magnetograms, and discuss nonlinear force-free (NLFF) solutions based on two different sets of the free model parameters. The two time series differ from each other concerning their force-free and solenoidal quality. Both force- and divergence-freeness are required for a consistent NLFF solution. Full satisfaction of the solenoidal property is inherent in the definition of relative magnetic helicity in order to insure gauge-independence. We apply two different magnetic helicity computation methods (Thalmann et al. 2011; Valori et al. 2012) to both NLFF time series and find that the output is highly dependent on the level to which the NLFF magnetic fields satisfy the divergence-free condition, with the computed magnetic energy being less sensitive than the relative helicity. Proxies for the non-potentiality and eruptivity derived from both quantities are also shown to depend strongly on the solenoidal property of the NLFF fields. As a reference for future applications, we provide quantitative thresholds for the force- and divergence-freeness, for the assurance of reliable computation of magnetic energy and helicity, and of their related eruptivity proxies. △ Less

Submitted 2 July, 2019; originally announced July 2019.

Comments: 8 pages, 4 Figures, 1 Table, accepted for publication in ApJ Letters

Relative magnetic field line helicity

Authors: K. Moraitis, E. Pariat, G. Valori, K. Dalmasse

Abstract: 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 analys… ▽ More 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. △ Less

Submitted 27 February, 2019; originally announced February 2019.

Comments: accepted by A&A

Journal ref: A&A 624, A51 (2019)

Time variations of the non-potential and volume-threading magnetic helicities

Authors: Luis Linan, Étienne Pariat, Kostas Moraitis, Gherardo Valori, James E. Leake

Abstract: 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 g… ▽ More 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. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 24 pages, 14 figures, accepted for publication in the Astrophysical Journal

Threshold of non-potential magnetic helicity ratios at the onset of solar eruptions

Authors: Francesco P. Zuccarello, Etienne Pariat, Gherardo Valori, Luis Linan

Abstract: 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 o… ▽ More 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. △ Less

Submitted 10 July, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 11 pages, 6 figures, in press in the Astrophysical Journal

Computation of Relative Magnetic Helicity in Spherical Coordinates

Authors: K. Moraitis, É. Pariat, A. Savcheva, G. Valori

Abstract: Magnetic helicity is a quantity of great importance in solar studies because it is conserved in ideal magneto-hydrodynamics. While many methods to compute magnetic helicity in Cartesian finite volumes exist, in spherical coordinates, the natural coordinate system for solar applications, helicity is only treated approximately. We present here a method to properly compute relative magnetic helicity… ▽ More Magnetic helicity is a quantity of great importance in solar studies because it is conserved in ideal magneto-hydrodynamics. While many methods to compute magnetic helicity in Cartesian finite volumes exist, in spherical coordinates, the natural coordinate system for solar applications, helicity is only treated approximately. We present here a method to properly compute relative magnetic helicity in spherical geometry. The volumes considered are finite, of shell or wedge shape, and the three-dimensional magnetic field is considered fully known throughout the studied domain. Testing of the method with well-known, semi-analytic, force-free magnetic-field models reveals that it has excellent accuracy. Further application to a set of nonlinear force-free reconstructions of the magnetic field of solar active regions, and comparison with an approximate method used in the past, indicates that the proposed methodology can be significantly more accurate, thus making our method a promising tool in helicity studies that employ the spherical geometry. Additionally, the range of applicability of the approximate method is determined and discussed. △ Less

Submitted 8 June, 2018; originally announced June 2018.

Comments: accepted for publication in Solar Physics

Studying the transfer of magnetic helicity in solar active regions with the connectivity-based helicity flux density method

Authors: K. Dalmasse, E. Pariat, G. Valori, J. Jing, P. Démoulin

Abstract: 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 transfe… ▽ More 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. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: 15 pages, 7 figures ; Accepted for publication in ApJ

Testing predictors of eruptivity using parametric flux emergence simulations

Authors: Chloé Guennou, Etienne Pariat, Nicole Vilmer, James E. Leake

Abstract: 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, o… ▽ More 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. △ Less

Submitted 15 June, 2017; originally announced June 2017.

Comments: 46 pages, 16 figures, accepted for publication in Space Weather and Space Climate on June, 8th

Magnetic Helicity Estimations in Models and Observations of the Solar Magnetic Field. Part III: Twist Number Method

Authors: Y. Guo, E. Pariat, G. Valori, S. Anfinogentov, F. Chen, M. Georgoulis, Y. Liu, K. Moraitis, J. K. Thalmann, S. Yang

Abstract: 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émoulin 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-sep… ▽ More 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émoulin 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. △ Less

Submitted 7 April, 2017; originally announced April 2017.

Comments: To be published in ApJ

Flux rope, hyperbolic flux tube, and late EUV phases in a non-eruptive circular-ribbon flare

Authors: S. Masson, E. Pariat, G. Valori, N. Deng, C. Liu, H. Wang, H. Reid

Abstract: 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… ▽ More 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~Å 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. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: Astronomy & Astrophysics, in press

Journal ref: A&A 604, A76 (2017)

Relative magnetic helicity as a diagnostic of solar eruptivity

Authors: E. Pariat, J. E. Leake, G. Valori, M. G. Linton, F. P. Zuccarello, K. Dalmasse

Abstract: 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 disc… ▽ More 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. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Journal ref: A&A 601, A125 (2017)

Observational Evidence of Magnetic Reconnection for Brightenings and Transition Region Arcades in IRIS observations

Authors: Jie Zhao, Brigitte Schmieder, Hui Li, Etienne Pariat, Xiaoshuai Zhu, Li Feng, Michalina Grubecka

Abstract: 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… ▽ More 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. △ Less

Submitted 29 January, 2017; originally announced January 2017.

Analysis and modelling of recurrent solar flares observed with Hinode/EIS on March 9, 2012

Authors: V. Polito, G. Del Zanna, G. Valori, E. Pariat, H. E. Mason, J. Dudik, M. Janvier

Abstract: 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 occurre… ▽ More 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. △ Less

Submitted 11 December, 2016; originally announced December 2016.

Journal ref: A&A 601, A39 (2017)

Magnetic helicity estimations in models and observations of the solar magnetic field. Part I: Finite volume methods

Authors: Gherardo Valori, Etienne Pariat, Sergey Anfinogentov, Feng Chen, Manolis K. Georgoulis, Yang Guo, Yang Liu, Kostas Moraitis, Julia K. Thalmann, Shangbin Yang

Abstract: 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 t… ▽ More 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. △ Less

Submitted 7 October, 2016; originally announced October 2016.

Comments: 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/

Blowout Jets and Impulsive Eruptive Flare in a Bald-Patch Topology

Authors: R. Chandra, C. H. Mandrini, B. Schmieder, B. Joshi, G. D. Cristiani, H. Cremades, E. Pariat, F. A. Nuevo, A. K. Srivastava, W. Uddin

Abstract: 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… ▽ More 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. △ Less

Submitted 6 October, 2016; originally announced October 2016.

Journal ref: A&A 598, A41 (2017)

A model for straight and helical solar jets: II. Parametric study of the plasma beta

Authors: E. Pariat, K. Dalmasse, C. R. DeVore, S. K. Antiochos, J. T. Karpen

Abstract: 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 coro… ▽ More 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 $β$ on the generation and properties of solar-like jets. The parametric study validates our model of jets for plasma $β$ 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 $β\le 1$. This study introduces the new result that the plasma $β$ 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. △ Less

Submitted 28 September, 2016; originally announced September 2016.

Comments: Accepted in Astronomy and Astrophysics

Journal ref: A&A 596, A36 (2016)

Solar Coronal Jets: Observations, Theory, and Modeling

Authors: N. E. Raouafi, S. Patsourakos, E. Pariat, P. R. Young, A. C. Sterling, A. Savcheva, M. Shimojo, F. Moreno-Insertis, C. R. DeVore, V. Archontis, T. Török, H. Mason, W. Curdt, K. Meyer, K. Dalmasse, Y. Matsui

Abstract: 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 explos… ▽ More 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. △ Less

Submitted 7 July, 2016; originally announced July 2016.

Comments: 53 pages, 24 figures

Journal ref: Space Science Reviews, 04 July 2016

Reconnection-Driven Coronal-Hole Jets with Gravity and Solar Wind

Authors: J. T. Karpen, C. R. DeVore, S. K. Antiochos, E. Pariat

Abstract: 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… ▽ More 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. △ Less

Submitted 29 June, 2016; originally announced June 2016.

Comments: 29 pages including 12 figures

Evolution of Flare Ribbons, Electric Currents and Quasi-separatrix Layers During an X-class Flare

Authors: M. Janvier, A. Savcheva, E. Pariat, S. Tassev, S. Millholland, V. Bommier, P. McCauley, S. McKillop, F. Dougan

Abstract: 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 mod… ▽ More 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. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 17 pages, 14 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 591, A141 (2016)

Hooked flare ribbons and flux-rope related QSL footprints

Authors: Jie Zhao, Stuart A. Gilchrist, Guillaume Aulanier, Brigitte Schmieder, Etienne Pariat, Hui Li

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2016; originally announced March 2016.

Comments: Accepted for publication in ApJ

A Circular-ribbon Solar Flare Following an Asymmetric Filament Eruption

Authors: Chang Liu, Na Deng, Rui Liu, Jeongwoo Lee, Etienne Pariat, Thomas Wiegelmann, Yang Liu, Lucia Kleint, Haimin Wang

Abstract: 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 forc… ▽ More 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. △ Less

Submitted 28 September, 2015; originally announced September 2015.

Comments: 7 pages, 5 figures, accepted to ApJ Letters

The origin of net electric currents in solar active regions

Authors: K. Dalmasse, G. Aulanier, P. Démoulin, B. Kliem, T. Török, E. Pariat

Abstract: 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… ▽ More 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. △ Less

Submitted 17 July, 2015; originally announced July 2015.

Comments: 14 pages and 11 figures (Accepted in The Astrophysical Journal)

Testing magnetic helicity conservation in a solar-like active event

Authors: E. Pariat, G. Valori, P. Démoulin, K. Dalmasse

Abstract: 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 bare… ▽ More 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. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Comments: Accepted in Astronomy & Astrophysics 16 Pages; 10 figures

Journal ref: A&A 580, A128 (2015)

The Relation between Solar Eruption Topologies and Observed Flare Features I: Flare Ribbons

Authors: A. Savcheva, E. Pariat, S. McKillop, P. McCauley, E. Hanson, Y. Su, E. Werner, E. E. DeLuca

Abstract: 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 f… ▽ More 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. △ Less

Submitted 10 June, 2015; originally announced June 2015.

Coronal magnetic reconnection driven by CME expansion -- the 2011 June 7 event

Authors: L. van Driel-Gesztelyi, D. Baker, T. Torok, E. Pariat, L. M. Green, D. R. Williams, J. Carlyle, G. Valori, P. Demoulin, B. Kliem, D. M. Long, S. A. Matthews, J. -M. Malherbe

Abstract: 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… ▽ More 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. △ Less

Submitted 12 June, 2014; originally announced June 2014.

Comments: 12 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 788, Issue 1, 85, 12 pp. (2014)

Temporal Evolution of the Magnetic Topology of the NOAA Active Region 11158

Authors: Jie Zhao, Hui Li, Etienne Pariat, Brigitte Schmieder, Yang Guo, Thomas Wiegelmann

Abstract: 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-separat… ▽ More 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. △ Less

Submitted 19 April, 2014; originally announced April 2014.

Comments: Accepted by ApJ

Electric current in flares ribbons: observations and 3D standard model

Authors: Miho Janvier, G. Aulanier, V. Bommier, B. Schmieder, P. Démoulin, E. Pariat

Abstract: 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 th… ▽ More 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. △ Less

Submitted 23 April, 2014; v1 submitted 9 February, 2014; originally announced February 2014.

Comments: 12 pages, 7 figures

Twist Accumulation and Topology Structure of a Solar Magnetic Flux Rope

Authors: Y. Guo, M. D. Ding, X. Cheng, J. Zhao, E. Pariat

Abstract: 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… ▽ More 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. △ Less

Submitted 8 November, 2013; originally announced November 2013.

Comments: 32 pages, 9 figures, accepted for publication in ApJ

First observational application of a connectivity--based helicity flux density

Authors: K. Dalmasse, E. Pariat, G. Valori, P. Démoulin, L. M. Green

Abstract: 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 connectivit… ▽ More 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. △ Less

Submitted 10 July, 2013; originally announced July 2013.

Comments: 4 pages, 3 figures; published online in A&A 555, L6 (2013)

Photospheric Injection of Magnetic Helicity: Connectivity--based Flux Density Method

Authors: K. Dalmasse, E. Pariat, P. Démoulin, G. Aulanier

Abstract: 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 magn… ▽ More 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. △ Less

Submitted 10 July, 2013; originally announced July 2013.

Comments: Solar Physics, June 2013 (this is the version of the author, a definitive version is now available in the online journal)

The standard flare model in three dimensions III. Slip-running reconnection properties

Authors: Miho Janvier, Guillaume Aulanier, Etienne Pariat, Pascal Demoulin

Abstract: 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 r… ▽ More 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. △ Less

Submitted 17 May, 2013; originally announced May 2013.

Comments: 15 pages, 13 figures

Accuracy of magnetic energy computations

Authors: G. Valori, P. Demoulin, E. Pariat, S. Masson

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2013; originally announced March 2013.

The standard flare model in three dimensions. II. Upper limit on solar flare energy

Authors: G. Aulanier, P. Demoulin, C. J. Schrijver, M. Janvier, E. Pariat, B. Schmieder

Abstract: 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 ab… ▽ More 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. △ Less

Submitted 10 December, 2012; originally announced December 2012.

Comments: 8 pages, 4 figures

X-ray and UV investigation into the magnetic connectivity of a solar flare

Authors: Hamish A. S. Reid, Nicole Vilmer, Guillaume Aulanier, Etienne Pariat

Abstract: 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 Å. The flare appears to have two distinct phases of emission separate… ▽ More 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 Å. 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 Å 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. △ Less

Submitted 10 October, 2012; originally announced October 2012.

Comments: 8 pages, 8 figures, published in Astronomy and Astrophysics

Interchange Slip-Running Reconnection and Sweeping SEP Beams

Authors: S. Masson, G. Aulanier, E. Pariat, K. -L. Klein

Abstract: 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… ▽ More 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. △ Less

Submitted 26 September, 2011; originally announced September 2011.

Homologous Flares and Magnetic Field Topology in Active Region NOAA 10501 on 20 November 2003

Authors: R. Chandra, B. Schmieder, C. H. Mandrini, P. Démoulin, E. Pariat, T. Torok, W. Uddin

Abstract: 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… ▽ More 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. △ Less

Submitted 4 November, 2010; originally announced November 2010.

Comments: 24 pages, 10 figures, Solar Physics (accepted)

How can a Negative Magnetic Helicity Active Region Generate a Positive Helicity Magnetic Cloud ?

Authors: R. Chandra, E. Pariat, B. Schmieder, C. H. Mandrini, W. Uddin

Abstract: 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… ▽ More 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. △ Less

Submitted 6 October, 2009; originally announced October 2009.

Comments: 26 pages, 11 figures, accepted in Solar Physics

STEREO/SECCHI Stereoscopic Observations Constraining the Initiation of Polar Coronal Jets

Authors: S. Patsourakos, E. Pariat, A. Vourlidas, S. K. Antiochos, J. P. Wuelser

Abstract: 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 a… ▽ More 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. △ Less

Submitted 30 April, 2008; originally announced April 2008.

Comments: ApJL, in press