Within the coronae of stars, abundances of those elements with low first ionization potential (FIP) often differ from their photospheric values. The coronae of the Sun and solar-type stars mostly show enhancements of low- FIP elements (the FIP effect) while more active stars such as M-dwarfs have coronae generally characterized by the inverse-FIP effect (I-FIP). Here we observe patches of I-FIP effect solar plasma in AR 12673, a highly complex beta/gamma/delta active region. We argue that the umbrae of coalescing sunspots and more specifically strong light bridges within the umbrae, are preferential locations for observing I-FIP effect plasma. Furthermore, the magnetic complexity of the active region and major episodes of fast flux emergence also lead to repetitive and intense flares. The induced evaporation of the chromospheric plasma in flare ribbons crossing umbrae enables the observation of four localized patches of I-FIP effect plasma in the corona of AR 12673. These observations can be interpreted in the context of the ponderomotive force fractionation model which predicts that plasma with I-FIP effect composition is created by the refraction of waves coming from below the chromosphere. We propose that the waves generating the I-FIP effect plasma in solar active regions are generated by sub-photospheric reconnection of coalescing flux systems. Although we only glimpse signatures of I-FIP effect fractionation produced by this interaction in patches on the Sun, on highly active M-stars it may be the dominant process.

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.

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

Freely-propagating global waves in the solar atmosphere are commonly observed using Extreme UltraViolet passbands (EUV or "EIT waves"), and less regularly in H-alpha (Moreton-Ramsey waves). Despite decades of research, joint observations of EUV and Moreton-Ramsey waves remain rare, complicating efforts to quantify the connection between these phenomena. We present observations of four homologous global waves originating from the same active region between 28-30 March 2014 and observed using both EUV and H-alpha data. Each global EUV wave was observed by the Solar Dynamics Observatory, with the associated Moreton-Ramsey waves identified using the Global Oscillations Network Group (GONG) network. All of the global waves exhibit high initial velocity (e.g., 842-1388 km s$^{-1}$ in the 193A passband) and strong deceleration (e.g., -1437 - -782 m s$^{-2}$ in the 193A passband) in each of the EUV passbands studied, with the EUV wave kinematics exceeding those of the Moreton-Ramsey wave. The density compression ratio of each global wave was estimated using both differential emission measure and intensity variation techniques, with both indicating that the observed waves were weakly shocked with a fast magnetosonic Mach number slightly greater than one. This suggests that, according to current models, the global coronal waves were not strong enough to produce Moreton-Ramsey waves, indicating an alternative explanation for these observations. Instead, we conclude that the evolution of the global waves was restricted by the surrounding coronal magnetic field, in each case producing a downward-angled wavefront propagating towards the north solar pole which perturbed the chromosphere and was observed as a Moreton-Ramsey wave.

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

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

Understanding elemental abundance variations in the solar corona provides an insight into how matter and energy flow from the chromosphere into the heliosphere. Observed variations depend on the first ionization potential (FIP) of the main elements of the Sun's atmosphere. High-FIP elements (>10 eV) maintain photospheric abundances in the corona, whereas low-FIP elements have enhanced abundances. Conversely, inverse FIP (IFIP) refers to the enhancement of high-FIP or depletion of low-FIP elements. We use spatially resolved spectroscopic observations, specifically the Ar XIV/Ca XIV intensity ratio, from Hinode's Extreme-ultraviolet Imaging Spectrometer to investigate the distribution and evolution of plasma composition within two confined flares in a newly emerging, highly sheared active region. During the decay phase of the first flare, patches above the flare ribbons evolve from the FIP to the IFIP effect, while the flaring loop tops show a stronger FIP effect. The patch and loop compositions then evolve toward the pre-flare basal state. We propose an explanation of how flaring in strands of highly sheared emerging magnetic fields can lead to flare-modulated IFIP plasma composition over coalescing umbrae which are crossed by flare ribbons. Subsurface reconnection between the coalescing umbrae leads to the depletion of low-FIP elements as a result of an increased wave flux from below. This material is evaporated when the flare ribbons cross the umbrae. Our results are consistent with the ponderomotive fractionation model (Laming2015) for the creation of IFIP-biased plasma.

