Relative Magnetic Helicity Research Articles

Relative field-line helicity in bounded domains

Models for astrophysical plasmas often have magnetic field lines that leave the boundary rather than closing within the computational domain. Thus, the relative magnetic helicity is frequently used in place of the usual magnetic helicity, so as to restore gauge invariance. We show how to decompose the relative helicity into a relative field-line helicity that is an ideal-magnetohydrodynamic invariant for each individual magnetic field line, and vanishes along any field line where the original field matches the reference field. Physically, this relative field-line helicity is a magnetic flux, whose specific definition depends on the gauge of the reference vector potential on the boundary. We propose a particular ‘minimal’ gauge that depends only on the reference field and minimises this boundary contribution, so as to reveal topological information about the original magnetic field. We illustrate the effect of different gauge choices using the Low–Lou and Titov–Démoulin models of solar active regions. Our numerical code to compute appropriate vector potentials and relative field-line helicity in Cartesian domains is open source and freely available.

Time Variations of the Nonpotential and Volume-threading Magnetic Helicities

Abstract Relative magnetic helicity is a gauge-invariant quantity suitable for the study of the magnetic helicity content of heliospheric plasmas. Relative magnetic helicity can be decomposed uniquely into two gauge-invariant quantities, the magnetic helicity of the nonpotential 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 nonpotential 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 nonpotential and volume-threading helicities. The validity of the analytical formulae derived are confirmed with analysis of 3D magnetohydrodynamics (MHD) simulations of solar coronal dynamics. Both the analytical investigation and the numerical application show that, unlike magnetic helicity, the nonpotential and the volume-threading helicities are not conserved quantities, even in the ideal MHD regime. A term corresponding to the transformation between the nonpotential 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 nonpotential and volume-threading helicities and the observational analysis of helicity-based solar-eruptivity proxies.

Dependence of Coronal Mass Ejection Properties on Their Solar Source Active Region Characteristics and Associated Flare Reconnection Flux

Abstract The near-Sun kinematics of coronal mass ejections (CMEs) determine the severity and arrival time of associated geomagnetic storms. We investigate the relationship between the deprojected speed and kinetic energy of CMEs and magnetic measures of their solar sources, reconnection flux of associated eruptive events, and intrinsic flux-rope characteristics. Our data covers the period 2010–2014 in solar cycle 24. Using vector magnetograms of source active regions, we estimate the size and nonpotentiality. We compute the total magnetic reconnection flux at the source regions of CMEs using the post-eruption arcade method. By forward modeling the CMEs, we find their deprojected geometric parameters and constrain their kinematics and magnetic properties. Based on an analysis of this database, we report that the correlation between CME speed and their source active region size and global nonpotentiality is weak, but not negligible. We find the near-Sun velocity and kinetic energy of CMEs to be well correlated with the associated magnetic reconnection flux. We establish a statistically significant empirical relationship between the CME speed and reconnection flux that may be utilized for prediction purposes. Furthermore, we find CME kinematics to be related with the axial magnetic field intensity and relative magnetic helicity of their intrinsic flux ropes. The amount of coronal magnetic helicity shed by CMEs is found to be well correlated with their near-Sun speeds. The kinetic energy of CMEs is well correlated with their intrinsic magnetic energy density. Our results constrain processes related to the origin and propagation of CMEs and may lead to better empirical forecasting of their arrival and geoeffectiveness.

Threshold of Non-potential Magnetic Helicity Ratios at the Onset of Solar Eruptions

Abstract The relative magnetic helicity is a quantity that is often used to describe the level of entanglement of non-isolated magnetic fields, such as the magnetic field of solar active regions. The aim of this paper is to investigate how different kinds of photospheric boundary flows accumulate relative magnetic helicity in the corona and if and how well magnetic-helicity-related quantities identify the 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, and 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 a threshold is not observed for other quantities as, for example, those related to the magnetic energy.

