Ilia I. Roussev

In this paper, we performed magnetohydrodynamics numerical experiments to look into the dynamic behaviour of the current sheet (CS) between the coronal mass ejection (CME) and the associated solar flare, especially the CS oscillation and plasmoid motions in coronal conditions. During the evolution, the disrupting magnetic configuration becomes asymmetric first in the buffer region at the bottom of the CME bubble. The Rayleigh−Taylor instability in the buffer region and the deflected motion of the plasma driven by the termination shock at the bottom of the CME bubble cause the buffer region to oscillate around the y-axis. The local oscillation propagates downwards through the CS, prompting an overall CS oscillation. As the buffer region grows, the oscillation period becomes longer, increasing from about 30 s to about 16 min. Meanwhile, there is another separated oscillation with a period between 0.25 and 1.5 min in the cusp region of the flare generated by velocity shearing. The tearing mode instability yields formations of plasmoids inside the CS. The motions of all the plasmoids observed in the experiment accelerate, which implies that the large-scale CME/flare CS itself in the true eruptive event is filled with the diffusion region according the the standard theory of magnetic reconnection.

Context. It is widely accepted that photospheric shearing motions play an important role in triggering the initiation of coronal mass ejections (CMEs). Even so, there are events for which the source signatures are difficult to locate, while the CMEs can be clearly observed in coronagraph data. These events are therefore called ‘stealth’ CMEs. They are of particular interest to space weather forecasters, since eruptions are usually discarded from arrival predictions if they appear to be backsided, which means not presenting any clear low-coronal signatures on the visible solar disc. Such assumptions are not valid for stealth CMEs since they can originate from the front side of the Sun and be Earth-directed, but they remain undetected and can therefore trigger unpredicted geomagnetic storms. Aims. We numerically model and investigate the effects of shearing motion variations onto the resulting eruptions and we focus in particular on obtaining a stealth CME in the trailing current sheet of a previous ejection. Methods. We used the 2.5D magnetohydrodynamics package of the code MPI-AMRVAC to numerically simulate consecutive CMEs by imposing shearing motions onto the inner boundary, which represents, in our case, the low corona. The initial magnetic configuration consists of a triple arcade structure embedded into a bimodal solar wind, and the sheared polarity inversion line is found in the southern loop system. The mesh was continuously adapted through a refinement method that applies to current carrying structures, allowing us to easily track the CMEs in high resolution, without resolving the grid in the entire domain. We also compared the obtained eruptions with the observed directions of propagation, determined using a forward modelling reconstruction technique based on a graduated cylindrical shell geometry, of an initial multiple coronal mass ejection (MCME) event that occurred in September 2009. We further analysed the simulated ejections by tracking the centre of their flux ropes in latitude and their total speed. Radial Poynting flux computation was employed as well to follow the evolution of electromagnetic energy introduced into the system. Results. Changes within 1% in the shearing speed result in three different scenarios for the second CME, although the preceding eruption seems insusceptible to such small variations. Depending on the applied shearing speed, we thus obtain a failed eruption, a stealth, or a CME driven by the imposed shear, as the second ejection. The dynamics of all eruptions are compared with the observed directions of propagation of an MCME event and a good correlation is achieved. The Poynting flux analysis reveals the temporal variation of the important steps of eruptions. Conclusions. For the first time, a stealth CME is simulated in the aftermath of a first eruption, originating from an asymmetric streamer configuration, through changes in the applied shearing speed, indicating it is not necessary for a closed streamer to exist high in the corona for such an event to occur. We also emphasise the high sensitivity of the corona to small changes in motions at the photosphere, or in our simulations, at the low corona.

In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona.