The evolution of coronal mass ejections in the inner heliosphere : Implementing the spheromak model with EUHFORIA

Aims. We introduce a new model for coronal mass ejections (CMEs) that has been implemented in the magnetohydrodynamics (MHD) inner heliosphere model EUHFORIA. Utilising a linear force-free spheromak (LFFS) solution, the model provides an intrinsic magnetic field structure for the CME. As a result, the new model has the potential to predict the magnetic components of CMEs at Earth. In this paper, we present the implementation of the new model and show the capability of the new model. Methods. We present initial validation runs for the new magnetised CME model by considering the same set of events as used in the initial validation run of EUHFORIA that employed the Cone model. In particular, we have focused on modelling the CME that was responsible for creating the largest geomagnetic disturbance (Dst index). Two scenarios are discussed: one where a single magnetised CME is launched and another in which we launch all five Earth-directed CMEs that were observed during the considered time period. Four out of the five CMEs were modelled using the Cone model. Results. In the first run, where the propagation of a single magnetized CME is considered, we find that the magnetic field components at Earth are well reproduced as compared to in-situ spacecraft data. Considering a virtual spacecraft that is separated approximately seven heliographic degrees from the position of Earth, we note that the centre of the magnetic cloud is missing Earth and a considerably larger magnetic field strength can be found when shifting to that location. For the second run, launching four Cone CMEs and one LFFS CME, we notice that the simulated magnetised CME is arriving at the same time as in the corresponding full Cone model run. We find that to achieve this, the speed of the CME needs to be reduced in order to compensate for the expansion of the CME due to the addition of the magnetic field inside the CME. The reduced initial speed of the CME and the added magnetic field structure give rise to a very similar propagation of the CME with approximately the same arrival time at 1 au. In contrast to the Cone model, however, the magnetised CME is able to predict the magnetic field components at Earth. However, due to the interaction between the Cone model CMEs and the magnetised CME, the magnetic field amplitude is significantly lower than for the run using a single magnetised CME. Conclusions. We have presented the LFFS model that is able to simulate and predict the magnetic field components and the propagation of magnetised CMEs in the inner heliosphere and at Earth. We note that shifting towards a virtual spacecraft in the neighbourhood of Earth can give rise to much stronger magnetic field components. This gives the option of adding a grid of virtual spacecrafts to give a range of values for the magnetic field components.Peer reviewe

The effect of drifts on the decay phase of SEP events

Aims. We study the effect of the magnetic gradient and curvature drifts on the pitch-angle dependent transport of solar energetic particles (SEPs) in the heliosphere, focussing on similar to 3-36 MeV protons. By considering observers located at different positions in the heliosphere, we investigate how drifts may alter the measured intensity-time profiles and energy spectra. We focus on the decay phase of solar energetic proton events in which a temporal invariant spectrum and disappearing spatial intensity gradients are often observed; a phenomenon known as the "reservoir effect" or the "SEP flood". We study the effects of drifts by propagating particles both in nominal and non-nominal solar wind conditions.Methods. We used a three-dimensional (3D) particle transport model, solving the focused transport equation extended with the effect of particle drifts in the spatial term. Nominal Parker solar wind configurations of different speeds and a magnetohydrodynamic (MHD) generated solar wind containing a corotating interaction region (CIR) were considered. The latter configuration gives rise to a magnetic bottle structure, with one bottleneck at the Sun and the other at the CIR. We inject protons from a fixed source at 0.1 AU, the inner boundary of the MHD model.Results. When the drift induced particle net-flux is zero, the modelled intensity-time profiles obtained at different radial distances along an IMF line show the same intensity fall-off after the prompt phase of the particle event, which is in accordance with the SEP flood phenomenon. However, observers magnetically connected close to the edges of the particle injection site can experience, as a result of drifts, a sudden drop in the intensities occurring at different times for different energies such that no SEP flood phenomenon is established. In the magnetic bottle structure, this effect is enhanced due to the presence of magnetic field gradients strengthening the nominal particle drifts. Moreover, anisotropies can be large for observers that only receive particles through drifts, illustrating the importance of pitch-angle dependent 3D particle modelling. We observe that interplanetary cross-field diffusion can mitigate the effects of particle drifts.Conclusions. Particle drifts can substantially modify the decay phase of SEP events, especially if the solar wind contains compression regions or shock waves where the drifts are enhanced. This is, for example, the case for our CIR solar wind configuration generated with a 3D MHD model, where the effect of drifts is strong. A similar decay rate in different energy channels and for different observers requires the mitigation of the effect of drifts. One way to accomplish this is through interplanetary cross-field diffusion, suggesting thus a way to determine a minimum value for the cross-field diffusion strength.Peer reviewe

