We observe plasma flows in cool loops using the Slit-Jaw Imager (SJI) onboard the Interface Region Imaging Spectrometer (IRIS). Huang et al. (2015) observed unusually broadened Si IV 1403 angstrom line profiles at the footpoints of such loops that were attributed to signatures of explosive events (EEs). We have chosen one such uni-directional flowing cool loop system observed by IRIS where one of the footpoints is associated with significantly broadened Si IV line profiles. The line profile broadening indirectly indicates the occurrence of numerous EEs below the transition region (TR), while it directly infers a large velocity enhancement /perturbation further causing the plasma flows in the observed loop system. The observed features are implemented in a model atmosphere in which a low-lying bi-polar magnetic field system is perturbed in the chromosphere by a velocity pulse with a maximum amplitude of 200 km/s. The data-driven 2-D numerical simulation shows that the plasma motions evolve in a similar manner as observed by IRIS in the form of flowing plasma filling the skeleton of a cool loop system. We compare the spatio-temporal evolution of the cool loop system in the framework of our model with the observations, and conclude that their formation is mostly associated with the velocity response of the transient energy release above their footpoints in the chromosphere/TR. Our observations and modeling results suggest that the velocity responses most likely associated to the EEs could be one of the main candidates for the dynamics and energetics of the flowing cool loop systems in the lower solar atmosphere.

We numerically simulate emergence of a magnetic plasma blob into a solar coronal hole. This blob may be associated with granulation and therefore it has a weak magnetic field. Two-dimensional simulations are performed using the MAGNUS code which solves magnetohydrodynamic equations, taking into account magnetic resistivity and thermal conduction. As a result of the interaction of the emerging blob with the ambient plasma, the magnetic lines experience reconnection with the blob getting flattened and deformed with time. Additionally, this process launches a vertical outflow of hot plasma and the chromosphere in its response increases its temperature. We perform parametric studies by varying the magnitude of the magnetic field of the blob and observing the net heating of the chromosphere. These studies are inspired by realistic simulations of granulation made with the use of two-fluid JOANNA code. In these simulations a number of magnetic blobs are detected in the convection zone and in the photosphere. From the numerical results, we conclude that as a result of granulation operating in a solar quiet region the emerging blob may trigger very complex dynamics in the upper regions of the solar atmosphere, and the associated outflows may be a source of heating of the chromosphere and possibly the solar corona.

With the use of our JOANNA code, which solves radiative equations for ion + electron and neutral fluids, we perform realistic 2.5D numerical simulations of plasma outflows associated with the solar granulation. These outflows exhibit physical quantities consistent to the order of magnitude with the observational findings for mass and energy losses in the upper chromosphere, transition region and inner corona, and they may originate the fast solar wind.

Using high-resolution numerical simulations we investigate the plasma heating driven by periodic two-fluid acoustic waves that originate at the bottom of the photosphere and propagate into the gravitationally stratified and partially ionized solar atmosphere. We consider ions+electrons and neutrals as separate fluids that interact between themselves via collision forces. The latter play an important role in the chromosphere, leading to significant damping of short-period waves. Long-period waves do not essentially alter the photospheric temperatures, but they exhibit the capability of depositing a part of their energy in the chromosphere. This results in up about a five times increase of ion temperature that takes place there on a time-scale of a few minutes. The most effective heating corresponds to waveperiods within the range of about 30-200 s with a peak value located at 80 s. However, we conclude that for the amplitude of the driver chosen to be equal to 0.1 km s$^{-1}$, this heating is too low to balance the radiative losses in the chromosphere.

We investigate the Interface Region Imaging Spectrograph (IRIS) observations of the quiet-Sun (QS) to understand the propagation of acoustic waves in transition region (TR) from photosphere. We selected a few IRIS spectral lines, which include the photospheric (Mn~{\sc i} 2801.25~{\AA}), chromospheric (Mg~{\sc ii} k 2796.35~{\AA}) and TR (C~{\sc ii} 1334.53~{\AA}), to investigate the acoustic wave propagation.The wavelet cross-spectrum reveals significant coherence (about 70\% locations) between photosphere and chromosphere. Few minutes oscillations (i.e., period range from 1.6 to 4.0 minutes) successfully propagate into chromosphere from photosphere, which is confirmed by dominance of positive phase lags. However, in higher period regime (i.e., greater than $\approx$ 4.5 minutes), the downward propagation dominates is evident by negative phase lags. The broad spectrum of waves (i.e., 2.5-6.0 minutes) propagates freely upwards from chromosphere to TR. We find that only about 45\% locations (out of 70\%) show correlation between chromosphere and TR. Our results indicate that roots of 3 minutes oscillations observed within chromosphere/TR are located in photosphere. Observations also demonstrate that 5 minute oscillations propagate downward from chromosphere. \textbf{However, some locations within QS also show successful propagation of 5 minute oscillations as revealed by positive phase lags, which might be the result of magnetic field}. In addition, our results clearly show that a significant power, within period ranging from 2.5 to 6.0 minutes, of solar chromosphere is freely transmitted into TR triggering atmospheric oscillations. Theoretical implications of our observational results are discussed.

