Issue 
A&A
Volume 636, April 2020



Article Number  A41  
Number of page(s)  16  
Section  The Sun and the Heliosphere  
DOI  https://doi.org/10.1051/00046361/202037548  
Published online  13 April 2020 
Energy and helicity fluxes in linetied eruptive simulations
^{1}
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 place Jules Janssen, 92195 Meudon, France
email: luis.linan@obspm.fr
^{2}
Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK
Received:
22
January
2020
Accepted:
26
February
2020
Context. Conservation properties of magnetic helicity and energy in the quasiideal and lowβ solar corona make these two quantities relevant for the study of solar active regions and eruptions.
Aims. 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 volumethreading and currentcarrying helicities. In this study, we perform a coupled analysis of their behaviors in a set of parametric 3D magnetohydrodynamic (MHD) simulations of solarlike eruptions.
Methods. We present the general formulations for the timevarying 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 boundarydriving motions. Thus, we investigated the impact of different active regions surface flows on the development of the energy and helicityrelated quantities.
Results. 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 volumethreading helicity can be governed by boundary fluxes or helicity transfer, depending on the simulation.
Conclusions. The eruption takes place for the same value of the ratio of the currentcarrying helicity to the total helicity in all simulations. However, our study highlights that this threshold can be reached in different ways, with different helicityrelated processes dominating for different photospheric flows. This means that the details of the preeruptive dynamics do not influence the eruptiononset helicityrelated threshold. Nevertheless, the helicityflux dynamics may be more or less efficient in changing the time required to reach the onset of the eruption.
Key words: magnetic fields / Sun: photosphere / Sun: corona / Sun: flares / magnetohydrodynamics (MHD) / Sun: activity
© L. Linan et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Magnetic helicity is a volumeintegrated ideal magnetohydrodynamic (MHD) invariant describing the level of twist and entanglement of the magnetic field lines. Initially introduced by Elsasser (1956), magnetic helicity is a conserved quantity within the ideal MHD paradigm (Woltjer 1958). However, the strict definition of Elsasser is gauge invariant only for magnetically bounded system, a condition that is not satisfied in most cases, such as the solar atmosphere. This led to the introduction by Berger & Field (1984) and Finn & Antonsen (1985) of the relative magnetic helicity, a gaugeinvariant quantity suitable for use in solar physics and more generally for natural plasmas.
Using a numerical simulation, Pariat et al. (2015) confirmed the hypothesis introduced by Taylor (1974) that even in presence of nonideal processes, the dissipation of magnetic helicity is negligible. Relative magnetic helicity cannot be dissipated or created within the corona, therefore it can only be transported or annihilated. This conservation property has several major consequences, one of which might be that coronal mass ejections (CMEs) are the consequence of the evacuation of an excess of helicity (Rust 1994; Low 1996).
In recent years, magnetic helicity has been at the heart of many studies dealing with various topics such as the generation of solar eruptions (e.g., Kusano et al. 2004; Longcope & Beveridge 2007; Priest et al. 2016), magnetic reconnection (e.g., Linton et al. 2001; Linton & Antiochos 2002; Del Sordo et al. 2010), solar filaments (e.g., Antiochos 2013; Knizhnik et al. 2015; Zhao et al. 2015), and solar and stellar dynamos (e.g., Brandenburg & Subramanian 2005; Simon 2012).
Magnetic energy is another relevant quantity in MHD with which eruptivity in the solar corona is studied because most of solar events are driven magnetically (e.g., Schrijver & Zwaan 2008). From the point of view of the energetic budget, magnetic energy is the only source of energy that can generate the powerful events that are observed in the solar atmosphere, such as coronal mass ejections, flares, and solar jets (Forbes 2000). Magnetic energy can be decomposed into a currentcarrying energy, known as free energy, and a potential energy (cf. Sect. 3.1). Solar flares and CMEs are characterized by a rapid change of the coronal magnetic field that does not change the radial component of the photospheric field. Because the potential field is determined by the radial magnetic field at the photosphere, only the free energy can therefore be converted into kinetic and thermal energies during fast coronal events (Aulanier et al. 2009; Karpen et al. 2012). The potential energy thus represents the lowest energy state of the magnetic field in the solar corona (e.g., Priest 2014).
An analysis of magnetic energies combined with a study of magnetic helicities appears a powerful tool for characterizing active regions and their evolutions toward eruptive events. However, measuring these quantities from observational data remains challenging. One possibility is to estimate the accumulation of magnetic helicity and energy in the solar corona by integrating their fluxes across the solar photosphere over time (Kusano et al. 2002; Nindos et al. 2003; Yamamoto et al. 2005; Yamamoto & Sakurai 2009). This method cannot trace the coronal evolution and requires highcadence timeseries magnetograms as well as the velocity fields on the photosphere. Because no direct observation of the photospheric velocity is available, it is obtained by inferring the magnetic field on the solar surface. Despite the progress made in deducing the velocity field (Kusano et al. 2002; Welsch et al. 2004; Longcope 2004) as well as further improvement on flux estimations (Pariat et al. 2005; Chae 2007; Liu & Schuck 2012, 2013; Dalmasse et al. 2014, 2018), the computation of magnetic energy and helicity fluxes remains very sensitive to the method that is used and to the quality of the observations. A different approach is to compute energy and helicity in coronal volumes. Because magnetic energy and helicity are volume integrals, properly computing them with this method requires the full 3D knowledge of the magnetic field in the volume that is studied. Currently, only 2D measurement on the solar surface are provided, therefore a 3D extrapolation of the magnetic field is a necessary step. The diverse methods based on the volumeintegration approach for estimating the magnetic relative helicity were benchmarked in Valori et al. (2016). Different solar active regions have previously been investigated (Valori et al. 2013; Moraitis et al. 2014, 2019; Guo et al. 2017; Polito et al. 2017; Temmer et al. 2017; James et al. 2018; Thalmann et al. 2019).
In parallel, the properties of both helicity and energy are still being studied in solarlike parametric simulations. Berger (2003) introduced the decomposition of the relative magnetic helicity into two gaugeinvariant components: a currentcarrying helicity related to the currentcarrying magnetic field, and a complementary volumethreading helicity. Pariat et al. (2017) followed and estimated these quantities in a set of seven simulations of the formation of solar active regions (Leake et al. 2013, 2014). The different simulations led to either stable or eruptive configurations. The authors found that it is possible to distinguish the two configurations by studying the ratio of the currentcarrying helicity to the relative helicity. The ratio before the eruption indeed presents high values only in the eruptive case. To better understand the properties of the relative helicity decomposition, Linan et al. (2018) provided the first analytical formulae of the timevariation of nonpotential and volumethreading helicity. They also computed and followed them in two simulations of Leake et al. (2013, 2014) and in a simulation of the generation of a coronal jet (Pariat et al. 2005). They found that the currentcarrying helicity does indeed not evolve as a result of boundary fluxes, but builds up through its exchange with the volumethreading helicity. The evolutions of the currentcarrying and the volumethreading helicities are partially controlled by a transfer term that reflects the exchange between these two types of helicity. This exchange term dominates the dynamics of the currentcarrying helicity at different instants of the simulations. This means that this helicity does not only evolve as a result of boundary fluxes. The eruption phases of these simulation follow the same dynamics: the currentcarrying helicity is first transformed into the volumethreading helicity, and then the latter is ejected from the domain by boundary fluxes. Moreover, the transfer term is expressed as a volume integral: consequently, these two helicities are not classically conserved quantities in the sense that they cannot be independently expressed as a flux through the boundaries, even in ideal MHD, unlike the relative magnetic helicity. This finding strengthens the knowledge of the properties of nonpotential and volumethreading helicity that was first studied by Moraitis et al. (2014).
Zuccarello et al. (2015) presented 3D parametric resistive MHD simulations of solar coronal eruptions. Simulations are distinguished by the different motions (linetied) applied on the photosphere with similar but distinct flux cancelation drivers. Their eruptions were driven by the torus instability (Aulanier et al. 2009; Démoulin & Aulanier 2010) and occurred at a precise time identified by a series a relaxation runs. Recently, these models were used to investigate the increase in the downward component of the Lorentz force density around an polarityinversion line in comparison with the photospheric observation (Barczynski et al. 2019). From these simulation, Zuccarello et al. (2018) studied the impact of the different boundary driving flows on the helicity and energy injection. They found that the helicity ratio of the currentcarrying helicity to the relative helicity is clearly associated with the eruption trigger because the eruptions within the different runs took place exactly when the ratio reached the very same threshold value.
Recently, the first preliminary observational tests confirmed the idea that the helicity ratio is a good predictor of eruptivity. Based on 3D extrapolation and using different nonlinear forcefree models, the time evolution of the helicity ratio has been investigated in three active regions: AR 12673 in Moraitis et al. (2019), the most active of the cycle 24; and AR 11158 and AR 12192 in Thalmann et al. (2019), two extensively studied active regions that generated eruptive and confined flares. However, complementary studies are still needed to understand how the different magnetic topologies observed in the solar corona are linked to the dynamics of the helicity ratio.
In the present study, we apply the helicity decomposition to the analysis of the simulations of Zuccarello et al. (2015, 2018) to investigate the timevariations of the different types of magnetic energy and helicity. In particular, we are interested in the way that the different boundary motion influence the helicity and energy dynamics. Zuccarello et al. (2018) showed that the different boundary motions lead to different efficiency in injecting helicity and energy in the domain. In the present work, we aim to explain the physical processes that are responsible for these differences: are they related to boundary fluxes, dissipation, or volume evolution? We also examine whether the transfer term between the two helicity components plays a major role in the helicity budgets, as has been observed in Linan et al. (2018).
Additionally, we study the dynamics of the helicities in comparison with the dynamics of their energy counterparts, for instance, currentcarrying helicity and free energy, and volumethreading helicity and potential energy. Our goal is to highlight the differences and the similarities in the helicity and energy buildup. This study aims to improve our knowledge on magnetic helicities and energies, and it is a necessary step to better understand the full topological and energetic complexity of solar active regions.
Our paper is divided into different sections that are organized as follows. First, we present the timevariation of nonpotential and volumethreading helicities (see Sect. 2). In the same way, we then introduce the different components of the magnetic energy and also their timevariation written for the specific case of resistive MHD (see Sect. 3). After presenting the simulations (see Sect. 4), we present the time evolution of the different quantities (see Sect. 5). Using our set of simulations, we investigate the role of the transfer term between the helicity components in their evolutions (see Sect. 6). While we investigate the difference in helicity dynamics in Sect. 7, we focus on the similarities between magnetic energy and helicity fluxes in Sect. 8. In the conclusion, we discuss the effect of the different boundarydriving motions on the energy and helicity injection and on the eruptivity helicity ratio.
2. Nonpotential and volumethreading helicities
In the fixed volume V bounded by the surface S, the magnetic helicity is defined as
with A the vector potential of the studied magnetic field B, i.e ∇ × A = B. In practice, this scalar description of the geometrical properties of magnetic field lines is general only if the magnetic field is tangential to the surface, that is, if V is a magnetically bounded volume. The magnetic helicity is gauge invariant if and only if this condition is respected. For the study of natural plasmas, especially in solar physics, the magnetic field does not satisfy this condition, the solar photosphere being subject to significant flux. Berger & Field (1984) introduced the relative magnetic helicity, a gaugeinvariant quantity, based on a reference field. Throughout the paper, we use the potential reference field B_{p} that has the same normal distribution of B throughout the surface S and satisfies
where n is the outwardpointing unit vector normal on S. The potential field can thus be defined by a scalar function, such as ∇ϕ = B_{p}, and ϕ is the solution of the Laplace equation,
When A_{p} is the vector potential of the potential field B_{p} = ∇×A_{p}, the relative magnetic helicity provided by Finn & Antonsen (1985) is defined as
In this form, the relative magnetic helicity is independently invariant to gauge transformation of both A and A_{p}. The difference between the potential field and the magnetic field can be written as a nonpotential magnetic field, B_{j} = B − B_{p}, associated with the vector A_{j}, defined as A_{j} = A − A_{p}, such as ∇ × A_{j} = B_{j}. When we use this unique decomposition of B and following the work of Berger (2003), Eq. (4) can be divided into two gaugeinvariant quantities:
where H_{j} is the currentcarrying magnetic helicity associated with only the currentcarrying component of the magnetic field B_{j}, and H_{pj} is the volumethreading helicity involving both B and B_{p}. By construction, both H_{j} and H_{pj} are gauge invariant because by virtue of Eq. (3), B_{j} has a vanishing normal component on the surface.
In resistive MHD, where E = −v × B + η∇×B (η being the resistivity, which is here assumed to be constant), Linan et al. (2018) established the following equation for the time evolution of the currentcarrying magnetic helicity H_{j}:
with
From this decomposition, Linan et al. (2018) obtained an equation for the timevariation that is composed only of gaugeinvariant terms:
with
All these terms initially provided by Linan et al. (2018) are recalled here because we analyze them and comment on them in the next sections.
Similarly, the time evolution of the volumethreading magnetic helicity H_{pj} can be decomposed as
with
The timevariation of H_{pj} can also be constructed with gaugeinvariant terms only:
with
These decompositions were obtained without any particular hypothesis on the gauge that is used. In particular, we are free to use the Coulomb gauge for A and A_{p}. With this choice, the volume terms dH_{j}/dt_{Bp, var} and dH_{pj}/dt_{Bp, var} both vanish. Thus dH_{j}/dt_{Own} and dH_{pj}/dt_{Own} only contain boundaryflux contributions. The transfer term dH_{j}/dt_{Transf} expresses the exchange between the helicities H_{j} and H_{pj} without any consequence on the evolution of the total relative helicity H_{V}. Furthermore, this quantity being a volume term, the timevariations of H_{j} and H_{pj} cannot be expressed solely through boundary fluxes. Therefore H_{j} and H_{pj} are not conserved quantities in resistive or ideal MHD.
3. Magnetic energy
3.1. Free and potential energies
With the decomposition of the magnetic field into currentcarrying and potential components, B = B_{p} + B_{j}, the total magnetic energy E_{v} can be also decomposed as
with
In this decomposition, E_{j} is the energy of the currentcarrying magnetic field, also known as free energy, and E_{p} the energy of the solenoidal magnetic field. Because the potential field shares the same surface distribution as the total magnetic field, the surface integral vanishes in Eq. (30). Numerically, the discretization of the mesh grid unavoidably induces a finite level of nonsolenoidality (∇ ⋅ B ≠ 0), and consequently, the last term in Eq. (30) is not exactly null (cf. Valori et al. 2013). However, considering a solenoidal field, Eq. (30) can be simplified into
This decomposition is similar to the decomposition of the helicity obtained in Eq. (5). However, here, the potential energy E_{p}, unlike the volumethreading helicity, only depends on the potential field without a dependence on the nonpotential field.
3.2. Timevariation of the total magnetic energy
We aim to determine the timevariation of the total magnetic energy E_{v} in a fixed volume V,
In the resistive MHD, we use the Faraday law, ∂B/∂t = −∇×E, and we obtain
Using the Gauss divergence theorem, we can decompose the first term of Eq. (34):
Here, the surface term corresponds to the surface integral of the poynting vector and can be divided into two terms:
Finally, assuming for simplicity that the resistivity is constant in space, the variation of the total magnetic energy can be decomposed as
with
We performed a similar time variation for E_{v} as Linan et al. (2018) did for helicities as we summarized in Sect. 2. We find two fluxes: F_{Vn, E} a shearing term associated with horizontal motion, and an emerging term F_{Bn, E} that is related to the emergence. Longcope et al. (2007) decomposed the shearing term at the photospheric level into two contributions by differentiating the motion between the different flux patches and the spin motion of isolated flux patches. As with H_{j} and H_{pj}, the timevariation of E_{v} cannot be expressed through boundary fluxes, and thus E_{v} is not a conserved quantity. Even in ideal MHD, when dE_{v}/dt_{Diss} is null, a volume term dE_{v}/dt_{var} survives. Unlike dH_{j}/dt and dH_{pj}/dt, the timevariation of the total energy E_{v} depends on neither A nor A_{p}: each term of Eq. (37) is gauge invariant.
3.3. Timevariation of the potential and free magnetic energies
Similarly to the analysis of dE_{v}/dt in the previous section, it is possible to obtain the timevariation of E_{p}. Using the scalar potential ϕ of B_{p} such as ∇ϕ = B_{p} and the Gauss divergence theorem, we write
with
For a purely solenoidal potential field, dE_{p}/dt_{ns} is null and thus the timevariation of E_{p} is written in a simple way as a single surface term depending on the variation of the potential field. A decrease in potential energy can therefore be associated in particular with a cancelation at the polarityinversion line (PIL) or a dispersion of the potential magnetic field. Unlike E_{v}, E_{p} is a conserved quantity in resistive MHD.
For the study of the evolution of nonpotential energy E_{j}, the easiest way is to consider only the difference between the decompositions obtained from Eqs. (37) and (42):
Still, this identity is truly accurate only in the case of purely solenoidal magnetic fields, that is, ∇ ⋅ B_{j} = 0 and when B_{j}_{S} = 0 (cf. Sect. 3.1). Hereafter, we therefore introduce a generic formulation where terms that explicitly account for nonsolenoidal errors are retained,
The last two terms on the right are very similar to dE_{p}/dt and dE_{v}/dt. They can thus be decomposed in the same way,
and
By definition, the curl of the potential field is null, and thus the second volume integral on the righthand side formally vanishes. Finally, by grouping all the terms, the timevariation of the nonpotential magnetic energy can be written as
with
As expected, the decomposition we obtain is similar to the one obtained with dE_{v}/dt (see Eq. (37)). The dissipation term dE_{j}/dt_{ns} is null if the magnetic field is solenoidal. However, we compute it in order to quantify the effect of purely numerical errors.
4. Linetied eruptive simulations
In order to comparatively analyze the evolution of helicities and energies and to study the timevariation of the energy, we used magnetic field data produced by parametric 3D MHD simulations of eruptive events of the solar corona that were initially presented in Zuccarello et al. (2015).
For this set of simulations, the OHMMPI code (Aulanier et al. 2005) solves MHD equations in the system’s nondimensional units for a volume covering the domain x ∈ [ − 10, 10], y ∈ [ − 10, 10] and z ∈ [0, 30]. The employed mesh is nonuniform and the smallest cell is centered at x = y = z = 0. In order to facilitate the computation of the energies and helicities, our study was performed on a subdomain excluding the z = 0 plane, interpolated into a uniform Cartesian grid composed of 333 cells in the x and y direction, and 500 cells in the z directions. The resulting analyzed volume is x ∈ [ − 10, 10], y ∈ [ − 10, 10] and z ∈ [0.006, 30]. As a result of the interpolation on a uniform grid (whose cell sizes are 0.06, to be compared to the original grid, whose cell sizes range from 0.006 to 0.32), the magnetic field we obtained has a lower solenoidality than the initial grid. This reduces the accuracy of the magnetic helicity and energy computations (Valori et al. 2016).
The system is delimited by a set of boundaries subject to “open” boundary conditions (except at z = 0). In the analyzed datasets, the bottom boundary is at z = 0.006, one mesh point above the surface corresponding to the photospheric level where linetied boundary conditions were imposed in the original numerical experiments. All the physical MHD quantities, such as the magnetic field, can leave the simulation domain through lateral and top boundaries during the evolution of the system.
Four parametric simulations were performed, all starting with a common phase that is referred to as the “shearing phase” in Zuccarello et al. (2015, 2018). Initially, the magnetic field is potential and generated by two unbalanced and asymmetric subphotospheric polarities. During the shearing phase, asymmetric vortices centered around the local maxima of B_{z}(z = 0) slowly evolve the initial potential magnetic field into a currentcarrying magnetic field. This shearing flow motion induces a magnetic shear close to the PIL and at the end of this phase, creates a currentcarrying magnetic field arcade surrounded by a quasipotential background field. During this entire phase, the distribution of B_{z} at the bottom boundary remains unchanged. This phase lasts from the time t = 10t_{A} until t = 100t_{A} in the system time coordinate, where t_{A} is the reference Alfvén time used in Zuccarello et al. (2015).
Then the four parametric simulations differ by the motion pattern that is imposed at the bottom boundary (cf. Fig. 1). We refer to this phase as the “preeruption” phase. In a first motion profile, labeled “convergence” (Fig. 1, left panel), the velocity flows only follow the horizontal direction x and are only applied close to the PIL. This creates a cancellation of the magnetic flux around the PIL but only slightly affects the periphery of the active regions. Unlike the previous case, for the run labeled “stretching” (Fig. 1, middle left panel), these horizontal motions are also applied at the periphery of the active region. For the other two runs, labeled “dispersion central” and “dispersion peripheral” (Fig. 1, middle right and right panels), the motions spread in all directions from the center. The only difference is in the portion of the active region that is subjected to these motions. In the dispersion peripheral run, only the periphery of the active region is concerned, while in dispersion central, the dispersion also occurs in the center of the polarity where the magnetic field is strongest.
Fig. 1. Applied boundarydriving motions for the four different numerical experiments. White represents the positive polarity (B_{z}(z = 0.006) > 0) and black the negative polarity (B_{z}(z = 0.006) < 0). Orange and cyan arrows indicate the distribution of the velocity flows we applied to the negative and positive polarity, respectively. 
The four runs all present a cancellation of magnetic flux at the PIL that is permitted by a finite photospheric diffusion. The shearedarcade configuration at the end of the shearing phase evolves into a baldpatch topology, and the magnetic reconnection process leads to the formation of a flux rope. The system then evolves until it reaches the instant where it becomes unstable and erupts (cf. Aulanier et al. 2009, for a description of the eruption process). The onset of the eruption, that is, the time t_{1} of the onset of the instability, is accurately determined by a series of relaxation runs for each simulation (Zuccarello et al. 2015). It occurs at t_{1} = 196, 214, 220 and 164t_{A} for the convergence, stretching, dispersion peripheral, and dispersion central runs, respectively.
To ensure the numerical stability of the code, a finite resistivity η and a pseudoviscosity ν are necessary (Zuccarello et al. 2015). During the common shearing phase until t < 100t_{A} and during the preeruption phase, the coronal diffusivity is η_{cor} = 4, 8 × 10^{−4} and the pseudoviscosity is fixed to ν′ = 25. After the eruption, during a phase referred to as eruption phase, η_{cor} is 2, 1 × 10^{−3} and ν′ is 41, 7. To allow later flux cancellation at the PIL, a photospheric resistivity η_{phot} = η_{cor} = 4, 8 × 10^{−4} is imposed only during the preeruption phase. The photospheric resistivity is set to zero in the shearing phase and before the eruption. The change in resistivity follows a rampdown time profile during the time t_{1} − 5t_{A} < t < t_{1} + 5t_{A}. This transitional period is called “eruption onset phase” and is represented as the yellow band in all the figures. We note that the time t_{1} corresponds to the middle of the rampdown time profile, therefore the boundary flows are null only at t > t_{1} + 5t_{A}. In this paper, we removed the first mesh point in the zdirection, which corresponds to the bottom boundary level. We therefore have a uniform resistivity throughout the domain in order to facilitate the calculation of the socalled “nonideal” terms.
5. Energy and helicity evolutions
In this section we first introduce the method for numerically computing the energies and helicities in our set of simulations. Then, we discuss the computation of the timevariations.
5.1. Energy and helicity estimations
In order to compute the different helicities and energies at each time t_{A}, we used the method of Valori et al. (2012). The datacubes of the magnetic field B, of the plasmavelocity field v, and of the plasma thermodynamic quantities allowed us to compute all the quantities that appear in Eqs. (8), (19), (48), and (42).
First the scalar potential ϕ(t) of the potential magnetic field B_{p}(t) was obtained from a numerical solution of the Laplace Eq. (3). The numerical methods we employed to solve this equation required an uniform grid and thus led to the interpolation of the initial grid, as mentioned in Sect. 4. The potential vectors A(t) and A_{p}(t) were computed according to Eq. (14) in Valori et al. (2012) and follow the DeVoreCoulomb gauge defined in Pariat et al. (2015):
This choice of gauge was complemented by the following relationship inherent in the integration method:
where ⊥ means the normal component. It corresponds to the 1D integration of magnetic fields starting at the top boundary of the domain at height z_{top}. Finally, we obtained the helicities and energies from Eqs. (5) and (30). In Fig. 2 we plot the time evolution of these quantities for the different runs. In order to facilitate the comparison between the different runs, the time is plotted with a modified time variable t − t_{1} in each figure, where t_{1} is the onset time defined in Sect. (4) and is different for each of the four simulations.
Fig. 2. Evolution of the different magnetic helicities (top panels), from left to right: relative magnetic helicity (H_{v}, Eq. (4)), volumethreading helicity (H_{pj}, Eq. (7)), and currentcarrying helicity (H_{j}, Eq. (6)). Time evolution of the different magnetic energies (bottom panel), from left to right: total magnetic energy (E_{v}, Eq. (32)), potential energy (E_{p}, Eq. (31)), and free energy (E_{j}, Eq. (31)). The different simulations are dispersion central (red line), dispersion peripheral (green line), stretching (yellow line) and convergence (blue line). The yellow vertical band corresponds to the onset phase of the eruption. 
Unlike Zuccarello et al. (2018), we are interested here in the evolution of the quantities after the common shearing phase, including the eruption phase (which was not studied by Zuccarello et al. 2018). We also recall that the domain studied here is slightly different from the one studied in Zuccarello et al. (2018) (cf. Sect. 4). Figure 2 shows that the dynamics of energies and helicities are qualitatively similar from one simulation to the other during the three different phases. The different boundaryforcing mainly affects the magnitude of the different quantities, but not the quality of their dynamical behaviors.
We also note that the dynamic of helicities and energies changes during the eruption. For the currentcarrying helicity, H_{j}, and the free energy, E_{j}, we observe an overall increase during the preeruption phase, followed by a decrease during the eruption phase. The quantities related to the potential magnetic fields, E_{p} and H_{pj}, both decrease in the preeruption phase. In the eruption phase E_{p} remains constant while H_{pj} continues to decrease, although at a lower rate than in the preeruption phase.
Because both H_{j} and H_{pj} decrease, the relative helicity H_{v} also decreases in the eruption phase. The total energy, E_{v}, has a dynamics similar to H_{v}: the system loses energy throughout the simulation, but at a different rate before and after the eruption.
Overall, the behavior of the helicities is similar to their energy counterparts, for example, H_{v} to E_{v}, H_{j} to E_{j}, H_{pj} to E_{p}. This is particularly visible for the currentcarrying helicity and the free energy: when H_{j} increases (decreases), E_{j} also increases (decreases). In the next sections we focus more on the reasons of these trends by studying the fluxes of the different quantities.
5.2. Timevariation estimation
After all the vectors were calculated, we computed the instantaneous timevariations of energies and helicities obtained from Eqs. (8), (19), (48), and (42). The surface integrals were calculated systematically as the sum of the contributions from the six boundaries. A study focusing on the contribution of the lower boundary alone is conducted in Sect. 9.3.
Linan et al. (2018) validated the timevariations equations established for the volumethreading helicity, H_{pj}, and the currentcarrying helicity, H_{j}. The accuracy of the computation is related to difference factors such as the discretization and the remapping of the data, spatial and temporal, as well as to nonexplicit numerical diffusive terms that are not accounted for in our analytical resistive MHD model. In our study, we present a complementary test by comparing the timevariation of the relative helicity, dH_{v}/dt computed from Eq. (23) in Pariat et al. (2015), with the sum of the timevariations of nonpotential and volumethreading helicities, dH_{j}/dt + dH_{pj}/dt from (8) to (19). In this way, both sides are computed with the same temporal accuracy.
The result of this comparison is presented in the right panel of Fig. 3. In this figure we plot dH_{v}/dt and dH_{j}/dt + dH_{pj}/dt for the dispersion peripheral simulation. For the other runs, the difference is on the same order of magnitude and varies in a similar way. We therefore do not plot this here. The difference is very low, with an average deviation smaller than 0.1%. This confirms the robustness of our calculation method as well as the validity of our analytical equations.
Fig. 3. Left panel: time evolution of the instantaneous timevariation, dE_{v}/dt (dashed black curves, Eq. (37)), and of the sum, dE_{j}/dt + dE_{p}/dt (continuous blue curves, Eqs. (8) and (19)) for the dispersion peripheral run. Right panel: time evolution of the instantaneous timevariation, dH_{v}/dt (dashed black curves, Eq. (23) in Pariat et al. 2015, and of the sum, dH_{j}/dt + dH_{pj}/dt (continuous blue curves, Eqs. (8) and (19)) for the dispersion peripheral run. The yellow bands correspond to the eruption onset phase. 
In the same way, we compare in Fig. 3 (left panel) the timevariation of the total magnetic energy, dE_{v}/dt, computed from Eq. (37), with the sum of the timevariations of free and potential energies, dE_{j}/dt + dE_{p}/dt, from Eqs. (42) to (48). Here, the difference is not negligible, with an average relative difference of 27%. The cause of this difference is likely mainly the nonsolenoidality of the magnetic field. As mentioned Sect. 3, with a finite level of nonsolenoidality, the equality of Eq. (45) is not fully accurate because E_{v} ≠ E_{j,s} + E_{p} (see Sect. 3.1).
In order to estimate the artificial nonsolenoidal energy contributions, Valori et al. (2013) introduced the following decomposition:
where E_{j,s} and E_{p, s} are the energies of the currentcarrying and potential solenoidal magnetic field. E_{j,ns} and E_{p, ns} are the nonsolenoidal components, whereas E_{mix} corresponds to all the remaining cross terms. For a solenoidal field we have E_{j,ns} = E_{p, ns} = E_{mix} = 0, E_{j,ns} = E_{j}, E_{p, ns} = E_{p}, and therefore dE_{j}/dt + dE_{p}/dt = dE_{v}/dt.
The finite nonsolenoidality (∇ ⋅ B ≠ 0) affects the precision of the helicity and energy computations, as studied by Valori et al. (2016). Following Valori et al. (2013, 2016), Thalmann et al. (2019), we used the energy criteria E_{div}/E_{v} to quantify the nonsolenoidality effect. The divergencebased energy is defined as
For the four simulation analyzed in this paper we find an average of E_{div}/E_{v} ≃ 2%. According to Valori et al. (2016), this value for the average of nonsolenoidality leads to a precision of ≤6% for our helicity computations, which is much lower than the 27% discrepancy found in Fig. 3. One possible cause of nonsolenoidality is our interpolation of the original data from a highly nonuniform mesh onto a uniform grid, which can increase the nondivergence of the magnetic field. Different tests have been made to degrade and also improve the interpolation to give a rough estimate of the effect. Finer interpolations, however, have required considering only a fraction of the whole numerical domain to keep the number of grid points manageable. The outcome of these tests is that neither presented results that differed significantly from our baseline interpolation. In particular, E_{div}/E_{v} always remained on the order of 2%, which means that this error therefore seems intrinsic to the numerical models. The level of interpolation chosen in our study therefore is a good compromise between the required computed power and the quality of our data. In addition, it is worth noting that various terms in our equations depend on time variations, so that the accuracy of our flux computations can also be limited by a relative coarseness of the time outputs of the available data. Testing for this would require recalculating the simulations with a higher cadence for its outputs, which is beyond the scope of this paper.
6. Helicity transfer
As mentioned in the introduction (see Sect. 1), Linan et al. (2018) showed for two different eruptive simulations that the exchange between H_{pj} and H_{j} is controlled by the gaugeinvariant transfer term dH_{j}/dt_{Transf} (see Eq. (10)), which therefore plays a key role in evolving these helicities. In order to confirm these different results in the particular case of our linetied simulations, we plot in Fig. 4 the gaugeinvariant terms of dH_{j}/dt (see Eq. (17)), and dH_{pj}/dt (see Eq. (28)) for the convergence and the dispersion peripheral runs.
Fig. 4. Time evolution of the helicity variation rates, dH_{j}/dt and dH_{pj}/dt (dashed black curves; Eqs. (8) and (19)), of the helicity transfer term, dH_{j}/dt_{Transf} and dH_{pj}/dt_{Transf} (solid red curves; Eqs. (10) and (21)), of the own terms, dH_{j}/dt_{Own} and dH_{pj}/dt_{Own} (solid blue curves; Eqs. (18) and (29)), and of the dissipation terms, dH_{j}/dt_{Diss} and dH_{pj}/dt_{Diss} (solid green curves; Eqs. (9) and (20)) for the convergence simulation (top panels) and for the dispersion peripheral simulation (bottom panels). Left and right panels: evolution of the current carrying helicity, H_{j}, and volumethreading helicity H_{pj}, respectively. 
In the convergence case (see Fig. 4, top left panel), the conversion of helicity from H_{pj} to H_{j} and the boundary flux have similar amplitudes during the preeruption phase. Both contributions are positive, thus H_{j} increases (e.g., dH_{j}/dt is positive). The transfer term, dH_{j}/dt_{Transf}, dominates the evolution of H_{j} when it is close in time to the eruption. Meanwhile, H_{pj} decreases mainly because of the strong helicity transfer (cf. Fig. 4, top right panel), that is, dH_{j}/dt_{Transf} is the dominating term of dH_{pj}/dt. In comparison, the flux of H_{pj} related to the own term, dH_{pj}/dt_{Own}, is very low. This means that the injection of H_{pj} is not enough to compensate for its conversion to H_{j}.
Moreover, as with the flux emergence simulations presented in Linan et al. (2018), the eruption is accompanied by a sharp decrease of the transfer term. However, unlike in Linan et al. (2018), the sign of dH_{j}/dt_{Transf} does not change here during the onset phase of the eruption. During the eruption phase, the dissipation terms, dH_{j}/dt_{Diss} and dH_{pj}/dt_{Diss}, dominate the variation of the helicities. This is related to the increase in resistivity η imposed in the numerical experiment during that period.
In the dispersion peripheral simulation (see Fig. 4, bottom panels), the variations of H_{j} and H_{pj} during the preeruption phase are noticeably dominated by the boundary fluxes, dH_{j}/dt_{Own} and dH_{pj}/dt_{Own}. The transfer terms are significantly less intense than the injection of H_{j} and H_{pj}, except close to the eruption. During the eruption phase, the dynamics is mostly dominated by the resistive dissipation terms, dH_{j}/dt_{Diss} and dH_{pj}/dt_{Diss}.
We thus observe that while the trends of H_{j} and H_{pj} are similar for the convergence and the dispersion peripheral simulations, as discussed in Sect. 5, these variations are in fact due to noticeably different dynamics of the helicity fluxes. For instance, the decrease of H_{pj} in the preeruption phase is due to an intense conversion of helicity (high values of dH_{pj}/dt_{Transf}) for the dispersion peripheral simulation, whereas in the convergence run, a similar evolution of H_{pj} is explained by an intense negative boundary flux, dH_{pj}/dt_{Own} during that period.
Finally, as has been noted in Linan et al. (2018), for both the simulations analyzed here but also for the two others, the transfer terms cannot be neglected. We confirm that the estimations of boundary fluxes are not sufficient to follow the dynamics of H_{j} and H_{pj}. However, unlike with the flux emergence and solar jet simulations studied in Linan et al. (2018), the precise mechanism of the buildup of H_{j} and H_{pj} does depend on the simulations during the preeruption phase. Studying this dependence is the goal of the next section.
7. Distinguishing between simulations in terms of helicity dynamics
In the previous section, we described that the magnitude of the different gaugeinvariant helicity variation terms could be significantly different in the dispersion peripheral and in the convergence simulations. This demonstrates that even if the general trends of H_{pj} and H_{j} are similar (see Sect. 5), their dynamics can be significantly different.
In order to estimate the effect of the different boundary driver, Fig. 5 displays the different gaugeinvariant variation terms for the four different numerical simulations: the boundary fluxes dH_{j}/dt_{Own} and dH_{pj}/dt_{Own}; the transfer term, dH_{j}/dt_{Transf}; and the dissipations terms dH_{j}/dt_{Diss} and dH_{pj}/dt_{Diss}.
Fig. 5. Time evolution of the different gauge invariant terms of dH_{j}/dt (top panels), from left to right: dissipation term (dH_{j}/dt_{Diss}, Eq. (9)), own term (dH_{j}/dt_{Own}, Eq. (18)), and helicity transfer term (dH_{j}/dt_{Transf}, Eq. (10)). Time evolution of the different gauge invariant terms of dH_{pj}/dt (bottom panels), from left to right: dissipation term (dH_{pj}/dt_{Diss}, Eq. (20)), and own term (dH_{pj}/dt_{Own}, Eq. (29)). Each curve color corresponds to a particular simulation: dispersion central (red line), dispersion peripheral (green line), stretching (yellow line), and convergence (blue line). The yellow band corresponds to the onset phase of the eruption. 
The dissipations terms (cf. Fig. 5, left panels) are not significantly different from one simulation to the other. The sudden increase in absolute values of the dissipations terms, observed during the eruption onset phase, is related to the imposed numerical increase in resistivity. The variations in magnitude, particularly in the preeruption phase, are minor compared to the variations in other gaugeinvariant terms.
The boundary flux term dH_{j}/dt_{Own} is also very similar from one simulation to another, except for the dispersion central run, for which more nonpotential helicity is markedly injected during the preeruption phase (see Fig. 5, middle top panel). Unlike dH_{j}/dt_{Own}, the boundary flux of H_{pj}, dH_{pj}/dt_{Own}, is strongly sensitive to the boundarydriving pattern (see Fig. 5, bottom left panel). Tthe sign and magnitude of dH_{pj}/dt_{Own} depend on the simulation. For the dispersion simulations, there is a significant injection of negative H_{pj}, whereas in the convergence and stretching case, the flux is significantly weaker, if not of the opposite sign.
Finally, Fig. 5 (top right panel) shows that the helicity transfer rate, dH_{j}/dt_{Transf}, is higher for the convergence and stretching simulations than for the dispersion cases in the preeruption phase. Unlike with the other cases where the transfer term is almost constant during an early period, in the dispersion central run dH_{j}/dt_{Transf} increases from the first moments of the simulation.
In summary, we observed that during the preeruption phase, the increase of H_{j} and reciprocally the decrease of H_{pj} (cf. Fig. 5) are not explained by the same physical process in the different simulations. We observe three significantly different dynamics:

The convergence and stretching simulations present a similar dynamics for their fluxes of helicity. They are characterized by a relatively weak boundary flux of H_{pj} counterbalanced by a strong transfer from H_{pj} to H_{j}. The own term of H_{j} is positive, involving an injection of currentcarrying helicity. Its magnitude is almost identical in these two runs.

For the dispersion peripheral run, H_{pj} decreases mostly because of the boundary flux, unlike with the previous cases. In comparison to the boundary flux, the transfer from H_{pj} to H_{j} is less important. The flux of H_{j} is similar to the convergence and stretching runs.

The dispersion central shares some similarities with the dispersion peripheral run regarding the variations of H_{pj}. However, this simulation is characterized by a high boundary flux of H_{j} that is distinct and significantly higher than the three other cases.
Finally, simulations with the largest injection of helicities due to their own terms (whether H_{j} or H_{pj}) have the lowest magnitude of the transfer term. Inversely, a strong exchange between H_{j} and H_{pj} is accompanied by lower fluxes through the surfaces. Both lead to a similar trend for the relative helicity H_{v}. This shows that the boundary fluxes of H_{j} or H_{pj} as well as the volume term, dH_{j}/dt_{Transf}, are directly related to the morphology and the evolution of the magnetic field at the bottom boundary. In Sect. 9.3 we discuss that a specific boundarydriven pattern may influence the different physical mechanisms of the evolution of the magnetic helicities.
8. Distinguishing between simulations in terms of energy dynamics
As shown in Sect. 5, the evolutions of H_{j} and E_{j} are very similar. Likewise, H_{pj} and E_{pj} evolve in the same way during the preeruption phase. The main difference appears after the eruption, where H_{pj} still decreases while E_{pj} remains constant. However, the similar overall behaviors of H_{pj} and H_{j} hide very different physical mechanisms, depending on the simulation, as shown in the previous section. We identified three types of evolution for the dynamics of the helicities. In this section we focus on the following questions: how do E_{pj} and E_{j} evolve? Does the dynamics of the energy fluxes also distinguish between the different simulations, as the helicity dynamics do?
For this purpose, we present in Fig. 6 all the flux that appear in the decomposition of dE_{p}/dt (see Eq. (42)) for the convergence (Fig. 6, left panel) and the dispersion central simulations (Fig. 6, middle panel). In both simulations, the nonideal term is almost null because of a very low nonsolenoidality of the potential field, for instance, ∇ ⋅ B_{p} ≃ 0. The evolution of the potential energy therefore depends only on F_{ϕ, Bz}, which results from the evolution of the magnetic field at the boundaries. During the preeruption phase, the magnitude of F_{ϕ, Bz} as the relative change of the magnetic field at the boundary becomes weaker. Then, during the eruption onset phase, F_{ϕ, Bz} decreases strongly before becoming null during the eruption phase. The magnetic field is indeed kept fixed at the bottom boundary during that period.
Fig. 6. Left and middle panel: time evolution of the potential energy variation term (dashed black line; dE_{p}/dt; Eq. (42)) and the different terms constituting the instantaneous timevariation of E_{p} (Eq. (42)): F_{ϕ, Bz} (blue line; Eq. (43)), and dE_{p}/dt_{ns} (orange line; Eq. (44)). Left and middle panels: evolution for the convergence and dispersion central simulation, respectively. Right panel: time evolution of F_{ϕ, Bz} for the four simulations dispersion central (red line), dispersion peripheral (green line), stretching (yellow line), and convergence (blue line). The yellow band corresponds to the onset phase of the eruption. 
From comparing F_{ϕ, Bz} in the different simulations, we note that F_{ϕ, Bz} presents the same evolution for all simulations except for the dispersion central (see Fig. 6, right panel). This indicates that except for the dispersion central simulation, the differences in boundarydriven motions do not affect the injection of the potential energy (cf. the discussion in Sect. 9.3). However, this run has the same functional form as the other, is only more efficient, and therefore quicker, in achieving the eruption.
Regarding E_{p}, only the dispersion central run presents a different behavior. The same conclusion was obtained for the fluxes of H_{j}, for instance, dH_{j}/dt_{Own} (see Fig. 5, middle top panel), but not for dH_{pj}/dt_{Own} (cf. Fig. 5l). First, F_{ϕ, Bz} is negative for the entire simulation, while the sign dH_{pj}/dt_{Own} depends on the simulation. This confirms that there is no direct link between the dynamics of E_{p} and the injection of H_{pj}.
In Fig. 7 we observe the different terms of dE_{j}/dt (Eq. (48)) for the four simulations. Unlike dH_{j}/dt and dH_{pj}/dt, the trends and dominant terms of dE_{j}/dt are similar in the simulations. Only the magnitude of each term may differ. Before the eruption, dE_{j}/dt is dominated by the emergence term F_{Bn, Ej} despite a significant magnitude of the dissipation term, dE_{j}/dt_{Diss}. Then, during the eruption, F_{Bn, Ej} becomes null as a result of the interruption of the boundarydriving flows. During the eruption phase, dE_{j}/dt is negative and dominated by dE_{j}/dt_{Diss}. The free energy mainly decreases because it is dissipated and not ejected out of the volume.
Fig. 7. Time evolution of the freeenergy variation rate (dashed black line; dE_{j}/dt; Eq. (48)) and the different terms constituting the instantaneous timevariation of E_{j} (Eq. (48)): dE_{j}/dt_{Diss} (blue line; Eq. (49)), dE_{j}/dt_{Var} (orange line; Eq. (50)), F_{Bn, Ej} (green line; Eq. (51)), F_{Vn, Ej} (red line; Eq. (52)), and dE_{j}/dt_{ns} (purple line; Eq. (54)). Each panel corresponds to a different simulation: convergence (top left), stretching (top right), dispersion peripheral (bottom left), and dispersion central (bottom right). The yellow band corresponds to the onset phase of the eruption. 
The dissipation term, dE_{j}/dt_{Diss} does not vary much between the simulations (see Fig. 8, bottom right panel). Similarly, the differences of F_{Vn, Ej} and dE_{j}/dt_{Var} are small between the runs during the preeruption phase. Only F_{Bn, Ej} (see Fig. 7, top left panel) presents significant differences in the simulations that affect the evolution of the free energy, E_{j}. The dispersion central simulation presents a distinctive trend. The magnitude of F_{Bn, Ej} starts very high and then decreases to values similar to the other runs during the onset phase of the eruption.
Fig. 8. Time evolution of the different gaugeinvariant terms of dE_{j}/dt: F_{Bn, Ej} (top right panel, Eq. (51)), F_{Vn, Ej} (top left panel, Eq. (52)), dE_{j}/dt_{Var} (bottom left panel, Eq. (50)), and dE_{j}/dt_{Diss} (bottom right panel, Eq. (49)). The different colors present one simulation: dispersion central (red line), dispersion peripheral (green line), stretching (yellow line), and convergence (blue line). The yellow band corresponds to the onset phase of the eruption. 
Unlike the evolution of the helicities, H_{j} and H_{pj}, only the dispersion central simulation stands out from the other runs. This simulation is characterized by a higher decrease of the potential field (see Fig. 6, right panel) and by a higher initial injection of E_{j} caused by F_{Bn, Ej} (see Fig. 8, top left panel). Before the eruption, another difference with the helicities is that the variations of the trend of E_{j} and E_{p} are purely related to the boundary fluxes. Finally, one key outcome of our study is that the dynamics of the energy fluxes do not distinguish between the simulations, unlike the helicity fluxes.
This shows that even if volume helicities and energies follow similar trends (cf. Fig. 2), the physical mechanisms that drive their dynamics are very different. First, the evolution of free energy, E_{j}, and potential energy, E_{p}, are independent, while the currentcarrying helicity, H_{j}, evolves in a correlated way with the volumethreading helicity, H_{pj}. Additionally, different boundary forcing only affects the magnitude of the energy fluxes. The dynamics of the helicity is more complex and varies drastically from one simulation to another. One group of simulations (dispersion central and peripheral) is dominated by the flux through the surfaces, while a second group (convergence and stretching) is controlled by volume exchange within the domain. We conclude that energy, helicity, and their decompositions have distinct properties whose analysis should be complementary for the study of the eruptivity of active regions.
9. Discussion
9.1. Summary
In Sect. 2 we introduced the formulation of the magnetic relative helicity, H_{v}, as well as the formulation of its decomposition into the currentcarrying helicity, H_{j}, and the volumethreading helicity, H_{pj}. We also recalled the analytical equations of their timevariations obtained in Linan et al. (2018). Similarly, we introduced the decomposition of magnetic energy, E_{v}, into the sum of the potential energy, E_{p}, and the free energy, E_{j}. Then, we obtained the timevariation of E_{v} (see Eq. (37)), E_{p} (see Eq. (42)), and E_{j} (see Eq. (48)) by analytically deriving their time derivative (see Sect. 3). These formulae are valid for any gauge choices and in the presence of finite level of nonsolenoidality for the magnetic field.
Our numerical study of timevariations of energies and helicities is based on a series of four eruptive numerical MHD simulations of solar active regions (see Sect. 4) that have been investigated in Zuccarello et al. (2015). The evolution of each simulation is characterized by different boundary forcing (linetied) until the eruption (see Sect. 4). After the same shearing phase, four driving photospheric flows were considered: convergence, stretching, and peripheral and central dispersion flows. In this study we were particularly interested in the fluxes of energies and helicities during the flux rope formation, during the eruption onset phase, when the torus instability occurs, and during a short time interval after the eruption onset, called the eruption phase.
Initially, the magnetic energy, E_{v}, decreases as a result of the decrease in the potential magnetic field, E_{p}, despite the increase in free energy, E_{j}. At the same time, the decrease in volumethreading helicity, H_{pj}, compensates for the injection of the currentcarrying helicity, H_{j}, which leads to a quasiconstant evolution of the relative helicity H_{v} (see Sect. 5). The relative helicity was mostly injected during the earlier shearing phase.
The fluxes of free and potential energies, E_{j} and E_{p} (see Sect. 8) showed that the effectiveness of the buildup of free energy within the domain is purely related to the magnitude of one surface term. Similarly, the evolution of potential energy is fully linked with its boundary fluxes.
We also used these simulations to investigate the importance of the exchange between H_{j} and H_{pj}, which was previously highlighted by Linan et al. (2018). The exchange of helicity between H_{j} and H_{pj} is controlled by a gaugeinvariant term, dH_{j}/dt_{Transf} (cf. Eq. (10)). As in Linan et al. (2018), we observed that this term plays a key role in the dynamics of both H_{j} and H_{pj}, in particular during the buildup phases where the transfer terms, dH_{pj}/dt_{Transf}, dominate the evolution of dH_{pj}/dt for the runs convergence and stretching. In the dispersion simulations, the evolutions of dH_{j}/dt and dH_{pj}/dt are dominated by their boundary fluxes, dH_{pj}/dt_{Own} and dH_{pj}/dt_{Own} (cf. Eqs. (18) and (29)), even though the magnitude of the transfer term remains significant. This means that neither H_{j} nor H_{pj} evolve as a result of boundary fluxes alone. This conclusion is consistent with the results of Linan et al. (2018) that were obtained for different numerical experiments. Specifically, H_{j} and H_{pj} cannot be estimated in observed active solar active regions by time integration of its flux through the solar photosphere, but rather with a volumeintegration method (Valori et al. 2016). This approach requires a 3D reconstruction of the coronal magnetic field from the 2D photospheric measurement with coronal field extrapolation techniques (Wiegelmann & Sakurai 2012; Wiegelmann et al. 2014). A more detailed discussion of the effect of the transfer term on the estimation of H_{j} and H_{pj} can be found in the conclusion of Linan et al. (2018).
The key outcome of the study is the observation that the dynamics of the transfer and fluxes of H_{j} and H_{pj} depend on the simulation and thus on the imposed driving motions, even though the variations in H_{j} and H_{pj} seemed relatively independent of the simulation setup (see Sect. 6). Despite the four boundary forcings, the different simulations remain very similar in terms of the magnetic field topology. Nonetheless, the dominant terms of dH_{j}/dt and dH_{j}/dt are not the same from one simulation to another.
We highlighted three distinct types of dynamics of the evolution of the helicities in the simulations. In the convergence and stretching simulations, H_{pj} does not evolve as a result of boundary fluxes, but because of its conversion into H_{j}. The opposite is observed during the two dispersion simulations, for which the evolution of H_{pj} is mainly related to its fluxes through the boundary, with a weak transfer term. The dispersion central simulation stands out from the others because its boundary flux of the currentcarrying helicity, H_{j}, is significantly higher.
Thus we were able to process several photospheric forcings to approach the diversity of active regions that are observed at the solar photosphere. We have come to the conclusion that the evolution of helicity during the formation of a flux rope is a complex process whose origin can be related to fluxes through the surface as well as to volume contributions.
9.2. Buildup of the helicity ratio
Zuccarello et al. (2018) have shown that the trigger of the eruptions is related to a threshold in the helicity ratio H_{j}/H_{v}: this ratio reaches the same value, H_{j}/H_{v}_{thresh}, for all simulations at the onset of the eruptions. In our runs, this threshold is 0.29 ± 0.01. However, as discussed in Sect. 7 of Zuccarello et al. (2018), the exact value of this threshold needs to be taken with care because relative helicity is not simply an additive quantity. We investigated how this helicity ratio is built up and eventually reached by studying the specific dynamics of H_{j} and H_{pj} (see Sects. 6 and 7). Despite the different boundary forcings, the simulations are very similar, so that it might have been thought that the fluxes of H_{j} and H_{pj} would also be similar. However, the key outcome of our study is that the terms that dominate the evolution of dH_{j}/dt or dH_{pj}/dt sensitively depend on the simulation even if the overall trends are the same (H_{j} increases and H_{pj} decreases, see Sect. 5).
Three very distinct ways to reach the helicity eruptivity threshold were found. We observed that the eruption was triggered at a specific value of H_{j}/H_{v} independently of the dynamics of H_{j} and H_{pj} to reach this threshold. Our work suggests that active regions could reach an eruptive state, either through strong increases of helicity fluxes or through magnetic configurations that induce strong helicity transfer.
The different ways to reach the helicity eruptivity threshold are not all equally effective. The eruption occurs more or less quickly after the end of the shearing phase. The dispersion central run is the most rapid simulation. Then we find the stretching and convergence runs, and last the dispersion peripheral case. The dispersion central simulation stands out from the others because it presents higher helicity fluxes (due to dH_{j}/dt_{Own} and dH_{pj}/dt_{Own}) and energy fluxes (due to F_{ϕ, Bz} and F_{Bn, Ej}) than the other cases.
Finally, using a set of resistive MHD simulations, we provided new knowledge of the energy and helicity properties. In particular, our analytical and numerical work emphasizes recent studies (Pariat et al. 2017; Zuccarello et al. 2018; Moraitis et al. 2019; Thalmann et al. 2019) that demonstrated how promising the helicity ratio is as a marker of eruptivity. Further studies are still needed, whether to analytically establish the link between the helicity ratio and torus instability or to properly estimate it from data that are measured in the solar atmosphere.
From direct observational data, the evolution of the ratio H_{j}/H_{v} was also analyzed in three active regions with different eruptive profiles. Moraitis et al. (2019) investigated the most active region of cycle 24 (AR 12673), while Thalmann et al. (2019) focused on an eruptive and a confined flare (AR 11158 and AR 12192). These recent observational analyses seem to qualitatively confirm that the H_{j}/H_{v} ratio is tightly related to the eruptivity of solar active regions.
9.3. Effect of the different flows on the helicity and energy injections
A key result of our study is that the specific driving flows that are applied at the bottom boundary are of significant importance on the dynamics of magnetic helicities and energies. They influence the magnitude and the sign of the own terms for the helicity as well as those of the main fluxes of dE_{j}/dt and dE_{p}/dt. Moreover, as mentioned in the previous section, even if the way to reach the helicity eruptivity threshold matters less than reaching the threshold, the spatial velocity and magnetic distributions at the boundary affect the time that is required to reach the threshold.
To reach an understanding of the effect of the linetied forcing, we here briefly discuss the distribution of different quantities at the bottom boundary. A more complete study is beyond the scope of the present paper but will be performed, however.
Our goal is to present quantities that might eventually be obtained from observed photospheric magnetograms. Figure 9 shows four quantities that are related to the energy and helicity fluxes in a (x − y) view at z = 0.006 for all simulations at the same modified time to the eruption (t_{1} − t_{A} = −58): (v ⋅ A_{j})B_{z}, the integrand of F_{Vn, Aj}, for the injection of H_{j} (see the first row in Fig. 9); (v ⋅ A_{j})B_{z}, the integrand of F_{Vn, Ap}, for the injection of H_{pj} (see the second row in Fig. 9); (∂ϕ/∂t)B_{p}, the integrand of F_{ϕ, Bp} (see the third row in Fig. 9); and (v ⋅ B_{j})B_{z}, the integrand of F_{Bn, Ej} (see the fourth row in Fig. 9). The modified time was taken arbitrarily in the preeruption phase.
Fig. 9. From top to bottom: dimensionless magnitude of (v ⋅ A_{j})B_{z}, (v ⋅ A_{p})B_{z}, (∂ϕ/∂t)(B_{p}), and (v ⋅ B_{j})B_{z} viewed in the (x, y) plane at z = 0.006, at the relative time of t − t_{1} = −58. Isocontours of B_{z} (dashed line for negative values, solid line for positive values) correspond to values of B_{z} = − 4.5, −2.0, 0, 2.0, and 4.5. Each column in the panels presents one simulation, from left to right: convergence, stretching, dispersion peripheral, and dispersion central. 
We recall that the injections of H_{j} or H_{pj} related to their own terms cannot be simply reduced to F_{Vn, Aj} and F_{Vn, Ap}. Other terms such as dH_{pj}/dt_{Transf} must be considered. The second issue is that neither F_{Vn, Aj} nor F_{Vn, Ap} are gaugeinvariant quantities. They are still good indicators of the spatial distribution of the helicity injection, however, because the same gauge was used for all four runs. Moreover, the emergence term of the relative helicity is one quantity that has traditionally been investigated at the solar photosphere (Liu & Schuck 2012; Bi et al. 2018). The distributions of F_{Vn, Aj} and F_{Bn, Ej} are very similar (see the first and last row in Fig. 9). This reveals that the injection of free energy appears to be directly connected to the injection of currentcarrying helicity.
The area where F_{Vn, Aj} and F_{Bn, Ej} are intense is along the PIL, except for the dispersion central simulation. In the convergence and stretching runs, the shear is favored because the angle between the velocity and the magnetic field at the PIL is high. For the dispersion peripheral case, the velocity is perpendicular to the PIL. Thus the free energy increases more slowly than in the other cases, and more time is needed to build up the flux rope. We also note a low contribution of F_{Vn, Aj} related to the dispersion of the magnetic field at the periphery of the polarities in the stretching and dispersion peripheral runs.
As highlighted in Sect. 7, the dynamics in the dispersion central simulation stands out from the others. Unlike the other cases, the intense regions of F_{Vn, Aj} and F_{Bn, Ej} are located close to the center of the polarities. The difference between this simulation and the other three appears clearly for F_{ϕ, Bp} (see the third row in Fig. 9). In Sect. 8 we observed that F_{ϕ, Bp} was distinctly higher for the dispersion central than in the other cases. The magnitude of F_{ϕ, Bp} is higher for the dispersion central simulation because the center of the polarity, where the magnetic field is the most intense, is displaced. The positive polarity, which possess a higher magnetic flux than the negative polarity, contributes more to F_{ϕ, Bp}. The three other cases are very similar because the regions subjected to the flow do not influence F_{ϕ, Bp}.
The distribution of the F_{Vn, Ap} term (see the second row in Fig. 9) is markedly different from the other quantities presented in Fig. 9. Two different behaviors appear. For the convergence and stretching runs, like F_{Vn, Aj}, the shear along the PIL is the main contribution of (v ⋅ A_{j})B_{z}. For the dispersion runs, F_{Vn, Ap} is distributed in a symmetric quadrupole at the location where the flow is applied. This leads to an almost null total flux, F_{Vn, Ap}.
Finally, the best way to reach the helicity threshold for these parametric simulations is to facilitate the dispersion of the intense magnetic field. For instance, in the dispersion central case, this corresponds to distributing the flow from the center of each polarity. For the other cases with lower dispersion, the efficient ways to build up the energies and helicities are mostly linked with fluxcancellation converging motions that contain a shearingflow component parallel to the PIL. The convergence and stretching runs have more intense contributions there.
The observations made above offer only a limited insight on the effect of the different flows on the evolution of energy and helicity. We saw that even if the four simulations look similar overall, some differences related to the applied motions clearly appear in the dynamic of the helicity and energy. New investigations are required to provide a better understanding of the energy and helicity dynamics. Our study especially provides new information for interpreting the injection of helicity and energy in observed active regions.
With HMI vector magnetic field data, the evolution of energy and helicity flux has previously been studied in different active regions (Liu & Schuck 2012). This is commonly based on the decomposition of the flux into two components: a shear component provided by photospheric tangential flow, and a vertical component linked with normal flows due to emergence (Bi et al. 2018). At the same time, new methods for properly measuring helicities in the solar corona are still being developed (Dalmasse et al. 2013, 2014, 2018; Valori et al. 2016; Guo et al. 2017; Moraitis et al. 2018; Gosain & Brandenburg 2019). For instance, new analytic expressions for the helicity transport allow us to estimate the injection of relative magnetic helicity into the solar atmosphere over an entire solar cycle (Hawkes & Yeates 2019; Pipin et al. 2019). However, the quality of the helicity and energy estimations from observational data greatly depends on the accuracy of the magnetic field measurements. In particular, it is still difficult to measure the tangential component of the magnetic field. The instrumentation on board the new Solar Orbiter mission, for example, the PHI magnetogram, will help us to make significant progress on the magnetic field measurement and consequently on the helicity and energy computations. Future results that will benefit from our study will provide new insight for a better understanding of the solar coronal activity.
Acknowledgments
The authors thank the anonymous reviewer for careful reading of the manuscript and helpful comments. L.L., E.P., K.M. acknowledge support of the French Agence Nationale pour la Recherche through the HELISOL project ANR15CE310001. L.L., E.P. acknowledge the support of the french Programme National SoleilTerre. GV acknowledges the support from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 824135 and of the STFC grant number ST/T000317/1. The numerical simulation used in this work was performed on the HPC center MesoPSL financed by the project Equip@Meso (reference ANR10EQPX2901) and the Region IledeFrance. This article profited from discussions during the meetings of the ISSI International Team Magnetic Helicity in Astrophysical Plasmas.
References
 Antiochos, S. K. 2013, ApJ, 772, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Aulanier, G., Démoulin, P., & Grappin, R. 2005, A&A, 430, 1067 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2009, ApJ, 708, 314 [Google Scholar]
 Barczynski, K., Aulanier, G., Masson, S., & Wheatland, M. S. 2019, ApJ, 877, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Berger, M. 2003, Fluid Mech. Astrophys. Geophys., 20030424, 23 [Google Scholar]
 Berger, M. A., & Field, G. B. 1984, J. Fluid Mech., 147, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Bi, Y., Liu, Y. D., Liu, Y., et al. 2018, ApJ, 865, 139 [Google Scholar]
 Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Chae, J. 2007, Adv. Space Res., 39, 1700 [NASA ADS] [CrossRef] [Google Scholar]
 Dalmasse, K., Pariat, E., Valori, G., Démoulin, P., & Green, L. 2013, A&A, 555, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dalmasse, K., Pariat, E., Démoulin, P., & Aulanier, G. 2014, Sol. Phys., 289, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Dalmasse, K., Pariat, E., Valori, G., Jing, J., & Démoulin, P. 2018, ApJ, 852, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Del Sordo, F., Candelaresi, S., & Brandenburg, A. 2010, Phys. Rev. E, 81, 036401 [NASA ADS] [CrossRef] [Google Scholar]
 Démoulin, P., & Aulanier, G. 2010, ApJ, 718, 1388 [NASA ADS] [CrossRef] [Google Scholar]
 Elsasser, W. M. 1956, Rev. Mod. Phys., 28, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Finn, J. M., & Antonsen, Jr., T. 1985, Plasma Phys. Control. Fusion, 9, 111 [Google Scholar]
 Forbes, T. G. 2000, J. Geophys. Res., 105, 23153 [NASA ADS] [CrossRef] [Google Scholar]
 Gosain, S., & Brandenburg, A. 2019, ApJ, 882, 80 [NASA ADS] [CrossRef] [Google Scholar]
 Guo, Y., Pariat, E., Valori, G., et al. 2017, ApJ, 840, 40 [NASA ADS] [CrossRef] [Google Scholar]
 Hawkes, G., & Yeates, A. R. 2019, A&A, 631, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 James, A. W., Valori, G., Green, L. M., et al. 2018, ApJ, 855, L16 [NASA ADS] [CrossRef] [Google Scholar]
 Karpen, J. T., Antiochos, S. K., & DeVore, C. R. 2012, ApJ, 760, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Knizhnik, K. J., Antiochos, S. K., & DeVore, C. R. 2015, ApJ, 809, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Kusano, K., Maeshiro, T., Yokoyama, T., & Sakurai, T. 2002, ApJ, 577, 501 [NASA ADS] [CrossRef] [Google Scholar]
 Kusano, K., Maeshiro, T., Yokoyama, T., & Sakurai, T. 2004, ApJ, 610, 537 [NASA ADS] [CrossRef] [Google Scholar]
 Leake, J. E., Linton, M. G., & Török, T. 2013, ApJ, 778, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Leake, J. E., Linton, M. G., & Antiochos, S. K. 2014, ApJ, 787, 46 [NASA ADS] [CrossRef] [Google Scholar]
 Linan, L., Pariat, É., Moraitis, K., Valori, G., & Leake, J. 2018, ApJ, 865, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Linton, M., & Antiochos, S. 2002, ApJ, 581, 703 [NASA ADS] [CrossRef] [Google Scholar]
 Linton, M., Dahlburg, R., & Antiochos, S. 2001, ApJ, 553, 905 [Google Scholar]
 Liu, Y., & Schuck, P. 2012, ApJ, 761, 105 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, Y., & Schuck, P. 2013, Sol. Phys., 283, 283 [NASA ADS] [CrossRef] [Google Scholar]
 Longcope, D. 2004, ApJ, 612, 1181 [Google Scholar]
 Longcope, D., & Beveridge, C. 2007, ApJ, 669, 621 [Google Scholar]
 Longcope, D., Ravindra, B., & Barnes, G. 2007, ApJ, 668, 571 [NASA ADS] [CrossRef] [Google Scholar]
 Low, B. 1996, Sol. Phys., 167, 217 [Google Scholar]
 Moraitis, K., Tziotziou, K., Georgoulis, M., & Archontis, V. 2014, Sol. Phys., 289, 4453 [Google Scholar]
 Moraitis, K., Pariat, É., Savcheva, A., & Valori, G. 2018, Sol. Phys., 293, 92 [NASA ADS] [CrossRef] [Google Scholar]
 Moraitis, K., Sun, X., Pariat, É., & Linan, L. 2019, A&A, 628, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nindos, A., Zhang, J., & Zhang, H. 2003, ApJ, 594, 1033 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pariat, E., Démoulin, P., & Berger, M. 2005, A&A, 439, 1191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pariat, E., Valori, G., Démoulin, P., & Dalmasse, K. 2015, A&A, 580, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pariat, E., Leake, J., Valori, G., et al. 2017, A&A, 601, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pipin, V. V., Pevtsov, A. A., Liu, Y., & Kosovichev, A. G. 2019, ApJ, 877, L36 [Google Scholar]
 Polito, V., Del Zanna, G., Valori, G., et al. 2017, A&A, 601, A39 [Google Scholar]
 Priest, E. 2014, Magnetohydrodynamics of the Sun (Cambridge: Cambridge University Press) [Google Scholar]
 Priest, E. R., Longcope, D., & Janvier, M. 2016, Sol. Phys., 291, 2017 [NASA ADS] [CrossRef] [Google Scholar]
 Rust, D. 1994, Geophys. Res. Lett., 21, 241 [NASA ADS] [CrossRef] [Google Scholar]
 Schrijver, C. J., & Zwaan, C. 2008, Solar and Stellar Magnetic Activity (Cambridge: Cambridge University Press), 34 [Google Scholar]
 Simon, C. 2012, PhD Thesis [Google Scholar]
 Taylor, J. B. 1974, Phys. Rev. Lett., 33, 1139 [NASA ADS] [CrossRef] [Google Scholar]
 Temmer, M., Thalmann, J. K., Dissauer, K., et al. 2017, Sol. Phys., 292, 93 [Google Scholar]
 Thalmann, J. K., Moraitis, K., Linan, L., et al. 2019, ApJ, 887, 64 [NASA ADS] [CrossRef] [Google Scholar]
 Valori, G., Démoulin, P., & Pariat, E. 2012, Sol. Phys., 278, 347 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valori, G., Demoulin, P., Pariat, E., & Masson, S. 2013, A&A, 553, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Valori, G., Pariat, E., Anfinogentov, S., et al. 2016, Space Sci. Rev., 201, 147 [NASA ADS] [CrossRef] [Google Scholar]
 Welsch, B., Fisher, G., Abbett, W., & Regnier, S. 2004, ApJ, 610, 1148 [NASA ADS] [CrossRef] [Google Scholar]
 Wiegelmann, T., & Sakurai, T. 2012, Liv. Rev. Sol. Phys., 9, 5 [Google Scholar]
 Wiegelmann, T., Thalmann, J. K., & Solanki, S. K. 2014, A&ARv, 22, 78 [Google Scholar]
 Woltjer, L. 1958, Proc. Natl. Acad. Sci. USA, 44, 833 [NASA ADS] [CrossRef] [Google Scholar]
 Yamamoto, T. T., & Sakurai, T. 2009, ApJ, 698, 928 [Google Scholar]
 Yamamoto, T. T., Kusano, K., Maeshiro, T., Yokoyama, T., & Sakurai, T. 2005, ApJ, 624, 1072 [NASA ADS] [CrossRef] [Google Scholar]
 Zhao, L., DeVore, C., Antiochos, S., & Zurbuchen, T. 2015, ApJ, 805, 61 [Google Scholar]
 Zuccarello, F., Aulanier, G., & Gilchrist, S. 2015, ApJ, 814, 126 [NASA ADS] [CrossRef] [Google Scholar]
 Zuccarello, F. P., Pariat, E., Valori, G., & Linan, L. 2018, ApJ, 863, 41 [Google Scholar]
All Figures
Fig. 1. Applied boundarydriving motions for the four different numerical experiments. White represents the positive polarity (B_{z}(z = 0.006) > 0) and black the negative polarity (B_{z}(z = 0.006) < 0). Orange and cyan arrows indicate the distribution of the velocity flows we applied to the negative and positive polarity, respectively. 

In the text 
Fig. 2. Evolution of the different magnetic helicities (top panels), from left to right: relative magnetic helicity (H_{v}, Eq. (4)), volumethreading helicity (H_{pj}, Eq. (7)), and currentcarrying helicity (H_{j}, Eq. (6)). Time evolution of the different magnetic energies (bottom panel), from left to right: total magnetic energy (E_{v}, Eq. (32)), potential energy (E_{p}, Eq. (31)), and free energy (E_{j}, Eq. (31)). The different simulations are dispersion central (red line), dispersion peripheral (green line), stretching (yellow line) and convergence (blue line). The yellow vertical band corresponds to the onset phase of the eruption. 

In the text 
Fig. 3. Left panel: time evolution of the instantaneous timevariation, dE_{v}/dt (dashed black curves, Eq. (37)), and of the sum, dE_{j}/dt + dE_{p}/dt (continuous blue curves, Eqs. (8) and (19)) for the dispersion peripheral run. Right panel: time evolution of the instantaneous timevariation, dH_{v}/dt (dashed black curves, Eq. (23) in Pariat et al. 2015, and of the sum, dH_{j}/dt + dH_{pj}/dt (continuous blue curves, Eqs. (8) and (19)) for the dispersion peripheral run. The yellow bands correspond to the eruption onset phase. 

In the text 
Fig. 4. Time evolution of the helicity variation rates, dH_{j}/dt and dH_{pj}/dt (dashed black curves; Eqs. (8) and (19)), of the helicity transfer term, dH_{j}/dt_{Transf} and dH_{pj}/dt_{Transf} (solid red curves; Eqs. (10) and (21)), of the own terms, dH_{j}/dt_{Own} and dH_{pj}/dt_{Own} (solid blue curves; Eqs. (18) and (29)), and of the dissipation terms, dH_{j}/dt_{Diss} and dH_{pj}/dt_{Diss} (solid green curves; Eqs. (9) and (20)) for the convergence simulation (top panels) and for the dispersion peripheral simulation (bottom panels). Left and right panels: evolution of the current carrying helicity, H_{j}, and volumethreading helicity H_{pj}, respectively. 

In the text 
Fig. 5. Time evolution of the different gauge invariant terms of dH_{j}/dt (top panels), from left to right: dissipation term (dH_{j}/dt_{Diss}, Eq. (9)), own term (dH_{j}/dt_{Own}, Eq. (18)), and helicity transfer term (dH_{j}/dt_{Transf}, Eq. (10)). Time evolution of the different gauge invariant terms of dH_{pj}/dt (bottom panels), from left to right: dissipation term (dH_{pj}/dt_{Diss}, Eq. (20)), and own term (dH_{pj}/dt_{Own}, Eq. (29)). Each curve color corresponds to a particular simulation: dispersion central (red line), dispersion peripheral (green line), stretching (yellow line), and convergence (blue line). The yellow band corresponds to the onset phase of the eruption. 

In the text 
Fig. 6. Left and middle panel: time evolution of the potential energy variation term (dashed black line; dE_{p}/dt; Eq. (42)) and the different terms constituting the instantaneous timevariation of E_{p} (Eq. (42)): F_{ϕ, Bz} (blue line; Eq. (43)), and dE_{p}/dt_{ns} (orange line; Eq. (44)). Left and middle panels: evolution for the convergence and dispersion central simulation, respectively. Right panel: time evolution of F_{ϕ, Bz} for the four simulations dispersion central (red line), dispersion peripheral (green line), stretching (yellow line), and convergence (blue line). The yellow band corresponds to the onset phase of the eruption. 

In the text 
Fig. 7. Time evolution of the freeenergy variation rate (dashed black line; dE_{j}/dt; Eq. (48)) and the different terms constituting the instantaneous timevariation of E_{j} (Eq. (48)): dE_{j}/dt_{Diss} (blue line; Eq. (49)), dE_{j}/dt_{Var} (orange line; Eq. (50)), F_{Bn, Ej} (green line; Eq. (51)), F_{Vn, Ej} (red line; Eq. (52)), and dE_{j}/dt_{ns} (purple line; Eq. (54)). Each panel corresponds to a different simulation: convergence (top left), stretching (top right), dispersion peripheral (bottom left), and dispersion central (bottom right). The yellow band corresponds to the onset phase of the eruption. 

In the text 
Fig. 8. Time evolution of the different gaugeinvariant terms of dE_{j}/dt: F_{Bn, Ej} (top right panel, Eq. (51)), F_{Vn, Ej} (top left panel, Eq. (52)), dE_{j}/dt_{Var} (bottom left panel, Eq. (50)), and dE_{j}/dt_{Diss} (bottom right panel, Eq. (49)). The different colors present one simulation: dispersion central (red line), dispersion peripheral (green line), stretching (yellow line), and convergence (blue line). The yellow band corresponds to the onset phase of the eruption. 

In the text 
Fig. 9. From top to bottom: dimensionless magnitude of (v ⋅ A_{j})B_{z}, (v ⋅ A_{p})B_{z}, (∂ϕ/∂t)(B_{p}), and (v ⋅ B_{j})B_{z} viewed in the (x, y) plane at z = 0.006, at the relative time of t − t_{1} = −58. Isocontours of B_{z} (dashed line for negative values, solid line for positive values) correspond to values of B_{z} = − 4.5, −2.0, 0, 2.0, and 4.5. Each column in the panels presents one simulation, from left to right: convergence, stretching, dispersion peripheral, and dispersion central. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.