Observation-based modelling of magnetised coronal mass ejections with EUHFORIA

Context. Coronal mass ejections (CMEs) are the primary source of strong space weather disturbances at Earth. Their geo-effectiveness is largely determined by their dynamic pressure and internal magnetic fields, for which reliable predictions at Earth are not possible with traditional cone CME models. Aims. We study two well-observed Earth-directed CMEs using the EUropean Heliospheric FORecasting Information Asset (EUH-FORIA) model, testing for the first time the predictive capabilities of a linear force-free spheromak CME model initialised using parameters derived from remote-sensing observations. Methods. Using observation-based CME input parameters, we performed magnetohydrodynamic simulations of the events with EU-HFORIA, using the cone and spheromak CME models. Results. Simulations show that spheromak CMEs propagate faster than cone CMEs when initialised with the same kinematic parameters. We interpret these differences as the result of different Lorentz forces acting within cone and spheromak CMEs, which lead to different CME expansions in the heliosphere. Such discrepancies can be mitigated by initialising spheromak CMEs with a reduced speed corresponding to the radial speed only. Results at Earth provide evidence that the spheromak model improves the predictions of B (B-z) by up to 12-60 (22-40) percentage points compared to a cone model. Considering virtual spacecraft located within +/- 10 degrees around Earth, B (Bz) predictions reach 45-70% (58-78%) of the observed peak values. The spheromak model shows inaccurate predictions of the magnetic field parameters at Earth for CMEs propagating away from the Sun-Earth line. Conclusions. The spheromak model successfully predicts the CME properties and arrival time in the case of strictly Earth-directed events, while modelling CMEs propagating away from the Sun-Earth line requires extra care due to limitations related to the assumed spherical shape. The spatial variability of modelling results and the typical uncertainties in the reconstructed CME direction advocate the need to consider predictions at Earth and at virtual spacecraft located around it.Peer reviewe

A type II solar radio burst without a coronal mass ejection

The Sun produces the most powerful explosions in the solar system, solar flares, that can also be accompanied by large eruptions of magnetised plasma, coronal mass ejections (CMEs). These processes can accelerate electron beams up to relativistic energies through magnetic reconnection processes during solar flares and CME-driven shocks. Energetic electron beams can in turn generate radio bursts through the plasma emission mechanism. CME shocks, in particular, are usually associated with type II solar radio bursts. However, on a few occasions, type II bursts have been reported to occur either in the absence of CMEs or shown to be more likely related with the flaring process. It is currently an open question how a shock generating type II bursts forms without the occurrence of a CME eruption. Here, we aim to determine the physical mechanism responsible for a type II burst which occurs in the absence a CME. By using radio imaging from the Nan{\c c}ay Radioheliograph, combined with observations from the Solar Dynamics Observatory and the Solar Terrestrial Relations Observatory spacecraft, we investigate the origin of a type II radio burst that appears to have no temporal association with a white-light CME. We identify a typical type II radio burst with band-split structure that is associated with a C-class solar flare. The type II burst source is located above the flaring active region and ahead of disturbed coronal loops observed in extreme ultraviolet images. The type II is also preceded by type III radio bursts, some of which are in fact J-bursts indicating that accelerated electron beams do not all escape along open field lines. The type II sources show single-frequency movement towards the flaring active region. The type II is located ahead of a faint extreme-ultraviolet (EUV) front propagating through the corona.Comment: 10 pages, 8 figure

Shock-accelerated electrons during the fast expansion of a coronal mass ejection

