The relationship between the peak velocities of high-speed solar wind streams near Earth and the areas of their solar source regions, i.e., coronal holes, has been known since the 1970s, but it is still physically not well understood. We perform 3D magnetohydrodynamic (MHD) simulations using the European Heliospheric Forecasting Information Asset (EUHFORIA) code to show that this empirical relationship forms during the propagation phase of high-speed streams from the Sun to Earth. For this purpose, we neglect the acceleration phase of high-speed streams, and project the areas of coronal holes to a sphere at 0.1 au. We then vary only the areas and latitudes of the coronal holes. The velocity, temperature, and density in the cross section of the corresponding highspeed streams at 0.1 au are set to constant, homogeneous values. Finally, we propagate the associated high-speed streams through the inner heliosphere using the EUHFORIA code. The simulated high-speed stream peak velocities at Earth reveal a linear dependence on the area of their source coronal holes. The slopes of the relationship decrease with increasing latitudes of the coronal holes, and the peak velocities saturate at a value of about 730 km/s, similar to the observations. These findings imply that the empirical relationship between the coronal hole areas and high-speed stream peak velocities does not describe the acceleration phase of high-speed streams, but is a result of the high-speed stream propagation from the Sun to Earth.

In collision-poor plasmas from space, three distinct ion-ion instabilities can be driven by the proton beams streaming along the background magnetic field: left-hand resonant, non-resonant, and right-hand resonant instabilities. These instabilities are in general investigated considering only idealized proton beams with Maxwellian velocity distributions, and ignoring the implications of suprathermal populations, usually reproduced by the Kappa power-laws. Moreover, the existing theories minimize the kinetic effects of electrons, assuming them isotropic and Maxwellian distributed. In an attempt to overcome these limitations, in the present paper we present the results of an extended investigation of ion-ion instabilities, which show that their dispersion and stability properties (e.g. growth rates, wave frequencies, and the unstable wave numbers) are highly sensitive to the influence of suprathermal populations and anisotropic electrons. These results offer valuable explanations for the origin of the enhanced low-frequency fluctuations, frequently observed in space plasmas and associated with proton beams.

Relative drifts between different species or particle populations are characteristic to solar plasma outflows, e.g., in the fast streams of the solar winds, coronal mass ejections and interplanetary shocks. This paper characterizes the dispersion and stability of the low-beta alpha/proton drifts in the absence of any intrinsic thermal anisotropies, which are usually invoked in order to stimulate various instabilities. The dispersion relations derived here describe the full spectrum of instabilities and their variations with the angle of propagation and plasma parameters. The results unveil a potential competition between instabilities of the electromagnetic proton cyclotron and alpha cyclotron modes. For conditions specific to a low-beta solar wind, e.g., at low heliocentric distances in the outer corona, the instability operates on the alpha cyclotron branch. The growth rates of the alpha cyclotron mode are systematically stimulated by the (parallel) plasma beta and/or the alpha-proton temperature ratio. One can therefore expect that this instability develops even in the absence of temperature anisotropies, with potential to contribute to a self-consistent regulation of the observed drift of alpha particles.

Context: Eruptive events such as coronal mass ejections (CMEs) and flares accelerate particles and generate shock waves which can arrive at Earth and can disturb the magnetosphere. Understanding the association between CMEs and CME-driven shocks is therefore highly important for space weather studies. Aims: We present a study of the CME/flare event associated with two type II bursts observed on September 27, 2012. The aim of the study is to understand the relationship between the observed CME and the two distinct shock wave signatures. Methods: The multi-wavelength study of the eruptive event (CME/flare) was complemented with radio triangulation of the associated radio emission and modelling of the CME and the shock wave employing MHD simulations. Results: We found that, although temporal association between the type II bursts and the CME is good, the low-frequency type II(LF-type II) burst occurs significantly higher in the corona than the CME and its relationship to the CME is not straightforward. The analysis of the EIT wave (coronal bright front) shows the fastest wave component to be in the southeast quadrant of the Sun. This is also the quadrant in which the source positions of the LF-type II were found to be located, probably resulting from the interaction between the shock wave and a streamer. Conclusions: The relationship between the CME/flare event and the shock wave signatures is discussed using the temporal association,as well as the spatial information of the radio emission. Further, we discuss the importance and possible effects of the frequently non-radial propagation of the shock wave.

The heat transport in the solar wind mainly involves the electron strahl population beaming along the local magnetic field. The rapid non-adiabatic decrease of the heat flux with the collisionless expansion of the solar wind is believed to be self-consistently controlled by the electron beam-plasma (or heat-flux) instabilities. However, multiple heat-flux instabilities (HFIs) of different nature are proposed in the literature, always claiming for relevance in the solar wind conditions and puzzling over their role in the self-regulation of electron strahl. Present paper describes by comparison the full spectrum of electromagnetic and electrostatic heat flux instabilities, as prescribed by the kinetic theory for high beta conditions ($\beta_e \gg 0.1$) and different beaming velocities (or drifts) of the strahl ($U_s$). The parametric study performed here reveals the existence of alternative (complementary) regimes of HFIs, and provides a detailed characterization of their dominance and interplay. For instance, for $\beta_e =2$ and drifts lower than the thermal speed of the strahl ($U_s / \alpha_s < 1$) the most probable to develop is the (parallel) whistler HFI (WHFI), while more energetic beams are susceptible to the oblique WHFI (for $U_s/\alpha_s \gtrsim 1$), or to the electrostatic instabilities of the electron acoustic and electron beam modes (for $U_s/\alpha_s > \sqrt{2}$). These results demonstrate that only a realistic parameterization combined with a selective spectral analysis may offer plausible explanations for the nature and origin of wave fluctuations reported by the observations in association with different types of electron strahls, e.g., in the slow or fast winds, streaming interaction regions and interplanetary shocks.

In this article, we present the application of the weighted horizontal gradient of magnetic field (WGM) flare prediction method to 3-dimensional (3D) extrapolated magnetic configurations of 13 flaring solar active regions (ARs). The main aim is to identify an optimal height range, if any, in the interface region between the photosphere and lower corona, where the flare onset time prediction capability of WGM is best exploited. The optimal height is where flare prediction, by means of the WGM method, is achieved earlier than at the photospheric level. 3D magnetic structures, based on potential and non-linear force-free field extrapolations, are constructed to study a vertical range from the photosphere up to the low corona with a 45 km step size. The WGM method is applied as a function of height to all 13 flaring AR cases that are subject to certain selection criteria. We found that applying the WGM method between 1000 and 1800 km above the solar surface would improve the prediction of the flare onset time by around 2-8 hours.Certain caveats and an outlook for future work along these lines are also discussed.

A new class of one-dimensional solar wind models is developed within the general polytropic, single-fluid hydrodynamic framework. The particular case of quasi-adiabatic radial expansion with a localized heating source is considered. We consider analytical solutions with continuous Mach number over the entire radial domain while allowing for jumps in the flow velocity, density, and temperature, provided that there exists an external source of energy in the vicinity of the critical point which supports such jumps in physical quantities. This is substantially distinct from both the standard Parker solar wind model and the original nozzle solutions, where such discontinuous solutions are not permissible. We obtain novel sample analytic solutions of the governing equations corresponding to both slow and fast wind.

We study the effects of heat flows and velocity shear on the parallel firehose instability in weakly collisional plasma flow. For this purpose we apply an anisotropic 16-moments MHD fluid closure model that takes into account the pressure and temperature anisotropy, as well as the effect of anisotropic heat flux. The linear stability analysis of the firehose modes is carried out in the incompressible limit, where the MHD flow is parallel to the background magnetic field, while the velocity is sheared in the direction transverse to the flow direction. It seems that an increase of the velocity shear parameter leads to higher growth rates of the firehose instability. The increase of the instability growth rate is most profound for perturbations with oblique wave-numbers $k_{\perp}/k_{\parallel} < 1$. The heat flux parameter introduces an asymmetry of the instability growth in the shear plane: perturbations with wave-vectors with a component in the direction of the velocity shear grow significantly stronger as compared to those with components in the opposite direction. We discuss the implications of the presented study on the observable features of the solar wind and possible measurements of local parameters of the solar wind based on the stability constraints set by the firehose instability.

Methods: Stealth CMEs represent a particular class of solar eruptions that are clearly distinguished in coronagraph observations, but they don't have a clear source signature. A particular type of stealth CMEs occurs in the trailing current sheet of a previous ejection, therefore, we used the 2.5D MHD package of the code MPI-AMRVAC to numerically simulate consecutive CMEs by imposing shearing motions onto the inner boundary. 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. We then compared the obtained eruptions with the observed directions of propagation 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. For the first time, a stealth CME is simulated in the aftermath of a first eruption, through changes in the applied shearing speed.

It is generally accepted that extreme space weather events tend to be related to strong flares and fast halo coronal mass ejections CMEs. In the present paper, we carefully identify the chain of events from the Sun to the Earth induced by all 12 X-class flares that occurred in 2002. In this small sample, we find an unusual high rate (58\%) of solar sources with a longitude larger than 74 degrees. Yet, all 12 X-class flares are associated with at least one CME. The fast halo CMEs (50\% ) are related to interplanetary CMEs (ICMEs) at L1 and weak Dst minimum values ($> -51\;$nT); while 5 (41\%) of the 12 X-class flares are related to solar proton events (SPE). We conclude that: (i) All twelve analyzed solar events, even those associated with fast halo CMEs originating from the central disk region, and those ICMEs and SPEs were not very geo-effective. This unexpected result demonstrates that the suggested events in the chain (fast halo CME, X-class flares, central disk region, ICME, SPE) are not infallible proxies for geo-effectiveness. (ii) The low value of integrated and normalized southward component of the IMF ($B^*_z$) may explain the low geo-effectiveness for this small sample. In fact, $B^*_z$ is well correlated to the weak Dst and low auroral electrojet (AE) activity. Hence, the only space weather impact at Earth in 2002 we can explain is based on $B^*_z$ at L1.