Aims. We perform a detailed observational analysis of network jets to understand their kinematics, rotational motion and underlying triggering mechanism(s). We have analyzed the quiet-Sun (QS) data. Methods. IRIS high resolution imaging and spectral observations (SJI: Si iv 1400.0 \AA, Raster: Si iv 1393.75 \AA) are used to analyze the omnipresent rotating network jets in the transition-region (TR). In addition, we have also used Atmospheric Imaging Assembly (AIA) onboard Solar Dynamic Observation (SDO) observations. Results. The statistical analysis of fifty-one network jets is performed to understand various their mean properties, e.g., apparent speed (140.16+/-39.41 km/s), length (3.16+/-1.18 Mm), lifetimes (105.49+/-51.75 s). The Si iv 1393.75 \AA line has secondary component along with its main Gaussian, which is formed due to the high-speed plasma flows (i.e., network jets). The variation of Doppler velocity across these jets (i.e., blue shift on one edge and red shift on the other) signify the presence of inherited rotational motion. The statistical analysis predicts that the mean rotational velocity (i.e., \delV) is 49.56 km/s. The network jets have high angular velocity in comparison to the other class of solar jets. Conclusions. The signature of network jets are inherited in TR spectral lines in terms of the secondary component of the Si iv 1393.75 \AA line. The rotational motion of network jets is omnipresent, which is reported firstly for this class of jet-like features. The magnetic reconnection seems to be the most favorable mechanism for the formation of these network jets.

We aim to study the formation and evolution of solar spicules by means of numerical simulations of the solar atmosphere. With the use of newly developed JOANNA code, we numerically solve two-fluid (for ions + electrons and neutrals) equations in 2D Cartesian geometry. We follow the evolution of a spicule triggered by the time-dependent signal in ion and neutral components of gas pressure launched in the upper chromosphere. We use the potential magnetic field, which evolves self-consistently, but mainly plays a passive role in the dynamics. Our numerical results reveal that the signal is steepened into a shock that propagates upward into the corona. The chromospheric cold and dense plasma lags behind this shock and rises into the corona with a mean speed of 20-25 km s$^{-1}$. The formed spicule exhibits the upflow/downfall of plasma during its total lifetime of around 3-4 minutes, and it follows the typical characteristics of a classical spicule, which is modeled by magnetohydrodynamics. The simulated spicule consists of a dense and cold core that is dominated by neutrals. The general dynamics of ion and neutral spicules are very similar to each other. Minor differences in those dynamics result in different widths of both spicules with increasing rarefaction of the ion spicule in time.

We perform numerical simulations of impulsively generated magnetic swirls in an isolated flux tube which is rooted in the solar photosphere. These swirls are triggered by an initial pulse in a horizontal component of the velocity. The initial pulse is launched either: (a) centrally, within the localized magnetic flux tube; or (b) off-central, in the ambient medium. The evolution and dynamics of the flux tube is described by three-dimensional, ideal magnetohydrodynamic equations. These equations are numerically solved to reveal that in case (a) dipole-like swirls associated with the fast magnetoacoustic kink and $m=1$ Alfv\'en waves are generated. In case (b), the fast magnetoacoustic kink and $m=0$ Alfv\'en modes are excited. In both these cases, the excited fast magnetoacoustic kink and Alfv\'en waves consist of similar flow pattern and magnetic shells are also generated with clockwise and counter-clockwise rotating plasma within them, which can be the proxy of dipole-shaped chromospheric swirls. The complex dynamics of vortices and wave perturbations reveals the channelling of sufficient amount of energy to fulfill energy losses in the chromosphere ($\sim$ 10$^{4}$ W m$^{-1}$) and in the corona ($\sim$ 10$^{2}$ W m$^{-1}$). Some of these numerical findings are reminiscent of signatures in recent observational data.