Quiescent solar prominences are observed to exist within the solar atmosphere for up to several solar rotations. Their eruption is commonly preceded by a slow increase in height that can last from hours to days. This increase in the prominence height is believed to be due to their host magnetic flux rope transitioning through a series of neighbouring quasi-equilibria before the main loss-of-equilibrium that drives the eruption. Recent work suggests that the removal of prominence mass from a stable, quiescent flux rope is one possible cause for this change in height. However, these conclusions are drawn from observations and are subject to interpretation. Here we present a simple model to quantify the effect of "mass-draining" during the pre-eruptive height-evolution of a solar flux rope. The flux rope is modeled as a line current suspended within a background potential magnetic field. We first show that the inclusion of mass, up to $10^{12}$~kg, can modify the height at which the line current experiences loss-of-equilibrium by up to 14\%. Next, we show that the rapid removal of mass prior to the loss-of-equilibrium can allow the height of the flux rope to increase sharply and without upper bound as it approaches its loss-of-equilibrium point. This indicates that the critical height for the loss-of-equilibrium can occur at a range of heights depending explicitly on the amount and evolution of mass within the flux rope. Finally, we demonstrate that for the same amount of drained mass, the effect on the height of the flux rope is up to two order of magnitude larger for quiescent than for active region prominences.

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.

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.

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.

Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the coronae of stars. Understanding the plasma processes involved in CME initiation has applications to space weather forecasting and laboratory plasma experiments. James et al. (Sol. Phys. 292, 71, 2017) used EUV observations to conclude that a magnetic flux rope formed in the solar corona above NOAA Active Region 11504 before it erupted on 14 June 2012 (SOL2012-06-14). In this work, we use data from the Solar Dynamics Observatory to model the coronal magnetic field of the active region one hour prior to eruption using a nonlinear force-free field extrapolation, and find a flux rope reaching a maximum height of 150 Mm above the photosphere. Estimations of the average twist of the strongly asymmetric extrapolated flux rope are between 1.35 and 1.88 turns, depending on the choice of axis, although the erupting structure was not observed to kink. The decay index near the apex of the axis of the extrapolated flux rope is comparable to typical critical values required for the onset of the torus instability, so we suggest that the torus instability drove the eruption.

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.

Determining the energy transport mechanisms in flares remains a central goal in solar flares physics that is still not adequately answered by the 'standard flare model'. In particular, the relative roles of particles and/or waves as transport mechanisms, the contributions of low energy protons and ions to the overall flare budget, and the limits of low energy non-thermal electron distribution are questions that still cannot be adequately reconciled with current instrumentation. In this 'White Paper' submitted in response to the call for inputs to the Next Generation Solar Physics Mission review process initiated by JAXA, NASA and ESA in 2016, we outline the open questions in this area and possible instrumentation that could provide the required observations to help answer these and other flare-related questions.

Sunquakes are created by the hydrodynamic response of the lower atmosphere to a sudden deposition of energy and momentum. In this study we investigate a sunquake that occurred in NOAA active region 11675 on 17 February 2013. Observations of the corona, chromosphere and photosphere are brought together for the first time with a non-linear force-free model of the active region's magnetic field in order to probe the magnetic environment in which the sunquake was initiated. We find that the sunquake was associated with the destabilization of a flux rope and an associated M-class GOES flare. Active region 11675 was in its emergence phase at the time of the sunquake and photospheric motions caused by the emergence heavily modified the flux rope and its associated quasi-separatrix layers, eventually triggering the flux rope's instability. The flux rope was surrounded by an extended envelope of field lines rooted in a small area at the approximate position of the sunquake. We argue that the configuration of the envelope, by interacting with the expanding flux rope, created a "magnetic lens" that may have focussed energy in one particular location the photosphere, creating the necessary conditions for the initiation of the sunquake.

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

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