Coronal mass ejections (CMEs) are the primary sources of intense disturbances at Earth, where their geo-effectiveness is largely determined by their dynamic pressure and internal magnetic field, which can be significantly altered during interactions with other CMEs in interplanetary space. We analyse three successive CMEs that erupted from the Sun during September 4-6, 2017, investigating the role of CME-CME interactions as source of the associated intense geomagnetic storm (Dst_min=-142 nT on September 7). To quantify the impact of interactions on the (geo-)effectiveness of individual CMEs, we perform global heliospheric simulations with the EUHFORIA model, using observation-based initial parameters with the additional purpose of validating the predictive capabilities of the model for complex CME events. The simulations show that around 0.45 AU, the shock driven by the September 6 CME started compressing a preceding magnetic ejecta formed by the merging of two CMEs launched on September 4, significantly amplifying its Bz until a maximum factor of 2.8 around 0.9 AU. The following gradual conversion of magnetic energy into kinetic and thermal components reduced the Bz amplification until its almost complete disappearance around 1.8 AU. We conclude that a key factor at the origin of the intense storm triggered by the September 4-6, 2017 CMEs was their arrival at Earth during the phase of maximum Bz amplification. Our analysis highlights how the amplification of the magnetic field of individual CMEs in space-time due to interaction processes can be characterised by a growth, a maximum, and a decay phase, suggesting that the time interval between the CME eruptions and their relative speeds are critical factors in determining the resulting impact of complex CMEs at various heliocentric distances (helio-effectiveness).

We investigate the evolution of both normal and inverse CMEs ejected at different initial velocities, and observe the effect of the background wind density and their magnetic polarity on their evolution up to 1 AU. We performed 2.5D simulations by solving the magnetohydrodynamic equations on a radially stretched grid, employing a block-based adaptive mesh refinement scheme based on a density threshold to achieve high resolution following the evolution of the magnetic clouds and the leading bow shocks. All the simulations discussed in the present paper were performed using the same initial grid and numerical methods. The polarity of the internal magnetic field of the CME has a substantial effect on its propagation velocity and on its deformation and erosion during its evolution towards Earth. We quantified the effects of the polarity of the internal magnetic field of the CMEs and of the density of the background solar wind on the arrival times of the shock front and the magnetic cloud. We determined the positions and propagation velocities of the magnetic clouds and thus also the stand-off distance of the leading shock fronts (i.e. the thickness of the magnetic sheath region) and the deformation and erosion of the magnetic clouds during their evolution from the Sun to the Earth. Inverse CMEs were found to be faster than normal CMEs ejected in the same initial conditions, but with smaller stand-off distances. They also have a higher magnetic cloud length, opening angle, and mass. Synthetic satellite time series showed that the shock magnitude is not affected by the polarity of the CME. However, the density peak of the magnetic cloud is dependent on the polarity and, in case of inverse CMEs, also on the background wind density. The magnitude of the z-component of the magnetic field was not influenced by either the polarity or the wind density.

The evolution of the magnetic field and plasma quantities inside a coronal mass ejection (CME) with distance are known from statistical studies using data from 1 au monitors, planetary missions, Helios, and Ulysses. This does not cover the innermost heliosphere, below 0.29 au, where no data are yet publicly available. Here, we describe the evolution of the properties of simulated CMEs in the inner heliosphere using two different initiation mechanisms. We compare the radial evolution of these properties with that found from statistical studies based on observations in the inner heliosphere by Helios and MESSENGER. We find that the evolution of the radial size and magnetic field strength is nearly indistinguishable for twisted flux rope from that of writhed CMEs. The evolution of these properties is also consistent with past studies, primarily with recent statistical studies using in situ measurements and with studies using remote observations of CMEs.

In the absence of efficient collisions, deviations from thermal equilibrium of plasma particle distributions are controlled by the self-generated instabilities. The whistler instability is a notorious example, usually responsible for the regulation of electron temperature anisotropy $A = T_{\perp}/T_\parallel>$ (with $\perp, \parallel$ respective to the magnetic field direction) observed in space plasmas, e.g., solar wind and planetary magnetospheres. Suprathermal electrons present in these environments change the plasma dispersion and stability properties, with expected consequences on the kinetic instabilities and the resulting fluctuations, which, in turn, scatter the electrons and reduce their anisotropy. In order to capture these mutual effects we use a quasilinear kinetic approach and PIC simulations, which provide a comprehensive characterization of the whistler instability under the influence of suprathermal electrons. Analysis is performed for a large variety of plasma conditions, ranging from low-beta plasmas encountered in outer corona or planetary magnetospheres to a high-beta solar wind characteristic to large heliospheric distances. Enhanced by the suprathermal electrons, whistler fluctuations stimulate the relaxation of temperature anisotropy, and this influence of suprathermals increases with plasma beta parameter.

In collision-poor plasmas from space, e.g., solar wind or stellar outflows, the heat-flux carried by the strahl or beaming electrons is expected to be regulated by the self-generated instabilities. Recently, simultaneous field and particle observations have indeed revealed enhanced whistler-like fluctuations in the presence of counter-beaming populations of electrons, connecting these fluctuations to the whistler heat-flux instability (WHFI). This instability is predicted only for limited conditions of electron beam-plasmas, and was not captured in numerical simulations yet. In this letter we report the first simulations of WHFI in particle-in-cell (PIC) setups, realistic for the solar wind conditions, and without temperature gradients or anisotropies to trigger the instability in the initiation phase. The velocity distributions have a complex reaction to the enhanced whistler fluctuations conditioning the instability saturation by a decrease of the relative drifts combined with induced (effective) temperature anisotropies (heating the core electrons and pitch-angle and energy scattering the strahl). These results are in good agreement with a recent quasilinear approach, and support therefore a largely accepted belief that WHFI saturates at moderate amplitudes. In anti-sunward direction the strahl becomes skewed with a pitch-angle distribution decreasing in width as electron energy increases, that seems to be characteristic to self-generated whistlers and not to small-scale turbulence.

In order to address the growing need for more accurate space weather predictions, a new model named EUHFORIA (EUropean Heliospheric FORecasting Information Asset) was recently developed (Pomoell and Poedts, 2018). We present first results of the performance assessment for the solar wind modeling with EUHFORIA and identify possible limitations of its present setup. Using the basic EUHFORIA 1.0.4. model setup with the default input parameters, we modeled background solar wind (no coronal mass ejections) and compared the obtained results with ACE, in situ measurements. For the need of statistical study we developed a technique of combining daily EUHFORIA runs into continuous time series. The combined time series were derived for the years 2008 (low solar activity) and 2012 (high solar activity) from which in situ speed and density profiles were extracted. We find for the low activity phase a better match between model results and observations compared to the considered high activity time interval. The quality of the modeled solar wind parameters is found to be rather variable. Therefore, to better understand the obtained results we also qualitatively inspected characteristics of coronal holes, sources of the studied fast streams. We discuss how different characteristics of the coronal holes and input parameters to EUHFORIA influence the modeled fast solar wind, and suggest possibilities for the improvements of the model.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

The adopted 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 \textsc{catch} (Collection of Analysis Tools for Coronal Holes) to extract CH areas and compare them to the open flux areas modelled 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_{\rm ss}$, of the outer boundary of the PFSS, known as the source surface, and the height, $R_{\rm i}$, of the inner boundary of the SCS are important parameters affecting the modelled CH areas. We investigate the impact the two model parameters can have in the modelled results. We vary $R_{\rm ss}$ within the interval [1.4, 3.2]$R_{\rm \odot}$ with a step of 0.1$R_{\rm \odot}$, and $R_{\rm i}$ within the interval [1.3, 2.8]$R_{\rm \odot}$ with the same step, and the condition that $R_{\rm i}<R_{\rm ss}$. This way we have a set of 184 initial parameters to the model and we assess the model results for all these possible height pairs. We conclude that the default heights used so far fail in modelling accurately CH areas and lower heights need to be considered.

We analyse in this work the propagation and geoeffectiveness of four successive coronal mass ejections (CMEs) that erupted from the Sun during 21--23 May 2013 and that were detected in interplanetary space by the Wind and/or STEREO-A spacecraft. All these CMEs featured critical aspects for understanding so-called "problem space weather storms" at Earth. In the first three events a limb CMEs resulted in moderately geoeffective in-situ structures at their target location in terms of the disturbance storm time (Dst) index (either measured or estimated). The fourth CME, which also caused a moderate geomagnetic response, erupted from close to the disc centre as seen from Earth, but it was not visible in coronagraph images from the spacecraft along the Sun--Earth line and appeared narrow and faint from off-angle viewpoints. Making the correct connection between CMEs at the Sun and their in-situ counterparts is often difficult for problem storms. We investigate these four CMEs using multiwavelength and multipoint remote-sensing observations (extreme ultraviolet, white light, and radio), aided by 3D heliospheric modelling, in order to follow their propagation in the corona and in interplanetary space and to assess their impact at 1 AU. Finally, we emphasise the difficulties in forecasting moderate space weather effects provoked by problematic and ambiguous events and the importance of multispacecraft data for observing and modelling problem storms.

Context. Solar outflows are a considerable source of free energy which accumulates in multiple forms like beaming (or drifting) components and/or temperature anisotropies. However, kinetic anisotropies of plasma particles do not grow indefinitely and particle-particle collisions are not efficient enough to explain the observed limits of these anisotropies. Instead, the self-generated wave instabilities can efficiently act to constrain kinetic anisotropies, but the existing approaches are simplified and do not provide satisfactory explanations. Thus, small deviations from isotropy shown by the electron temperature ($T$) in fast solar winds are not explained yet. Aims. This paper provides an advanced quasilinear description of the whistler instability driven by the anisotropic electrons in conditions typical for the fast solar winds. The enhanced whistler-like fluctuations may constrain the upper limits of temperature anisotropy $A \equiv T_\perp /T_\parallel > 1$, where $\perp, \parallel$ are defined with respect to the magnetic field direction. Methods. Studied are the self-generated whistler instabilities, cumulatively driven by the temperature anisotropy and the relative (counter)drift of the electron populations, e.g., core and halo electrons. Recent studies have shown that quasi-stable states are not bounded by the linear instability thresholds but an extended quasilinear approach is necessary to describe them in this case. Results. Marginal conditions of stability are obtained from a quasilinear theory of the cumulative whistler instability, and approach the quasi-stable states of electron populations reported by the observations.The instability saturation is determined by the relaxation of both the temperature anisotropy and the relative drift of electron populations.

