One of the very common in situ signatures of interplanetary coronal mass ejections (ICMEs), as well as other interplanetary transients, are Forbush decreases (FDs), i.e. short-term reductions in the galactic cosmic ray (GCR) flux. A two-step FD is often regarded as a textbook example, which presumably owes its specific morphology to the fact that the measuring instrument passed through the ICME head-on, encountering first the shock front (if developed), then the sheath and finally the CME magnetic structure. The interaction of GCRs and the shock/sheath region, as well as the CME magnetic structure, occurs all the way from Sun to Earth, therefore, FDs are expected to reflect the evolutionary properties of CMEs and their sheaths. We apply modelling to different ICME regions in order to obtain a generic two-step FD profile, which qualitatively agrees with our current observation-based understanding of FDs. We next adapt the models for energy dependence to enable comparison with different GCR measurement instruments (as they measure in different particle energy ranges). We test these modelling efforts against a set of multi-spacecraft observations of the same event, using the Forbush decrease model for the expanding flux rope (ForbMod). We find a reasonable agreement of the ForbMod model for the GCR depression in the CME magnetic structure with multi-spacecraft measurements, indicating that modelled FDs reflect well the CME evolution.

Forbush decreases (FDs) are depletions in the galactic cosmic ray (GCR) count rate that last typically for about a week and can be caused by coronal mass ejections (CMEs) or corotating interacting regions (CIRs). Fast CMEs that drive shocks cause large FDs that often show a two-step decrease where the first step is attributed to the shock/sheath region, while the second step is attributed to the closed magnetic structure. Since the difference in size of shock and sheath region is significant, and since there are observed effects that can be related to shocks and not necessarily to the sheath region we expect that the physical mechanisms governing the interaction with GCRs in these two regions are different. We therefore aim to analyse interaction of GCRs with heliospheric shocks only. We approximate the shock by a structure where the magnetic field linearly changes with position within this structure. We assume protons of different energy, different pitch angle and different incoming direction. We also vary the shock parameters such as the magnetic field strength and orientation, as well as the shock thickness. The results demonstrate that protons with higher energies are less likely to be reflected. Also, thicker shocks and shocks with stronger field reflect protons more efficiently.

Interplanetary evolution of eleven magnetic clouds (MCs) recorded by at least two radially aligned spacecraft is studied. The in situ magnetic field measurements are fitted to a cylindrically symmetric Gold-Hoyle force-free uniform-twist flux-rope configuration. The analysis reveals that in a statistical sense the expansion of studied MCs is compatible with self-similar behavior. However, individual events expose a large scatter of expansion rates, ranging from very weak to very strong expansion. Individually, only four events show an expansion rate compatible with the isotropic self-similar expansion. The results indicate that the expansion has to be much stronger when MCs are still close to the Sun than in the studied 0.47 - 4.8 AU distance range. The evolution of the magnetic field strength shows a large deviation from the behavior expected for the case of an isotropic self-similar expansion. In the statistical sense, as well as in most of the individual events, the inferred magnetic field decreases much slower than expected. Only three events show a behavior compatible with a self-similar expansion. There is also a discrepancy between the magnetic field decrease and the increase of the MC size, indicating that magnetic reconnection and geometrical deformations play a significant role in the MC evolution. About half of the events show a decay of the electric current as expected for the self-similar expansion. Statistically, the inferred axial magnetic flux is broadly consistent with it remaining constant. However, events characterized by large magnetic flux show a clear tendency of decreasing flux.

We present an analytical diffusion-expansion Forbush decrease (FD) model ForbMod which is based on the widely used approach of the initially empty, closed magnetic structure (i.e. flux rope) which fills up slowly with particles by perpendicular diffusion. The model is restricted to explain only the depression caused by the magnetic structure of the interplanetary coronal mass ejection (ICME). We use remote CME observations and a 3D reconstruction method (the Graduated Cylindrical Shell method) to constrain initial boundary conditions of the FD model and take into account CME evolutionary properties by incorporating flux rope expansion. Several flux rope expansion modes are considered, which can lead to different FD characteristics. In general, the model is qualitatively in agreement with observations, whereas quantitative agreement depends on the diffusion coefficient and the expansion properties (interplay of the diffusion and the expansion). A case study was performed to explain the FD observed 2014 May 30. The observed FD was fitted quite well by ForbMod for all expansion modes using only the diffusion coefficient as a free parameter, where the diffusion parameter was found to correspond to expected range of values. Our study shows that in general the model is able to explain the global properties of FD caused by FR and can thus be used to help understand the underlying physics in case studies.