Coronal mass ejections (CMEs) are one of the primary manifestations of solar activity and can drive severe space weather effects. Therefore, it is vital to work towards being able to predict their occurrence. However, many aspects of CME formation and eruption remain unclear, including whether magnetic flux ropes are present before the onset of eruption and the key mechanisms that cause CMEs to occur. In this work, the pre-eruptive coronal configuration of an active region that produced an interplanetary CME with a clear magnetic flux rope structure at 1 AU is studied. A forward-S sigmoid appears in extreme-ultraviolet (EUV) data two hours before the onset of the eruption (SOL2012-06-14), which is interpreted as a signature of a right-handed flux rope that formed prior to the eruption. Flare ribbons and EUV dimmings are used to infer the locations of the flux rope footpoints. These locations, together with observations of the global magnetic flux distribution, indicate that an interaction between newly emerged magnetic flux and pre-existing sunspot field in the days prior to the eruption may have enabled the coronal flux rope to form via tether-cutting-like reconnection. Composition analysis suggests that the flux rope had a coronal plasma composition, supporting our interpretation that the flux rope formed via magnetic reconnection in the corona. Once formed, the flux rope remained stable for 2 hours before erupting as a CME.

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.

"EIT waves" are freely-propagating global pulses in the low corona which are strongly associated with the initial evolution of coronal mass ejections (CMEs). They are thought to be large-amplitude, fast-mode magnetohydrodynamic waves initially driven by the rapid expansion of a CME in the low corona. An "EIT wave" was observed on 6 July 2012 to impact an adjacent trans-equatorial loop system which then exhibited a decaying oscillation as it returned to rest. Observations of the loop oscillations were used to estimate the magnetic field strength of the loop system by studying the decaying oscillation of the loop, measuring the propagation of ubiquitous transverse waves in the loop and extrapolating the magnetic field from observed magnetograms. Observations from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory (SDO/AIA) and the Coronal Multi-channel Polarimeter (CoMP) were used to study the event. An Empirical Mode Decomposition analysis was used to characterise the oscillation of the loop system in CoMP Doppler velocity and line width and in AIA intensity. The loop system was shown to oscillate in the 2nd harmonic mode rather than at the fundamental frequency, with the seismological analysis returning an estimated magnetic field strength of ~5.5+/-1.5 G. This compares to the magnetic field strength estimates of ~1-9 G and ~3-9 G found using the measurements of transverse wave propagation and magnetic field extrapolation respectively.

A key aim in space weather research is to be able to use remote-sensing observations of the solar atmosphere to extend the lead time of predicting the geoeffectiveness of a coronal mass ejection (CME). In order to achieve this, the magnetic structure of the CME as it leaves the Sun must be known. In this article we address this issue by developing a method to determine the intrinsic flux rope type of a CME solely from solar disk observations. We use several well known proxies for the magnetic helicity sign, the axis orientation, and the axial magnetic field direction to predict the magnetic structure of the interplanetary flux rope. We present two case studies: the 2 June 2011 and the 14 June 2012 CMEs. Both of these events erupted from an active region and, despite having clear in situ counterparts, their eruption characteristics were relatively complex. The first event was associated with an active region filament that erupted in two stages, while for the other event the eruption originated from a relatively high coronal altitude and the source region did not feature the presence of a filament. Our magnetic helicity sign proxies include the analysis of magnetic tongues, soft X-ray and/or EUV sigmoids, coronal arcade skew, filament emission and absorption threads, and filament rotation. Since the inclination of the post-eruption arcades was not clear, we use the tilt of the polarity inversion line to determine the flux rope axis orientation, and coronal dimmings to determine the flux rope footpoints and, therefore, the direction of the axial magnetic field. The comparison of the estimated intrinsic flux rope structure to in situ observations at the Lagrangian point L1 indicated a good agreement with the predictions. Our results highlight the flux rope type determination techniques that are particularly useful for active region eruptions, where most geoeffective CMEs originate.

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

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