Computation of Relative Magnetic Helicity in Spherical Coordinates

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.

Evolution of relative magnetic helicity

Context. For a better understanding of the dynamics of the solar corona, it is important to analyse the evolution of the helicity of the magnetic field. Since the helicity cannot be directly determined by observations, we have recently proposed a method to calculate the relative magnetic helicity in a finite volume for a given magnetic field, which however required the flux to be balanced separately on all the sides of the considered volume. Aims. We developed a scheme to obtain the vector potential in a volume without the above restriction at the boundary. We studied the dissipation and escape of relative magnetic helicity from an active region. Methods. In order to allow finite magnetic fluxes through the boundaries, a Coulomb gauge was constructed that allows for global magnetic flux balance. The property of sinusoidal function was used to obtain the vector potentials at the 12 edges of the considered rectangular volume extending above an active region. We tested and verified our method in a theoretical fore-free magnetic field model. Results. We applied the new method to the former calculation data and found a difference of less than 1.2%. We also applied our method to the magnetic field above active region NOAA 11429 obtained by a new photospheric-data-driven magnetohydrodynamics (MHD) model code GOEMHD3. We analysed the magnetic helicity evolution in the solar corona using our new method. We find that the normalized magnetic helicity (H∕Φ2) is equal to −0.038 when fast magnetic reconnection is triggered. This value is comparable to the previous value (−0.029) in the MHD simulations when magnetic reconnection happened and the observed normalized magnetic helicity (−0.036) from the eruption of newly emerging active regions. We find that only 8% of the accumulated magnetic helicity is dissipated after it is injected through the bottom boundary. This is in accordance with the Woltjer conjecture. Only 2% of the magnetic helicity injected from the bottom boundary escapes through the corona. This is consistent with the observation of magnetic clouds, which could take magnetic helicity into the interplanetary space. In the case considered here, several halo coronal mass ejections (CMEs) and two X-class solar flares originate from this active region.

Optimization of Photospheric Electric Field Estimates for Accurate Retrieval of Total Magnetic Energy Injection

Estimates of the photospheric magnetic, electric and plasma velocity fields are essential for studying the dynamics of the solar atmosphere, for example through the derivative quantities of Poynting and relative helicity flux and by using of the fields to obtain the lower boundary condition for data-driven coronal simulations. In this paper we study the performance of a data processing and electric field inversion approach that requires only high-resolution and high-cadence line-of-sight or vector magnetograms -- which we obtain from Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO). The approach does not require any photospheric velocity estimates, and the lacking velocity information is compensated using ad hoc assumptions. We show that the free parameters of these assumptions can be optimized to reproduce the time evolution of the total magnetic energy injection through the photosphere in NOAA AR 11158, when compared to the recent estimates for this active region. However, we find that the relative magnetic helicity injection is reproduced poorly reaching at best a modest underestimation. We discuss also the effect of some of the data processing details on the results, including the masking of the noise-dominated pixels and the tracking method of the active region, both of which have not received much attention in the literature so far. In most cases the effect of these details is small, but when the optimization of the free parameters of the ad hoc assumptions is considered a consistent use of the noise mask is required. The results found in this paper imply that the data processing and electric field inversion approach that uses only the photospheric magnetic field information offers a flexible and straightforward way to obtain photospheric magnetic and electric field estimates suitable for practical applications such as coronal modeling studies.

A Helicity-Based Method to Infer the CME Magnetic Field Magnitude in Sun and Geospace: Generalization and Extension to Sun-Like and M-Dwarf Stars and Implications for Exoplanet Habitability