Coronal Mass Ejections (CMEs) are the primary source of strong space weather disturbances at Earth. Their geoeffectiveness 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. We study two Earth-directed CMEs using the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) model, testing the predictive capabilities of a linear force-free spheromak CME model initialised using parameters derived from remote-sensing observations. Using observation-based CME input parameters, we perform MHD simulations of the events with the cone and spheromak CME models. Results show that spheromak CMEs propagate faster than cone CMEs when initialised with the same kinematic parameters. We interpret these differences as due to different Lorentz forces acting within cone and spheromak CMEs, leading to different CME expansions. Such discrepancies can be mitigated by initialising spheromak CMEs with a reduced speed corresponding to the radial speed only. Results at Earth evidence that the spheromak model improves the predictions of B(Bz) up to 12-60(22-40) percentage points compared to a cone model. Considering virtual spacecraft located around Earth, B(Bz) predictions reach 45-70%(58-78%) of the observed peak values. The spheromak model predicts inaccurate magnetic field parameters at Earth for CMEs propagating away from the Sun-Earth line, while it successfully predicts the CME properties and arrival time in the case of strictly Earth-directed events. 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 around it.

As is well known, CME propagation often results in the fragmentation of the solar atmosphere on smaller regions of density (magnetic field) enhancement (depletion). It is expected that this type of fragmentation may have radio signatures. The general aim of the present paper is to perform a comparative analysis of type III solar and narrow-band type-III-like radio burst properties before and during CME events, respectively. The main goal is to analyze radio observational signatures of the dynamical processes in solar corona. In particular, we aim to perform a comparison of local plasma parameters without and with CME propagation, based on the analysis of decameter radio emission data. We performed a comparison of usual type III bursts before the CME with narrow-band type-III-like bursts, which are observationally detectable on top of the background type IV radio bursts associated with CME propagation. We focused on the analysis of in total 429 type III and 129 narrow-band type-III-like bursts. We studied their main characteristic parameters such as frequency drift rate, duration, and instantaneous frequency bandwidth using standard statistical methods. Furthermore, we inferred local plasma parameters using known definitions of frequency drift, duration, and instantaneous frequency bandwidth.The analysis reveals that the physical parameters of coronal plasma before CMEs considerably differ from those during the propagation of CMEs (the observational periods 2 and 4 with type IV radio bursts associated with CMEs). Local density radial profiles and the characteristic spatial scales of radio emission sources vary with radial distance more drastically during the CME propagation compared to the cases of quasistatic solar atmosphere without CME(s) (observational periods 1 and 3). The results of the work enable us to distinguish different regimes of plasma state in the solar corona.

The hot beaming (or strahl) electrons responsible for the main electron heat-flux in the solar wind are believed to be self-regulated by the electromagnetic beaming instabilities, also known as the heat-flux instabilities. Here we report the first quasi-linear theoretical approach of the whistler unstable branch able to characterize the long-term saturation of the instability as well as the relaxation of the electron velocity distributions. The instability saturation is not solely determined by the drift velocities, which undergo only a minor relaxation, but mainly from a concurrent interaction of electrons with whistlers that induces (opposite) temperature anisotropies of the core and beam populations and reduces the effective anisotropy. These results might be able to (i) explain the low intensity of the whistler heat-flux fluctuations in the solar wind (although other explanations remain possible and need further investigation), and (ii) confirm a reduced effectiveness of these fluctuations in the relaxation and isotropization of the electron strahl and in the regulation of the electron heat-flux.

In the solar wind an equipartition of kinetic energy densities can be easily established between thermal and suprathermal electrons and the instability conditions are markedly altered by the interplay of these two populations. The new thresholds derived here for the periodic branch of firehose instability shape the limits of temperature anisotropy reported by the observations for both electron populations. This instability constraint is particularly important for the suprathermal electrons which, by comparison to thermal populations, are even less controlled by the particle-particle collisions. An extended quasilinear approach of this instability confirms predictions from linear theory and unveil the mutual effects of thermal and suprathermal electrons in the relaxation of their temperature anisotropies and the saturation of growing fluctuations.

In collision-poor plasmas from space, e.g., solar wind, terrestrial magnetospheres, kinetic instabilities are expected to play a major role in constraining the temperature anisotropy of plasma particles, but a definitive answer can be given only after ascertaining their properties in these environments. Present study describes the full spectrum of electron firehose instabilities in the presence of suprathermal electron populations which are ubiquitous in space plasmas. Suprathermal electrons stimulate both the periodic and aperiodic branches, remarkable being the effects shown by the aperiodic mode propagating obliquely to the ambient magnetic field which markedly exceeds the growth rates of the parallel (periodic) branch reported recently in Lazar et al., (2017a, MNRAS 464, 564). Derived exclusively in terms of the plasma parameters, the anisotropy thresholds of this instability are also lowered in the presence of suprathermal electrons, predicting an enhanced effectiveness in the solar wind conditions. These results may also be relevant in various other astrophysical contexts where the firehose instabilities involve, e.g., solar flares, sites of magnetic field reconnection, accretion flows or plasma jets leading to shocks and co-rotating interactions in heliosphere, interstellar medium and galaxy clusters.

We studied the dynamics of all prominence tornadoes detected by the Solar Dynamics Observatory/Atmospheric Imaging Assembly from 2011 January 01 to December 31. In total, 361 events were identified during the whole year, but only 166 tornadoes were traced until the end of their lifetime. Out of 166 tornadoes, 80 (48%) triggered CMEs in hosting prominences, 83 (50%) caused failed coronal mass ejections (CMEs) or strong internal motion in the prominences, and only 3 (2%) finished their lifetimes without any observed activity. Therefore, almost all prominence tornadoes lead to the destabilization of their hosting prominences and half of them trigger CMEs. Consequently, prominence tornadoes may be used as precursors for CMEs and hence for space weather predictions.

In the solar wind electron velocity distributions reveal two counter-moving populations which may induce electromagnetic (EM) beaming instabilities known as heat flux instabilities. Depending on plasma parameters two distinct branches of whistler and firehose instabilities can be excited. These instabilities are invoked in many scenarios, but their interplay is still poorly understood. An exact numerical analysis is performed to resolve the linear Vlasov-Maxwell dispersion and characterize these two instabilities, e.g., growth rates, wave frequencies and thresholds, enabling to identify their dominance for conditions typically experienced in space plasmas. Of particular interest are the effects of suprathermal Kappa-distributed electrons which are ubiquitous in these environments. The dominance of whistler or firehose instability is highly conditioned by the beam-core relative velocity, core plasma beta and the abundance of suprathermal electrons. Derived in terms of relative drift velocity the instability thresholds show an inverse correlation with the core plasma beta for the whistler modes, and a direct correlation with the core plasma beta for the firehose instability. Suprathermal electrons reduce the effective (beaming) anisotropy inhibiting the firehose modes while the whistler instability is stimulated.

In space plasmas kinetic instabilities are driven by the beaming (drifting) components and/or the temperature anisotropy of charged particles. The heat-flux instabilities are known in the literature as electromagnetic modes destabilized by the electron beams (or strahls) aligned to the interplanetary magnetic field. A new kinetic approach is proposed here in order to provide a realistic characterization of heat-flux instabilities under the influence of electrons with temperature anisotropy. Numerical analysis is based on the kinetic Vlasov-Maxwell theory for two electron counter-streaming (core and beam) populations with temperature anisotropies, and stationary, isotropic protons. The main properties of electromagnetic heat-flux instabilities are found to be markedly changed by the temperature anisotropy of electron beam $A_b = T_\perp / T_\parallel \ne 1$, leading to stimulation of either the whistler branch if $A_b > 1$, or the firehose branch for $A_b<1$. For a high temperature anisotropy whistlers switch from heat-flux to a standard regime, when their instability is inhibited by the beam.

Results of solar observations at 20 and 25 MHz by the UTR-2 (\textit{Ukrainian T-shaped Radio telescope of the second modification}) radio telescope in the interferometric session from 27 May to 2 June 2014 are presented. In such a case the different baselines 225, 450, and 675 m between sections of East--West and North--South arms of the radio telescope UTR-2 were used. On 29 May 2014, strong sporadic radio emission consisting of Type III, a Type II and a Type IV bursts was observed. On other days there was no solar radio activity in the decameter range. We discuss the results of observations of such the quiet Sun. Fluxes and sizes of the Sun in East--West and North--South directions were measured. The average fluxes were 1050--1100 Jy and 1480--1570 Jy at 20 and 25 MHz, respectively. Angular sizes of the quiet Sun in equatorial and polar directions were 55' and 49' at 20 MHz and 50' and 42' at 25 MHz. Brightness temperatures of radio emission were Tb = 5.1 10^5 K and Tb = 5.7 10^5 K at 20 and 25 MHz, respectively.

A set of 23 observations of coronal jet events that occurred in coronal bright points has been analyzed. The focus was on the temporal evolution of the mean brightness before and during coronal jet events. In the absolute majority of the cases either single or recurrent coronal jets were preceded by slight precursor disturbances observed in the mean intensity curves. The key conclusion is that we were able to detect quasi-periodical oscillations with characteristic periods from sub-minute up to 3-4 min values in the bright point brightness which precede the jets. Our basic claim is that along with the conventionally accepted scenario of bright point evolution through new magnetic flux emergence and its reconnection with the initial structure of the bright point and the coronal hole, certain MHD oscillatory and wave-like motions can be excited and these can take an important place in the observed dynamics. These quasi-oscillatory phenomena might play the role of links between different epochs of the coronal jet ignition and evolution. They can be an indication of the MHD wave excitation processes due to the system entropy variations, density variations or shear flows. It is very likely a sharp outflow velocity transverse gradients at the edges between the open and closed field line regions. We suppose that magnetic reconnections can be the source of MHD waves due to impulsive generation or rapid temperature variations, and shear flow driven nonmodel MHD wave evolution (self-heating and/or overreflection mechanisms).

