A new generation of coronagraphs to study the solar wind and CMEs are being developed and launched. These coronagraphs will heavily rely on multi-channel observations where visible light (VL) and UV-EUV observations provide new plasma diagnostics. One of these instruments, Metis on board ESA-Solar Orbiter, will simultaneously observe VL and the UV Lyman-$\alpha$ line. The number of neutral Hydrogen atoms (a small fraction of coronal protons) is a key parameter for deriving plasma properties such as temperature from the observed Lyman-$\alpha$ line intensity. However, these measurements are significantly affected if non-equilibrium ionisation effects occur. The aim of this work is to determine if non-equilibrium ionisation effects are relevant in CMEs and in particular when and in which regions of the CME. We use a magneto-hydrodynamic simulation of a magnetic flux rope ejection to generate a CME. From this we then reconstruct the ionisation state of Hydrogen atoms in the CME by evaluating both the advection of neutral and ionised Hydrogen atoms and the ionisation and recombination rates in the MHD simulation. We find that the equilibrium ionisation assumption holds mostly in the core of the CME. In contrast non-equilibrium ionisation effects are significant at the CME front, where we find about 100 times more neutral Hydrogen atoms than prescribed by ionisation equilibrium conditions, even if this neutral Hydrogen excess might be difficult to identify due to projection effects. This work provides key information for the development of a new generation of diagnostic techniques that aim at combining visible light and Lyman-$\alpha$ line emissions. The results show that non-ionisation equilibrium effects need to be considered when we analyse CME fronts. To incorrectly assume equilibrium ionisation in these regions would lead to a systematic underestimate of plasma temperatures.

Phase mixing of Alfv\'en waves has been studied extensively as a possible coronal heating mechanism but without the full thermodynamic consequences considered self-consistently. It has been argued that in some cases, the thermodynamic feedback of the heating could substantially affect the transverse density gradient and even inhibit the phase mixing process. In this paper, we use MHD simulations with the appropriate thermodynamical terms included to quantify the evaporation following heating by phase mixing of Alfv\'en waves in a coronal loop and the effect of this evaporation on the transverse density profile. The numerical simulations were performed using the Lare2D code. We set up a 2D loop model consisting of a field-aligned thermodynamic equilibrium and a cross-field (background) heating profile. A continuous, sinusoidal, high-frequency Alfv\'en wave driver was implemented. As the Alfv\'en waves propagate along the field, they undergo phase mixing due to the cross-field density gradient in the coronal part of the loop. We investigated the presence of field-aligned flows, heating from the dissipation of the phase-mixed Alfv\'en waves, and the subsequent evaporation from the lower atmosphere. We find that phase mixing of Alfv\'en waves leads to modest heating in the shell regions of the loop and evaporation of chromospheric material into the corona with upflows of the order of only 5-20 m/s. Although the evaporation leads to a mass increase in the shell regions of the loop, the effect on the density gradient and, hence, on the phase mixing process, is insignificant. This paper self-consistently investigates the effect of chromospheric evaporation on the cross-field density gradient and the phase mixing process in a coronal loop. We found that the effects in our particular setup (small amplitude, high frequency waves) are too small to significantly change the density gradient.

One of the main goals of solar physics is the timely identification of eruptive active regions. Space missions such as Solar Orbiter or future Space Weather forecasting missions would largely benefit from this achievement. Our aim is to produce a relatively simple technique that can provide real time indications or predictions that an active region will produce an eruption. We expand on the theoretical work of \citet{Pagano2019fp} that was able to distinguish eruptive from non-eruptive active regions. From this we introduce a new operational metric that uses a combination of observed line-of-sight magnetograms, 3D data-driven simulations and the projection of the 3D simulations forward in time. Results show that the new metric correctly distinguishes active regions as eruptive when observable signatures of eruption have been identified and as non-eruptive when there are no observable signatures of eruption. After successfully distinguishing eruptive from non-eruptive active regions we illustrate how this metric may be used in a "real-time" operational sense were three levels of warning are categorised. These categories are: high risk (red), medium risk (amber) and low risk (green) of eruption. Through considering individual cases we find that the separation into eruptive and non-eruptive active regions is more robust the longer the time series of observed magnetograms used to simulate the build up of magnetic stress and free magnetic energy within the active region. Finally, we conclude that this proof of concept study delivers promising results where the ability to categorise the risk of an eruption is a major achievement.