We investigate numerically Alfv\'en waves propagating along an axisymmetric and non-isothermal solar flux tube embedded in the solar atmosphere. The tube magnetic field is current-free and diverges with height, and the waves are excited by a periodic driver along the tube magnetic field lines. The main results are that the two wave variables, the velocity and magnetic field perturbations in the azimuthal direction, behave differently as a result of gradients of physical parameters along the tube. To explain these differences in the wave behavior, the time evolution of the wave variables and the resulting cutoff period for each wave variable are calculated, and used to determine regions in the solar chromosphere where strong wave reflection may occur.

Recently Wi\'sniewska et al. demonstrated observationally how the acoustic cutoff frequency varies with height in the solar atmosphere including the upper photosphere and the lower and middle chromosphere, and showed that the observational results cannot be accounted for by the existing theoretical formulas for the acoustic cutoff. In order to reproduce the observed variation of the cutoff with atmospheric height, numerical simulations of impulsively generated acoustic waves in the solar atmosphere are performed, and the spectral analysis of temporal wave profiles is used to compute numerically changes of the acoustic cutoff with height. Comparison of the numerical results with the observational data shows good agreement, which clearly indicates that the obtained results may be used to determine the structure of the background solar atmosphere.

We present analytical and numerical models of a normal-polarity quiescent prominence that are based on the model of Pikelner (Solar Phys. 1971, 17, 44 ). We derive the general analytical expressions for the two-dimensional equilibrium plasma quantities such as the mass density and a gas pressure, and we specify magnetic-field components for the prominence, which corresponds to a dense and cold plasma residing in the dip of curved magnetic-field lines. With the adaptation of these expressions, we solve numerically the 2D, nonlinear, ideal MHD equations for a Pikelner's model of a prominence that is initially perturbed by reducing the gas pressure at the dip of magnetic-field lines. Our findings reveal that as a result of pressure perturbations the prominence plasma starts evolving in time and this leads to the antisymmetric magnetoacoustic--gravity oscillations as well as to the mass-density growth at the magnetic dip, and the magnetic-field lines subside there. This growth depends on the depth of magnetic dip. For a shallower dip, less plasma is condensed and vice-versa. We conjecture that the observed long-period magnetoacoustic-gravity oscillations in various prominence systems are in general the consequence of the internal pressure perturbations of the plasma residing in equilibrium at the prominence dip.

Without any doubt solar flaring loops possess a multi-thread internal structure that is poorly resolved and there are no means to observe heating episodes and thermodynamic evolution of the individual threads. These limitations cause fundamental problems in numerical modelling of flaring loops, such as selection of a structure and a number of threads, and an implementation of a proper model of the energy deposition process. A set of 1D hydrodynamic and 2D magnetohydrodynamic models of a flaring loop are developed to compare energy redistribution and plasma dynamics in the course of a prototypical solar flare. Basic parameters of the modeled loop are set according to the progenitor M1.8 flare recorded in the AR10126 on September 20, 2002 between 09:21 UT and 09:50 UT. The non-ideal 1D models include thermal conduction and radiative losses of the optically thin plasma as energy loss mechanisms, while the non-ideal 2D models take into account viscosity and thermal conduction as energy loss mechanisms only. The 2D models have a continuous distribution of the parameters of the plasma across the loop, and are powered by varying in time and space along and across the loop heating flux. We show that such 2D models are a borderline case of a multi-thread internal structure of the flaring loop, with a filling factor equal to one. Despite the assumptions used in applied 2D models, their overall success in replicating the observations suggests that they can be adopted as a correct approximation of the observed flaring structures.