We investigate whether flux cancellation is responsible for the formation of a very massive filament resulting in the spectacular 2011 June 7 eruption. We analyse and quantify the amount of flux cancellation that occurs in NOAA AR 11226 and its two neighbouring ARs (11227 & 11233) using line-of-sight magnetograms from the Heliospheric Magnetic Imager. During a 3.6-day period building up to the filament eruption, 1.7 x 10^21 Mx, 21% of AR 11226's maximum magnetic flux, was cancelled along the polarity inversion line (PIL) where the filament formed. If the flux cancellation continued at the same rate up until the eruption then up to 2.8 x 10^21 Mx (34% of the AR flux) may have been built into the magnetic configuration that contains the filament plasma. The large flux cancellation rate is due to an unusual motion of the positive polarity sunspot, which splits, with the largest section moving rapidly towards the PIL. This motion compresses the negative polarity and leads to the formation of an orphan penumbra where one end of the filament is rooted. Dense plasma threads above the orphan penumbra build into the filament, extending its length, and presumably injecting material into it. We conclude that the exceptionally strong flux cancellation in AR 11226 played a significant role in the formation of its unusually massive filament. In addition, the presence and coherent evolution of bald patches in the vector magnetic field along the PIL suggests that the magnetic field configuration supporting the filament material is that of a flux rope.

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

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.

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

We study a filament eruption, two-ribbon flare, and coronal mass ejection (CME) that occurred in Active Region NOAA 10898 on 6 July 2006. The filament was located South of a strong sunspot that dominated the region. In the evolution leading up to the eruption, and for some time after it, a counter-clockwise rotation of the sunspot of about 30 degrees was observed. We suggest that the rotation triggered the eruption by progressively expanding the magnetic field above the filament. To test this scenario, we study the effect of twisting the initially potential field overlying a pre-existing flux-rope, using three-dimensional zero-$\beta$ MHD simulations. We first consider a relatively simple and symmetric system, and then study a more complex and asymmetric magnetic configuration, whose photospheric flux distribution and coronal structure are guided by the observations and a potential field extrapolation. In both cases, we find that the twisting leads to the expansion of the overlying field. As a consequence of the progressively reduced magnetic tension, the flux-rope quasi-statically adapts to the changed environmental field, rising slowly. Once the tension is sufficiently reduced, a distinct second phase of evolution occurs where the flux-rope enters an unstable regime characterized by a strong acceleration. Our simulations thus suggest a new mechanism for the triggering of eruptions in the vicinity of rotating sunspots.

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.

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.

Extrapolations of solar photospheric vector magnetograms into three-dimensional magnetic fields in the chromosphere and corona are usually done under the assumption that the fields are force-free. The field calculations can be improved by preprocessing the photospheric magnetograms. We compare two preprocessing methods presently in use, namely the methods of Wiegelmann et al. (2006) and Fuhrmann et al. (2007). The two preprocessing methods were applied to a recently observed vector magnetogram. We examine the changes in the magnetogram effected by the two preprocessing algorithms. Furthermore, the original magnetogram and the two preprocessed magnetograms were each used as input data for nonlinear force-free field extrapolations by means of two different methods, and we analyze the resulting fields. Both preprocessing methods managed to significantly decrease the magnetic forces and magnetic torques that act through the magnetogram area and that can cause incompatibilities with the assumption of force-freeness in the solution domain. Both methods also reduced the amount of small-scale irregularities in the observed photospheric field, which can sharply worsen the quality of the solutions. For the chosen parameter set, the Wiegelmann et al. method led to greater changes in strong-field areas, leaving weak-field areas mostly unchanged, and thus providing an approximation of the magnetic field vector in the chromosphere, while the Fuhrmann et al. method weakly changed the whole magnetogram, thereby better preserving patterns present in the original magnetogram. Both preprocessing methods raised the magnetic energy content of the extrapolated fields to values above the minimum energy, corresponding to the potential field. Also, the fields calculated from the preprocessed magnetograms fulfill the solenoidal condition better than those calculated without preprocessing.