We study the properties of 115 coronal holes in the time-range from 2010/08 to 2017/03, the peak velocities of the corresponding high-speed streams as measured in the ecliptic at 1 AU, and the corresponding changes of the Kp index as marker of their geo-effectiveness. We find that the peak velocities of high-speed streams depend strongly on both the ar- eas and the co-latitudes of their solar source coronal holes with regard to the heliospheric latitude of the satellites. Therefore, the co-latitude of their source coronal hole is an im- portant parameter for the prediction of the high-speed stream properties near the Earth. We derive the largest solar wind peak velocities normalized to the coronal hole areas for coronal holes located near the solar equator, and that they linearly decrease with increas- ing latitudes of the coronal holes. For coronal holes located at latitudes & 60{\deg}, they turn statistically to zero, indicating that the associated high-speed streams have a high chance to miss the Earth. Similar, the Kp index per coronal hole area is highest for the coronal holes located near the solar equator and strongly decreases with increasing latitudes of the coronal holes. We interpret these results as an effect of the three-dimensional propaga- tion of high-speed streams in the heliosphere, i.e., high-speed streams arising from coro- nal holes near the solar equator propagate in direction towards and directly hit the Earth, whereas solar wind streams arising from coronal holes at higher solar latitudes only graze or even miss the Earth.

The drag-based model (DBM) for heliospheric propagation of coronal mass ejections (CMEs) is a widely used analytical model which can predict CME arrival time and speed at a given heliospheric location. It is based on the assumption that the propagation of CMEs in interplanetary space is solely under the influence of magnetohydrodynamical drag, where CME propagation is determined based on CME initial properties as well as the properties of the ambient solar wind. We present an upgraded version, covering ensemble modelling to produce a distribution of possible ICME arrival times and speeds, the drag-based ensemble model (DBEM). Multiple runs using uncertainty ranges for the input values can be performed in almost real-time, within a few minutes. This allows us to define the most likely ICME arrival times and speeds, quantify prediction uncertainties and determine forecast confidence. The performance of the DBEM is evaluated and compared to that of ensemble WSA-ENLIL+Cone model (ENLIL) using the same sample of events. It is found that the mean error is $ME=-9.7$ hours, mean absolute error $MAE=14.3$ hours and root mean square error $RMSE=16.7$ hours, which is somewhat higher than, but comparable to ENLIL errors ($ME=-6.1$ hours, $MAE=12.8$ hours and $RMSE=14.4$ hours). Overall, DBEM and ENLIL show a similar performance. Furthermore, we find that in both models fast CMEs are predicted to arrive earlier than observed, most probably owing to the physical limitations of models, but possibly also related to an overestimation of the CME initial speed for fast CMEs.

The propagation of 15 interplanetary coronal mass ejections (ICMEs) from Earth's orbit (1 AU) to Mars (~ 1.5 AU) has been studied with their propagation speed estimated from both measurements and simulations. The enhancement of magnetic fields related to ICMEs and their shock fronts cause the so-called Forbush decrease, which can be de- tected as a reduction of galactic cosmic rays measured on-ground. We have used galactic cosmic ray (GCR) data from in-situ measurements at Earth, from both STEREO A and B as well as GCR measurements by the Radiation Assessment Detector (RAD) instrument onboard Mars Science Laboratory (MSL) on the surface of Mars. A set of ICME events has been selected during the periods when Earth (or STEREO A or B) and Mars locations were nearly aligned on the same side of the Sun in the ecliptic plane (so-called opposition phase). Such lineups allow us to estimate the ICMEs' transit times between 1 and 1.5 AU by estimating the delay time of the corresponding Forbush decreases measured at each location. We investigate the evolution of their propagation speeds before and after passing Earth's orbit and find that the deceleration of ICMEs due to their interaction with the ambient solar wind may continue beyond 1 AU. We also find a substantial variance of the speed evolution among different events revealing the dynamic and diverse nature of eruptive solar events. Furthermore, the results are compared to simulation data obtained from two CME propagation models, namely the Drag-Based Model and ENLIL plus cone model.