Using the Hinode/EIS 2$"$ spectroscopic observations, we study the intensity, velocity, and FWHM variations of the strongest Fe XII 195.12 \AA\ line along the jet to find the signature of Alfv\'en waves. We simulate numerically the impulsively generated Alfv\'en waves within the vertical Harris current-sheet, forming the jet plasma flows, and mimicking their observational signatures. Using the FLASH code and the atmospheric model with embedded weakly expanding magnetic field configuration within a vertical Harris current-sheet, we solve the two and half-dimensional (2.5-D) ideal magnetohydrodynamic (MHD) equations to study the evolution of Alfv\'en waves and vertical flows forming the plasma jet. At a height of $\sim 5~\mathrm{Mm}$ from the base of the jet, the red-shifted velocity component of Fe XII 195.12 \AA\ line attains its maximum ($5~\mathrm{km\,s}^{-1}$) which converts into a blue-shifted one between the altitude of $5-10~\mathrm{Mm}$. The spectral intensity continously increases up to $10~\mathrm{Mm}$, while FWHM still exhibits the low values with almost constant trend. This indicates that the reconnection point within the jet's magnetic field topology lies in the corona $5-10~\mathrm{Mm}$ from its footpoint anchored in the Sun's surface. Beyond this height, FWHM shows a growing trend. This may be the signature of Alfv\'en waves that impulsively evolve due to reconnection and propagate along the jet. From our numerical data, we evaluate space- and time- averaged Alfv\'en waves velocity amplitudes at different heights in the jet's current-sheet, which contribute to the non-thermal motions and spectral line broadening. The synthetic width of Fe XII $195.12~\mathrm{\AA}$ line exhibits similar trend of increment as in the observational data, possibly proving the existence of impulsively generated (by reconnection) Alfv\'en waves which propagate along the jet.

We perform numerical simulations of impulsively generated Alfv\'en waves in an isolated photospheric flux tube, and explore the propagation of these waves along such magnetic structure that extends from the photosphere, where these waves are triggered, to the solar corona, and analyze resulting magnetic shells. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical model and specifying the curved magnetic field lines that constitute the magnetic flux tube which is rooted in the solar photosphere. The evolution of the solar atmosphere is described by 3D, ideal magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal, based on the physical properties of the multi-shell magnetic twisters and the amount of energy and momentum associated with them, that these multi-shell magnetic twisters may be responsible for the observed heating of the lower solar corona and for the formation of solar wind. Moreover, it is likely that the existence of these twisters can be verified by high-resolution observations.

Aims: Propagation and energy transfer of torsional Alfv\'en waves in solar magnetic flux tubes of axial symmetry is studied. Methods: An analytical model of a solar magnetic flux tube of axial symmetry is developed by specifying a magnetic flux and deriving general analytical formulae for the equilibrium mass density and a gas pressure. The main advantage of this model is that it can be easily adopted to any axisymmetric magnetic structure. The model is used to simulate numerically the propagation of nonlinear Alfv\'en waves in such 2D flux tubes of axial symmetry embedded in the solar atmosphere. The waves are excited by a localized pulse in the azimuthal component of velocity and launched at the top of the solar photosphere, and they propagate through the solar chromosphere, transition region, and into the solar corona. Results: The results of our numerical simulations reveal a complex scenario of twisted magnetic field lines and flows associated with torsional Alfv\'en waves as well as energy transfer to the magnetoacoustic waves that are triggered by the Alfv\'en waves and are akin to the vertical jet flows. Alfv\'en waves experience about 5 % amplitude reflection at the transition region. Magnetic (velocity) field perturbations experience attenuation (growth) with height is agreement with analytical findings. Kinetic energy of magnetoacoustic waves consists of 25 % of the total energy of Alfv\'en waves. The energy transfer may lead to localized mass transport in the form of vertical jets, as well as to localized heating as slow magnetoacoustic waves are prone to dissipation in the inner corona.

Aims. We construct a new analytical model of a solar coronal loop that is embedded in a gravitationally stratified and magnetically confined atmosphere. On the basis of this analytical model, we devise a numerical model of solar coronal loops. We adopt it to perform the numerical simulations of its vertical kink oscillations excited by an external driver. Methods. Our model of the solar atmosphere is constructed by adopting a realistic temperature distribution and specifying the curved magnetic field lines that constitute a coronal loop. This loop is described by 2D, ideal magnetohydro- dynamic equations that are numerically solved by the FLASH code. Results. The vertical kink oscillations are excited by a periodic driver in the vertical component of velocity, acting at the top of the photosphere. For this forced driver with its amplitude 3 km/s, the excited oscillations exhibit about 1.2 km/s amplitude in their velocity and the loop apex oscillates with its amplitude in displacement of about 100 km. Conclusions. The newly devised analytical model of the coronal loops is utilized for the numerical simulations of the vertical kink oscillations, which match well with the recent observations of decay-less kink oscillations excited in solar loops. The model will have further implications on the study of waves and plasma dynamics in coronal loops, revealing physics of energy and mass transport mechanisms in the localized solar atmosphere.