Publisher Copyright: © D. E. Morosan et al. 2022.Context. Some of of the most prominent sources for energetic particles in our Solar System are huge eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs), which usually drive shocks that accelerate charged particles up to relativistic energies. In particular, energetic electron beams can generate radio bursts through the plasma emission mechanism. The main types of bursts associated with CME shocks are type II and herringbone bursts. However, it is currently unknown where early accelerated electrons that produce metric type II bursts and herringbones propagate and when they escape the solar atmosphere. Aims. Here, we investigate the acceleration location, escape, and propagation directions of electron beams during the early evolution of a strongly expanding CME-driven shock wave associated with herrinbgone bursts. Methods. We used ground-based radio observations from the Nançay Radioheliograph combined with space-based extreme-ultraviolet and white-light observations from the Solar Dynamics Observatory and and the Solar Terrestrial Relations Observatory. We produced a three-dimensional (3D) representation of the electron acceleration locations which, combined with results from magneto-hydrodynamic (MHD) models of the solar corona, was used to investigate the origin of the herringbone bursts observed. Results. Multiple herringbone bursts are found close to the CME flank in plane-of-sky images. Some of these herringbone bursts have unusual inverted J shapes and opposite drifting herringbones also show opposite senses of circular polarisation. By using a 3D approach combined with the radio properties of the observed bursts, we find evidence that the first radio emission in the CME eruption most likely originates from electrons that initially propagate in regions of low Alfvén speeds and along closed magnetic field lines forming a coronal streamer. The radio emission appears to propagate in the same direction as a coronal wave in three dimensions. Conclusions. The CME appears to inevitably expand into a coronal streamer where it meets ideal conditions to generate a fast shock which, in turn, can accelerate electrons. However, at low coronal heights, the streamer consists of exclusively closed field lines indicating that the early accelerated electron beams do not escape. This is in contrast with electrons which, in later stages, escape the corona so that they are detected by spacecraft.Peer reviewe

Three-dimensional reconstruction of multiple particle acceleration regions during a coronal mass ejection

Context. Some of the most prominent sources for particle acceleration in our Solar System are large eruptions of magnetised plasma from the Sun called coronal mass ejections (CMEs). These accelerated particles can generate radio emission through various mechanisms. Aims. CMEs are often accompanied by a variety of solar radio bursts with different shapes and characteristics in dynamic spectra. Radio bursts directly associated with CMEs often show movement in the direction of CME expansion. Here, we aim to determine the emission mechanism of multiple moving radio bursts that accompanied a flare and CME that took place on 14 June 2012. Methods. We used radio imaging from the Nancay Radioheliograph, combined with observations from the Solar Dynamics Observatory and Solar Terrestrial Relations Observatory spacecraft, to analyse these moving radio bursts in order to determine their emission mechanism and three-dimensional (3D) location with respect to the expanding CME. Results. In using a 3D representation of the particle acceleration locations in relation to the overlying coronal magnetic field and the CME propagation, for the first time, we provide evidence that these moving radio bursts originate near the CME flanks and that some are possible signatures of shock-accelerated electrons following the fast CME expansion in the low corona. Conclusions. The moving radio bursts, as well as other stationary bursts observed during the eruption, occur simultaneously with a type IV continuum in dynamic spectra, which is not usually associated with emission at the CME flanks. Our results show that moving radio bursts that could traditionally be classified as moving type IVs can represent shock signatures associated with CME flanks or plasma emission inside the CME behind its flanks, which are closely related to the lateral expansion of the CME in the low corona. In addition, the acceleration of electrons generating this radio emission appears to be favoured at the CME flanks, where the CME encounters coronal streamers and open field regions.Peer reviewe

Reconstructing Coronal Hole Areas With EUHFORIA and Adapted WSA Model : Optimizing the Model Parameters

The adopted Wang-Sheeley-Arge (WSA) model embedded in EUHFORIA (EUropean Heliospheric FORecasting Information Asset) is compared to EUV observations. According to the standard paradigm, coronal holes are sources of open flux; thus, we use remote sensing EUV observations and CATCH (Collection of Analysis Tools for Coronal Holes) to extract CH areas and compare them to the open flux areas modeled by EUHFORIA. From the adopted WSA model we employ only the Potential Field Source Surface (PFSS) model for the inner corona and the Schatten Current Sheet (SCS) model for the outer (PFSS+SCS). The height, R-ss, of the outer boundary of the PFSS, known as the source surface, and the height, R-i, of the inner boundary of the SCS are important parameters affecting the modeled CH areas. We investigate the impact the two model parameters can have in the modeled results. We vary R-ss within the interval [1.4, 3.2]R-circle dot with a step of 0.1R(circle dot), and R-i within the interval [1.3, 2.8]R-circle dot with the same step, and the condition that R-iPeer reviewe