Active regions are the source of the majority of magnetic flux rope ejections that become Coronal Mass Ejections (CMEs). To identify in advance which active regions will produce an ejection is key for both space weather prediction tools and future science missions such as Solar Orbiter. The aim of this study is to develop a new technique to identify which active regions are more likely to generate magnetic flux rope ejections. The new technique will aim to: (i) produce timely space weather warnings and (ii) open the way to a qualified selection of observational targets for space-borne instruments. We use a data-driven Non-linear Force-Free Field (NLFFF) model to describe the 3D evolution of the magnetic field of a set of active regions. We determine a metric to distinguish eruptive from non-eruptive active regions based on the Lorentz force. Furthermore, using a subset of the observed magnetograms, we run a series of simulations to test whether the time evolution of the metric can be predicted. The identified metric successfully differentiates active regions observed to produce eruptions from the non-eruptive ones in our data sample. A meaningful prediction of the metric can be made between 6 to 16 hours in advance. This initial study presents an interesting first step in the prediction of CME onset using only LOS magnetogram observations combined with NLFFF modelling. Future studies will address how to generalise the model such that it can be used in a more operational sense and for a variety of simulation approaches.

The determination of the speed of Coronal Mass Ejections (CMEs) is usually done by tracking brighter features (such as the CME front and core) in visible light coronagraphic images and by deriving unidimensional profiles of the CME speed as a function of altitude or time. Nevertheless, CMEs are usually characterized by the presence of significant density inhomogeneities propagating outward with different radial and latitudinal projected speeds, resulting in a complex evolution eventually forming the Interplanetary CME. In this work, we demonstrate for the first time how coronagraphic image sequences can be analyzed with cross-correlation technique to derive 2D maps of the almost instantaneous plasma speed distribution within the body of CMEs. The technique is first tested with the analysis of synthetic data, and then applied to real observations. Results from this work allow to characterize the distribution and time evolution of kinetic energy inside CMEs, as well as the mechanical energy (combined with the kinetic and potential energy) partition between the core and front of the CME. In the future, CMEs will be observed by two channels (VL and UV Ly-$\alpha$) coronagraphs, such as Metis on-board ESA Solar Orbiter mission as well as Ly-$\alpha$ Solar Telescope (LST) on-board Chinese Advanced Space-based Solar Observatory (ASO-S) mission. Our results will help the analysis of these future observations, helping in particular to take into account the 2D distribution of Ly-$\alpha$ Doppler dimming effect.

MHD waves are ubiquitous in the solar corona where the highly structured magnetic fields provide efficient wave guides for their propagation. While MHD waves have been observed originating from lower layers of the solar atmosphere, recent studies have shown that some can be generated in-situ by the collision of dense counter propagating flows. In this theoretical study, we analyse the mechanism that triggers the propagation of kink and sausage modes in the solar corona, following the collision of counter-propagating flows and how the properties of the flows affect the properties of the generated waves. In order to study in detail such a mechanism we run a series of ideal 2D and 3D MHD simulations where we vary the properties of the counter-propagating flows and by means of a simple technique to estimate the amplitudes of the kink and sausage modes, we investigate their role on the generation and propagation of the MHD waves. We find that the amplitude of the waves is largely dependent on the kinetic energy of the flows and that the onset of kink or sausage modes depends on the asymmetries between the colliding blobs. Moreover, the initial wavelength of the MHD waves is associated with the magnetic configuration resulting from the collision of the flows. We also find that genuine 3D systems respond with smaller wave amplitudes. In this study, we present a parameter space description of the mechanism that leads to the generation of MHD waves from the collision of flows in the corona. Future observations of these waves can be used to understand the properties of the plasma and magnetic field of the solar corona.