We investigate the statistics of 288 low-latitude coronal holes extracted from SDO/AIA-193 filtergrams over the time range 2011/01/01 to 2013/12/31. We analyse the distribution of characteristic coronal hole properties, such as the areas, mean AIA-193 intensities, and mean magnetic field densities, the local distribution of the SDO/AIA-193 intensity and the magnetic field within the coronal holes, and the distribution of magnetic flux tubes in coronal holes. We find that the mean magnetic field density of all coronal holes under study is 3.0 +- 1.6 G, and the percentage of unbalanced magnetic flux is 49 +- 16 %. The mean magnetic field density, the mean unsigned magnetic field density, and the percentage of unbalanced magnetic flux of coronal holes depend strongly pairwise on each other, with correlation coefficients cc > 0.92. Furthermore, we find that the unbalanced magnetic flux of the coronal holes is predominantly concentrated in magnetic flux tubes: 38 % (81 %) of the unbalanced magnetic flux of coronal holes arises from only 1 % (10 %) of the coronal hole area, clustered in magnetic flux tubes with field strengths > 50 G (10 G). The average magnetic field density and the unbalanced magnetic flux derived from the magnetic flux tubes correlate with the mean magnetic field density and the unbalanced magnetic flux of the overall coronal hole (cc > 0.93). These findings give evidence that the overall magnetic characteristics of coronal holes are governed by the characteristics of the magnetic flux tubes.

Predicting transit times of Coronal Mass Ejections (CMEs) from their initial parameters is a very important subject, not only from the scientific perspective, but also because CMEs represent a hazard for human technology. We used a neural network to analyse transit times for 153 events with only two input parameters: initial velocity of the CME, $v$, and Central Meridian Distance, CMD, of its associated flare. We found that transit time dependence on $v$ is showing a typical drag-like pattern in the solar wind. The results show that the speed at which acceleration by drag changes to deceleration is $v\approx$500 km s$^{-1}$. Transit times are also found to be shorter for CMEs associated with flares on the western hemisphere than those originating on the eastern side of the Sun. We attribute this difference to the eastward deflection of CMEs on their path to 1 AU. The average error of the NN prediction in comparison to observations is $\approx$12 hours which is comparable to other studies on the same subject.

The so-called drag-based model (DBM) simulates analytically the propagation of coronal mass ejections (CMEs) in interplanetary space and allows the prediction of their arrival times and impact speeds at any point in the heliosphere ("target"). The DBM is based on the assumption that beyond a distance of about 20 solar radii from the Sun, the dominant force acting on CMEs is the "aerodynamic" drag force. In the standard form of DBM, the user provisionally chooses values for the model input parameters, by which the kinematics of the CME over the entire Sun--"target" distance range is defined. The choice of model input parameters is usually based on several previously undertaken statistical studies. In other words, the model is used by ad hoc implementation of statistics-based values of the input parameters, which are not necessarily appropriate for the CME under study. Furthermore, such a procedure lacks quantitative information on how well the simulation reproduces the coronagraphically observed kinematics of the CME, and thus does not provide an estimate of the reliability of the arrival prediction. In this paper we advance the DBM by adopting it in a form that employs the CME observations over a given distance range to evaluate the most suitable model input parameters for a given CME by means of the least-squares fitting. Furthermore, the new version of the model automatically responds to any significant change of the conditions in the ambient medium (solar wind speed, density, CME--CME interactions, etc.) by changing the model input parameters according to changes in the CME kinematics. The advanced DBM is shaped in a form that can be readily employed in an operational system for real-time space-weather forecasting by promptly adjusting to a successively expanding observational dataset, thus providing a successively improving prediction of the CME arrival.

Forecasting the in situ properties of coronal mass ejections (CMEs) from remote images is expected to strongly enhance predictions of space weather, and is of general interest for studying the interaction of CMEs with planetary environments. We study the feasibility of using a single heliospheric imager (HI) instrument, imaging the solar wind density from the Sun to 1 AU, for connecting remote images to in situ observations of CMEs. We compare the predictions of speed and arrival time for 22 CMEs (in 2008-2012) to the corresponding interplanetary coronal mass ejection (ICME) parameters at in situ observatories (STEREO PLASTIC/IMPACT, Wind SWE/MFI). The list consists of front- and backsided, slow and fast CMEs (up to $2700 \: km \: s^{-1}$). We track the CMEs to $34.9 \pm 7.1$ degrees elongation from the Sun with J-maps constructed using the SATPLOT tool, resulting in prediction lead times of $-26.4 \pm 15.3$ hours. The geometrical models we use assume different CME front shapes (Fixed-$\Phi$, Harmonic Mean, Self-Similar Expansion), and constant CME speed and direction. We find no significant superiority in the predictive capability of any of the three methods. The absolute difference between predicted and observed ICME arrival times is $8.1 \pm 6.3$ hours ($rms$ value of 10.9h). Speeds are consistent to within $284 \pm 288 \: km \: s^{-1}$. Empirical corrections to the predictions enhance their performance for the arrival times to $6.1 \pm 5.0$ hours ($rms$ value of 7.9h), and for the speeds to $53 \pm 50 \: km \: s^{-1}$. These results are important for Solar Orbiter and a space weather mission positioned away from the Sun-Earth line.