We examine the nonlinear magnetofrictional extrapolation scheme using the solar active region model by Titov and D\'emoulin as test field. This model consists of an arched, line-tied current channel held in force-free equilibrium by the potential field of a bipolar flux distribution in the bottom boundary. A modified version, having a parabolic current density profile, is employed here. We find that the equilibrium is reconstructed with very high accuracy in a representative range of parameter space, using only the vector field in the bottom boundary as input. Structural features formed in the interface between the flux rope and the surrounding arcade-"hyperbolic flux tube" and "bald patch separatrix surface"-are reliably reproduced, as are the flux rope twist and the energy and helicity of the configuration. This demonstrates that force-free fields containing these basic structural elements of solar active regions can be obtained by extrapolation. The influence of the chosen initial condition on the accuracy of reconstruction is also addressed, confirming that the initial field that best matches the external potential field of the model quite naturally leads to the best reconstruction. Extrapolating the magnetogram of a Titov-D\'emoulin equilibrium in the unstable range of parameter space yields a sequence of two opposing evolutionary phases which clearly indicate the unstable nature of the configuration: a partial buildup of the flux rope with rising free energy is followed by destruction of the rope, losing most of the free energy.

Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and powered by the partial dissipation of electrical currents in the solar atmosphere. A critical unanswered question is whether the currents involved are induced by the motion of pre-existing atmospheric magnetic flux subject to surface plasma flows, or whether these currents are associated with the emergence of flux from within the solar convective zone. We address this problem by applying state-of-the-art nonlinear force-free field (NLFFF) modeling to the highest resolution and quality vector-magnetographic data observed by the recently launched Hinode satellite on NOAA Active Region 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with 4 different codes and a variety of boundary conditions. We find that the model fields differ markedly in geometry, energy content, and force-freeness. We discuss the relative merits of these models in a general critique of present abilities to model the coronal magnetic field based on surface vector field measurements. For our application in particular, we find a fair agreement of the best-fit model field with the observed coronal configuration, and argue (1) that strong electrical currents emerge together with magnetic flux preceding the flare, (2) that these currents are carried in an ensemble of thin strands, (3) that the global pattern of these currents and of field lines are compatible with a large-scale twisted flux rope topology, and (4) that the ~10^32 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.

CONTEXT: As the coronal magnetic field can usually not be measured directly, it has to be extrapolated from photospheric measurements into the corona. AIMS: We test the quality of a non-linear force-free coronal magnetic field extrapolation code with the help of a known analytical solution. METHODS: The non-linear force-free equations are numerically solved with the help of an optimization principle. The method minimizes an integral over the force-free and solenoidal condition. As boundary condition we use either the magnetic field components on all six sides of the computational box in Case I or only on the bottom boundary in Case II. We check the quality of the reconstruction by computing how well force-freeness and divergence-freeness are fulfilled and by comparing the numerical solution with the analytical solution. The comparison is done with magnetic field line plots and several quantitative measures, like the vector correlation, Cauchy Schwarz, normalized vector error, mean vector error and magnetic energy. RESULTS: For Case I the reconstructed magnetic field shows good agreement with the original magnetic field topology, whereas in Case II there are considerable deviations from the exact solution. This is corroborated by the quantitative measures, which are significantly better for Case I. CONCLUSIONS: Despite the strong nonlinearity of the considered force-free equilibrium, the optimization method of extrapolation is able to reconstruct it; however, the quality of reconstruction depends significantly on the consistency of the input data, which is given only if the known solution is provided also at the lateral and top boundaries, and on the presence or absence of flux concentrations near the boundaries of the magnetogram.