Patsourakos et al. (Astrophys. J. 817, 14, 2016) and Patsourakos and Georgoulis (Astron. Astrophys. 595, A121, 2016) introduced a method to infer the axial magnetic field in flux-rope coronal mass ejections (CMEs) in the solar corona and farther away in the interplanetary medium. The method, based on the conservation principle of magnetic helicity, uses the relative magnetic helicity of the solar source region as input estimates, along with the radius and length of the corresponding CME flux rope. The method was initially applied to cylindrical force-free flux ropes, with encouraging results. We hereby extend our framework along two distinct lines. First, we generalize our formalism to several possible flux-rope configurations (linear and nonlinear force-free, non-force-free, spheromak, and torus) to investigate the dependence of the resulting CME axial magnetic field on input parameters and the employed flux-rope configuration. Second, we generalize our framework to both Sun-like and active M-dwarf stars hosting superflares. In a qualitative sense, we find that Earth may not experience severe atmosphere-eroding magnetospheric compression even for eruptive solar superflares with energies ~ 10^4 times higher than those of the largest Geostationary Operational Environmental Satellite (GOES) X-class flares currently observed. In addition, the two recently discovered exoplanets with the highest Earth-similarity index, Kepler 438b and Proxima b, seem to lie in the prohibitive zone of atmospheric erosion due to interplanetary CMEs (ICMEs), except when they possess planetary magnetic fields that are much higher than that of Earth.

Relative magnetic helicity as a diagnostic of solar eruptivity

Context. 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. Aims. 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. Methods. From the magnetic field generated by the three-dimensional magnetohydrodynamical simulations, we compute and analyze 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. Results. 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. Conclusions. 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.

Multi-region relaxed magnetohydrodynamics in plasmas with slowly changing boundaries—Resonant response of a plasma slab

The adiabatic limit of a recently proposed dynamical extension of Taylor relaxation, multi-region relaxed magnetohydrodynamics (MRxMHD), is summarized, with special attention to the appropriate definition of a relative magnetic helicity. The formalism is illustrated using a simple two-region, sheared-magnetic-field model similar to the Hahm–Kulsrud–Taylor (HKT) rippled-boundary slab model. In MRxMHD, a linear Grad–Shafranov equation applies, even at finite ripple amplitude. The adiabatic switching on of boundary ripple excites a shielding current sheet opposing reconnection at a resonant surface. The perturbed magnetic field as a function of ripple amplitude is calculated by invoking the conservation of magnetic helicity in the two regions separated by the current sheet. At low ripple amplitude, “half islands” appear on each side of the current sheet, locking the rotational transform at the resonant value. Beyond a critical amplitude, these islands disappear and the rotational transform develops a discontinuity across the current sheet.

Sunspot rotation

Context. Solar eruptions and high flare activity often accompany the rapid rotation of sunspots. The study of sunspot rotation and the mechanisms driving this motion are therefore key to our understanding of how the solar atmosphere attains the conditions necessary for large energy release. Aims. We aim to demonstrate and investigate the rotation of sunspots in a 3D numerical experiment of the emergence of a magnetic flux tube as it rises through the solar interior and emerges into the atmosphere. Furthermore, we seek to show that the sub-photospheric twist stored in the interior is injected into the solar atmosphere by means of a definitive rotation of the sunspots. Methods. A numerical experiment is performed to solve the 3D resistive magnetohydrodynamic (MHD) equations using a Lagrangian-Remap code. We track the emergence of a toroidal flux tube as it rises through the solar interior and emerges into the atmosphere investigating various quantities related to both the magnetic field and plasma. Results. Through detailed analysis of the numerical experiment, we find clear evidence that the photospheric footprints or sunspots of the flux tube undergo a rotation. Significant vertical vortical motions are found to develop within the two polarity sources after the field emerges. These rotational motions are found to leave the interior portion of the field untwisted and twist up the atmospheric portion of the field. This is shown by our analysis of the relative magnetic helicity as a significant portion of the interior helicity is transported to the atmosphere. In addition, there is a substantial transport of magnetic energy to the atmosphere. Rotation angles are also calculated by tracing selected fieldlines; the fieldlines threading through the sunspot are found to rotate through angles of up to 353 degrees over the course of the experiment.