We apply a novel pre-flare tracking of sunspot groups towards improving the estimation of flare onset time by focusing on the evolution of the 3D magnetic field construction of AR 11429. The 3D magnetic structure is based on potential field extrapolation encompassing a vertical range from the photosphere through the chromosphere and transition region into the low corona. The basis of our proxy measure of activity prediction is the so-called weighted horizontal gradient of magnetic field (WG_M) defined between spots of opposite polarities close to the polarity inversion line of an active region. The temporal variation of the distance of the barycenter of the opposite polarities is also found to possess potentially important diagnostic information about the flare onset time estimation as function of height similar to its counterpart introduced initially in an application at the photosphere only in Korsos et al. (2015). We apply the photospheric pre-flare behavioural patterns of sunspot groups to the evolution of their associated 3D-constructed AR 11429 as function of height. We found that at a certain height in the lower solar atmosphere the onset time may be estimated much earlier than at the photosphere or at any other heights. Therefore, we present a tool and recipe that may potentially identify the optimum height for flare prognostic in the solar atmosphere allowing to improve our flare prediction capability and capacity.

In this study we present a statistical analysis of 53 fast Earth-directed halo CMEs observed by the SOHO/LASCO instrument during the period Jan. 2009-Sep. 2015, and we use this CME sample to test the capabilities of a Sun-to-Earth prediction scheme for CME geoeffectiveness. First, we investigate the CME association with other solar activity features by means of multi-instrument observations of the solar magnetic and plasma properties. Second, using coronagraphic images to derive the CME kinematical properties at 0.1 AU, we propagate the events to 1 AU by means of the WSA-ENLIL+Cone model. Simulation results at Earth are compared with in-situ observations at L1. By applying the pressure balance condition at the magnetopause and a solar wind-Kp index coupling function, we estimate the expected magnetospheric compression and geomagnetic activity level, and compare them with global data records. The analysis indicates that 82% of the CMEs arrived at Earth in the next 4 days. Almost the totality of them compressed the magnetopause below geosynchronous orbits and triggered a geomagnetic storm. Complex sunspot-rich active regions associated with energetic flares result the most favourable configurations from which geoeffective CMEs originate. The analysis of related SEP events shows that 74% of the CMEs associated with major SEPs were geoeffective. Moreover, the SEP production is enhanced in the case of fast and interacting CMEs. In this work we present a first attempt at applying a Sun-to-Earth geoeffectiveness prediction scheme - based on 3D simulations and solar wind-geomagnetic activity coupling functions - to a statistical set of potentially geoeffective halo CMEs. The results of the prediction scheme are in good agreement with geomagnetic activity data records, although further studies performing a fine-tuning of such scheme are needed.

The present paper discusses results of a statistical study of the characteristics of coronal hole (CH) rotation in order to find connections to the internal rotation of the Sun. The goal is to measure CH rotation rates and study their distribution over latitude and their area sizes. In addition, the CH rotation rates are compared with the solar photospheric and inner layer rotational profiles. We study coronal holes observed within $\pm 60$ latitude and longitude degrees from the solar disc centre during the time span from the 1 January 2013 to 20 April 2015, which includes the extended peak of solar cycle 24.We used data created by the Spatial Possibilistic Clustering Algorithm (SPoCA), which provides the exact location and characterisation of solar coronal holes using SDO=AIA 193 {\AA} channel images. The CH rotation rates are measured with four-hour cadence data to track variable positions of the CH geometric centre. North-south asymmetry was found in the distribution of coronal holes: about 60 percent were observed in the northern hemisphere and 40 percent were observed in the southern hemisphere. The smallest and largest CHs were present only at high latitudes. The average sidereal rotation rate for 540 examined CHs is $13:86 (\pm 0:05)$ degrees/d. Conclusions. The latitudinal characteristics of CH rotation do not match any known photospheric rotation profile. The CH angular velocities exceed the photospheric angular velocities at latitudes higher than 35-40 degrees. According to our results, the CH rotation profile perfectly coincides with tachocline and the lower layers of convection zone at around 0.71 $R_{\odot}$; this indicates that CHs may be linked to the solar global magnetic field, which originates in the tachocline region.

This letter presents the results of an advanced parametrization of the solar wind electron temperature anisotropy and the instabilities resulting from the interplay of the (bi-)Maxwellian core and (bi-)Kappa halo populations in the slow solar wind. A large set of observational data (from the Ulysses, Helios and Cluster missions) is used to parametrize these components and establish their correlations. The instabilities are significantly stimulated in the presence of suprathermals, and the instability thresholds shape the limits of the temperature anisotropy for both the core and halo populations re-stating the incontestable role that the selfgenerated instabilities can play in constraining the electron anisotropy. These results confirm a particular implication of the suprathermal electrons which are less dense but hotter than thermal electrons.

Context. Recent studies on Kappa distribution functions invoked in space plasma applications have emphasized two alternative approaches which may assume the temperature parameter either dependent or independent of the power-index $\kappa$. Each of them can obtain justification in different scenarios involving Kappa-distributed plasmas, but direct evidences supporting any of these two alternatives with measurements from laboratory or natural plasmas are not available yet. Aims. This paper aims to provide more facts on this intriguing issue from direct fitting measurements of suprathermal electron populations present in the solar wind, as well as from their destabilizing effects predicted by these two alternating approaches. Methods. Two fitting models are contrasted, namely, the global Kappa and the dual Maxwellian-Kappa models, which are currently invoked in theory and observations. The destabilizing effects of suprathermal electrons are characterized on the basis of a kinetic approach which accounts for the microscopic details of the velocity distribution. Results. In order to be relevant, the model is chosen to accurately reproduce the observed distributions and this is achieved by a dual Maxwellian-Kappa distribution function. A statistical survey indicates a $\kappa$-dependent temperature of the suprathermal (halo) electrons for any heliocentric distance. Only for this approach the instabilities driven by the temperature anisotropy are found to be systematically stimulated by the abundance of suprathermal populations, i.e., lowering the values of $\kappa$-index.

On 2010 August 14, a wide-angled coronal mass ejection (CME) was observed. This solar eruption originated from a destabilized filament that connected two active regions and the unwinding of this filament gave the eruption an untwisting motion that drew the attention of many observers. In addition to the erupting filament and the associated CME, several other low-coronal signatures that typically indicate the occurrence of a solar eruption were associated to this event. However, contrary to what is expected, the fast CME ($\mathrm{v}>900~\mathrm{km}~\mathrm{s}^{-1}$) was accompanied by only a weak C4.4 flare. We investigate the various eruption signatures that were observed for this event and focus on the kinematic evolution of the filament in order to determine its eruption mechanism. Had this solar eruption occurred just a few days earlier, it could have been a significant event for space weather. The risk to underestimate the strength of this eruption based solely on the C4.4 flare illustrates the need to include all eruption signatures in event analyses in order to obtain a complete picture of a solar eruption and assess its possible space weather impact.

Context. The dynamics of the flaring loops in active region (AR) 11429 are studied. The observed dynamics consist of several evolution stages of the flaring loop system during both the ascending and descending phases of the registered M-class flare. The dynamical properties can also be classified by different types of magnetic reconnection, related plasma ejection and aperiodic flows, quasi-periodic oscillatory motions, and rapid temperature and density changes, among others. The focus of the present paper is on a specific time interval during the ascending (pre-flare) phase. Aims. The goal is to understand the quasi-periodic behavior in both space and time of the magnetic loop structures during the considered time interval. Methods.We have studied the characteristic location, motion, and periodicity properties of the flaring loops by examining space-time diagrams and intensity variation analysis along the coronal magnetic loops using AIA intensity and HMI magnetogram images (from the Solar Dynamics Observatory(SDO)). Results. We detected bright plasma blobs along the coronal loop during the ascending phase of the solar flare, the intensity variations of which clearly show quasi-periodic behavior. We also determined the periods of these oscillations. Conclusions. Two different interpretations are presented for the observed dynamics. Firstly, the oscillations are interpreted as the manifestation of non-fundamental harmonics of longitudinal standing acoustic oscillations driven by the thermodynamically nonequilibrium background (with time variable density and temperature). The second possible interpretation we provide is that the observed bright blobs could be a signature of a strongly twisted coronal loop that is kink unstable.

A variety of nonthermal characteristics like kinetic, e.g., temperature, anisotropies and suprathermal populations (enhancing the high energy tails of the velocity distributions) are revealed by the in-situ observations in the solar wind indicating quasistationary states of plasma particles out of thermal equilibrium. Large deviations from isotropy generate kinetic instabilities and growing fluctuating fields which should be more efficient than collisions in limiting the anisotropy (below the instability threshold) and explain the anisotropy limits reported by the observations. The present paper aims to decode the principal instabilities driven by the temperature anisotropy of electrons and protons in the solar wind, and contrast the instability thresholds with the bounds observed at 1~AU for the temperature anisotropy. The instabilities are characterized using linear kinetic theory to identify the appropriate (fastest) instability in the relaxation of temperature anisotropies $A_{e,p} = T_{e,p,\perp}/ T_{e,p,\parallel} \ne 1$. The analysis focuses on the electromagnetic instabilities driven by the anisotropic protons ($A_p \lessgtr 1$) and invokes for the first time a dynamical model to capture the interplay with the anisotropic electrons by correlating the effects of these two species of plasma particles, dominant in the solar wind.

In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ion-neutral interactions in the partially ionized solar chromosphere, we have developed a new multi-fluid computational model, accounting for ionization and recombination reactions in gravitationally stratified magnetized collisional media. The two-fluid model used in our 2D numerical simulations treats neutrals as a separate fluid and considers charged species (electrons and ions) within the resistive MHD approach with Coulomb collisions and anisotropic heat flux determined by Braginskii's transport coefficients. The electromagnetic fields are evolved according to the full Maxwell equations and the solenoidality of the magnetic field is enforced with a hyperbolic divergence cleaning scheme. The initial density and temperature profiles are similar to VAL III chromospheric model in which dynamical, thermal and chemical equilibrium are considered to ensure comparison to existing MHD models and avoid artificial numerical heating. In this initial setup we include simple homogeneous flux tube magnetic field configuration and an external photospheric velocity driver to simulate the propagation of MHD waves in the partially ionized reactive chromosphere. In particular, we investigate the loss of chemical equilibrium and the plasma heating related to the steepening of fast magnetosonic wave fronts in the gravitationally stratified medium.

