Search

We investigate which scalar quantity or quantities can best predict the loss of equilibrium and subsequent eruption of magnetic flux ropes in the solar corona. Our models are initialized with a potential magnetic arcade, which is then evolved by means of two effects on the lower boundary: firstly a gradual shearing of the arcade, modelling differential rotation on the solar surface, and secondly supergranular diffusion. These result in flux cancellation at the polarity inversion line and the formation of a twisted flux rope. We use three model setups: full magnetohydrodynamics (MHD) in cartesian coordinates, and the magnetofrictional model in both cartesian and polar coordinates. The flux ropes are translationally-invariant, allowing for very fast computational times and thus a comprehensive parameter study, comprising hundreds of simulations and thousands of eruptions. Similar flux rope behavior is observed using either magnetofriction or MHD, and there are several scalar criteria that could be used as proxies for eruptivity. The most consistent predictor of eruptions in either model is the squared current in the axial direction of the rope, normalised by the relative helicity, although a variation on the previously proposed 'eruptivity index' is also found to perform well in both the magnetofrictional and MHD simulations., Comment: 20 pages, 11 figures, 1 table. Accepted to the Astrophysical Journal on 13th August 2023

We apply the magneto-frictional approach to investigate which quantity or quantities can best predict the loss of equilibrium of a translationally-invariant magnetic flux rope. The flux rope is produced self-consistently by flux cancellation combined with gradual footpoint shearing of a coronal arcade which is open at the outer boundary. This models the magnetic field in decaying active regions on the Sun. Such a model permits two types of eruption: episodic small events caused by shearing and relaxation of the overlying arcade, and major eruptions of the main low-lying coronal flux rope. Through a parameter study, we find that the major eruptions are best predicted not by individual quantities but by thresholds in the ratios of squared rope current to either magnetic energy or relative magnetic helicity. We show how to appropriately define the latter quantity for translationally-invariant magnetic fields, along with a related eruptivity index that has recently been introduced for three-dimensional magnetic fields. In contrast to previous configurations studied, we find that the eruptivity index has only a weak predictive skill, and in fact is lower prior to eruption, rather than higher. This is because the overlying background magnetic field has the same direction as the arcade itself. Thus we propose that there are a whole class of solar eruptions that cannot be predicted by a high eruptivity index., Comment: 25 pages, 11 figures. To be published in 'Frontiers in Astronomy and Space Sciences', accepted for publication 8th March 2022

Given a known radial magnetic field distribution on the Sun's photospheric surface, there exist well-established methods for computing a potential magnetic field in the corona above. Such potential fields are routinely used as input to solar wind models, and to initialize magneto-frictional or full magnetohydrodynamic simulations of the coronal and heliospheric magnetic fields. We describe an improved magnetic field model which calculates a magneto-frictional equilibrium with an imposed solar wind profile (which can be Parker's solar wind solution, or any reasonable equivalent). These `outflow fields' appear to approximate the real coronal magnetic field more closely than a potential field, take a similar time to compute, and avoid the need to impose an artificial source surface. Thus they provide a practical alternative to the potential field model for initializing time-evolving simulations or modeling the heliospheric magnetic field. We give an open-source Python implementation in spherical coordinates and apply the model to data from Solar Cycle 24. The outflow tends to increase the open magnetic flux compared to the potential field model, reducing the well known discrepancy with in situ observations., Comment: 16 Pages, 6 Figures. To be published in the Astrophysical Journal (accepted 1/10/21)