THE INFLUENCE OF SPATIAL RESOLUTION ON NONLINEAR FORCE-FREE MODELING

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.

Energy and helicity injection in solar quiet regions

Aims. We investigate the free magnetic energy and relative magnetic helicity injection in solar quiet regions.Methods. We use the DAVE4VM method to infer the photospheric velocity field and calculate the free magnetic energy and relative magnetic helicity injection rates in 16 quiet-Sun vector magnetograms sequences.Results. We find that there is no dominant sense of helicity injection in quiet-Sun regions, and that both helicity and energy injections are mostly due to surface shuffling motions that dominate the respective emergence by factors slightly larger than two. We, furthermore, estimate the helicity and energy rates per network unit area as well as the respective budgets over a complete solar cycle.Conclusions. Derived helicity and energy budgets over the entire solar cycle are similar to respective budgets derived in a recent work from the instantaneous helicity and free magnetic energy budgets and higher than previously reported values that relied on similar approaches to this analysis. Free-energy budgets, mostly generated like helicity at the network, are high enough to power the dynamics of fine-scale structures residing at the network, such as mottles and spicules, while corresponding estimates of helicity budgets are provided, pending future verification from high-resolution magneto-hydrodynamic simulations and/or observations.

Comparison between the eruptive X2.2 flare on 2011 February 15 and confined X3.1 flare on 2014 October 24

We compare two contrasting X-class flares in terms of magnetic free energy, relative magnetic helicity and decay index of the active regions (ARs) in which they occurred. The events in question are the eruptive X2.2 flare from AR 11158 accompanied by a halo coronal mass ejection (CME) and the confined X3.1 flare from AR 12192 with no associated CME. These two flares exhibit similar behavior of free magnetic energy and helicity buildup for a few days preceding them. A major difference between the two flares is found to lie in the time-dependent change of magnetic helicity of the ARs that hosted them. AR 11158 shows a significant decrease in magnetic helicity starting ∼4 hours prior to the flare, but no apparent decrease in helicity is observed in AR 12192. By examining the magnetic helicity injection rates in terms of sign, we confirmed that the drastic decrease in magnetic helicity before the eruptive X2.2 flare was not caused by the injection of reversed helicity through the photosphere but rather the CME-related change in the coronal magnetic field. Another major difference we find is that AR 11158 had a significantly larger decay index and therefore weaker overlying field than AR 12192. These results suggest that the coronal magnetic helicity and the decay index of the overlying field can provide a clue about the occurrence of CMEs.

Testing magnetic helicity conservation in a solar-like active event

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.

Evolution of field line helicity during magnetic reconnection

We investigate the evolution of field line helicity for magnetic fields that connect two boundaries without null points, with emphasis on localized finite-B magnetic reconnection. Total (relative) magnetic helicity is already recognized as an important topological constraint on magnetohydrodynamic processes. Field line helicity offers further advantages because it preserves all topological information and can distinguish between different magnetic fields with the same total helicity. Magnetic reconnection changes field connectivity and field line helicity reflects these changes; the goal of this paper is to characterize that evolution. We start by deriving the evolution equation for field line helicity and examining its terms, also obtaining a simplified form for cases where dynamics are localized within the domain. The main result, which we support using kinematic examples, is that during localized reconnection in a complex magnetic field, the evolution of field line helicity is dominated by a work-like term that is evaluated at the field line endpoints, namely, the scalar product of the generalized field line velocity and the vector potential. Furthermore, the flux integral of this term over certain areas is very small compared to the integral of the unsigned quantity, which indicates that changes of field line helicity happen in a well-organized pairwise manner. It follows that reconnection is very efficient at redistributing helicity in complex magnetic fields despite having little effect on the total helicity.

JOINT INVERSE CASCADE OF MAGNETIC ENERGY AND MAGNETIC HELICITY IN MHD TURBULENCE