Whilst there are observational indications that transverse MHD waves carry enough energy to maintain the thermal structure of the solar corona, it is not clear whether such energy can be efficiently and effectively converted into heating. Phase-mixing of Alfv\'en waves is considered a candidate mechanism, as it can develop transverse gradient where magnetic energy can be converted into thermal energy. However, phase-mixing is a process that crucially depends on the amplitude and period of the transverse oscillations, and only recently have we obtained a complete measurement of the power spectrum for transverse oscillations in the corona. We aim to investigate the heating generated by phase-mixing of transverse oscillations triggered by buffeting of a coronal loop that follows from the observed coronal power spectrum as well as the impact of these persistent oscillations on the structure of coronal loops. We consider a 3D MHD model of an active region coronal loop and we perturb its footpoints with a 2D horizontal driver that represents a random buffeting motion of the loop footpoints. Our driver is composed of 1000 pulses superimposed to generate the observed power spectrum. We find that the heating supply from the observed power spectrum in the solar corona through phase mixing is not sufficient to maintain the million degree active region solar corona. We also find that the development of Kelvin-Helmholtz instabilities could be a common phenomenon in coronal loops that could affect their apparent life time. This study concludes that is unlikely that phase-mixing of Alfv\'en waves resulting from an observed power spectrum of transverse coronal loop oscillations can heat the active region solar corona. However, transverse waves could play an important role in the development of small scale structures.

Transverse MHD waves permeate the solar atmosphere and are a candidate for coronal heating. However, the origin of these waves is still unclear. In this work, we analyse coordinated observations from \textit{Hinode}/SOT and \textit{IRIS} of a prominence/coronal rain loop-like structure at the limb of the Sun. Cool and dense downflows and upflows are observed along the structure. A collision between a downward and an upward flow with an estimated energy flux of $10^{7}-10^{8}$~erg~cm$^{-2}$~s$^{-1}$ is observed to generate oscillatory transverse perturbations of the strands with an estimated $\approx40~$km~s$^{-1}$ total amplitude, and a short-lived brightening event with the plasma temperature increasing to at least $10^{5}~$K. We interpret this response as sausage and kink transverse MHD waves based on 2D MHD simulations of plasma flow collision. The lengths, density and velocity differences between the colliding clumps and the strength of the magnetic field are major parameters defining the response to the collision. The presence of asymmetry between the clumps (angle of impact surface and/or offset of flowing axis) is crucial to generate a kink mode. Using the observed values we successfully reproduce the observed transverse perturbations and brightening, and show adiabatic heating to coronal temperatures. The numerical modelling indicates that the plasma $\beta$ in this loop-like structure is confined between $0.09$ and $0.36$. These results suggest that such collisions from counter-streaming flows can be a source of in-situ transverse MHD waves, and that for cool and dense prominence conditions such waves could have significant amplitudes.

Kink oscillations of a coronal loop are observed and studied in detail because they provide a unique probe into the structure of coronal loops through MHD seismology and a potential test of coronal heating through the phase-mixing of Alfv\'en waves. In particular, recent observations show that standing oscillations of loops often involve also higher harmonics, beside the fundamental mode. The damping of these kink oscillations is explained by mode coupling with Alfv\'en waves. We investigate the consequences for wave-based coronal heating of higher harmonics and what coronal heating observational signatures we may use to infer the presence of higher harmonic kink oscillations. We perform a set of non-ideal MHD simulations where the damping of the kink oscillation of a flux tube via mode coupling is modelled. Our MHD simulation parameters are based on the seismological inversion of an observation for which the first three harmonics are detected. We study the phase-mixing of Alfv\'en waves that leads to the deposition of heat in the system, and we apply the seismological inversion techniques to the MHD simulation output. We find that the heating due to phase-mixing of the Alfv\'en waves triggered by the damping of the kink oscillation is relatively small, however we can illustrate i) how the heating location drifts due to the subsequent damping of lower order harmonics. We also address the role of the higher order harmonics and the width of the boundary shell in the energy deposition. We conclude that the coronal heating due to phase-mixing seems not to provide enough energy to maintain the thermal structure of the solar corona even when multi-harmonics oscillations are included, and these oscillations play an inhibiting role in the development of smaller scale structures.