We study a filament eruption, two-ribbon flare, and coronal mass ejection (CME) that occurred in Active Region NOAA 10898 on 6 July 2006. The filament was located South of a strong sunspot that dominated the region. In the evolution leading up to the eruption, and for some time after it, a counter-clockwise rotation of the sunspot of about 30 degrees was observed. We suggest that the rotation triggered the eruption by progressively expanding the magnetic field above the filament. To test this scenario, we study the effect of twisting the initially potential field overlying a pre-existing flux-rope, using three-dimensional zero-$\beta$ MHD simulations. We first consider a relatively simple and symmetric system, and then study a more complex and asymmetric magnetic configuration, whose photospheric flux distribution and coronal structure are guided by the observations and a potential field extrapolation. In both cases, we find that the twisting leads to the expansion of the overlying field. As a consequence of the progressively reduced magnetic tension, the flux-rope quasi-statically adapts to the changed environmental field, rising slowly. Once the tension is sufficiently reduced, a distinct second phase of evolution occurs where the flux-rope enters an unstable regime characterized by a strong acceleration. Our simulations thus suggest a new mechanism for the triggering of eruptions in the vicinity of rotating sunspots.

The Type-II solar radio burst recorded on 13 June 2010 by the radio spectrograph of the Hiraiso Solar Observatory was employed to estimate the magnetic-field strength in the solar corona. The burst was characterized by a well pronounced band-splitting, which we used to estimate the density jump at the shock and Alfven Mach number using the Rankine-Hugoniot relations. The plasma frequency of the Type-II bursts is converted into height [R] in solar radii using the appropriate density model, then we estimated the shock speed [Vs], coronal Alfven velocity [Va], and the magnetic-field strength at different heights. The relative bandwidth of the band-split is found to be in the range 0.2 -- 0.25, corresponding to the density jump of X = 1.44 -- 1.56, and the Alfven Mach number of MA = 1.35 -- 1.45. The inferred mean shock speed was on the order of V ~ 667 km/s. From the dependencies V(R) and MA(R) we found that Alfven speed slightly decreases at R ~ 1.3 -- 1.5. The magnetic-field strength decreases from a value between 2.7 and 1.7 G at R ~ 1.3 -- 1.5 Rs depending on the coronal-density model employed. We find that our results are in good agreement with the empirical scaling by Dulk and McLean (Solar Phys. 57, 279, 1978) and Gopalswamy et al. (Astrophys. J. 744, 72, 2012). Our result shows that Type-II band splitting method is an important tool for inferring the coronal magnetic field, especially when independent measurements were made from white light observations.

Numerical simulations of magnetosonic wave formation driven by an expanding cylindrical piston are performed to get better physical insight into the initiation and evolution of large-scale coronal waves. Several very basic initial configurations are employed to analyze intrinsic characteristics of the MHD wave formation that do not depend on specific properties of the environment. It turns out that these simple initial configurations result in piston/wave morphologies and kinematics that reproduce common characteristics of coronal waves. In the initial stage the wave and the expanding source-region cannot be clearly resolved. During the acceleration stage of the source-region inflation, the wave is driven by the piston expansion, so its amplitude and phase-speed increase, whereas the wavefront profile steepens. At a given point, a discontinuity forms in the wavefront profile. The time/distance required for the shock formation is shorter for a more impulsive source-region expansion. After the piston stops, the wave amplitude and phase-speed start decreasing. During the expansion, most of the source region becomes strongly rarified, which reproduces the coronal dimming left behind the eruption. On the other hand, the density increases at the source-region boundary, and stays enhanced even after the expansion stops, which might explain stationary brightenings that are sometimes observed at the edges of the erupted coronal structure. In addition, in the rear of the wave a weak density depletion develops, trailing the wave, which is sometimes observed as weak transient coronal dimming. Finally, we find a well defined relationship between the impulsiveness of the source-region expansion and the wave amplitude and phase speed. The results for the cylindrical piston are also compared with the outcome for a planar wave, to find out how different geometries affect the evolution of the wave.