Estimating the Magnetic Structure of an Erupting CME Flux Rope From AR12158 Using Data-Driven Modeling

We investigate here the magnetic properties of a large-scale magnetic flux rope related to a coronal mass ejection (CME) that erupted from the Sun on September 12, 2014 and produced a well-defined flux rope in interplanetary space on September 14-15, 2014. We apply a fully data-driven and time-dependent magnetofrictional method (TMFM) using Solar Dynamics Observatory (SDO) magnetograms as the lower boundary condition. The simulation self-consistently produces a coherent flux rope and its ejection from the simulation domain. This paper describes the identification of the flux rope from the simulation data and defining its key parameters (e.g., twist and magnetic flux). We define the axial magnetic flux of the flux rope and the magnetic field time series from at the apex and at different distances from the apex of the flux rope. Our analysis shows that TMFM yields axial magnetic flux values that are in agreement with several observational proxies. The extracted magnetic field time series do not match well with in-situ components in direct comparison presumably due to interplanetary evolution and northward propagation of the CME. The study emphasizes also that magnetic field time-series are strongly dependent on how the flux rope is intercepted which presents a challenge for space weather forecasting.Peer reviewe

The impact of coronal hole characteristics and solar cycle activity in reconstructing coronal holes with EUHFORIA

Modelling with high accuracy the open magnetic field and the fast solar wind in the heliosphere is essential for space weather forecasting purposes. Primary sources of open magnetic field flux are Coronal Holes (CH), uni-polar regions that appear as dark patches in the solar corona when observed in X-ray and extreme-ultraviolet (EUV) images due to having significantly lower density and temperature to their surroundings. Therefore, when assessing how well the open magnetic field and the fast solar wind are modelled one can look at how well the model performs on one of its fundamental functions, that of reconstructing coronal hole areas. In this study we investigate how the CH morphology (i.e. latitudinal position of the centre of mass, area, intensity, elongation) and the solar variability, from high to low activity periods, can affect the results. We also investigated the possibility that the model is reconstructing CHs that are systematically shifted with respect to their observed position. The study is applied on 15 CHs exhibiting different latitudinal position and geometry. We compare the modelled CH areas with boundaries obtained by remote sensing EUV observations using the CATCH tool (Collection of Analysis Tools for Coronal Holes). We found no apparent effect of the CH characteristics on the modelling capabilities. In addition, solar cycle activity seems not to have any effect either. However, we emphasize that our sample is small and this outcome highlights the need for an extended research.Peer reviewe

Probing the Effect of Cadence on the Estimates of Photospheric Energy and Helicity Injections in Eruptive Active Region NOAA AR 11158

We study how the input-data cadence affects the photospheric energy and helicity injection estimates in eruptive NOAA Active Region 11158. We sample the novel 2.25-minute vector magnetogram and Dopplergram data from the Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) spacecraft to create input datasets of variable cadences ranging from 2.25 minutes to 24 hours. We employ state-of-the-art data processing, velocity, and electric-field inversion methods for deriving estimates of the energy and helicity injections from these datasets. We find that the electric-field inversion methods that reproduce the observed magnetic-field evolution through the use of Faraday's law are more stable against variable cadence: the PDFI (PTD-Doppler-FLCT-Ideal, where PTD refers to Poloidal-Toroidal Decomposition, and FLCT to Fourier Local Correlation Tracking) electric-field inversion method produces consistent injection estimates for cadences from 2.25 minutes up to two hours, implying that the photospheric processes acting on time scales below two hours contribute little to the injections, or that they are below the sensitivity of the input data and the PDFI method. On other hand, the electric-field estimate derived from the output of DAVE4VM (Differential Affine Velocity Estimator for Vector Magnetograms), which does not fulfill Faraday's law exactly, produces significant variations in the energy and helicity injection estimates in the 2.25 minutes - two hours cadence range. We also present a third, novel DAVE4VM-based electric-field estimate, which corrects the poor inductivity of the raw DAVE4VM estimate. This method is less sensitive to the changes of cadence, but it still faces significant issues for the lowest of considered cadences (two hours). We find several potential problems in both PDFI- and DAVE4VM-based injection estimates and conclude that the quality of both should be surveyed further in controlled environments.Peer reviewe