We perform numerical simulations of impulsively generated Alfv\'en waves in an isolated solar arcade, which is gravitationally stratified and magnetically confined. We study numerically the propagation of Alfv\'en waves along such magnetic structure that extends from the lower chromosphere, where the waves are generated, to the solar corona, and analyze influence of the arcade size and width of the initial pulses on the wave propagation and reflection. Our model of the solar atmosphere is constructed by adopting the temperature distribution based on the semi-empirical VAL-C model and specifying the curved magnetic field lines that constitute the asymmetric magnetic arcade. The propagation and reflection of Alfv\'en waves in this arcade is described by 2.5D magnetohydrodynamic equations that are numerically solved by the FLASH code. Our numerical simulations reveal that the Alfv\'en wave amplitude decreases as a result of a partial reflection of Alfv\'en waves in the solar transition region, and that the waves which are not reflected leak through the transition region and reach the solar corona. We also find the decrement of the attenuation time of Alfv\'en waves for wider initial pulses. Moreover, our results show that the propagation of Alfv\'en waves in the arcade is affected by spatial dependence of the Alfv\'en speed, which leads to phase-mixing that is stronger for more curved and larger magnetic arcades. We discuss processes that affect the Alfv\'en wave propagation in an asymmetric solar arcade and conclude that besides phase-mixing in the magnetic field configuration, plasma properties of the arcade and size of the initial pulse as well as structure of the solar transition region all play a vital role in the Alfv\'en wave propagation.

Using one-arcsecond-slit scan observations from the Hinode/EUV Imaging Spectrometer (EIS) on 05 February 2007, we find the plasma outflows in the open and expanding coronal funnels at the eastern boundary of AR 10940. The Doppler velocity map of Fe XII 195.120 A shows that the diffuse close-loop system to be mostly red-shifted. The open arches (funnels) at the eastern boundary of AR exhibit blue-shifts with a maximum speed of about 10-15 km/s. This implies outflowing plasma through these magnetic structures. In support of these observations, we perform a 2D numerical simulation of the expanding coronal funnels by solving the set of ideal MHD equations in appropriate VAL-III C initial temperature conditions using the FLASH code. We implement a rarefied and hotter region at the footpoint of the model funnel, which results in the evolution of slow plasma perturbations propagating outward in the form of plasma flows. We conclude that the heating, which may result from magnetic reconnection, can trigger the observed plasma outflows in such coronal funnels. This can transport mass into the higher corona, giving rise to the formation of the nascent solar wind.

We present results of 3-D numerical simulations of a fast magnetic twister excited above a foot-point of the potential solar coronal arcade that is embedded in the solar atmosphere with the initial VAL-IIIC temperature profile, which is smoothly extended into the solar corona. With the use of the FLASH code, we solve 3-D ideal magnetohydrodynamic equations by specifying a twist in the azimuthal component of magnetic field in the solar chromosphere. The imposed perturbation generates torsional Alfv\'en waves as well as plasma swirls that reach the other foot-point of the arcade and partially reflect back from the transition region. The two vortex channels are evident in the generated twisted flux-tube with a fragmentation near its apex that results from the initial twist as well as from the morphology of the tube. The numerical results are compared to observational data of plasma motions in a solar prominence. The comparison shows that the numerical results and the data qualitatively agree even though the observed plasma motions occur over comparatively large spatio-temporal scales in the prominence.

We present simulation results of the impulsively generated linear and non-linear Alfv\'en waves in the weakly curved coronal magnetic flux-tubes (coronal funnels) and discuss their implications for the coronal heating and solar wind acceleration. We solve numerically the time-dependent magnetohydrodynamic (MHD) equations to obtain the temporal signatures of the small (linear) and large-amplitude (non-linear) Alfv\'en waves in the model atmosphere of expanding open magnetic field configuration (e.g., coronal funnels) by considering a realistic temperature distribution. We compute the maximum transversal velocity of both linear and non-linear Alfv\'en waves at different heights in the coronal funnel, and study their response in the solar corona during the time of their propagation. We infer that the pulse-driven non-linear Alfv\'en waves may carry sufficient wave energy fluxes to heat the coronal funnels and also to power the solar wind that originates in these funnels. Our study of linear Alfv\'en waves show that they can contribute only to the plasma dynamics and heating of the funnel-like magnetic flux-tubes associated with the polar coronal holes.