Active regions (ARs) are the main sources of variety in solar dynamic events. Automated detection and identification tools need to be developed for solar features for a deeper understanding of the solar cycle. Of particular interest here are the dynamical properties of the ARs, regardless of their internal structure and sunspot distribution. We studied the oscillatory dynamics of two ARs: NOAA 11327 and NOAA 11726 using two different methods of pattern recognition. We developed a novel method of automated AR border detection and compared it to an existing method for the proof-of-concept. The first method uses least-squares fitting on the smallest ellipse enclosing the AR, while the second method applies regression on the convex hull.} After processing the data, we found that the axes and the inclination angle of the ellipse and the convex hull oscillate in time. These oscillations are interpreted as the second harmonic of the standing long-period kink oscillations (with the node at the apex) of the magnetic flux tube connecting the two main sunspots of the ARs. In both ARs we have estimated the distribution of the phase speed magnitude along the magnetic tubes (along the two main spots) by interpreting the obtained oscillation of the inclination angle as the standing second harmonic kink mode. After comparing the obtained results for fast and slow kink modes, we conclude that both of these modes are good candidates to explain the observed oscillations of the AR inclination angles, as in the high plasma $\beta$ regime the phase speeds of these modes are comparable and on the order of the Alfv\'{e}n speed. Based on the properties of the observed oscillations, we detected the appropriate depth of the sunspot patterns, which coincides with estimations made by helioseismic methods. The latter analysis can be used as a basis for developing a magneto-seismological tool for ARs.

Many natural processes exhibit power-law behavior. The power-law exponent is linked to the underlying physical process and therefore its precise value is of interest. With respect to the energy content of nanoflares, for example, a power-law exponent steeper than 2 is believed to be a necessary condition to solve the enigmatic coronal heating problem. Studying power-law distributions over several orders of magnitudes requires sufficient data and appropriate methodology. In this paper we demonstrate the shortcomings of some popular methods in solar physics that are applied to data of typical sample sizes. We use synthetic data to study the effect of the sample size on the performance of different estimation methods and show that vast amounts of data are needed to obtain a reliable result with graphical methods (where the power-law exponent is estimated by a linear fit on a log-transformed histogram of the data). We revisit published results on power laws for the angular width of solar coronal mass ejections and the radiative losses of nanoflares. We demonstrate the benefits of the maximum likelihood estimator and advocate its use.

Solar activity undergoes a variation over time scales of several months known as Rieger-type periodicity, which usually occurs near maxima of sunspot cycles. An early analysis showed that the periodicity appears only in some cycles, and is absent in other cycles. But the appearance/absence during different cycles has not been explained. We performed a wavelet analysis of sunspot data from the Greenwich Royal Observatory and the Royal Observatory of Belgium during cycles 14-24. We found that the Rieger-type periods occur in all cycles, but they are cycle-dependent: shorter periods occur during stronger cycles. Our analysis revealed a periodicity of 185-195 days during the weak cycles 14-15 and 24, and a periodicity of 155-165 days during the stronger cycles 16-23. We derived the dispersion relation of the spherical harmonics of the magnetic Rossby waves in the presence of differential rotation and a toroidal magnetic field in the dynamo layer near the base of the convection zone. This showed that the harmonic of fast Rossby waves with m=1 and n=4, where m (n) indicate the toroidal (poloidal) wavenumbers, respectively, perfectly fit with the observed periodicity. The variation of the toroidal field strength from weaker to stronger cycles may lead to the different periods found in those cycles, which explains the observed enigmatic feature of the Rieger-type periodicity. Finally, we used the observed periodicity to estimate the dynamo field strength during cycles 14-24. Our estimations suggest a field strength of 40 kG for the stronger cycles, and 20 kG for the weaker cycles.

In collision-poor plasmas from space, e.g., the solar wind and planetary magnetospheres, the kinetic anisotropy of the plasma particles is expected to be regulated by the kinetic instabilities. Driven by an excess of ion (proton) temperature perpendicular to the magnetic field $(~T_\perp >T_\parallel)$, the electromagnetic ion-cyclotron (EMIC) instability is fast enough to constrain the proton anisotropy, but the observations do not conform to the instability thresholds predicted by the standard theory for bi-Maxwellian models of the plasma particles. This paper presents an extended investigation of the EMIC instability in the presence of suprathermal electrons which are ubiquitous in these environments. The analysis is based on the kinetic (Vlasov-Maxwell) theory assuming that both species, protons and electrons, may be anisotropic, and the EMIC unstable solutions are derived numerically providing an accurate description for conditions typically encountered in space plasmas. The effects of suprathermal populations are triggered by the electron anisotropy and the temperature contrast between electrons and protons. For certain conditions the anisotropy thresholds exceed the limits of the proton anisotropy measured in the solar wind considerably restraining the unstable regimes of the EMIC modes.

The indefinite increase of temperature predicted by the solar wind expansion in the direction parallel to the interplanetary magnetic field is already notorious for not being confirmed by the observations. In hot and dilute plasmas from space particle-particle collisions are not efficient in constraining large deviations from isotropy, but the resulting firehose instability provides in this case plausible limitations for the temperature anisotropy of the thermal (core) populations of both the electron and proton species. The present paper takes into discussion the suprathermal (halo) electrons, which are ubiquitous in the solar wind. Less dense but hotter than the core, suprathermals may be highly anisotropic and susceptible to the firehose instability. The main features of the instability are here derived from a first-order theory for conditions specific to the suprathermal electrons in the solar wind and terrestrial magnetospheres. Unveiled here, new regimes of the electron firehose instability may be exclusively controlled by the suprathermals. The instability is found to be systematically stimulated by the suprathermal electrons, with thresholds that approach the limits of the temperature anisotropy reported by the observations. These results represent new and valuable evidences for the implication of the firehose instability in the relaxation of the temperature anisotropy in space plasmas.

Estimating the temperature of the solar wind particles and their anisotropies is particularly important for understanding the origin of these deviations from thermal equilibrium as well as their effects. In the absence of energetic events the velocity distribution of electrons reveal a dual structure with a thermal (Maxwellian) core and a suprathermal (Kappa) halo. This paper presents a detailed observational analysis of these two components, providing estimations of their temperatures and temperature anisotropies and decoding any potential interdependence that their properties may indicate. The data set used in this study includes more than 120$\,$000 the events detected by three missions in the ecliptic within an extended range of heliocentric distances from 0.3 to over 4~AU. The anti-correlation found for the core and halo temperatures is consistent with the radial evolution of the Kappa model, clarifying an apparent contradiction in previous observational analysis and providing valuable clues about the temperature of the Kappa-distributed populations. However, these two components manifest a clear tendency to deviate from isotropy in the same direction, that seems to confirm the existence of mechanisms with similar effects on both components, e.g., the solar wind expansion, or the particle heating by the fluctuations. On the other hand, the existence of plasma states with anti-correlated anisotropies of the core and halo populations suggests a dynamic interplay of these components, mediated, most probably, by the anisotropy-driven instabilities.

We present a linear stability analysis of the perturbation modes in anisotropic MHD flows with velocity shear and strong magnetic field. Collisionless or weakly collisional plasma is described within the 16-momentum MHD fluid closure model, that takes into account not only the effect of pressure anisotropy, but also the effect of anisotropic heat fluxes. In this model the low frequency acoustic wave is revealed into a standard acoustic mode and higher frequency fast thermo-acoustic and lower frequency slow thermo-acoustic waves. It is shown that thermo-acoustic waves become unstable and grow exponentially when the heat flux parameter exceeds some critical value. It seems that velocity shear makes thermo-acoustic waves overstable even at subcritical heat flux parameters. Thus, when the effect of heat fluxes is not profound acoustic waves will grow due to the velocity shear, while at supercritical heat fluxes the flow reveals compressible thermal instability. Anisotropic thermal instability should be also important in astrophysical environments, where it will limit the maximal value of magnetic field that a low density ionized anisotropic flow can sustain.

The characteristics of the type II bursts with herringbone structure observed both by ground based radio telescopes (UTR-2, URAN-2) and spaceborn spectrometers (STEREO A-B) are discussed. The burst was recorded on 7 June, 2011 in the frequency band 3--33~MHz. It was characterized by extremely rich fine structure. The statistical analysis of more than 300 herringbone sub-bursts constituting the burst was performed separately for the positively (reverse) and negatively (forward) drifting sub-bursts. The sense and the degree of circular polarization of the herringbone sub-bursts were measured in the respectively wide frequency band (16--32~MHz). A second order fine frequency structure of the herringbone sub-bursts was firstly observed and processed. Using STEREO COR1 (A,B) and SOHO LASCO C2 images the direction and radial speed of the CME responsible for the studied type II burst were determined. The possible location of the type II burst source on the flank of the shock was found.

We study the detailed dynamics of a solar prominence tornado using time series of 171, 304, 193 and 211 {\AA} spectral lines obtained by Solar Dynamics Observatory/ Atmospheric Imaging Assembly during 2012 November 7-8. The tornado first appeared at 08:00 UT, November 07, near the surface, gradually rose upwards with the mean speed of $\sim$ 1.5 km s$^{-1}$ and persisted over 30 hr. Time-distance plots show two patterns of quasi-periodic transverse displacements of the tornado axis with periods of 40 and 50 minute at different phases of the tornado evolution. The first pattern occurred during the rising phase and can be explained by the upward motion of the twisted tornado. The second pattern occurred during the later stage of evolution when the tornado already stopped rising and could be caused either by MHD kink waves in the tornado or by the rotation of two tornado threads around a common axis. The later hypothesis is supported by the fact that the tornado sometimes showed a double structure during the quasi-periodic phase. 211 and 193 {\AA} spectral lines show a coronal cavity above the prominence/tornado, which started expansion at $\sim$ 13:00 UT and continuously rose above the solar limb. The tornado finally became unstable and erupted together with the corresponding prominence as coronal mass ejection (CME) at 15:00 UT, November 08. The final stage of the evolution of the cavity and the tornado-related prominence resembles the magnetic breakout model. On the other hand, the kink instability may destabilize the twisted tornado, and consequently prominence tornadoes can be used as precursors for CMEs.