In recent years, space weather research has focused on developing modelling techniques to predict the arrival time and properties of coronal mass ejections (CMEs) at the Earth. The aim of this paper is to propose a new modelling technique suitable for the next generation of Space Weather predictive tools that is both efficient and accurate. The aim of the new approach is to provide interplanetary space weather forecasting models with accurate time dependent boundary conditions of erupting magnetic flux ropes in the upper solar corona. To produce boundary conditions, we couple two different modelling techniques, MHD simulations and a quasi-static non-potential evolution model. Both are applied on a spatial domain that covers the entire solar surface. The non-potential model uses a time series of observed synoptic magnetograms to drive the non-potential quasi-static evolution of the coronal magnetic field. This allows us to follow the formation and loss of equilibrium of magnetic flux ropes. Following this a MHD simulation captures the dynamic evolution of the erupting flux rope. The present paper focuses on the MHD simulations that follow the ejection of magnetic flux ropes to 4$R_\odot$. We first propose a technique for specifying the pre-eruptive plasma properties in the corona. Next, time dependent MHD simulations describe the ejection of two magnetic flux ropes, that produce time dependent boundary conditions for the magnetic field and plasma at 4$R_{\odot}$. In the present paper, we show that the dual use of quasi-static non-potential magnetic field simulations and full time dependent MHD simulations can produce realistic inhomogeneous boundary conditions for space weather forecasting tools. Before a fully operational model can be produced there are a number of technical and scientific challenges that still need to be addressed.

Phase-mixing of Alfv\'en waves in the solar corona has been identified as one possible candidate to explain coronal heating. While this scenario is supported by observations of ubiquitous oscillations in the corona carrying sufficient wave energy and by theoretical models that have described the concentration of energy in small scale structures, it is still unclear whether this wave energy can maintain the million degree solar corona. The aim of this work is to assess how much energy can be converted by a phase-mixing process triggered by the propagation of Alfv\'enic waves in a cylindric coronal structure, such as a coronal loop, and to estimate the impact on the coronal heating. To this end, we run 3D MHD simulations of a magnetised cylinder where the Alfv\'en speed varies through a boundary shell and a footpoint driver is set to trigger kink modes which mode couple to torsional Alfv\'en modes in the boundary shell. These Alfv\'en waves are expected to phase-mix and the system allows us to study the subsequent thermal energy deposition. We run a reference simulation to explain the main process and then we vary simulation parameters. When we take into consideration high values of magnetic resistivity and strong footpoint drivers, we find i) that phase-mixing leads to a temperature increase of the order of $10^5$ K or less, depending on the structure of the boundary shell, ii) that this energy is able to balance the radiative losses only in the localised region involved in the heating, iii) and how the boundary layer and the persistence of the driver influence the thermal structure of the system. Our conclusion is that due to the extreme physical parameters we adopted and the moderate impact on the heating of the system, it is unlikely that phase-mixing can contribute on a global scale to the heating of the solar corona.

Understanding the 3D structure of coronal mass ejections (CMEs) is crucial for understanding the nature and origin of solar eruptions. To derive information on the 3D structure of CMEs from white-light (total and polarized brightness) images, the polarization ratio technique is widely used. The soon-to-be-launched METIS coronagraph on board Solar Orbiter will use this technique to produce new polarimetric images. We determine the accuracy at which the position of the centre of mass, direction and speed of propagation, and the column density of the CME can be determined along the line of sight. We perform a 3D MHD simulation of a flux rope ejection where a CME is produced. From the simulation we (i) synthesize the corresponding METIS white-light (total and polarized brightness) images and (ii) apply the polarization ratio technique to these synthesized images and compare the results with the known density distribution from the MHD simulation. We find that the polarization ratio technique reproduces with high accuracy the position of the centre of mass along the line of sight. However, some errors are inherently associated with this determination. The polarization ratio technique also allows information to be derived on the real 3D direction of propagation of the CME. In addition, we find that the column density derived from white-light images is accurate and we propose an improved technique where the combined use of the polarization ratio technique and white-light images minimizes the error in the estimation of column densities. Our method allows us to thoroughly test the performance of the polarization ratio technique and allows a determination of the errors associated with it, which means that it can be used to quantify the results from the analysis of the forthcoming METIS observations in white light (total and polarized brightness).