With the adaptation of the FLASH code we simulate magnetohydrodynamic-gravity waves and vortices as well as their response in the magnetized three-dimensional (3D) solar atmosphere at different heights to understand the localized energy transport processes. In the solar atmosphere strongly structured by gravitational and magnetic forces, we launch a localized velocity pulse (in horizontal and vertical components) within a bottom layer of 3D solar atmosphere modelled by initial VAL-IIIC conditions, which triggers waves and vortices. The rotation direction of vortices depends on the orientation of an initial perturbation. The vertical driver generates magnetoacoustic-gravity waves which result in oscillations of the transition region, and it leads to the eddies with their symmetry axis oriented vertically. The horizontal pulse excites all magnetohydrodynamic-gravity waves and horizontally oriented eddies. These waves propagate upwards, penetrate the transition region, and enter the solar corona. In the high-beta plasma regions the magnetic field lines move with the plasma and the temporal evolution show that they swirl with eddies. We estimate the energy fluxes carried out by the waves in the magnetized solar atmosphere and conclude that such wave dynamics and vortices may be significant in transporting the energy to sufficiently balance the energy losses in the localized corona. Moreover, the structure of the transition region highly affects such energy transports, and causes the channelling of the propagating waves into the inner corona.

We report an observation of a small scale flux-tube that undergoes kinking and triggers the macrospicule and a jet on November 11, 2010 in the north polar corona. The small-scale flux-tube emerged well before the triggering of macrospicule and as the time progresses the two opposite halves of this omega shaped flux-tube bent transversely and approached towards each other. After $\sim$ 2 minutes, the two approaching halves of the kinked flux-tube touch each-other and internal reconnection as well as energy release takes place at the adjoining location and a macrospicule was launched which goes upto a height of 12 Mm. Plasma starts moving horizontally as well as vertically upward along with the onset of macrospicule and thereafter converts into a large-scale jet which goes up to $\sim$ 40 Mm in the solar atmosphere with a projected speed of $\sim$ 95 km s$^{-1}$. We perform 2-D numerical simulation by considering the VAL-C initial atmospheric conditions to understand the physical scenario of the observed macrospicule and associated jet. The simulation results show that reconnection generated velocity pulse in the lower solar atmosphere steepens into slow shock and the cool plasma is driven behind it in form of macrospicule. The horizontal surface waves also appeared with the shock fronts at different heights, which most likely drove and spread the large-scale jet associated with the macrospicule.

We observe a solar surge in NOAA AR11271 using SDO AIA 304 image data on 25 August, 2011. The surge rises vertically from its origin upto a height of 65 Mm with a terminal velocity of 100 km/s, and thereafter falls and fades gradually. The total life time of the surge was 20 min. We also measure the temperature and density distribution of the observed surge during its maximum rise, and found an average temperature and density of 2.0 MK and 4.1 x 109 cm-3, respectively. The temperature map shows the expansion and mixing of cool plasma lagging behind the hot coronal plasma along the surge. As SDO/HMI temporal image data does not show any detectable evidence of the significant photospheric magnetic field cancellation for the formation of the observed surge, we infer that it is probably driven by magnetic reconnection generated thermal energy in the lower chromosphere. The radiance (thus mass density) oscillations near the base of the surge are also evident, which may be the most likely signature of its formation by a reconnection generated pulse. In support of the present observational base-line of the triggering of the surge due to chromospheric heating, we devise a numerical model with conceivable implementation of VAL-C atmosphere and a thermal pulse as an initial trigger. We find that the pulse steepens into a slow shock at higher altitudes that triggers plasma perturbations exhibiting the observed features of the surge, e.g., terminal velocity, height, width, life-time, and heated fine structures near its base.

We investigate the excitation of magnetoacoustic-gravity waves generated from localized pulses in the gas pressure as well as in vertical component of velocity. These pulses are initially launched at the top of the solar photosphere that is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical and oblique to the gravitational force. We numerically model magnetoacoustic-gravity waves by implementing a realistic (VAL-C) model of solar temperature. We solve two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic-gravity waves, which are later reflected from the transition region due to the large temperature gradient. The atmospheric cavities in the lower solar atmosphere are found to be the ideal places that may act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics.