The formation and dynamics of coronal rain are currently not fully understood. Coronal rain is the fall of cool and dense blobs formed by thermal instability in the solar corona towards the solar surface with acceleration smaller than gravitational free fall. We aim to study the observational evidence of the formation of coronal rain and to trace the detailed dynamics of individual blobs. We used time series of the 171 \AA\, and 304 \AA\, spectral lines obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) above active region AR 11420 on February 22, 2012. Observations show that a coronal loop disappeared in the 171 \AA\ channel and appeared in the 304 \AA\ line$\text{}\text{}$ more than one hour later, which indicates a rapid cooling of the coronal loop from 1 MK to 0.05 MK. An energy estimation shows that the radiation is higher than the heat input, which indicates so-called catastrophic cooling. The cooling was accompanied by the formation of coronal rain in the form of falling cold plasma. We studied two different sequences of falling blobs. The first sequence includes three different blobs. The mean velocities of the blobs were estimated to be 50 km s$^{-1}$, 60 km s$^{-1}$ and 40 km s$^{-1}$. A polynomial fit shows the different values of the acceleration for different blobs, which are lower than free-fall in the solar corona. The first and second blob move along the same path, but with and without acceleration, respectively. We performed simple numerical simulations for two consecutive blobs, which show that the second blob moves in a medium that is modified by the passage of the first blob. Therefore, the second blob has a relatively high speed and no acceleration, as is shown by observations. The second sequence includes two different blobs with mean velocities of 100 km s$^{-1}$ and 90 km s$^{-1}$, respectively.

In plasmas where the thermal energy density exceeds the magnetic energy density ($\beta_\parallel > 1$), the aperiodic ordinary mode (O-mode) instability is driven by an excess of parallel temperature $A = T_\perp /T_\parallel < 1$ (where $\parallel$ and $\perp$ denote directions relative to the uniform magnetic field). When stimulated by parallel plasma streams the instability conditions extend to low beta states, i.e., $\beta_\parallel <1$, and recent studies have proven the existence of a new regime, where the anisotropy threshold decreases steeply with lowering $\beta_\parallel \to 0$ if the streaming velocity is sufficiently high. However, the occurrence of this instability is questionable especially in the low-beta plasmas, where the electrostatic two-stream instabilities are expected to develop much faster in the process of relaxation of the counterstreams. It is therefore proposed here to identify the instability conditions for the O-mode below those required for the onset of the electrostatic instability. An hierarchy of these two instabilities is established for both the low $\beta_\parallel <1$ and large $\beta_\parallel > 1$ plasmas. The conditions where the O-mode instability can operate efficiently are markedly constrained by the electrostatic instabilities especially in the low-beta plasmas.

Solar eruptions are usually associated with a variety of phenomena occurring in the low corona before, during, and after onset of eruption. Though easily visible in coronagraph observations, so-called stealth coronal mass ejections (CMEs) do not obviously exhibit any of these low-coronal signatures. The presence or absence of distinct low coronal signatures can be linked to different theoretical models to establish the mechanisms by which the eruption is initiated and driven. In this study, 40 CMEs without low coronal signatures, occurring in 2012, are identified. Their observational and kinematic properties are analyzed and compared to those of regular CMEs. Solar eruptions without clear on-disk or low coronal signatures can lead to unexpected space weather impacts, since many early warning signs for significant space weather activity are not present in these events. A better understanding of their initiation mechanism(s) will considerably improve the ability to predict such space weather events.

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.

Counterstreaming beams of electrons are ubiquitous in coronal mass ejections (CMEs) - although their existence is not unanimously accepted as a necessary and/or sufficient signature of these events. We continue the investigations of a high-latitude CME registered by the \emph{Ulysses} spacecraft on January 18\,--\,19, 2002 (Dumitrache, Popescu, and Oncica, Solar Phys. {\bf 272}, 137, 2011), by surveying the solar wind electron distributions associated with this event. The temporal-evolution of the pitch-angle distributions reveal populations of electrons distinguishable through their anisotropy, with clear signatures of i) electron strahls, ii) counter-streaming in the magnetic clouds and their precursors, and iii) unidirectional in the fast wind preceding the CME. The analysis of the counter-streams inside the CME allows us to elucidate the complexity of the magnetic-cloud structures embeded in the CME and to refine the borders of the event. Identifying such strahls in CMEs, which preserve properties of the low $\beta < 1$ coronal plasma, gives more support to the hypothesis that these populations are remnants of the hot coronal electrons that escape from the electrostatic potential of the Sun into the heliosphere.

Solar filaments are magnetic structures often observed in the solar atmosphere and consist of plasma that is cooler and denser than their surroundings. They are visible for days -- and even weeks -- which suggests that they are often in equilibrium with their environment before disappearing or erupting. Several eruption models have been proposed that aim to reveal what mechanism causes (or triggers) these solar eruptions. Validating these models through observations represents a fundamental step in our understanding of solar eruptions. We present an analysis of the observation of a filament eruption that agrees with the torus instability model. This model predicts that a magnetic flux rope embedded in an ambient field undergoes an eruption when the axis of the flux rope reaches a critical height that depends on the topology of the ambient field. We use the two vantage points of SDO and STEREO to reconstruct the three-dimensional shape of the filament, to follow its morphological evolution and to determine its height just before eruption. The magnetograms acquired by SDO/HMI are used to infer the topology of the ambient field and to derive the critical height for the onset of the torus instability. Our analysis shows that the torus instability is the trigger of the eruption. We also find that some pre-eruptive processes, such as magnetic reconnection during the observed flares and flux cancellation at the neutral line, facilitated the eruption by bringing the filament to a region where the magnetic field was more vulnerable to the torus instability.

We investigate electrostatic plasma instabilities of Farley-Buneman (FB) type driven by quasi-stationary neutral gas flows in the solar chromosphere. The role of these instabilities in the chromosphere is clarified. We find that the destabilizing ion thermal effect is highly reduced by the Coulomb collisions and can be ignored for the chromospheric FB-type instabilities. On the contrary, the destabilizing electron thermal effect is important and causes a significant reduction of the neutral drag velocity triggering the instability. The resulting threshold velocity is found as function of chromospheric height. Our results indicate that the FB type instabilities are still less efficient in the global chromospheric heating than the Joule dissipation of the currents driving these instabilities. This conclusion does not exclude the possibility that the FB type instabilities develop in the places where the cross-field currents overcome the threshold value and contribute to the heating locally. Typical length-scales of plasma density fluctuations produced by these instabilities are determined by the wavelengths of unstable modes, which are in the range $10-10^2$ cm in the lower chromosphere, and $10^2-10^3$ cm in the upper chromosphere. These results suggest that the decimetric radio waves undergoing scattering (scintillations) by these plasma irregularities can serve as a tool for remote probing of the solar chromosphere at different heights.

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 study the variation of the shape of the proton intensity-time profiles in simulated gradual Solar Energetic Particle (SEP) events with the relative observer's position in space with respect to the main direction of propagation of an interplanetary (IP) shock. Using a three-dimensional (3D) magnetohydrodynamic (MHD) code to simulate such a shock, we determine the evolution of the downstream-to-upstream ratios of the plasma variables at its front. Under the assumption of an existing relation between the normalized ratio in speed across the shock front and the injection rate of shock-accelerated particles, we model the transport of the particles and we obtain the proton flux profiles to be measured by a grid of 18 virtual observers located at 0.4 and 1.0 AU, with different latitudes and longitudes with respect to the shock nose. The differences among flux profiles are the result of the way each observer establishes a magnetic connection with the shock front, and we find that changes in the observer's latitude may result in intensity changes of up to one order of magnitude at both radial distances considered here. The peak intensity variation with the radial distance for the pair of observers located at the same angular position is also derived. This is the first time that the latitudinal dependence of the peak intensity with the observer's heliocentric radial distance has been quantified within the framework of gradual SEP event simulations.

The selfgenerated wave fluctuations are particularly interesting in the solar wind and magnetospheric plasmas, where Coulomb collisions are rare and cannot explain the observed states of quasi-equilibrium. Linear theory predicts that the firehose and the ordinary-mode instabilities can develop under the same conditions, confusing the role of these instabilities in conditioning the space-plasma properties. The hierarchy of these two instabilities is reconsidered here for nonstreaming plasmas with an electron temperature anisotropy $T_\parallel > T_\perp$, where $\parallel$ and $\perp$ denote directions with respect to the local mean magnetic field. In addition to the previous comparative analysis, here the entire 3D wave-vector spectrum of the competing instabilities is investigated, paying particular attention to the oblique firehose instability and the relatively poorly known ordinary-mode instability. Results show a dominance of the oblique firehose instability with a threshold lower than the parallel firehose instability and lower than the ordinary-mode instability. For larger anisotropies, the ordinary mode can grow faster, with maximum growth rates exceeding the ones of the oblique firehose instability. In contrast to previous studies that claimed a possible activity of the ordinary-mode in the small $\beta [< 1]$ regimes, here it is rigorously shown that only the large $\beta [> 1]$ regimes are susceptible to these instabilities.

On 15 May 2005, a huge interplanetary coronal mass ejection (ICME) was observed near Earth. It triggered one of the most intense geomagnetic storms of solar cycle 23 (Dst peak = -263 nT). This structure has been associated with the two-ribbon flare, filament eruption, and coronal mass ejection originating in active region 10759 (NOAA number). We analyze here the sequence of events, from solar wind measurements (at 1 AU) and back to the Sun, to understand the origin and evolution of this geoeffective ICME. From a detailed observational study of in situ magnetic field observations and plasma parameters in the interplanetary (IP) medium and the use of appropriate models we propose an alternative interpretation of the IP observations, different to those discussed in previous studies. In our view, the IP structure is formed by two extremely close consecutive magnetic clouds (MCs) that preserve their identity during their propagation through the interplanetary medium. Consequently, we identify two solar events in H{\alpha} and EUV which occurred in the source region of the MCs. The timing between solar and IP events, as well as the orientation of the MC axes and their associated solar arcades are in good agreement. Additionally, interplanetary radio type II observations allow the tracking of the multiple structures through inner heliosphere and pin down the interaction region to be located midway between the Sun and the Earth. The chain of observations from the photosphere to interplanetary space is in agreement with this scenario. Our analysis allows the detection of the solar sources of the transients and explains the extremely fast changes of the solar wind due to the transport of two attached (though nonmerging) MCs which affect the magnetosphere.