We investigate the interaction of three consecutive large-scale coronal waves with a polar coronal hole, simultaneously observed on-disk by the Solar TErrestrial Relations Observatory (STEREO)-A spacecraft and on the limb by the PRoject for On-Board Autonomy 2 (PROBA2) spacecraft on January 27, 2011. All three extreme-ultraviolet(EUV) waves originate from the same active region NOAA 11149 positioned at N30E15 in the STEREO-A field-of-view and on the limb in PROBA2. We derive for the three primary EUV waves start velocities in the range of ~310 km/s for the weakest up to ~500 km/s for the strongest event. Each large-scale wave is reflected at the border of the extended coronal hole at the southern polar region. The average velocities of the reflected waves are found to be smaller than the mean velocities of their associated direct waves. However, the kinematical study also reveals that in each case the end velocity of the primary wave matches the initial velocity of the reflected wave. In all three events the primary and reflected waves obey the Huygens-Fresnel principle, as the incident angle with ~10{\deg} to the normal is of the same size as the angle of reflection. The correlation between the speed and the strength of the primary EUV waves, the homologous appearance of both the primary and the reflected waves, and in particular the EUV wave reflections themselves implicate that the observed EUV transients are indeed nonlinear large-amplitude MHD waves.

We use high time cadence images acquired by the STEREO EUVI and COR instruments to study the evolution of coronal mass ejections (CMEs), from their initiation, through the impulsive acceleration to the propagation phase. For a set of 95 CMEs we derived detailed height, velocity and acceleration profiles and statistically analysed characteristic CME parameters: peak acceleration, peak velocity, acceleration duration, initiation height, height at peak velocity, height at peak acceleration and size of the CME source region. The CME peak accelerations derived range from 20 to 6800 m s^2 and are inversely correlated to the acceleration duration and to the height at peak acceleration. 74% of the events reach their peak acceleration at heights below 0.5 Rsun. CMEs which originate from compact sources low in the corona are more impulsive and reach higher peak accelerations at smaller heights. These findings can be explained by the Lorentz force, which drives the CME accelerations and decreases with height and CME size.

A study and classification of super-short structures (SSSs) recorded during metric type IV bursts is presented. The most important property of SSSs is their duration, at half power ranging from 4-50 ms, what is up to 10 times shorter than spikes at corresponding frequencies. The solar origin of the SSSs is confirmed by one-to-one correspondence between spectral recordings of Artemis-IV1 and high time resolution single frequency measurements of the TSRS2. We have divided the SSSs in the following categories: 1. Broad-Band SSSs: They were partitioned in two subcategories, the SSS-Pulses and Drifting SSSs; 2. Narrow-band: They appear either as Spike-Like SSSs or as Patch-Like SSSs; 3. Complex SSS: They consist of the absorption-emission segments and were morphologically subdivided into Rain-drop Bursts (narrow-band emission head and a broad-band absorption tail) and Blinkers.

Eruption of a coronal mass ejection (CME) drags and "opens" the coronal magnetic field, presumably leading to the formation of a large-scale current sheet and the field relaxation by magnetic reconnection. We analyze physical characteristics of ray-like coronal features formed in the aftermath of CMEs, to check if the interpretation of this phenomenon in terms of reconnecting current sheet is consistent with the observations. The study is focused on measurements of the ray width, density excess, and coronal velocity field as a function of the radial distance. The morphology of rays indicates that they occur as a consequence of Petschek-like reconnection in the large scale current sheet formed in the wake of CME. The hypothesis is supported by the flow pattern, often showing outflows along the ray, and sometimes also inflows into the ray. The inferred inflow velocities range from 3 to 30 km s$^{-1}$, consistent with the narrow opening-angle of rays, adding up to a few degrees. The density of rays is an order of magnitude larger than in the ambient corona. The density-excess measurements are compared with the results of the analytical model in which the Petschek-like reconnection geometry is applied to the vertical current sheet, taking into account the decrease of the external coronal density and magnetic field with height. The model results are consistent with the observations, revealing that the main cause of the density excess in rays is a transport of the dense plasma from lower to larger heights by the reconnection outflow.