This work is aimed at quantifying the uncertainties in the 3D reconstruction of the location of coronal mass ejections (CMEs) obtained with the polarization ratio technique. The method takes advantage of the different distributions along the line of sight (LOS) of total (tB) and polarized (pB) brightnesses to estimate the average location of the emitting plasma. To this end, we assumed two simple electron density distributions along the LOS (a constant density and Gaussian density profiles) for a plasma blob and synthesized the expected tB and pB for different distances $z$ of the blob from the plane of the sky (POS) and different projected altitudes $\rho$. Reconstructed locations of the blob along the LOS were thus compared with the real ones, allowing a precise determination of uncertainties in the method. Independently of the analytical density profile, when the blob is centered at a small distance from the POS (i.e. for limb CMEs) the distance from the POS starts to be significantly overestimated. Polarization ratio technique provides the LOS position of the center of mass of what we call folded density distribution, given by reflecting and summing in front of the POS the fraction of density profile located behind that plane. On the other hand, when the blob is far from the POS, but with very small projected altitudes (i.e. for halo CMEs, $\rho < 1.4$ R$_\odot$), the inferred distance from that plane is significantly underestimated. Better determination of the real blob position along the LOS is given for intermediate locations, and in particular when the blob is centered at an angle of $20^\circ$ from the POS. These result have important consequences not only for future 3D reconstruction of CMEs with polarization ratio technique, but also for the design of future coronagraphs aimed at providing a continuous monitoring of halo-CMEs for space weather prediction purposes.

Extreme ultraviolet (EUV) images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) are providing new insights into the early phase of CME evolution. Observations now show the ejection of magnetic flux ropes from the solar corona and how they evolve into CMEs. These observations are difficult to interpret in terms of basic physical mechanisms and quantities. To fully understand CMEs we need to compare equivalent quantities derived from both observations and theoretical models. To this end we aim to produce synthesised AIA observations from simluations of a flux rope ejection. To carry this out we include the role of thermal conduction and radiative losses, both of which are important for determining the temperature distribution of the solar corona during a CME. We perform a simulation where a flux rope is ejected from the solar corona. From the density and temperature of the plasma in the simulation we synthesise AIA observations. The emission is then integrated along the line of sight using the instrumental response function of AIA. We sythesise observations of AIA in the channels at 304 A, 171 A, 335 A, and 94 A. The synthesised observations show a number of features similar to actual observations and in particular reproduce the general development of CMEs in the low corona as observed by AIA. In particular we reproduce an erupting and expanding arcade in the 304 A and 171 A channels with a high density core. The ejection of a flux rope reproduces many of the features found in the AIA observations. This work is therefore a step forward in bridging the gap between observations and models, and can lead to more direct interpretations of EUV observations in terms of flux rope ejections. We plan to improve the model in future studies in order to perform a more quantitative comparison.