The non-modal self-heating mechanism driven by the velocity shear in kinematically complex magnetohydrodynamic (MHD) plasma flows is considered. The study is based on the full set of MHD equations including dissipative terms. The equations are linearized and unstable modes in the flow are looked for. Two different cases are specified and studied: (a) the instability related to an exponential evolution of the wave vector; and (b) the parametric instability, which takes place when the components of the wave vector evolve in time periodically. By examining the dissipative terms, it is shown that the self-heating rate provided by viscous damping is of the same order of magnitude as that due to the magnetic resistivity. It is found that the heating efficiency of the exponential instability is higher than that of the parametric instability.

This study aims to provide a reference to different magnetic field models and reconstruction methods for interplanetary coronal mass ejections (ICMEs). In order to understand the differences in the outputs of those models and codes, we analyze 59 events from the Coordinated Data Analysis Workshop (CDAW) list, using four different magnetic field models and reconstruction techniques; force-free fitting (Goldstein,1983,Burlaga,1988,Lepping et al.,1990), magnetostatic reconstruction using a numerical solution to the Grad-Shafranov equation (Hu and Sonnerup, 2001), fitting to a self-similarly expanding cylindrical configuration (Marubashi and Lepping, 2007) and elliptical, non-force free fitting (Hidalgo,2003). The resulting parameters of the reconstructions for the 59 events are compared statistically, as well as in selected case studies. The ability of a method to fit or reconstruct an event is found to vary greatly: the Grad-Shafranov reconstruction is successful for most magnetic clouds (MCs) but for less than 10% of the non-MC ICMEs; the other three methods provide a successful fit for more than 65% of all events. The differences between the reconstruction and fitting methods are discussed, and suggestions are proposed as to how to reduce them. We find that the magnitude of the axial field is relatively consistent across models but not the orientation of the axis of the ejecta. We also find that there are a few cases for which different signs of the magnetic helicity are found for the same event when we do not fix the boundaries, illustrating that this simplest of parameters is not necessarily always well constrained by fitting and reconstruction models. Finally, we look at three unique cases in depth to provide a comprehensive idea of the different aspects of how the fitting and reconstruction codes work.

Space plasmas are collisionpoor and kinetic effects prevail leading to wave fluctuations, which transfer the energy to small scales: wave-particle interactions replace collisions and enhance dispersive effects heating particles and producing suprathermal populations observed at any heliospheric distance in the solar wind. At large distances collisions are not efficient, and the selfgenerated instabilities constrain the solar wind anisotropy including the thermal core and the suprathermal components. The generalized power-laws of Kappa-type are the best fitting model for the observed distributions of particles, and a convenient mathematical tool for modeling their dynamics. But the anisotropic Kappa models are not correlated with the observations leading, in general, to inconsistent effects. This review work aims to reconcile some of the existing Kappa models with the observations.

Magnetic clouds (MCs) are "magnetized plasma clouds" moving in the solar wind. MCs transport magnetic flux and helicity away from the Sun. These structures are not stationary but feature temporal evolution as they propagate in the solar wind. Simplified analytical models are frequently used for the description of MCs, and fit certain observational data well. The goal of the present study is to investigate numerically the validity of an analytical model which is widely used for the description of MCs, and to determine under which conditions this model's implied assumptions cease to be valid. A numerical approach is applied. Analytical solutions that have been derived in previous studies are implemented in a \textbf{3-D magnetohydrodynamic} simulation code as initial conditions. Initially, the analytical model represents the main observational features of the MCs. However, these characteristics prevail only if the structure moves with a velocity close to the velocity of the background flow. In this case an MC's evolution can quite accurately be described using an analytic, self-similar approach. The dynamics of the magnetic structures which move with a velocity significantly above or below that of the velocity of the solar wind is investigated in detail. Besides the standard case in which MCs only expand and propagate in the solar wind, the case of an MC rotating around its axis of symmetry is also considered, and the resulting influence on the MC's dynamics is studied.

