Accuracy of magnetic energy computations

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

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

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.Comment: 4 pages, 3 figures; published online in A&A 555, L6 (2013

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

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.Comment: Astronomy & Astrophysics, in pres

Magneto-frictional Modeling Of Coronal Nonlinear Force-free Fields. I. Testing With Analytic Solutions

We report our implementation of the magneto-frictional method in the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The method aims at applications where local adaptive mesh refinement (AMR) is essential to make follow-up dynamical modeling affordable. We quantify its performance in both domain-decomposed uniform grids and block-adaptive AMR computations, using all frequently employed force-free, divergence-free, and other vector comparison metrics. As test cases, we revisit the semi-analytic solution of Low and Lou in both Cartesian and spherical geometries, along with the topologically challenging Titov-Démoulin model. We compare different combinations of spatial and temporal discretizations, and find that the fourth-order central difference with a local Lax-Friedrichs dissipation term in a single-step marching scheme is an optimal combination. The initial condition is provided by the potential field, which is the potential field source surface model in spherical geometry. Various boundary conditions are adopted, ranging from fully prescribed cases where all boundaries are assigned with the semi-analytic models, to solar-like cases where only the magnetic field at the bottom is known. Our results demonstrate that all the metrics compare favorably to previous works in both Cartesian and spherical coordinates. Cases with several AMR levels perform in accordance with their effective resolutions. The magneto-frictional method in MPI-AMRVAC allows us to model a region of interest with high spatial resolution and large field of view simultaneously, as required by observation-constrained extrapolations using vector data provided with modern instruments. The applications of the magneto-frictional method to observations are shown in an accompanying paper

The 2013 February 17 sunquake in the context of the active region's magnetic field configuration

© 2017. The American Astronomical Society. All rights reserved. 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 2013 February 17. Observations of the corona, chromosphere, and photosphere are brought together for the first time with a nonlinear 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 on one particular location of the photosphere, creating the necessary conditions for the initiation of the sunquake

Flux cancellation and the evolution of the eruptive filament of 2011 June 7

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.Comment: 18 pages, 7 figures. Submitted to ApJ in December 2015, accepted in June 201

Localized Acceleration of Energetic Particles by a Weak Shock in the Solar Corona

Globally-propagating shocks in the solar corona have long been studied to quantify their involvement in the acceleration of energetic particles. However, this work has tended to focus on large events associated with strong solar flares and fast coronal mass ejections (CMEs), where the waves are sufficiently fast to easily accelerate particles to high energies. Here we present observations of particle acceleration associated with a global wave event which occurred on 1 October 2011. Using differential emission measure analysis, the global shock wave was found to be incredibly weak, with an Alfv\'en Mach number of ~1.008-1.013. Despite this, spatially-resolved type III radio emission was observed by the Nan\c{c}ay RadioHeliograph at distinct locations near the shock front, suggesting localised acceleration of energetic electrons. Further investigation using a magnetic field extrapolation identified a fan structure beneath a magnetic null located above the source active region, with the erupting CME contained within this topological feature. We propose that a reconfiguration of the coronal magnetic field driven by the erupting CME enabled the weak shock to accelerate particles along field lines initially contained within the fan and subsequently opened into the heliosphere, producing the observed type III emission. These results suggest that even weak global shocks in the solar corona can accelerate energetic particles via reconfiguration of the surrounding magnetic field

Testing non-linear force-free coronal magnetic field extrapolations with the Titov-Demoulin equilibrium

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.Comment: 6 pages, 2 figures, Research Not

On-disc observations of flux rope formation prior to its eruption

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 two hours before erupting as a CME