We study the impulsively generated non-linear Alfv\'en waves in the solar atmosphere, and describe their most likely role in the observed non-thermal broadening of some spectral lines in solar coronal holes. We solve numerically the time-dependent magnetohydrodynamic equations to find temporal signatures of large-amplitude Alfv\'en waves in the model atmosphere of open and expanding magnetic field configuration, with a realistic temperature distribution. We calculate the temporally and spatially averaged, instantaneous transversal velocity of non-linear Alfv\'en waves at different heights of the model atmosphere, and estimate its contribution to the unresolved non-thermal motions caused by the waves. We find that the pulse-driven nonlinear Alfv\'en waves with the amplitude $A_{\rm v}$=50 km s$^{-1}$ are the most likely candidates for the non-thermal broadening of Si VIII $\lambda$1445.75 \AA\ line profiles in the polar coronal hole as reported by Banerjee et al. (1998). We also demonstrate that the Alfv\'en waves driven by comparatively smaller velocity pulse with its amplitude $A_{\rm v}$=25 km s$^{-1}$ may contribute to the spectral line width of the same line at various heights in coronal hole without any significant broadening. The main conclusion of this paper is that non-linear Alfv\'en waves excited impulsively in the lower solar atmosphere are responsible for the observed spectral line broadening in polar coronal holes. This is an important result as it allows us to conclude that such large amplitude and pulse-driven Alfv\'en waves do indeed exist in solar coronal holes. The existence of these waves and their undamped growth may impart the required momentum to accelerate the solar wind.

Using multi-wavelength observations of Solar and Heliospheric Observatory (SoHO)/Michelson Doppler Imager (MDI), Transition Region and Coronal Explorer (TRACE) 171 \AA, and H$\alpha$ from Culgoora Solar Observatory at Narrabri, Australia, we present a unique observational signature of a propagating supersonic plasma blob before an M6.2 class solar flare in AR10808 on 9th September 2005. The blob was observed between 05:27 UT to 05:32 UT with almost a constant shape for the first 2-3 minutes, and thereafter it quickly vanished in the corona. The observed lower bound speed of the blob is estimated as $\sim$215 km s$^{-1}$ in its dynamical phase. The evidence of the blob with almost similar shape and velocity concurrent in H$\alpha$ and TRACE 171 \AA\ supports its formation by multi-temperature plasma. The energy release by a recurrent 3-D reconnection process via the separator dome below the magnetic null point, between the emerging flux and pre-existing field lines in the lower solar atmosphere, is found to be the driver of a radial velocity pulse outwards that accelerates this plasma blob in the solar atmosphere. In support of identification of the possible driver of the observed eruption, we solve the two-dimensional ideal magnetohydrodynamic equations numerically to simulate the observed supersonic plasma blob. The numerical modelling closely match the observed velocity, evolution of multi-temperature plasma, and quick vanishing of the blob found in the observations. Under typical coronal conditions, such blobs may also carry an energy flux of 7.0$\times10^{6}$ ergs cm$^{-2}$ s$^{-1}$ to re-balance the coronal losses above active regions.

The P1/P2 period ratio of transversal loop oscillations is currently used for the diagnostics of longitudinal structuring of coronal loops as its deviation from 2 is intrinsically connected to the density scale-height along coronal loops and/or the sub-resolution structure of the magnetic field. The same technique can be applied not only to coronal structures, but also to other oscillating magnetic structures. The oscillations in magnetic structures are described by differential equations whose coefficients depend on the longitudinal structure of the plasma. Using a variational principle written for the transversal component of the velocity vector, developed earlier by McEwan et al. (2008), we investigate how the different temperature of the environment compared to the temperature of the magnetic structure will influence the P1/P2 ratio for typical coronal and prominence conditions. The possible changes are translated into quantities that are used in the process of remote plasma diagnostics in the solar atmosphere.

We observe the motion of cool and hot plasma in a multi-stranded post flare loop system that evolved in the decay phase of a two ribbon M1.0 class flare in AR 11093 on 7 August 2010 using SDO/AIA 304 \AA\ and 171 \AA\ filters. The moving intensity feature and its reflected counterpart are observed in the loop system at multi-temperature. The observed hot counterpart of the plasma that probably envelopes the cool confined plasma, moves comparatively faster ($\sim$34 km s$^{-1}$) to the later (29 km s$^{-1}$) in form of the spreaded intensity feature. The propagating plasma and intensity reflect from the region of another footpoint of the loop. The subsonic speed of the moving plasma and associated intensity feature may be most likely evolved in the post flare loop system through impulsive flare heating processes. Complementing our observations of moving multi-temperature intensity features in the post flare loop system and its reflection, we also attempt to solve two-dimensional ideal magnetohydrodynamic equations numerically using the VAL-IIIC atmosphere as an initial condition to simulate the observed plasma dynamics. We consider a localized thermal pulse impulsively generated near one footpoint of the loop system during the flare processes, which is launched along the magnetic field lines at the solar chromosphere. The pulse steepens into a slow shock at higher altitudes while moving along this loop system, which triggers plasma perturbations that closely exhibit the observed plasma dynamics.