We show that oppositely directed fluxes of energy and magnetic helicity coexist in the inertial range in fully developed magnetohydrodynamic (MHD) turbulence with small-scale sources of magnetic helicity. Using a helical shell model of MHD turbulence, we study the high Reynolds number magnetohydrodynamic turbulence for helicity injection at a scale that is much smaller than the scale of energy injection. In a short range of scales larger than the forcing scale of magnetic helicity, a bottleneck-like effect appears, which results in a local reduction of the spectral slope. The slope changes in a domain with a high level of relative magnetic helicity, which determines that part of the magnetic energy related to the helical modes at a given scale. If the relative helicity approaches unity, the spectral slope tends to $-3/2$. We show that this energy pileup is caused by an inverse cascade of magnetic energy associated with the magnetic helicity. This negative energy flux is the contribution of the pure magnetic-to-magnetic energy transfer, which vanishes in the non-helical limit. In the context of astrophysical dynamos, our results indicate that a large-scale dynamo can be affected by the magnetic helicity generated at small scales. The kinetic helicity, in particular, is not involved in the process at all. An interesting finding is that an inverse cascade of magnetic energy can be provided by a small-scale source of magnetic helicity fluctuations without a mean injection of magnetic helicity.

Lower and upper bounds on the energy of braided magnetic fields

Let ℬ be a class of topologically equivalent regular magnetic fields occupying a cylindrical domain of axis Oz and produced by ideal MHD deformations of the uniform field , and let ℬff ⊂ ℬ be constituted of the fields of ℬ which are force-free. We derive topological lower and upper bounds on the energy of the fields in ℬ and in ℬff, respectively. We also establish new formulae for the field line helicity and the relative magnetic helicity of the fields in ℬ, and give a sufficient criterion of ideal MHD linear stability for the fields in ℬff.

A complete topological invariant for braided magnetic fields

A topological flux function is introduced to quantify the topology of magnetic braids: non-zero line-tied magnetic fields whose field lines all connect between two boundaries. This scalar function is an ideal invariant defined on a cross-section of the magnetic field, whose integral over the cross-section yields the relative magnetic helicity. Recognising that the topological flux function is an action in the Hamiltonian formulation of the field line equations, a simple formula for its differential is obtained. We use this to prove that the topological flux function uniquely characterises the field line mapping and hence the magnetic topology. A simple example is presented.

A COMBINED STUDY OF PHOTOSPHERIC MAGNETIC AND CURRENT HELICITIES AND SUBSURFACE KINETIC HELICITIES OF SOLAR ACTIVE REGIONS DURING 2006-2013

We compare the average photospheric current helicity $H_c$, photospheric twist parameter $\alpha$ (a well-known proxy for the full relative magnetic helicity), and subsurface kinetic helicity $H_k$ for 194 active regions observed between 2006-2013. We use 2440 Hinode photospheric vector magnetograms, and the corresponding subsurface fluid velocity data derived from GONG (2006-2012) and HMI (2010-2013) Dopplergrams. We find a significant hemispheric bias in all three parameters. The subsurface kinetic helicity is preferentially positive in the southern hemisphere and negative in the northern hemisphere. The photospheric current helicity and the $\alpha$ parameter have the same bias for strong fields ($|B|>1000$ G) and no significant bias for weak fields (100 G $<|B|<500$ G). We find no significant region-by-region correlation between the subsurface kinetic helicity and either the strong-field current helicity or $\alpha$. Subsurface fluid motions of a given handedness correspond to photospheric helicities of both signs in approximately equal numbers. However, common variations appear in annual averages of these quantities over all regions. Furthermore, in a subset of 77 regions we find significant correlations between the temporal profiles of the subsurface and photospheric helicities. In these cases, the sign of the linear correlation coefficient matches the sign relationship between the helicities, indicating that the photospheric magnetic field twist is sensitive to the twisting motions below the surface.