Magnetic clouds (MCs) are "magnetized plasma clouds" moving in the solar wind. MCs transport magnetic flux and helicity away from the Sun. These structures are not stationary but feature temporal evolution. Commonly, simplified MC models are considered. The goal of the present study is to investigate the dynamics of more general, radially expanding MCs. They are considered as cylindrically symmetric magnetic structures with low plasma {\beta}. In order to study MC`evolution the self-similar approach method and a numerical approach are used. It is shown that the forces are balanced in the considered self-similarly evolving, cylindrically symmetric magnetic structures. Explicit analytical expressions for magnetic field, plasma velocity, density and pressure within MCs are derived. These solutions are characterized by conserved values of magnetic flux and helicity. We also investigate the dynamics of self-similarly evolving MCs by means of the numerical code "Graale". In addition, their expansion in a medium with higher density and higher plasma {\beta} is studied. It is shown that the physical parameters of the MCs maintain their self-similar character throughout their evolution. Conclusions. A comparison of the different self-similar and numerical solutions allows us to conclude that the evolving MCs are quite adequately described by our self-similar solutions - they retain their self-similar, coherent nature for quite a long time and over large distances from the Sun.

The solar atmosphere is structured and inhomogeneous both horizontally and vertically. The omnipresence of coronal magnetic loops implies gradients of the equilibrium plasma quantities like the density, magnetic field and temperature. These gradients are responsible for the excitation of drift waves that grow both within the two-component fluid description (in the presence of collisions and without it) and within the two-component kinetic descriptions (due to purely kinetic effects). In the present work the effects of the density gradient in the direction perpendicular to the magnetic field vector are investigated within the kinetic theory, in both electrostatic and electromagnetic regimes. The electromagnetic regime implies the coupling of the gradient-driven drift wave with the Alfven wave. The growth rates for the two cases are calculated and compared. It is found that, in general, the electrostatic regime is characterized by stronger growth rates, as compared with the electromagnetic perturbations. Also discussed is the stochastic heating associated with the drift wave. The released amount of energy density due to this heating should be more dependent on the magnitude of the background magnetic field than on the coupling of the drift and Alfven waves. The stochastic heating is expected to be much higher in regions with a stronger magnetic field. On the whole, the energy release rate caused by the stochastic heating can be several orders of magnitude above the value presently accepted as necessary for a sustainable coronal heating.

In the solar corona, several mechanisms of the drift wave instability can make the mode growing up to amplitudes at which particle acceleration and stochastic heating by the drift wave take place. The stochastic heating, well known from laboratory plasma physics where it has been confirmed in numerous experiments, has been completely ignored in past studies of coronal heating. However, in the present study and in our very recent works it has been shown that the inhomogeneous coronal plasma is, in fact, a perfect environment for fast growing drift waves. As a matter of fact, the large growth rates are typically of the same order as the plasma frequency. The consequent heating rates may exceed the required values for a sustained coronal heating by several orders of magnitude. Some aspects of these phenomena are investigated here. In particular the analysis of the particle dynamics within the growing wave is compared with the corresponding fluid analysis. While both of them predict the stochastic heating, the threshold for the heating obtained from the single particle analysis is higher. The explanation for this effect is given.

The efficiency of nonmodal self-heating by acoustic wave perturbations is examined. Considering different kinds of kinematically complex velocity patterns we show that nonmodal instabilities arising in these inhomogeneous flows may lead to significant amplification of acoustic waves. Subsequently, the presence of viscous dissipation damps these amplified waves and causes the energy transfer back to the background flow in the form of heat; viz. closes the "self-heating" cycle and contributes to the net heating of the flow patterns and the chromospheric network as a whole. The acoustic self-heating depends only on the presence of kinematically complex flows and dissipation. It is argued that together with other mechanisms of nonlinear nature the self-heating \textit{may be} a probable additinal mechanism of nonmagnetic chromospheric heating in the Sun and other solar-type stars with slow rotation and extended convective regions.

The electrostatic instabilities driven by the gradients of the density, temperature and magnetic field, are discussed in their application to solar magnetic structures. Strongly growing modes are found for some typical plasma parameters. These instabilities i) imply the presence of electric fields that can accelerate the plasma particles in both perpendicular and parallel directions with respect to the magnetic field vector, and ii) can stochastically heat ions. The perpendicular acceleration is to the leading order determined by the $\bmath{E}\times \bmath{B}$-drift acting equally on both ions and electrons, while the parallel acceleration is most effective on electrons. The experimentally confirmed stochastic heating is shown to act mainly in the direction perpendicular to the magnetic field vector and acts stronger on heavier ions. The energy release rate and heating may exceed for several orders of magnitude the value accepted as necessary for a self-sustained heating in the solar corona. The energy source for both the acceleration and the heating is stored in the mentioned background gradients.

The heating of the plasma in the solar atmosphere is discussed within both frameworks of fluid and kinetic drift wave theory. We show that the basic ingredient necessary for the heating is the presence of density gradients in the direction perpendicular to the magnetic field vector. Such density gradients are a source of free energy for the excitation of drift waves. We use only well established basic theory, verified experimentally in laboratory plasmas. Two mechanisms of the energy exchange and heating are shown to take place simultaneously: one due to the Landau effect in the direction parallel to the magnetic field, and another one, stochastic heating, in the perpendicular direction. The stochastic heating i) is due to the electrostatic nature of the waves, ii) is more effective on ions than on electrons, iii) acts predominantly in the perpendicular direction, iv) heats heavy ions more efficiently than lighter ions, and v) may easily provide a drift wave heating rate that is orders of magnitude above the value that is presently believed to be sufficient for the coronal heating, i.e., $\simeq 6 \cdot 10^{-5} $J/(m$^3$s) for active regions and $\simeq 8 \cdot 10^{-6} $J/(m$^3$s) for coronal holes. This heating acts naturally through well known effects that are, however, beyond the current standard models and theories.

The main goal is to study the dynamics of the gravitationally stratified, field-free cavities in the solar atmosphere, located under small-scale, cylindrical magnetic canopies, in response to explosive events in the lower-lying regions (due to granulation, small-scale magnetic reconnection, etc.). We derive the two-dimensional Klein-Gordon equation for isothermal density perturbations in cylindrical coordinates. The equation is first solved by a standard normal mode analysis in order to obtain the free oscillation spectrum of the cavity. Then, the equation is solved in the case of impulsive forcing associated to a pressure pulse specified in the lower-lying regions. The normal mode analysis shows that the entire cylindrical cavity of granular dimensions tends to oscillate with frequencies of 5-8 mHz and also with the atmospheric cut-off frequency. Furthermore, the passage of a pressure pulse, excited in the convection zone, sets up a wake in the cavity oscillating with the same cut-off frequency. The wake oscillations can resonate with the free oscillation modes, which leads to an enhanced observed oscillation power. The resonant oscillations of these cavities explain the observed power halos near magnetic network cores and active regions.

The solar coronal heating problem refers to the question why the temperature of the Sun's corona is more than two orders of magnitude higher than that of its surface. Almost 70 years after the discovery, this puzzle is still one of the major challenges in astrophysics. The current basic paradigm of coronal heating is unable to explain all the observational features of the heating. Here we argue that a new paradigm is required to solve the puzzle in a self-consistent manner. The alternative approach is based on the kinetic theory of drift waves. We show, with qualitative and quantitative arguments, that the drift waves have the potential to satisfy all coronal heating requirements.

The Farley-Buneman instability is studied in the partially ionized plasma of the solar chromosphere taking into account the finite magnetization of the ions and Coulomb collisions. We obtain the threshold value for the relative velocity between ions and electrons necessary for the instability to develop. It is shown that Coulomb collisions play a destabilizing role in the sense that they enable the instability even in the regions where the ion magnetization is greater than unity. By applying these results to chromospheric conditions, we show that the Farley-Buneman instability can not be responsible for the quasi-steady heating of the solar chromosphere. However, in the presence of strong cross-field currents it can produce small-scale, $\sim 0.1-3$ m, density irregularities in the solar chromosphere. These irregularities can cause scintillations of radio waves with similar wave lengths and provide a tool for remote chromospheric sensing.

Well-known analytical results dealing with ion cyclotron and drift waves and which follow from the kinetic theory are used and the dispersion equation, which describes coupled two modes, is solved numerically. The numerical results obtained by using the values for the plasma density, magnetic field and temperature applicable to the solar corona clearly show the coupling and the instability (growing) of the two modes. The coupling happens at very short wavelengths, that are of the order of the ion gyro radius, and for characteristic scale lengths of the equilibrium density that are altitude dependent and may become of the order of only a few meters. The demonstrated instability of the two coupled modes (driven by the equilibrium density gradient) is obtained by using a rigorous kinetic theory model and for realistic parameter values. The physical mechanism which is behind the coupling is simple and is expected to take place throughout the solar atmosphere and the solar wind which contain a variety of very elongated density structures of various sizes. The mode grows on account of the density gradient, it is essentially an ion mode, and its further dissipation should result in an increased ion heating.

The overshooting convective motions in the solar photosphere are frequently proposed as the source for the excitation of Alfv\'en waves. However, the photosphere is a) very weakly ionized, and, b) the dynamics of the plasma particles in this region is heavily influenced by the plasma-neutral collisions. The purpose of this work is to check the consequences of these two facts on the above scenario and their effects on the electromagnetic waves. It is shown that the ions and electrons in the photosphere are both un-magnetized; their collision frequency with neutrals is much larger than the gyro-frequency. This implies that eventual Alfv\'en-type electromagnetic perturbations must involve the neutrals as well. This has the following serious consequences: i) in the presence of perturbations, the whole fluid (plasma + neutrals) moves; ii) the Alfv\'en velocity includes the total (plasma + neutrals) density and is thus considerably smaller compared to the collision-less case; iii) the perturbed velocity of a unit volume, which now includes both plasma and neutrals, becomes much smaller compared to the ideal (collision-less) case; and iv) the corresponding wave energy flux for the given parameters becomes much smaller compared to the ideal case.

Observations show the increase of high-frequency wave power near magnetic network cores and active regions in the solar lower atmosphere. This phenomenon can be explained by the interaction of acoustic waves with a magnetic field. We consider small-scale, bipolar, magnetic field canopy structure near the network cores and active regions overlying field-free cylindrical cavities of the photosphere. Solving the plasma equations we get the analytical dispersion relation of acoustic oscillations in the field-free cavity area. We found that the m = 1 mode, where m is azimuthal wave number, cannot be trapped under the canopy due to energy leakage upwards. However, higher ($m \geq 2$) harmonics can be easily trapped leading to the observed acoustic power halos under the canopy.

The propagation of compressional MHD waves is studied for an externally driven system. It is assumed that the combined action of the external sources and sinks of the entropy results in the harmonic oscillation of the entropy (and temperature) in the system. It is found that with the appropriate resonant conditions fast and slow waves get amplified due to the phenomenon of parametric resonance. Besides, it is shown that the considered waves are mutually coupled as a consequence of the nonequilibrium state of the background medium. The coupling is strongest when the plasma $\beta \approx 1$. The proposed formalism is sufficiently general and can be applied for many dynamical systems, both under terrestrial and astrophysical conditions.

The dynamics of linear perturbations is studied in magnetized plasma shear flows with a constant shearing rate and with gravity-induced stratification. The general set of linearized equations is derived and the two-dimensional case is considered in detail. The Boussinesq approximation is used in order to examine relatively small-scale perturbations of low-frequency modes: Gravito-Alfven waves (GAW) and Entropy Mode (EM) perturbations. It is shown that for flows with arbitrary shearing rate there exists a finite time interval of non-adiabatic evolution of the perturbations. The non-adiabatic behavior manifests itself in a twofold way, viz. by the over-reflection of the GAWs and by the generation of GAWs from EM perturbations. It is shown that these phenomena act as efficient transformers of the equilibrium flow energy into the energy of the perturbations for moderate and high shearing rate solar plasma flows. Efficient generation of GAW by EM takes place for shearing rates about an order of magnitude smaller than necessary for development of a shear instability. The latter fact could have important consequences for the problem of angular momentum redistribution within the Sun and solar-type stars.

The possibility of velocity shear-induced linear transformations of different magnetohydrodynamic waves in the solar wind is studied both analytically and numerically. A quantitative analysis of the wave transformation processes for all possible plasma-$\beta$ regimes is performed. By applying the obtained criteria for effective wave coupling to the solar wind parameters, we show that velocity shear-induced linear transformations of Alfv\'en waves into magneto-acoustic waves could effectively take place for the relatively low-frequency Alfv\'en waves in the energy containing interval. The obtained results are in a good qualitative agreement with the observed features of density perturbations in the solar wind.

The observed power spectrum of high-degree solar p-modes (l>200) shows discrepancies with the power spectrum predicted by the stochastic excitement and damping theory. In an attempt to explain these discrepancies, the present paper is concerned with the influence of the observed subsurface flows on the trapped acoustic modes (p-modes). The effect of these inhomogeneous background flows is investigated by means of a non-modal analysis and a multi-layer model. It is shown that the rotational and meridional components of the velocity field change the wavelengths of the oscillation modes which, in turn, results in modifications of the corresponding modal frequencies. The magnitudes of the frequency residuals depend on the spatial scales of the modes and on the gradients of the different components of the flow velocity. Together with other mechanisms (e.g. the scattering of modes by the large scale convection (Goldreich & Murray 1994), the non-modal effect of the variation of the frequencies in time may contribute: 1) to the observed widening of the corresponding peaks in the observed power spectrum with increasing angular degree; 2) to the partial dissipation of spectral power, and, as a result, 3) to the discrepancies between the predicted and the observed power spectrum of solar p-modes.

The recently suggested swing interaction between fast magnetosonic and Alfv\'en waves (2002) is generalized to inhomogeneous media. We show that the fast magnetosonic waves propagating across an applied non-uniform magnetic field can parametrically amplify the Alfv\'en waves propagating along the field through the periodical variation of the Alfv\'en speed. The resonant Alfv\'en waves have half the frequency and the perpendicular velocity polarization of the fast waves. The wavelengths of the resonant waves have different values across the magnetic field, due to the inhomogeneity in the Alfv\'en speed. Therefore, if the medium is bounded along the magnetic field, then the harmonics of the Alfv\'en waves, which satisfy the condition for onset of a standing pattern, have stronger growth rates. In these regions the fast magnetosonic waves can be strongly 'absorbed', their energy going in transversal Alfv\'en waves. We refer to this phenomenon as 'Swing Absorption'. This mechanism can be of importance in various astrophysical situations.

Recently it was pointed out that nonmodally (transiently and/or adiabatically) pre-amplified waves in shear flows, undergoing subsequent viscous damping, can ultimately heat the ambient flow. The key ingredient of this process is the ability of waves to grow, by extracting energy from the spatially inhomogeneous mean flow. In this paper we examine this mechanism in the context of the solar coronal plasma flows. "Self-heating" (SH) processes are examined when both viscous damping and magnetic resistivity are at work. We show that if the plasma viscosity is in the favorable range of values the asymptotic SH rate in these flows can be quite substantial.

The recently suggested swing interaction between fast magnetosonic and Alfven waves (Zaqarashvili and Roberts 2002) is generalized to inhomogeneous media. We show that the fast magnetosonic waves propagating across an applied non-uniform magnetic field can parametrically amplify the Alfven waves propagating along the field through the periodical variation of the Alfven speed. The resonant Alfven waves have half the frequency and the perpendicular velocity polarization of the fast waves. The wave lengths of the resonant waves have different values across the magnetic field, due to the inhomogeneity in the Alfven speed. Therefore, if the medium is bounded along the magnetic field, then the harmonics of the Alfven waves, which satisfy the condition for onset of a standing pattern, have stronger growth rates. In these regions the fast magnetosonic waves can be strongly 'absorbed', their energy going in transversal Alfven waves. We refer to this phenomenon as 'Swing Absorption'. This mechanism can be of importance in various astrophysical situations.