Context. We consider a localized pulse in the component of velocity, parallel to the ambient magnetic field lines, that is initially launched in the solar chromosphere. Aims. We aim to generalize the recent numerical model of spicule formation (Murawski & Zaqarashvili 2010) by implementing a VAL-C model of solar temperature. Methods. With the use of the code FLASH we solve two-dimensional ideal magnetohydrodynamic equations numerically to simulate the solar macrospicules. Results. Our numerical results reveal that the pulse located below the transition region triggers plasma perturbations, which exhibit many features of macrospicules. We also present an observational (SDO/AIA 304 A) case study of the macrospicule that approximately mimics the numerical simulations. Conclusions. In the frame of the model we devised, the solar macrospicules can be triggered by velocity pulses launched from the chromosphere.

Context. We observe a solar jet at north polar coronal hole (NPCH) using SDO AIA 304 {\deg}A image data on 3 August 2010. The jet rises obliquely above the solar limb and then retraces its propagation path to fall back. Aims. We numerically model this observed solar jet by implementing a realistic (VAL-C) model of solar temperature. Methods. We solve two-dimensional ideal magnetohydrodynamic equations numerically to simulate the observed solar jet. We consider a localized velocity pulse that is essentially parallel to the background magnetic field lines and initially launched at the top of the solar photosphere. The pulse steepens into a shock at higher altitudes, which triggers plasma perturbations that exhibit the observed features of the jet. The typical direction of the pulse also clearly exhibits the leading front of the moving jet. Results. Our numerical simulations reveal that a large amplitude initial velocity pulse launched at the top of the solar photosphere produces in general the observed properties of the jet, e.g., upward and backward average velocities, height, width, life-time, and ballistic nature. Conclusions. The close matching between the jet observations and numerical simulations provides first strong evidence for the formation of this jet by a single velocity pulse. The strong velocity pulse is most likely generated by the low- atmospheric reconnection in the polar region which results in triggering of the jet. The downflowing material of the jet most likely vanishes in the next upcoming velocity pulses from lower solar atmosphere, and therefore distinctly launched a single jet upward in the solar atmosphere is observed.

We aim to study excitation of the observed 5-min oscillations in the solar corona by localized pulses that are launched in the photosphere. We solve the full set of nonlinear one-dimensional Euler equations numerically for the velocity pulse propagating in the solar atmosphere that is determined by the realistic temperature profile. Numerical simulations show that an initial velocity pulse quickly steepens into a leading shock, while the nonlinear wake in the chromosphere leads to the formation of consecutive pulses. The time interval between arrivals of two neighboring pulses to a detection point in the corona is approximately 5 min. Therefore, the consecutive pulses may result in the 5-min oscillations that are observed in the solar corona. The 5-min oscillations observed in the solar corona can be explained in terms of consecutive shocks that result from impulsive triggers launched within the solar photosphere by granulation and/or reconnection.

We study the upward propagation of a localized velocity pulse that is initially launched below the transition region within the solar atmosphere. The pulse quickly steepens into a shock, which may lead to the formation of spicules. We solve two-dimensional time-dependent magnetohydrodynamic equations numerically to find spatial and temporal dynamics of spicules. The numerical simulations show that the strong initial pulse may lead to the quasi periodic rising of chromospheric material into the lower corona in the form of spicules. The periodicity results from the nonlinear wake that is formed behind the pulse in the stratified atmosphere. The superposition of raising and falling off plasma portions resembles the time sequence of single and double (sometimes even triple) spicules, which is consistent with observational findings. The two-dimensional rebound shock model may explain the observed speed, width, and heights of type I spicules, as well as observed multi-structural and bi-directional flows. The model also predicts the appearance of spicules with 3-5 min period due to the consecutive shocks.

We explore the effect of an inhomogeneous mass density field on frequencies and wave profiles of torsional Alfven oscillations in solar coronal loops. Dispersion relations for torsional oscillations are derived analytically in limits of weak and strong inhomogeneities. These analytical results are verified by numerical solutions, which are valid for a wide range of inhomogeneity strength. It is shown that the inhomogeneous mass density field leads to the reduction of a wave frequency of torsional oscillations, in comparison to that of estimated from mass density at the loop apex. This frequency reduction results from the decrease of an average Alfven speed as far as the inhomogeneous loop is denser at its footpoints. The derived dispersion relations and wave profiles are important for potential observations of torsional oscillations which result in periodic variations of spectral line widths. Torsional oscillations offer an additional powerful tool for a development of coronal seismology.