Our aim is to study the role of gravitational stratification on the propagation of CMEs. In particular, we assess how it influences the speed and shape of CMEs and under what conditions the flux rope ejection becomes a CME or when it is quenched. We ran a set of MHD simulations that adopt an eruptive initial magnetic configuration that has already been shown to be suitable for a flux rope ejection. We varied the temperature of the backgroud corona and the intensity of the initial magnetic field to tune the gravitational stratification and the amount of ejected magnetic flux. We used an automatic technique to track the expansion and the propagation of the magnetic flux rope in the MHD simulations. From the analysis of the parameter space, we evaluate the role of gravitational stratification on the CME speed and expansion. Our study shows that gravitational stratification plays a significant role in determining whether the flux rope ejection will turn into a full CME or whether the magnetic flux rope will stop in the corona. The CME speed is affected by the background corona where it travels faster when the corona is colder and when the initial magnetic field is more intense. The fastest CME we reproduce in our parameter space travels at 850 km/s. Moreover, the background gravitational stratification plays a role in the side expansion of the CME, and we find that when the background temperature is higher, the resulting shape of the CME is flattened more. Our study shows that although the initiation mechanisms of the CME are purely magnetic, the background coronal plasma plays a key role in the CME propagation, and full MHD models should be applied when one focusses especially on the production of a CME from a flux rope ejection.

We provide a guideline to interpret the UVCS emission lines (in particular O VI and Si XII) during shock wave propagation in the outer solar corona. We use a numerical MHD model performing a set of simulations of shock waves generated in the corona and from the result we compute the plasma emission for the O VI and Si XII including the effects of NEI. We analyze the radiative and spectral properties of our model with the support of a detailed radiation model including Doppler dimming and an analytical model for shocks, and, finally, we synthesize the expected O VI 1032A line profile. We explain several spectral features of the observations like the absence of discontinuities in the O VI emission during the shock passage, the brightening of Si XII emission and the width of the lines. We use our model also to give very simple and general predictions for the strength of the line wings due to the ions shock heating and on the line shape for Limb CMEs or Halo CMEs. The emission coming from post-shock region in the solar corona roughly agrees with the emission from a simple planar and adiabatic shock, but the effect of thermal conduction and the magnetic field may be important depending on the event parameters. Doppler dimming significantly influences the O VI emission while Si XII line brightens mainly because of the shock compression. Significant shock heating is responsible for the wide and faint component of the O VI line usually observed which may be taken as a shock signature in the solar corona.

We explore the importance of magnetic-field-oriented thermal conduction in the interaction of supernova remnant (SNR) shocks with radiative gas clouds and in determining the mass and energy exchange between the clouds and the hot surrounding medium. We perform 2.5D MHD simulations of a shock impacting on an isolated gas cloud, including anisotropic thermal conduction and radiative cooling; we consider the representative case of a Mach 50 shock impacting on a cloud ten-fold denser than the ambient medium. We consider different configurations of the ambient magnetic field and compare MHD models with or without the thermal conduction. The efficiency of the thermal conduction in the presence of magnetic field is, in general, reduced with respect to the unmagnetized case. The reduction factor strongly depends on the initial magnetic field orientation, and it is minimum when the magnetic field is initially aligned with the direction of shock propagation. The thermal conduction contributes to suppress hydrodynamic instabilities, reducing the mass mixing of the cloud and preserving the cloud from complete fragmentation. Depending on the magnetic field orientation, the heat conduction may determine a significant energy exchange between the cloud and the hot surrounding medium which, while remaining always at levels less than those in the unmagnetized case, leads to a progressive heating and evaporation of the cloud. This additional heating may contrast the radiative cooling of some parts of the cloud, preventing the onset of thermal instabilities.

Detailed hydrodynamic modeling explained several features of a fragment of the core of a Coronal Mass Ejection observed with SoHO/UVCS at 1.7 Ro on 12 December 1997, but some questions remained unsolved. We investigate the role of the magnetic fields in the thermal insulation and the expansion of an ejected fragment (cloud) traveling upwards in the outer corona. We perform MHD simulations including the effects of thermal conduction and radiative losses of a dense spherical or cylindrical cloud launched upwards in the outer corona, with various assumptions on the strength and topology of the ambient magnetic field; we also consider the case of a cylindrical cloud with an internal magnetic field component along its axis. We find that a weak ambient magnetic field (beta~20) with open topology provides both significant thermal insulation and large expansion. The cylindrical cloud expands more than the spherical one.