We present the Solar Bayesian Analysis Toolkit (SoBAT) which is a new easy to use tool for Bayesian analysis of observational data, including parameter inference and model comparison. SoBAT is aimed (but not limited) to be used for the analysis of solar observational data. We describe a new Interactive Data Language (IDL) code designed to facilitate the comparison of user-supplied model with data. Bayesian inference allows prior information to be taken into account. The use of Markov chain Monte Carlo (MCMC) sampling allows efficient exploration of large parameter spaces and provides reliable estimation of model parameters and their uncertainties. The Bayesian evidence for different models can be used for quantitative comparison. The code is tested to demonstrate its ability to accurately recover a variety of parameter probability distributions. Its application to practical problems is demonstrated using studies of the structure and oscillation of coronal loops.

Negative energy wave phenomena may appear in shear flows in the presence of a wave decay mechanism and external energy supply. We study the appearance of negative energy surface waves in a plasma cylinder in the incompressible limit. The cylinder is surrounded by an axial magnetic field and by a plasma of different density. Considering flow inside and viscosity outside the flux tube, we derive dispersion relations, and obtain analytical solutions for the phase speed and growth rate (increment) of the waves. It is found that the critical speed shear for the occurrence of the dissipative instability associated with negative energy waves (NEWs) and the threshold of Kelvin--Helmholtz instability (KHI) depend on the axial wavelength. The critical shear for the appearance of sausage NEW is lowest for the longest axial wavelengths, while for kink waves the minimum value of the critical shear is reached for the axial wavelength comparable to the diameter of the cylinder. The range between the critical speed of the dissipative instability and the KHI threshold is shown to depend on the difference of the Alfv\'{e}n speeds inside and outside of the cylinder. For all axial wavenumbers, NEW appears for the shear flow speeds lower than the KHI threshold. It is easier to excite NEW in an underdense cylinder than in an overdense one. The negative energy surface waves can be effectively generated for azimuthal number $m=0$ with a large axial wave number and for higher modes ($m>0$) with a small axial wave number.

Using full three-dimensional magnetohydrodynamic numerical simulations, we study the effects of magnetic field sigmoidity or helicity on the properties of the fundamental kink oscillation of solar coronal loops. Our model consists of a single denser coronal loop, embedded in a plasma with dipolar force-free magnetic field with a constant alpha-parameter. For the loop with no sigmoidity, we find that the numerically determined oscillation period of the fundamental kink mode matches the theoretical period calculated using WKB theory. In contrast, with increasing sigmoidity of the loop, the actual period is increasingly smaller than the one estimated by WKB theory. Translated through coronal seismology, increasing sigmoidity results in magnetic field estimates which are increasingly shifting towards higher values, and even surpassing the average value for the highest alpha value considered. Nevertheless, the estimated range of the coronal magnetic field value lies within the mimimal/maximal limits, proving the robustness coronal seismology. We propose that the discrepancy in the estimations of the absolute value of the force-free magnetic field could be exploited seismologically to determine the free energy of coronal loops, if averages of the internal magnetic field and density can be reliably estimated by other methods.

A light bridge is a prominent structure commonly observed within a sunspot. Its presence usually triggers a wealth of dynamics in a sunspot, and has a lasting impact on sunspot evolution. However, the fundamental structure of light bridges is still not well understood. In this study, we used the high-resolution spectropolarimetry data obtained by the Solar Optical Telescope onboard the Hinode satellite to analyze the magnetic and thermal structure of a light bridge at \AR. We also combined the high-cadence $1700\unit{\AA}$ channel data provided by the Atmospheric Imaging Assembly onboard the Solar Dynamic Observatory to study the dynamics on this bridge. We found that a pair of blue and red Doppler shift patches at two ends of this bridge, this pattern appears to be the convective motion directed by the horizontal component of the magnetic field aligned with the spine of this bridge. Paired upward and downward motions implies that the light bridge could have a two-legged or undulate magnetic field. Significant four minute oscillations in the emission intensity of the $1700\unit{\AA}$ bandpass were detected at two ends, which had overlap with the paired blue and red shift patches. The oscillatory signals at the light bridge and the penumbra were highly correlated with each other. Although they are separated in space at the photosphere, the periodicity seems to have a common origin from the underneath. Therefore, we infer that the light bridge and penumbra could share a common magnetic source and become fragmented at the photosphere by magneto-convection.

The relatively large-amplitude decaying regime of transverse oscillations of coronal loops has been known for two decades and interpreted in terms of MHD kink modes of cylindrical plasma waveguides. Recent observational analysis has revealed the so-called decay-less small-amplitude oscillations, with a multi-harmonic structure being detected. Several models have been proposed to explain them. In particular, decay-less oscillations have been described in terms of standing kink waves driven with continuous monoperiodic motions of loop footpoints, in terms of a simple oscillator model of forced oscillations due to harmonic external force, and as a self-oscillatory process due to the interaction of a loop with quasi-steady flows. However, an alternative mechanism is needed to explain the simultaneous excitation of several longitudinal harmonics of the oscillation. We study the mechanism of random excitation of decay-less transverse oscillations of coronal loops. With a spatially one-dimensional and time-dependent analytical model taking into account effects of the wave damping and kink speed variation along the loop, we consider transverse loop oscillations driven by random motions of footpoints. The footpoint motions are modelled by broad-band coloured noise. We have found the excitation of loop eigenmodes and analysed their frequency ratios as well as the spatial structure of the oscillations along the loop. The obtained results successfully reproduce the observed properties of decay-less oscillations. In particular, excitation of eigenmodes of a loop as a resonator can explain the observed quasi-monochromatic nature of decay-less oscillations and generation of multiple harmonics detected recently. We propose the mechanism that can interpret decay-less transverse oscillations of coronal loops in terms of kink waves randomly driven at the loop footpoints.

We investigate quasi-periodic pulsations (QPPs) of high-energy nonthermal emissions from an X9.3 flare (SOL2017-Sep-06T11:53), the most powerful flare since the beginning of solar cycle 24. The QPPs are identified as a series of regular and repeating peaks in the light curves in the gamma- and hard X-ray (HXR) channels recorded by the Konus-Wind, as well as the radio and microwave fluxes measured by the CALLISTO radio spectrograph during the impulsive phase. The periods are determined from the global wavelet and Fourier power spectra, as 24-30 s in the HXR and microwave channels which are associated with nonthermal electrons, and ~20 s in the gamma-ray band related to nonthermal ions. Both nonthermal electrons and ions may be accelerated by repetitive magnetic reconnection during the impulsive phase. However, we could not rule out other mechanisms such as the MHD oscillation in a sausage mode. The QPP detected in this study is useful for understanding the particle acceleration and dynamic process in solar flares and also bridging the gap between stellar and solar flares since the energy realm of the X9.3 solar flare is almost compared with a typical stellar flare.

A new type of chromospheric jets in a polar limb of a coronal hole is discovered in the Ca II filtergram of the Solar Optical Telescope on board the \textit{Hinode}. We identify 30 jets in the Ca II movie of duration of 53 min. The average speed at their maximum heights is found to be 132$\pm$44 km s$^{-1}$ ranging from 57 km s$^{-1}$ to 264 km s$^{-1}$ along the propagation direction. The average lifetime is 20$\pm$6 ranging from 11 seconds to 36 seconds. The speed and lifetime of the jets are located at end-tails of those parameters determined for type II spicules, hence implying a new type of jets. To confirm whether these jets are different from conventional spicules, we construct a time-height image averaged over horizontal region of 1$\arcsec$, and calculate lagged cross-correlations of intensity profiles at each height with the intensity at 2 Mm. From this, we obtain a cross-correlation map as a function of lag and height. We find that the correlation curve as a function of lag time is well fitted into three different Gaussian functions whose standard deviations of the lag time are 193 seconds, 42 seconds, and 17 seconds. The corresponding propagation speeds are calculated to be 9 km s$^{-1}$, 67 km s$^{-1}$, and 121 km s$^{-1}$, respectively. The kinematic properties of the former two components seem to correspond to the 3 minutes oscillations and type II spicules, while the latter component to the jets addressed in this study.

Quasi-periodic pulsations (QPPs) appear to be a common feature observed in the light curves of both solar and stellar flares. However, their quasi-periodic nature, along with the fact that they can be small in amplitude and short-lived, makes QPPs difficult to unequivocally detect. In this paper, we test the strengths and limitations of state-of-the-art methods for detecting QPPs using a series of hare-and-hounds exercises. The hare simulated a set of flares, both with and without QPPs of a variety of forms, while the hounds attempted to detect QPPs in blind tests. We use the results of these exercises to create a blueprint for anyone who wishes to detect QPPs in real solar and stellar data. We present eight clear recommendations to be kept in mind for future QPP detections, with the plethora of solar and stellar flare data from new and future satellites. These recommendations address the key pitfalls in QPP detection, including detrending, trimming data, accounting for colored noise, detecting stationary-period QPPs, detecting QPPs with nonstationary periods, and ensuring that detections are robust and false detections are minimized. We find that QPPs can be detected reliably and robustly by a variety of methods, which are clearly identified and described, if the appropriate care and due diligence are taken.

The ubiquity of recently discovered low-amplitude decayless kink oscillations of plasma loops allows for the seismological probing of the corona on a regular basis. In particular, in contrast to traditionally applied seismology which is based on the large-amplitude decaying kink oscillations excited by flares and eruptions, decayless oscillations can potentially provide the diagnostics necessary for their forecasting. We analysed decayless kink oscillations in several distinct loops belonging to active region NOAA 12107 on 10 July 2010 during its quiet time period, when it was observed on the West limb in EUV by the Atmospheric Imaging Assembly on-board Solar Dynamics Observatory. The oscillation periods were estimated with the use of the motion magnification technique. The lengths of the oscillating loops were determined within the assumption of its semicircular shape by measuring the position of their foot-points. The density contrast in the loops was estimated from the observed intensity contrast accounting for the unknown spatial scale of the background plasma. The combination of those measurements allows us to determine the distribution of kink and Alfv\'en speeds in the active region. Thus, we demonstrate the possibility to obtain seismological information about coronal active regions during the quiet periods of time.

Vortex shedding is an oscillating flow that is commonly observed in fluids due to the presence of a blunt body in a flowing medium. Numerical simulations have shown that the phenomenon of vortex shedding could also develop in the magnetohydrodynamic (MHD) domain. The dimensionless Strouhal number, the ratio of the blunt body diameter to the product of the period of vortex shedding and the speed of a flowing medium, is a robust indicator for vortex shedding, and, generally of the order of 0.2 for a wide range of Reynolds number. Using an observation from the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, we report a wavelike or oscillating plasma flow propagating upward against the Sun's gravitational force. A newly formed shrinking loop in the post-flare region possibly generates the oscillation of the upflow in the wake of the hot and dense loop through vortex shedding. The computed Strouhal number is consistent with the prediction from previous MHD simulations. Our observation suggests the possibility of vortex shedding in the solar corona.

Slow magnetoacoustic waves are omnipresent in both natural and laboratory plasma systems. The wave-induced misbalance between plasma cooling and heating processes causes the amplification or attenuation, and also dispersion, of slow magnetoacoustic waves. The wave dispersion could be attributed to the presence of characteristic time scales in the system, connected with the plasma heating or cooling due to the competition of the heating and cooling processes in the vicinity of the thermal equilibrium. We analysed linear slow magnetoacoustic waves in a plasma in a thermal equilibrium formed by a balance of optically thin radiative losses, field-align thermal conduction, and an unspecified heating. The dispersion is manifested by the dependence of the effective adiabatic index of the wave on the wave frequency, making the phase and group speeds frequency-dependent. The mutual effect of the wave amplification and dispersion is shown to result into the occurrence of an oscillatory pattern in an initially broadband slow wave, with the characteristic period determined by the thermal misbalance time scales, i.e. by the derivatives of the combined radiation loss and heating function with respect to the density and temperature, evaluated at the equilibrium. This effect is illustrated by estimating the characteristic period of the oscillatory pattern, appearing because of thermal misbalance in the plasma of the solar corona. It is found that by an order of magnitude the period is about the typical periods of slow magnetoacoustic oscillations detected in the corona.

Rapidly decaying slow magnetoacoustic waves are regularly observed in the solar coronal structures, offering a promising tool for a seismological diagnostics of the coronal plasma, including its thermodynamical properties. The effect of damping of standing slow magnetoacoustic oscillations in the solar coronal loops is investigated accounting for the field-aligned thermal conductivity and a wave-induced misbalance between radiative cooling and some unspecified heating rates. The non-adiabatic terms were allowed to be arbitrarily large, corresponding to the observed values. The thermal conductivity was taken in its classical form, and a power-law dependence of the heating function on the density and temperature was assumed. The analysis was conducted in the linear regime and in the infinite magnetic field approximation. The wave dynamics is found to be highly sensitive to the characteristic time scales of the thermal misbalance. Depending on certain values of the misbalance time scales three regimes of the wave evolution were identified, namely the regime of a suppressed damping, enhanced damping where the damping rate drops down to the observational values, and acoustic over-stability. The specific regime is determined by the dependences of the radiative cooling and heating functions on thermodynamical parameters of the plasma in the vicinity of the perturbed thermal equilibrium. The comparison of the observed and theoretically derived decay times and oscillation periods allows us to constrain the coronal heating function. For typical coronal parameters, the observed properties of standing slow magnetoacoustic oscillations could be readily reproduced with a reasonable choice of the heating function.

Recently, spiral wave patterns (SWPs) have been detected in 3 minute oscillations of sunspot umbrae, but the nature of this phenomenon has remained elusive. We present a theoretical model that interprets the observed SWPs as the superposition of two different azimuthal modes of slow magnetoacoustic waves driven below the surface in an untwisted and non-rotating magnetic cylinder. We apply this model to SWPs of the line-of-sight (LOS) velocity in a pore observed by the Fast Imaging Solar Spectrograph installed at the 1.6 m Goode Solar Telescope. One- and two-armed SWPs were identified in instantaneous amplitudes of LOS Doppler velocity maps of 3 minute oscillations. The associated oscillation periods are about 160 s, and the durations are about 5 minutes. In our theoretical model, the observed spiral structures are explained by the superposition of non-zero azimuthal modes driven 1600 km below the photosphere in the pore. The one-armed SWP is produced by the slow-body sausage (m = 0) and kink (m = 1) modes, and the two-armed SWP is formed by the slow-body sausage (m = 0) and fluting (m = 2) modes of the magnetic flux tube forming the pore.

Quasi-periodic pulsations (QPPs) are a common feature of solar flares, but previously there has been a lack of observational evidence to support any of the theoretical models that might explain the origin of QPPs. We aimed to determine if there are any relationships between the QPP period and other properties of the flaring region, using the sample of flares with QPPs from Pugh et al. (2017b). If any relationships exist then these can be compared with scaling laws for the theoretical QPP mechanisms. To obtain the flaring region properties we made use of the AIA 1600 and HMI data. The AIA 1600 images allow the flare ribbons to be seen while the HMI magnetograms allow the positive and negative magnetic polarity ribbons to be distinguished and the magnetic properties determined. The ribbon properties calculated in this study were the ribbon separation distance, area, total unsigned magnetic flux, and average magnetic field strength. Only the flares that occurred within \pm 60{\deg} of the solar disk centre were included, which meant a sample of 20 flares with 22 QPP signals. Positive correlations were found between the QPP period and the ribbon properties. The strongest correlations were with the separation distance and magnetic flux. Because these ribbon properties also correlate with the flare duration, and the relationship between the QPP period and flare duration may be influenced by observational bias, we also made use of simulated data to check if artificial correlations could be introduced. These simulations show that although QPPs cannot be detected for certain combinations of QPP period and flare duration, this does not introduce an apparent correlation. There is evidence of relationships between the QPP period and flare ribbon properties, and in the future the derived scaling laws between these properties can be compared to equivalent scaling laws for theoretical QPP mechanisms.

The Alfv\'{e}n speed and plasma beta in photospheric bright points observed by the Broadband Filter Imager (BFI) of the Solar Optical Telescope (SOT) on board the \textit{Hinode} satellite, are estimated seismologically. The diagnostics is based on the theory of slow magnetoacoustic waves in a non-isothermally stratified photosphere with a uniform vertical magnetic field. We identify and track bright points in a G-band movie by using the 3D region growing method, and align them with blue continuum images to derive their brightness temperatures. From the Fourier power spectra of 118 continuum light curves made in the bright points, we find that light curves of 91 bright points have oscillations with properties which are significantly different from oscillation in quiet regions, with the periods ranging 2.2--16.2~min. We find that the model gives a moderate value of the plasma beta when $\gamma$ lies at around 5/3. The calculated Alfv\'{e}n speed is 9.68$\pm$2.02~km~s$^{-1}$, ranging in 6.3--17.4~km~s$^{-1}$. The plasma beta is estimated to be of 0.93$\pm$0.36, ranging in 0.2--1.9.

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

Measurement of the solar wind speed near the Sun is important for understanding the acceleration mechanism of the solar wind. In this study, we determine 2D solar wind speeds from 6 to 26 solar radii by applying Fourier motion filters to \textit{SOHO}/LASCO C3 movies observed from 1999 to 2010. Our method successfully reproduces the original flow speeds in the artificially generated data as well as streamer blobs. We measure 2D solar wind speeds from 1-day to 1-year timescales and their variation in solar cycle 24. We find that the solar wind speeds at timescales longer than a month in the solar maximum period are relatively uniform in the azimuthal direction, while they are clearly bimodal in the minimum period, as expected from the \textit{Ulysses} observations and IPS reconstruction. The bimodal structure appears at around 2006, becomes most distinctive in 2009, and abruptly disappears in 2010. The radial evolution of the solar wind speeds resembles the Parker's solar wind solution.

Observational detection of quasi-periodic drifting fine structures in a type III radio burst associated with a solar flare SOL2015-04-16T11:22, with Low Frequency Array, is presented. Although similar modulations of the type III emission have been observed before and were associated with the plasma density fluctuations, the origin of those fluctuations was unknown. Analysis of the striae of the intensity variation in the dynamic spectrum allowed us to reveal two quasi-oscillatory components. The shorter component has the apparent wavelength of $\sim2$ Mm, phase speed of $\sim657$ km s$^{-1}$, which gives the oscillation period of $\sim3$ s, and the relative amplitude of $\sim0.35$%. The longer component has the wavelength of $\sim12$ Mm, and relative amplitude of $\sim5.1$%. The short frequency range of the detection does not allow us to estimate its phase speed. However, the properties of the shorter oscillatory component allowed us to interpret it as a fast magnetoacoustic wave guided by a plasma non-uniformity along the magnetic field outwards from the Sun. The assumption that the intensity of the radio emission is proportional to the amount of plasma in the emitting volume allowed us to show that the superposition of the plasma density modulation by a fast wave and a longer-wavelength oscillation of an unspecified nature could readily reproduce the fine structure of the observed dynamic spectrum. The observed parameters of the fast wave give the absolute value of the magnetic field in the emitting plasma of $\sim1.1$ G which is consistent with the radial magnetic field model.

We present analysis of the magnetic field in seven solar flare regions accompanied by the pulsations of hard X-ray (HXR) emission. These flares were studied by Kuznetsov et al. (2016) (Paper~I), and chosen here because of the availability of the vector magnetograms for their parent active regions (ARs) obtained with the SDO/HMI data. In Paper~I, based on the observations only, it was suggested that a magnetic flux rope (MFR) might play an important role in the process of generation of the HXR pulsations. The goal of the present paper is to test this hypothesis by using the extrapolation of magnetic field with the non-linear force-free field (NLFFF) method. Having done this, we found that before each flare indeed there was an MFR elongated along and above a magnetic polarity inversion line (MPIL) on the photosphere. In two flare regions the sources of the HXR pulsations were located at the footpoints of different magnetic field lines wrapping around the central axis, and constituting an MFR by themselves. In five other flares the parent field lines of the HXR pulsations were not a part of an MFR, but surrounded it in the form of an arcade of magnetic loops. These results show that, at least in the analyzed cases, the "single flare loop" models do not satisfy the observations and magnetic field modeling, while are consistent with the concept that the HXR pulsations are a consequence of successive episodes of energy release and electron acceleration in different magnetic flux tubes (loops) of a complex AR. An MFR could generate HXR pulsations by triggering episodes of magnetic reconnection in different loops in the course of its non-uniform evolution along an MPIL. However, since three events studied here were confined flares, actual eruptions may not be required to trigger sequential particle acceleration episodes in the magnetic systems containing an MFR.

Quasi-periodic pulsations (QPP) are common in solar flares and are now regularly observed in stellar flares. We present the detection of two different types of QPP signals in the thermal emission light curves of the X9.3 class solar flare SOL2017-09-06T12:02, which is the most powerful flare of Cycle 24. The period of the shorter-period QPP drifts from about 12 to 25 seconds during the flare. The observed properties of this QPP are consistent with a sausage oscillation of a plasma loop in the flaring active region. The period of the longer-period QPP is about 4 to 5 minutes. Its properties are compatible with standing slow magnetoacoustic oscillations, which are often detected in coronal loops. For both QPP signals, other mechanisms such as repetitive reconnection cannot be ruled out, however. The studied solar flare has an energy in the realm of observed stellar flares, and the fact that there is evidence of a short-period QPP signal typical of solar flares along with a long-period QPP signal more typical of stellar flares suggests that the different ranges of QPP periods typically observed in solar and stellar flares is likely due to observational constraints, and that similar physical processes may be occurring in solar and stellar flares.

Context. During their journey to perihelion, comets may appear in the field-of-view of space-borne optical instruments, showing in some cases a nicely developed plasma tail extending from their coma and exhibiting an oscillatory behaviour. Aims. The oscillations of cometary tails may be explained in terms of vortex shedding because of the interaction of the comet with the solar wind streams. Therefore, it is possible to exploit these oscillations in order to infer the value of the Strouhal number $St$, which quantifies the vortex shedding phenomenon, and the physical properties of the local medium. Methods. We used the Heliospheric Imager (HI) data of the Solar TErrestrial Relations Observatory (STEREO) mission to study the oscillations of the tails of the comets 2P/Encke and C/2012 S1 (ISON) during their perihelion in Nov 2013, determining the Strouhal numbers from the estimates of the halo size, the relative speed of the solar wind flow and the period of the oscillations. Results. We found that the estimated Strouhal numbers are very small, and the typical value of $St\sim0.2$ would be extrapolated for size of the halo larger than $\sim10^6$ km. Conclusions. Despite the vortex shedding phenomenon has not been unambiguously revealed, the findings suggest that some MHD instability process is responsible for the observed behaviour of cometary tails, which can be exploited for probing the physical conditions of the near-Sun region.

The effects of finite amplitudes on the transverse oscillations of a quiescent prominence represented by a magnetic rope are investigated in terms of the model proposed by Kolotkov et al. 2016. We consider a weakly nonlinear case governed by a quadratic nonlinearity, and also analyse the fully nonlinear equations of motion. We treat the prominence as a massive line current located above the photosphere and interacting with the magnetised dipped environment via the Lorentz force. In this concept the magnetic dip is produced by two external current sources located at the photosphere. Finite amplitude horizontal and vertical oscillations are found to be strongly coupled between each other. The coupling is more efficient for larger amplitudes and smaller attack angles between the direction of the driver and the horizontal axis. Spatial structure of oscillations is represented by Lissajous-like curves with the limit cycle of a hourglass shape, appearing in the resonant case, when the frequency of the vertical mode is twice the horizontal mode frequency. A metastable equilibrium of the prominence is revealed, which is stable for small amplitude displacements, and becomes horizontally unstable, when the amplitude exceeds a threshold value. The maximum oscillation amplitudes are also analytically derived and analysed. Typical oscillation periods are determined by the oscillation amplitude, prominence current, its mass and position above the photosphere, and the parameters of the magnetic dip. The main new effects of the finite amplitude are the coupling of the horizontally and vertically polarised transverse oscillations (i.e. the lack of a simple, elliptically polarised regime) and the presence of metastable equilibria of prominences.

Solar flare emission is detected in all EM bands and variations in flux density of solar energetic particles. Often the EM radiation generated in solar and stellar flares shows a pronounced oscillatory pattern, with characteristic periods ranging from a fraction of a second to several minutes. These oscillations are referred to as quasi-periodic pulsations (QPPs), to emphasise that they often contain apparent amplitude and period modulation. We review the current understanding of quasi-periodic pulsations in solar and stellar flares. In particular, we focus on the possible physical mechanisms, with an emphasis on the underlying physics that generates the resultant range of periodicities. These physical mechanisms include MHD oscillations, self-oscillatory mechanisms, oscillatory reconnection/reconnection reversal, wave-driven reconnection, two loop coalescence, MHD flow over-stability, the equivalent LCR-contour mechanism, and thermal-dynamical cycles. We also provide a histogram of all QPP events published in the literature at this time. The occurrence of QPPs puts additional constraints on the interpretation and understanding of the fundamental processes operating in flares, e.g. magnetic energy liberation and particle acceleration. Therefore, a full understanding of QPPs is essential in order to work towards an integrated model of solar and stellar flares.

We investigate the properties of a set of solar flares originating from a single active region (AR) that exhibit QPPs, and look for signs of the QPP periods relating to AR properties. The AR studied, best known as NOAA 12192, was unusually long-lived and produced 181 flares. Data from the GOES, EVE, Fermi, Vernov and NoRH observatories were used to determine if QPPs were present in the flares. For the soft X-ray GOES and EVE data, the time derivative of the signal was used. Power spectra of the time series data (without any form of detrending) were inspected, and flares with a peak above the 95% confidence level in the spectrum were labelled as having candidate QPPs. The confidence levels were determined taking account of uncertainties and the possible presence of red noise. AR properties were determined using HMI line of sight magnetograms. A total of 37 flares (20% of the sample) show good evidence of having QPPs, and some of the pulsations can be seen in data from multiple instruments and in different wavebands. The QPP periods show a weak correlation with the flare amplitude and duration, but this may be due to an observational bias. A stronger correlation was found between the QPP period and duration of the QPP signal, which can be partially but not entirely explained by observational constraints. No correlations were found with the AR area, bipole separation, or average magnetic field strength. The fact that a substantial fraction of the flare sample showed evidence of QPPs using a strict detection method with minimal processing of the data demonstrates that these QPPs are a real phenomenon, which cannot be explained by the presence of red noise or the superposition of multiple unrelated flares. The lack of correlation between the QPP periods and AR properties implies that the small-scale structure of the AR is important, and/or that different QPP mechanisms act in different cases.

This paper presents an observation of quasi-periodic rapidly-propagating waves observed in the AIA 171/193 \AA~ channels during the impulsive phase of an M1.9 flare occurred on 7 May 2012. The instant period was found to decrease from 240~s to 120~s, the speed of the wave fronts was in the range of ~664-1416 km/s. Almost simultaneously, quasi-periodic bursts with similar instant periods, ~70~s and ~140~s, occur in the microwave emission and in decimetric type IV, and type III radio bursts, and in the soft X-ray emission. The magnetic field configuration of the flare site was consistent with a breakout topology, i.e., a quadrupolar field along with a magnetic null point. The quasi-periodic rapidly-propagating wavefronts of the EUV emission are interpreted as a fast magnetoacoustic wave train. The observations suggest that the fast-mode waves are generated during the quasi-periodic magnetic reconnection in the cusp-region above the flare arcade loops. For the first time, we provide the evidence of a tadpole wavelet signature at about 70--140~s in decimetric (245/610~MHz) radio bursts, along with the direct observation of a coronal fast-mode wave train in EUV. In addition, at AIA 131/193 \AA\ we observed quasi-periodic EUV disturbances with the periods of 95~s and 240~s propagating downward at the apparent speed of 172-273 km/s. The nature of these downward propagating disturbances is not revealed, but they could be connected with magnetoacoustic waves or periodically shrinking loops.

We analyse nonlinear evolution of torsional Alfv\'en waves in a straight magnetic flux tube filled in with a low-$\beta$ plasma, and surrounded with a plasma of lower density. Such magnetic tubes model, in particular, a segment of a coronal loop or a polar plume. The wavelength is taken comparable to the tube radius. We perform numerical simulation of the wave propagation using ideal magnetohydrodynamics. We find that a torsional wave nonlinearly induces three kinds of compressive flows: the parallel flow at the Alfv\'en speed, that constitutes a bulk plasma motion along the magnetic field, the tube wave, and also transverse flows in the radial direction, associated with sausage fast magnetoacoustic modes. In addition, nonlinear torsional wave steepens and its propagation speed increases. The latter effect leads to the progressive distortion of torsional wave front, i.e. nonlinear phase mixing. Because of the intrinsic non-uniformity of the torsional wave amplitude across the tube radius, the nonlinear effects are more pronounced in regions with higher wave amplitude. They are always absent at the axes of the flux tube. In the case of a linear radial profile of the wave amplitude, the nonlinear effects are localised in an annulus region near the tube boundary. Thus, the parallel compressive flows driven by torsional Alfv\'en waves in the solar and stellar coronae, are essentially non-uniform in the perpendicular direction. The presence of additional sinks for the wave energy reduces the efficiency of the nonlinear parallel cascade in torsional Alfv\'en waves.

The robust detection of quasi-periodic pulsations (QPPs) in solar and stellar flares has been the topic of recent debate. In light of this, we have adapted a method described by Vaughan (2005) to aid with the search for QPPs in flare time series data. The method identifies statistically significant periodic signals in power spectra, and properly accounts for red noise as well as the uncertainties associated with the data. We show how the method can be further developed to be used with rebinned power spectra, allowing us to detect QPPs whose signal is spread over more than one frequency bin. An advantage of these methods is that there is no need to detrend the data prior to creating the power spectrum. Examples are given where the methods have been applied to synthetic data, as well as real flare time series data with candidate QPPs from the Nobeyama Radioheliograph. These show that, despite the transient nature of QPPs, peaks corresponding to the QPPs can be seen at a significant level in the power spectrum without any form of detrending or other processing of the original time series data, providing the background trends are not too steep.

This paper presents the multi-wavelength analysis of a 13 s quasi-periodic pulsation (QPP) observed in hard X-ray (12-300 keV) and microwave (4.9-34 GHz) emissions during a C-class flare occurred on 21 September 2015. AIA 304 and 171 \AA~ images show an emerging loop/flux tube (L1) moving radially outward, which interacts with preexisting structures within the active region. The QPP was observed during the expansion and rising motion of L1. The Nobeyama Radioheliograph (NoRH) microwave images in 17/34 GHz channels reveal a single radio source, which was co-spatial with a neighboring loop (L2). In addition, using AIA 304 \AA~ images, we detected intensity oscillations in the legs of loop L2 with a period of about 26 s. A similar oscillation period was observed in the GOES soft X-ray flux derivative. This oscillation period seems to increase with time. We suggest that the observed QPP is most likely generated by the interaction between loops L2 and L3 observed in the AIA hot channels (131 and 94 \AA). The merging speed of loops L2 and L3 was ~35 km/s. Loop L1 destroyed possibly by its interaction with preexisting structures in the active region and produced a cool jet with the speed ~106-118 km/s associated with a narrow CME (~770 km/s). Another mechanism of the QPP in terms of a sausage oscillation of the loop (L2) is also possible.

We present observations of a B7.9-class flare that occurred on the 24th January, 2015, using the Atmopsheric Imaging Assembly (AIA) of the Solar Dynamics Observatory (SDO), the EUV Imaging Spectrometer (EIS) and the X-Ray Telescope of Hinode. The flare triggers the eruption of a dense cool plasma blob as seen in AIA 171\AA,\, which is unable to completely break out and remains confined within a local bundle of active region loops. During this process, transverse oscillations of the threads are observed. The cool plasma is then observed to descend back to the chromosphere along each loop strand. At the same time, a larger diffuse co-spatial loop observed in the hot wavebands of SDO/AIA and Hinode/XRT is formed, exhibiting periodic intensity variations along its length. The formation and evolution of magnetohydrodynamic (MHD) waves depend upon the values of the local plasma parameters (e.g. density, temperature and magnetic field), which can hence be inferred by coronal seismology. In this study we aim to assess how the observed MHD modes are affected by the variation of density and temperature. We combined analysis of EUV/X-ray imaging and spectroscopy using SDO/AIA, Hinode/EIS and XRT. We show that the evolution of the detected waves is determined by the temporal variations of the local plasma parameters, caused by the flare heating and the consequent cooling. We apply coronal seismology to both waves obtaining estimations of the background plasma parameters.

We introduce a new method for the investigation of low-amplitude transverse oscillations of solar plasma non-uniformities, such as coronal loops, individual strands in coronal arcades, jets, prominence fibrils, polar plumes, and other contrast features, observed with imaging instruments. The method is based on the two-dimensional dual tree complex wavelet transform (DT$\mathbb{C}$WT). It allows us to magnify transverse, in the plane-of-the-sky, quasi-periodic motions of contrast features in image sequences. The tests performed on the artificial data cubes imitating exponentially decaying, multi-periodic and frequency-modulated kink oscillations of coronal loops showed the effectiveness, reliability and robustness of this technique. The algorithm was found to give linear scaling of the magnified amplitudes with the original amplitudes provided they are sufficiently small. Also, the magnification is independent of the oscillation period in a broad range of the periods. The application of this technique to SDO/AIA EUV data cubes of a non-flaring active region allowed for the improved detection of low-amplitude decay-less oscillations in the majority of loops.

Radio emission observations from the Learmonth and Bruny Island radio spectrographs are analysed to determine the nature of a train of discrete, periodic radio \lq sparks\rq (finite-bandwidth, short-duration isolated radio features) which precede a type II burst. We analyse extreme ultraviolet (EUV) imaging from SDO/AIA at multiple wavelengths and identify a series of quasi-periodic rapidly-propagating enhancements, which we interpret as a fast wave train, and link these to the detected radio features. The speeds and positions of the periodic rapidly propagating fast waves and the coronal mass ejection (CME) were recorded using running-difference images and time-distance analysis. From the frequency of the radio sparks the local electron density at the emission location was estimated for each. Using an empirical model for the scaling of density in the corona, the calculated electron density was used to obtain the height above the surface at which the emission occurs, and the propagation velocity of the emission location. The period of the radio sparks, $\delta t_\mathrm{r}$=1.78$\pm$0.04 min, matches the period of the fast wave train observed at 171 $\AA$, $\delta t_\mathrm{EUV}$=1.7 $\pm$ 0.2 min. The inferred speed of the emission location of the radio sparks, 630 km s$^{-1}$, is comparable to the measured speed of the CME leading edge, 500 km s$^{-1}$, and the speeds derived from the drifting of the type II lanes. The calculated height of the radio emission (obtained from the density) matches the observed location of the CME leading edge. From the above evidence we propose that the radio sparks are caused by the quasi-periodic fast waves, and the emission is generated as they catch up and interact with the leading edge of the CME.

Sausage modes that are axisymmetric fast magnetoacoustic oscillations of solar coronal loops are characterized by variation of the plasma density and magnetic field, and hence cause time variations of the electron plasma frequency and cyclotron frequency. The latter parameters determine the condition for the double plasma resonance (DPR), which is responsible for the appearance of zebra-pattern (ZP) structures in time spectra of solar type IV radio bursts. We perform numerical simulations of standing and propagating sausage oscillations in a coronal loop modeled as a straight, field-aligned plasma slab, and determine the time variation of the DPR layer locations. Instant values of the plasma density and magnetic field at the DPR layers allowed us to construct skeletons of the time variation of ZP stripes in radio spectra. In the presence of a sausage oscillation, the ZP structures are shown to have characteristic wiggles with the time period prescribed by the sausage oscillation. Standing and propagating sausage oscillations are found to have different signatures in ZP patterns. We conclude that ZP wiggles can be used for the detection of short-period sausage oscillations and the exploitation of their seismological potential.

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

We embark on a study of quasi-periodic pulsations (QPPs) in the decay phase of white-light stellar flares observed by Kepler. Out of the 1439 flares on 216 different stars detected in the short-cadence data using an automated search, 56 flares are found to have pronounced QPP-like signatures in the light curve, of which 11 have stable decaying oscillations. No correlation is found between the QPP period and the stellar temperature, radius, rotation period and surface gravity, suggesting that the QPPs are independent of global stellar parameters. Hence they are likely to be the result of processes occurring in the local environment. There is also no significant correlation between the QPP period and flare energy, however there is evidence that the period scales with the QPP decay time for the Gaussian damping scenario, but not to a significant degree for the exponentially damped case. This same scaling has been observed for MHD oscillations on the Sun, suggesting that they could be the cause of the QPPs in those flares. Scaling laws of the flare energy are also investigated, supporting previous reports of a strong correlation between the flare energy and stellar temperature/radius. A negative correlation between the flare energy and stellar surface gravity is also found.

Standing long-period (with the periods longer than several minutes) oscillations in large hot (with the temperature higher than 3 MK) coronal loops have been observed as the quasi-periodic modulation of the EUV and microwave intensity emission and the Doppler shift of coronal emission lines, and have been interpreted as standing slow magnetoacoustic (longitudinal) oscillations. Quasi-periodic pulsations of shorter periods, detected in thermal and non-thermal emissions in solar flares could be produced by a similar mechanism. We present theoretical modelling of the standing slow magnetoacoustic mode, showing that this mode of oscillation is highly sensitive to peculiarities of the radiative cooling and heating function. We generalised the theoretical model of standing slow magnetoacoustic oscillations in a hot plasma, including the effects of the radiative losses, and accounting for plasma heating. The heating mechanism is not specified and taken empirically to compensate the cooling by radiation and thermal-conduction. It is shown that the evolution of the oscillations is described by a generalised Burgers equation. Numerical solution of an initial value problem for the evolutionary equation demonstrates that different dependences of the radiative cooling and plasma heating on the temperature lead to different regimes of the oscillations, including growing, quasi-stationary and rapidly decaying. Our findings provide a theoretical foundation for probing the coronal heating function, and may explain the observations of decayless long-period quasi-periodic pulsations in flares. The hydrodynamic approach employed in this study should be considered with caution in the modelling of non-thermal emission associated with flares, as it misses potentially important non-hydrodynamic effects.

We present multi-wavelength analysis of a quasi-periodic pulsation (QPP) observed in the hard X-ray, radio, and extreme-ultraviolet (EUV) channels during an M1.9 flare occurred on 23-24 September 2011. The non-thermal hard X-ray emission in 25-50 keV observed by RHESSI shows five distinct impulsive peaks of decaying amplitude with a period of about three minutes. Similar QPP was observed in the microwave emission recorded by the Nobeyama Radioheliograph and Polarimeter in the 8.8, 15, 17 GHz channels. Interestingly, the 3-min QPP was also observed in the metric and decimetric radio frequencies (25-180, 245, 610 MHz) as repetitive type III bursts. Multi-wavelength observations from the SDO/AIA, Hinode/SOT, and STEREO/SECCHI suggest a fan-spine topology at the eruption site, associated with the formation of a quasi-circular ribbon during the flare. A small filament was observed below the fan-loops before the flare onset. The filament rose slowly and interacted with the ambient field. This behaviour was followed by an untwisting motion of the filament. Two different structures of the filament showed $\sim$3-min periodic alternate rotation in the clockwise and counterclockwise directions. The 3-min QPP was found to highly correlate with 3-min oscillations in a nearby sunspot. We suggest that the periodic reconnection (modulated either by sunspot slow-mode wave or by untwisting filament) at a magnetic null-point most likely causes the repetitive particle acceleration, generating the QPP observed in hard X-ray, microwave and type III radio bursts.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

Despite intensive studies of kink oscillations of coronal loops in the last decade, a large scale statistically significant investigation of the oscillation parameters has not been made using data from the Solar Dynamics Observatory (SDO). We carry out a statistical study of kink oscillations using Extreme Ultra-Violet (EUV) imaging data from a previously compiled catalogue. We analysed 58 kink oscillation events observed by the Atmospheric Imaging Assembly (AIA) onboard SDO during its first four years of operation (2010-2014). Parameters of the oscillations, including the initial apparent amplitude, period, length of the oscillating loop, and damping are studied for 120 individual loop oscillations. Analysis of the initial loop displacement and oscillation amplitude leads to the conclusion that the initial loop displacement prescribes the initial amplitude of oscillation in general. The period is found to scale with the loop length, and a linear fit of the data cloud gives a kink speed of Ck =(1330+/-50) km s-1 . The main body of the data corresponds to kink speeds in the range Ck =(800-3300) km s-1. Measurements of 52 exponential damping times were made, and it was noted that at least 22 of the damping profiles may be better approximated by a combination of non-exponential and exponential profiles, rather than a purely exponential damping envelope. There are an additional 10 cases where the profile appears to be purely non-exponential, and no damping time was measured. A scaling of the exponential damping time with the period is found, following the previously established linear scaling between these two parameters.

Fast magnetoacoustic waves guided along the magnetic field by plasma non-uniformities, in particular coronal loops, fibrils and plumes, are known to be highly dispersive, which leads to the formation of quasi-periodic wave trains excited by a broadband impulsive driver, e.g. a solar flare. We investigated effects of cylindrical geometry on the fast sausage wave train formation. We performed magnetohydrodynamic numerical simulations of fast magnetoacoustic perturbations of a sausage symmetry, propagating from a localised impulsive source along a field-aligned plasma cylinder with a smooth radial profile of the fast speed. The wave trains are found to have pronounced period modulation, with the longer instant period seen in the beginning of the wave train. The wave trains have also a pronounced amplitude modulation. Wavelet spectra of the wave trains have characteristic tadpole features, with the broadband large-amplitude heads preceding low-amplitude quasi-monochromatic tails. The mean period of the wave train is about the transverse fast magnetoacoustic transit time across the cylinder. The mean parallel wavelength is about the diameter of the waveguiding plasma cylinder. Instant periods are longer than the sausage wave cutoff period. The wave train characteristics depend on the fast magnetoacoustic speed in both the internal and external media, and the smoothness of the transverse profile of the equilibrium quantities, and also the spatial size of the initial perturbation. If the initial perturbation is localised at the axis of the cylinder, the wave trains contain higher radial harmonics that have shorter periods.

Flares that are orders of magnitude larger than the most energetic solar flares are routinely observed on Sun-like stars, raising the question of whether the same physical processes are responsible for both solar and stellar flares. In this letter we present a white-light stellar superflare on the star KIC9655129, observed by NASA's Kepler mission, with a rare multi-period quasi-periodic pulsation (QPP) pattern. Two significant periodic processes were detected using the wavelet and autocorrelation techniques, with periods of 78 +/- 12 min and 32 +/- 2 min. By comparing the phases and decay times of the two periodicities, the QPP signal was found to most likely be linear, suggesting that the two periodicities are independent, possibly corresponding either to different magnetohydrodynamic modes of the flaring region, or different spatial harmonics of the same mode. The presence of multiple periodicities is a good indication that the QPPs were caused by magnetohydrodynamic oscillations, and suggests that the physical processes in operation during stellar flares could be the same as those in solar flares.

We investigate the decayless regime of coronal kink oscillations recently discovered in the Solar Dynamics Observatory (SDO)/AIA data. In contrast to decaying kink oscillations that are excited by impulsive dynamical processes, this type of transverse oscillations is not connected to any external impulsive impact, such as a flare or CME, and does not show any significant decay. Moreover the amplitude of these decayless oscillations is typically lower than that of decaying oscillations. The aim of this research is to estimate the prevalence of this phenomenon and its characteristic signatures. We analysed 21 active regions (NOAA 11637--11657) observed in January 2013 in the 171 A channel of SDO/AIA. For each active region we inspected six hours of observations, constructing time-distance plots for the slits positioned across pronounced bright loops. The oscillatory patterns in time-distance plots were visually identified and the oscillation periods and amplitudes were measured. We also estimated the length of each oscillating loop. Low-amplitude decayless kink oscillations are found to be present in the majority of the analysed active regions. The oscillation periods lie in the range from 1.5 to 10~minutes. In two active regions with insufficient observation conditions we did not identify any oscillation patterns. The oscillation periods are found to increase with the length of the oscillating loop. The considered type of coronal oscillations is a common phenomenon in the corona. The established dependence of the oscillation period on the loop length is consistent with their interpretation in terms of standing kink waves.

Measurements of the magnetic field in the interplanetary medium, of the sunspots area, and of the heliospheric current sheet position, reveal a possible North-South asymmetry in the magnetic field of the Sun. We study the North-South asymmetry as inferred from measurements of the deflection of polar coronal hole jets when they propagate throughout the corona. Since the corona is an environment where the magnetic pressure is greater than the kinetic pressure, we can assume that magnetic field controls the dynamics of plasma. On average, jets during their propagation follow the magnetic field lines, highlighting its local direction. The average jet deflection is studied both in the plane perpendicular to the line of sight, and, for a reduced number of jets, in three dimensional space. The observed jet deflection is studied in terms of an axisymmetric magnetic field model comprising dipole We measured the position angles at 1 rs and at 2 rs of the 79 jets from the catalogue of Nistico et al 2009., based on the STEREO ultraviolet and white-light coronagraph observations during the solar minimum period March 2007-April 2008. We found that the propagation is not radial, in agreement with the deflection due to magnetic field lines. Moreover, the amount of the deflection is different between jets over the north and those from the south pole. Comparison of jet deflections and field line tracing shows that a ratio g2/g1 ~ -0.5 for the quadrupole and a ratio g3/g1 ~ 1.6-2.0 for the esapole can describe the field. The presence of a non-negligible quadrupole moment. We find that the magnetic deflection of jets is larger in the North than in the South of the order of 25-40%, with an asymmetry which is consistent with a southward deflection of the heliospheric current sheet of the order of 10 deg, consistent with that inferred from other, independent, datasets and instruments.

The last solar minimum was, by recent standards, unusually deep and long. We are now close to the maximum of the subsequent solar cycle, which is relatively weak. In this article we make comparisons between different global (unresolved) measures of the Sun's magnetic activity, to investigate how they are responding to this weak-activity epoch. We focus on helioseismic data, which are sensitive to conditions, including the characteristics of the magnetic field, in the solar interior. Also considered are measures of the magnetic field in the photosphere (sunspot number and sunspot area), the chromosphere and corona (10.7cm radio flux and 530.3nm green coronal index), and two measures of the Sun's magnetic activity closer to Earth (the interplanetary magnetic field and the galactic cosmic-ray intensity). Scaled versions of the activity proxies diverge from the helioseismic data around 2000, indicating a change in relationship between the proxies. The degree of divergence varies from proxy to proxy with sunspot area and 10.7cm flux showing only small deviations, while sunspot number, coronal index, and the two interplanetary proxies show much larger departures. In Cycle 24 the deviations in the solar proxies and the helioseismic data decrease, raising the possibility that the deviations observed in Cycle 23 are just symptomatic of a 22-year Hale cycle. However, the deviations in the helioseismic data and the interplanetary proxies increase in Cycle 24. Interestingly the divergence in the solar proxies and the helioseismic data are not reflected in the shorter-term variations (often referred to as quasi-biennial oscillations) observed on top of the dominant 11-year solar cycle. However, despite being highly correlated in Cycle 22, the short-term variations in the interplanetary proxies show very little correlation with the helioseismic data during Cycles 23 and 24.

The fields of solar radiophysics and solar system radio physics, or radio heliophysics, will benefit immensely from an instrument with the capabilities projected for SKA. Potential applications include interplanetary scintillation (IPS), radio-burst tracking, and solar spectral radio imaging with a superior sensitivity. These will provide breakthrough new insights and results in topics of fundamental importance, such as the physics of impulsive energy releases, magnetohydrodynamic oscillations and turbulence, the dynamics of post-eruptive processes, energetic particle acceleration, the structure of the solar wind and the development and evolution of solar wind transients at distances up to and beyond the orbit of the Earth. The combination of the high spectral, time and spatial resolution and the unprecedented sensitivity of the SKA will radically advance our understanding of basic physical processes operating in solar and heliospheric plasmas and provide a solid foundation for the forecasting of space weather events.

Two different mechanisms may act to induce quasi-periodic pulsations (QPP) in whole-disk observations of stellar flares. One mechanism may be magneto-hydromagnetic (MHD) forces and other processes acting on flare loops as seen in the Sun. The other mechanism may be forced local acoustic oscillations due to the high-energy particle impulse generated by the flare (known as `sunquakes' in the Sun). We analyze short-cadence Kepler data of 257 flares in 75 stars to search for QPP in the flare decay branch or post-flare oscillations which may be attributed to either of these two mechanisms. About 18 percent of stellar flares show a distinct bump in the flare decay branch of unknown origin. The bump does not seem to be a highly-damped global oscillation because the periods of the bumps derived from wavelet analysis do not correlate with any stellar parameter. We detected damped oscillations covering several cycles (QPP), in seven flares on five stars. The periods of these oscillations also do not correlate with any stellar parameter, suggesting that these may be a due to flare loop oscillations. We searched for forced global oscillations which might result after a strong flare. To this end, we investigated the behaviour of the amplitudes of solar-like oscillations in eight stars before and after a flare. However, no clear amplitude change could be detected. We also analyzed the amplitudes of the self-excited pulsations in two delta Scuti stars and one gamma Doradus star before and after a flare. Again, no clear amplitude changes were found. Our conclusions are that a new process needs to be found to explain the high incidence of bumps in stellar flare light curves, that flare loop oscillations may have been detected in a few stars and that no conclusive evidence exists as yet for flare induced global acoustic oscillations (starquakes).

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

We report decaying quasi-periodic intensity oscillations in the X-ray (6-12 keV) and extreme ultraviolet (EUV) channels (131, 94, 1600, 304 \AA) observed by the Fermi GBM (Gamma-ray Burst Monitor) and SDO/AIA, respectively, during a C-class flare. The estimated period of oscillation and decay time in the X-ray channel (6-12 keV) was about 202 s and 154 s, respectively. A similar oscillation period was detected at the footpoint of the arcade loops in the AIA 1600 and 304 \AA channels. Simultaneously, AIA hot channels (94 and 131 \AA) reveal propagating EUV disturbances bouncing back and forth between the footpoints of the arcade loops. The period of the oscillation and decay time were about 409 s and 1121 s, respectively. The characteristic phase speed of the wave is about 560 km/s for about 115 Mm loop length, which is roughly consistent with the sound speed at the temperature about 10-16 MK (480-608 km/s). These EUV oscillations are consistent with the SOHO/SUMER Doppler-shift oscillations interpreted as the global standing slow magnetoacoustic wave excited by a flare. The flare occurred at one of the footpoints of the arcade loops, where the magnetic topology was a 3D fan-spine with a null-point. Repetitive reconnection at this footpoint could cause the periodic acceleration of non-thermal electrons that propagated to the opposite footpoint along the arcade and precipitating there, causing the observed 202-s periodicity. Other possible interpretations, e.g. the second harmonics of the slow mode are also discussed.

Magnetohydrodynamic waves interact with structured plasmas and reveal the internal magnetic and thermal structures therein, thereby having seismological applications in the solar atmosphere. We investigate the evolution of fast magnetoacoustic pulses in randomly structured plasmas, in the context of large-scale propagating waves in the solar atmosphere. We perform one dimensional numerical simulations of fast wave pulses propagating perpendicular to a constant magnetic field in a low-$\beta$ plasma with a random density profile across the field. Both linear and nonlinear regimes are considered. We study how the evolution of the pulse amplitude and width depends on their initial values and the parameters of the random structuring. A randomly structured plasma acts as a dispersive medium for a fast magnetoacoustic pulse, causing amplitude attenuation and broadening of the pulse width. After the passage of the main pulse, secondary propagating and standing fast waves appear in the plasma. Width evolution of both linear and nonlinear pulses can be well approximated by linear functions; however, narrow pulses may have zero or negative broadening. This arises because a narrow pulse is prone to splitting, while a broad pulse usually deviates less from their initial Gaussian shape and form ripple structures on top of the main pulse. A linear pulse decays at almost a constant rate, while a nonlinear pulse decays exponentially. A pulse interacts most efficiently with a random medium which has a correlation length of about half of its initial pulse width. The development of a detailed model of a fast MHD pulse propagating in highly structured medium substantiates the interpretation of EIT waves as fast magnetoacoustic waves. Evolution of a fast pulse provides us with a novel method to diagnose the sub-resolution filamentation of the solar atmosphere.

The high spatial and time resolution data obtained with SDO/AIA for the sunspot in active region NOAA 11131 on 08 December 2010 were analysed with the time-distance plot technique and the pixelised wavelet filtering method. Oscillations in the 3 min band dominate in the umbra. The integrated spectrum of umbral oscillations contains distinct narrowband peaks at 1.9 min, 2.3 min, and 2.8 min. The power significantly varies in time, forming distinct oscillation trains. The oscillation power distribution over the sunspot in the horizontal plane reveals that the enhancements of the oscillation amplitude, or wave fronts, have a distinct structure consisting of an evolving two-armed spiral and a stationary circular patch at the spiral origin, situated near the umbra centre. This structure is seen from the temperature minimum to the corona. In time, the spiral rotates anti-clockwise. The wave front spirality is most pronounced during the maximum amplitude phases of the oscillations. In the low-amplitude phases the spiral breaks into arc-shaped patches. The 2D cross-correlation function shows that the oscillations at higher atmospheric levels occur later than at lower layers. The phase speed is estimated to be about 100 km/s. The fine spectral analysis shows that the central patch corresponds to the high-frequency oscillations, while the spiral arms highlight the lower-frequency oscillations in the 3-min band. The vertical and horizontal radial structure of the oscillations is consistent with the model that interprets umbral oscillations as slow magnetoacoustic waves filtered by the atmospheric temperature non-uniformity in the presence of the magnetic field inclination from the vertical. The mechanism for the polar-angle structure of the oscillations, in particular the spirality of the wave fronts, needs to be revealed.

Solar Optical Telescope onboard Hinode observed a sunspot (AR 11836) with two light bridges (LBs) on 31 Aug 2013. We analysed a 2-hour \ion{Ca}{2} H emission intensity data set and detected strong 5-min oscillation power on both LBs and in the inner penumbra. The time-distance plot reveals that 5-min oscillation phase does not vary significantly along the thin bridge, indicating that the oscillations are likely to originate from the underneath. The slit taken along the central axis of the wide light bridge exhibits a standing wave feature. However, at the centre of the wide bridge, the 5-min oscillation power is found to be stronger than at its sides. Moreover, the time-distance plot across the wide bridge exhibits a herringbone pattern that indicates a counter-stream of two running waves originated at the bridge sides. Thus, the 5-min oscillations on the wide bridge also resemble the properties of running penumbral waves. The 5-min oscillations are suppressed in the umbra, while the 3-min oscillations occupy all three cores of the sunspot's umbra, separated by the LBs. The 3-min oscillations were found to be in phase at both sides of the LBs. It may indicate that either LBs do not affect umbral oscillations, or umbral oscillations at different umbral cores share the same source. Also, it indicates that LBs are rather shallow objects situated in the upper part of the umbra. We found that umbral flashes follow the life cycles of umbral oscillations with much larger amplitudes. They cannot propagate across LBs. Umbral flashes dominate the 3-min oscillation power within each core, however, they do not disrupt the phase of umbral oscillation.

We used the Pixelised Wavelet Filtering (PWF) method to compute narrow-band power maps of SDO/AIA imaging datasets in the 1700 \AA{}, 1600 \AA{} and 304 \AA{} bandpasses that correspond to different heights. The cut-off frequency was defined as contours where the spectral power droped to the median level. It was measured as a function of the spatial location. We inferred the magnetic field inclination according to the MAG wave theory in the low-$\beta$ limit and compared it with the potential field extrapolation.{We analysed intensity oscillations in a symmetric sunspot AR11131(08 Dec 2010) and an asymmetric sunspot AR11330 (27 Oct 2011). We reconstructed the magnetic field inclination in the radial direction for the symmetric sunspot and in both radial and azimuthal directions for the asymmetric sunspot. We observed 3D variation of the main oscillation periods in sunspots. We found that shorter-period oscillations were mostly constrained in sunspot umbrae, while longer-period oscillations formed an annular shape enclosing the umbra. Longer periods are found to be distributed further away from the sunspot centre. Our results indicate that 3-min oscillation are generated in the chromosphere, possibly by the acoustic resonator model, while 5-min and longer-period oscillations seemed to originate in a level under the photosphere. The reconstructed field inclinations gives the values of the field inclination that are systematically larger than the values obtained by the potential field extrapolation. The inclined magnetic field line can account for 60-80% of cut-offfrequency lowering only.

Analysis of the microwave data, obtained in the 17 GHz channel of the Nobeyama Radioheliograph during the M1.6 flare on 4th Nov 2010, revealed the presence of 11.8-min oscillations of the emitting plasma density. The oscilla- tions decayed with the characteristic time of about 25-min. These oscillations are also well-seen in the variation of EUV emission intensity measured in the 335 A channel of SDO/AIA. The observed properties of the oscillations are consistent with the properties of so-called SUMER oscillations, observed in the EUV and soft X-ray bands usually as a periodic Doppler shift. The accepted interpretation of SUMER oscillations is a standing slow magnetoacoustic wave. Our analysis presents the first direct observation of the slow magnetoacoustic oscillations in the microwave emission of a solar flare.

Quasi-periodic wiggles of microwave zebra pattern structures with period range from about 0.5 s to 1.5 s are found in a X-class solar flare on 2006 December 13 at the 2.6-3.8 GHz with the Chinese Solar Broadband Radio Spectrometer (SBRS/Huairou). Periodogram and correlation analysis show that the wiggles have two-three significant periodicities and almost in phase between stripes at different frequency. The Alfven speed estimated from the zebra pattern structures is about 700 Km/s. We obtain the spatial size of the waveguiding plasma structure to be about 1 Mm with the detected period of about 1 s. It suggests the ZP wiggles can be associated with the fast mag- netoacoustic oscillations in the flaring active region. The lack of a significant phase shift between wiggles of different stripes suggests that the ZP wiggles are caused by a standing sausage oscillation.

Aims: We perform an advanced study of the fast magnetoacoustic sausage oscillations of coronal loops in the context of MHD coronal seismology to establish the dependence of the sausage mode period and cut-off wavenumber on the plasma-beta of the loop-filling plasma. A parametric study of the ratios for different harmonics of the mode is also carried out. Methods: Full magnetohydrodynamic numerical simulations were performed using Lare2d, simulating hot, dense loops in a magnetic slab environment. The symmetric Epstein profile and a simple step-function profile were both used to model the density structure of the simulated loops. Analytical expressions for the cut-off wavenumber and the harmonic ratio between the second longitudinal harmonic and the fundamental were also examined. Results: It was established that the period of the global sausage mode is only very weakly dependent on the value of the plasma-beta inside a coronal loop, which justifies the application of this model to hot flaring loops. The cut-off wavenumber kc for the global mode was found to be dependent on both internal and external values of the plasma-beta, again only weakly. By far the most important factor in this case was the value of the density contrast ratio between the loop and the surroundings. Finally, the deviation of the harmonic ratio P1/2P2 from the ideal non-dispersive case was shown to be considerable at low k, again strongly dependent on plasma density. Quantifying the behaviour of the cut-off wavenumber and the harmonic ratio has significant applications to the field of coronal seismology.

Recent observations have revealed that MHD waves and oscillations are ubiquitous in the solar atmosphere, with a wide range of periods. We give a brief review of some aspects of MHD waves and coronal seismology which have recently been the focus of intense debate or are newly emerging. In particular, we focus on four topics: (i) the current controversy surrounding propagating intensity perturbations along coronal loops, (ii) the interpretation of propagating transverse loop oscillations, (iii) the ongoing search for coronal (torsional) Alfven waves and (iv) the rapidly developing topic of quasi-periodic pulsations (QPP) in solar flares.

Context. The height structure of 3-min oscillations over sunspots is studied in the context of the recently discovered effect of height inversion: over the umbra, the spatial location of the maximum of chromospheric 3-min oscillation power corresponds to the relative decrease in the power of photospheric oscillations. Aims. We investigate whether the height inversion of the power of 3-min oscillations is a common feature of the spatial structure of the oscillations for the majority of sunspots. Methods. Spectrogram sequences of H-alpha 6563{\AA} and FeI 6569{\AA} over sunspots, acquired with very high cadency (about 2 s or better) are obtained. The distribution of the oscillation power of the line-of-sight velocity signal is studied by using methods of wavelet frequency filtration and Fourier analysis. Results. The effect of the height inversion is found in 9 of 11 analyzed active regions. The interpretation of this effect is possibly connected to both the decrease in the level of photosphere in sunspot umbrae and the magnetic field topology.

Three-minute oscillations over sunspot's umbra in AR 11131 were observed simultaneously in UV/EUV emission by SDO/AIA and in radio emission by Nobeyama Radioheliograph (NoRH). We use 24-hours series of SDO and 8-hours series of NoRH observations to study spectral, spatial and temporal variations of pulsations in the 5-9 mHz frequency range at different layers of the solar atmosphere. High spatial and temporal resolution of SDO/AIA in combination with long-duration observations allowed us to trace the variations of the cut-off frequency and spectrum of oscillations across the umbra. We found that higher frequency oscillations are more pronounced closer to the umbra's center, while the lower frequencies concentrate to the peripheral parts. We interpreted this discovery as a manifestation of variation of the magnetic field inclination across the umbra at the level of temperature-minimum. Possible implications of this interpretation for the diagnostics of sunspot atmospheres is discussed.

We analyse 3-min oscillations of microwave and EUV emission generated at different heights of a sunspot atmosphere, studying the amplitude and frequency modulation of the oscillations, and its relationship with the variation of the spatial structure of the oscillations. High-resolution data obtained with the Nobeyama Radioheliograph, TRACE and SDO/AIA are analysed with the use of the Pixelised Wavelet Filtering and wavelet skeleton techniques. 3-min oscillations in sunspots appear in the form of repetitive trains of the duration 8-20 min (13 min in average). The typical interval between the trains is 30-50 min. The oscillation trains are transient in frequency and power. We detected a repetitive frequency drifts of 3-min oscillations during the development of individual trains. Wavelet analysis shows three types of the frequency drift: positive, negative and fluctuations without drift. The start and end of the drifts coincide with the start time and end of the train. The comparative study of 3-min oscillations in the sequences of microwave and EUV images show the appearance of new sources of the oscillations in sunspots during the development of the trains. These structures can be interpreted as waveguides that channel upward propagating waves, responsible for 3-min oscillations. A possible explanation of the observed properties is the operation of two simultaneous factors: dispersive evolution of the upwardly-propagating wave pulses and the non-uniformity of the distribution of the oscillation power over the sunspot umbra with different wave sources corresponding to different magnetic flux tubes with different physical conditions and line-of-sight angles.

Solar transients and eruptive phenomena which are ubiquitous in the solar atmosphere, can shed new light to the understanding of the outstanding problems like coronal heating and the solar wind acceleration. Observations in the entire electromagnetic spectrum of such dynamical processes of large and small-scale transient/eruptive events, with highly dynamic magnetic field configuration, and energetic particles, provide crucial information about the plasma processes at mega-Kelvin temperature embedded in a complex magnetic field, and also energy build-up/energy-release processes, taking place in such events. One of the most important phenomenological aspects of solar eruptive phenomena is the induced magnetohydrodynamic (MHD) waves generated during these energetic processes, which carry a potential signature to probing the solar active regions. In this paper, we briefly review the recent trends of the transient (e.g., flares) induced quasi-periodic oscillations in the solar atmosphere and discuss their implications in diagnosing the solar active regions, providing the clue to understanding local plasma dynamics and heating.

Aims: We seek to illustrate the analysis problems posed by RHESSI spacecraft motion by studying persistent instrumental oscillations found in the lightcurves measured by RHESSI's X-ray detectors in the 6-12 keV and 12-25 keV energy range during the decay phase of the flares of 2004 November 4 and 6. Methods: The various motions of the RHESSI spacecraft which may contribute to the manifestation of oscillations are studied. The response of each detector in turn is also investigated. Results: We find that on 2004 November 6 the observed oscillations correspond to the nutation period of the RHESSI instrument. These oscillations are also of greatest amplitude for detector 5, while in the lightcurves of many other detectors the oscillations are small or undetectable. We also find that the variation in detector pointing is much larger during this flare than the counterexample of 2004 November 4. Conclusions: Sufficiently large nutation motions of the RHESSI spacecraft lead to clearly observable oscillations in count rates, posing a significant hazard for data analysis. This issue is particularly problematic for detector 5 due to its design characteristics. Dynamic correction of the RHESSI counts, accounting for the livetime, data gaps, and the transmission of the bi-grid collimator of each detector, is required to overcome this issue. These corrections should be applied to all future oscillation studies.

Electromagnetic (radio, visible-light, UV, EUV, X-ray and gamma-ray) emission generated by solar and stellar flares often contains pronounced quasi-periodic pulsations (QPP). Physical mechanisms responsible for the generation of long-period QPP (with the periods longer than one second) are likely to be associated with MHD processes. The observed modulation depths, periods and anharmonicity of QPP suggest that they can be linked with some kind of MHD auto-oscillations, e.g. an oscillatory regime of magnetic reconnection. Such regimes, of both spontaneous and induced nature, have been observed in resistive-MHD numerical simulations. The oscillations are essentially nonlinear and non-stationary. We demonstrate that a promising novel method for their analysis is the Empirical Mode Decomposition technique.

A phenomenological relationship between oscillations in a sunspot and quasi-periodic pulsations in flaring energy releases at an active region above the sunspot, is established. The analysis of the microwave emission recorded by the Nobeyama Radioheliograph at 17 GHz shows a gradual increase in the power of the 3-min oscillation train in the sunspot associated with AR 10756 before flares in this active region. The flaring light curves are found to be bursty with a period of 3 min. Our analysis of the spatial distribution of the 3-min oscillation power implies that the oscillations follow from sunspots along coronal loops towards the flaring site. It is proposed that quasi-periodic pulsations in the flaring energy releases can be triggered by 3-min slow magnetoacoustic waves leaking from sunspots.

A web-based, interactive system for the remote processing of imaging data sets (i.e., EUV, X-ray and microwave) and the automated interactive detection of wave and oscillatory phenomena in the solar atmosphere is presented.The system targets localised, but spatially resolved, phenomena, such as kink, sausage, and longitudinal propagating and standing waves. The system implements the methods of Periodmapping for pre-analysis, and Pixelised Wavelet Filtering for detailed analysis of the imaging data cubes. The system is implemented on the dedicated data processing server http://pwf.iszf.irk.ru, which is situated at the Institute of Solar-Terrestrial Physics, Irkutsk, Russia. The input data in the .sav, .fits or .txt formats can be submitted via the local and/or global network (the Internet). The output data can be in the png, jpeg and binary formats, on the user's request. The output data are periodmaps; narrowband amplitude, power, phase and correlation maps of the wave's sources at significant harmonics and in the chosen spectral intervals, and mpeg-movies of their evolution. The system was tested by the analysis throughout the EUV and microwave emission from the active region NOAA 10756 on 4 May 2005 observed with TRACE and the Nobeyama Radioheliograph. The similarity of the spatial localisation of three-minute propagating waves, near the footpoint of locally open magnetic field lines determined by the potential field extrapolation, in both the transition region and the corona was established. In the transition region the growth of the three-minute amplitude was found to be accompanied by the decrease in the line of sight angle to the wave propagation direction.

In the recent work of Nakariakov et al. (2004), it has been shown that the time dependences of density and velocity in a flaring loop contain pronounced quasi-harmonic oscillations associated with the 2nd harmonic of a standing slow magnetoacoustic wave. That model used a symmetric heating function (heat deposition was strictly at the apex). This left outstanding questions: A) is the generation of the 2nd harmonic a consequence of the fact that the heating function was symmetric? B) Would the generation of these oscillations occur if we break symmetry? C) What is the spectrum of these oscillations? Is it consistent with a 2nd spatial harmonic? The present work (and partly Tsiklauri et al. (2004b)) attempts to answer these important outstanding questions. Namely, we investigate the physical nature of these oscillations in greater detail: we study their spectrum (using periodogram technique) and how heat positioning affects the mode excitation. We found that excitation of such oscillations is practically independent of location of the heat deposition in the loop. Because of the change of the background temperature and density, the phase shift between the density and velocity perturbations is not exactly a quarter of the period, it varies along the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating. We also were able to model successfully SUMER oscillations observed in hot coronal loops.

Low-frequency longitudinal oscillations of a flaring coronal loop are studied numerically. In the recent work of Nakariakov et al., Astron. and Astrophys., 414, L25-L28 (2004) it has been shown that the time dependences of density and velocity in a flaring loop contain well-pronounced quasi-harmonic oscillations associated with a 2nd harmonic of a standing slow magnetoacoustic wave. In this work we investigate physical nature of these oscillations in greater detail, namely, we study their spectrum (using periodogram technique) and how does heat positioning affects the mode excitation. We found that excitation of such oscillations practically independent of the positioning of the heat deposition in the loop. Because of the change of the background temperature and density, the phase shift between the density and velocity perturbations is not exactly equal to the quarter of the period, it varies along the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating. KEYWORDS -- Sun: flares -- Sun: oscillations -- Sun: Corona -- Stars: flare -- Stars: oscillations -- Stars: coronae

A 1-D loop radiative hydrodynamic model that incorporates the effects of gravitational stratification, heat conduction, radiative losses, external heat input, presence of helium, and Braginskii viscosity is used to simulate elementary flare loops. The physical parameters for the input are taken from observations of the Bastille-Day flare of 2000 July 14. The present analysis shows that: (a) The obtained maximum values of the electron density can be considerably higher (4.2 X 10^{11} cm^{-3} or more) in the case of footpoint heating than in the case of apex heating (2.5 X 10^{11} cm^{-3}). (b) The average cooling time after the flare peak takes less time in the case of footpoint heating than in the case of apex heating. (c) The peak apex temperatures are significantly lower (by about 10 MK) for the case of footpoint heating than for apex heating (for the same average loop temperature of about 30 MK). This characteristic allows us to discriminate between different heating positioning. (d) In both cases (of apex and footpoint heating), the maximum obtained apex temperature T^{max} is practically independent of the heating duration sigma_{t}, but scales directly with the heating rate E_{H0}. (e) The maximum obtained densities at the loop apex, n_e^{max}, increase with the heating rate E_{H0} and heating duration sigma_{t} for both footpoint and apex heating. In Paper II we will use the outputs of these hydrodynamic simulations, which cover a wide range of the parameter space of heating rates and durations, as an input for forward-fitting of the multi-loop arcade of the Bastille-day flare. KEYWORDS: Sun: Flares -- Sun: Activity -- Sun: Corona

Evolution of a coronal loop in response to an impulsive energy release is numerically modelled. It is shown that the loop density evolution curves exhibit quasi-periodic perturbations with the periods given approximately by the ratio of the loop length to the average sound speed, associated with the second standing harmonics of an acoustic wave. The density perturbations have a maximum near the loop apex. The corresponding field-aligned flows have a node near the apex. We suggest that the quasi-periodic pulsations with periods in the range 10-300 s, frequently observed in flaring coronal loops in the radio, visible light and X-ray bands, may be produced by the second standing harmonic of the acoustic mode.

A recent study by Nakariakov et al. pointed out that the dispersion relation of MHD sausage mode oscillations has been incorrectly applied to coronal loops, neglecting the highly dispersive nature of the phase speed and the long-wavelength cutoff of the wave number. In the light of these new insights we revisit previous observations that have been interpreted in terms of MHD sausage mode oscillations in coronal loops and come to the following conclusions: (1) Fast sausage MHD mode oscillations require such a high electron density imposed by the wave number cutoff that they can only occur in flare loops; (2) In the previously reported radio observations ($\nu \approx 100$ MHz to 1 GHz) with periods of $P\approx 0.5-5$ s, the fast sausage MHD mode oscillation is likely to be confined to a small segment (corresponding to a high harmonic node) near the apex of the loop, rather than involving a global oscillation over the entire loop length. The recent microwave and soft X-ray observations of fast periods ($P\approx 6-17$ s) by Asai et al. and Melnikov et al., however, are consistent with fast sausage MHD oscillations at the fundamental harmonic.

We investigate the distribution of fluctuations in solar irradiance when integrated over the full disk, obtained using extreme ultraviolet/soft X-ray observations from the SOHO CELIAS/SEM instrument. This time series sums over both the contributions of single distinguishable flares, and of many other processes. By detrending we select events with timescales of less than a few hours such as waves, slow flows, and CMEs. The statistics generated by this range of phenomena can be characterised by power-law-tailed distributions. We show that (i) during the high-activity period 2000 Jan-June the tail exponent a(T)=1.5 +/- 0.1; (ii) during the low-activity period 1996 Jan-June a(T)=3.0 +/- 0.2; and (iii) in general a(T) decreases with increasing activity.

Phase mixing of a three dimensional magnetohydrodynamic (MHD) pulse is studied in the compressive, three-dimensional (without an ignorable coordinate) regime. It is shown that the efficiency of decay of an Alfvenic part of a compressible MHD pulse is related linearly to the degree of non-planeness of the pulse in the homogeneous transverse direction. In the developed stage of phase mixing (for large times), coupling to its compressive part does not alter the power-law decay of an Alfvenic part of a compressible MHD pulse. The same applies to the dependence upon the resistivity of the Alfvenic part of the pulse. All this implies that the dynamics of Alfven waves can still be qualitatively understood in terms of the previous 2.5D models. Thus, the phase mixing remains a relevant paradigm for the coronal heating applications in the realistic 3D geometry and compressive plasma.

We investigate the interaction of a plane, linearly polarized, Alfvenic pulse with a one-dimensional, perpendicular to the magnetic field, plasma density inhomogeneity in the strongly nonlinear regime. Our numerical study of the full MHD equations shows that: (i) Plasma density inhomogeneity substantially enhances (by about a factor of 2) the generation of longitudinal compressive waves. (ii) Attained maximal values of the generated transverse compressive perturbations are insensitive to the strength of the plasma density inhomogeneity, plasma beta and the initial amplitude of the Alfven wave. Typically, they reach about 40% of the initial Alfven wave amplitude. (iii) Attained maximal values of the generated relative density perturbations are within 20-40% for 0.5 <= beta <= 2.0. They depend upon plasma beta strongly; and scale almost linearly with the initial Alfven wave amplitude.

Kink modes of solar coronal structures, perturbing the loop in the direction along the line-of-sight (LOS), can be observed as emission intensity disturbances propagating along the loop provided the angle between the LOS and the structure is not ninety degrees. The effect is based upon the change of the column depth of the loop (along the LOS) by the wave. The observed amplitude of the emission intensity variations can be larger than the actual amplitude of the wave by a factor of two and there is an optimal angle maximizing the observed amplitude. For other angles this effect can also attenuate the observed wave amplitude. The observed amplitude depends upon the ratio of the wave length of kink perturbations to the width of the structure and on the angle between the LOS and the axis of the structure. Sausage modes are always affected negatively from the observational point of view, as the observed amplitude is always less than the actual one. This effect should be taken into account in the interpretation of wave phenomena observed in the corona with space-borne and ground-based imaging telescopes.

Interaction of impulsively generated MHD waves with a one-dimensional plasma inhomogeneity, transverse to the magnetic field, is considered in the three-dimensional regime. Because of the transverse inhomogeneity, MHD fluctuations, even if they do not include initially any density perturbation, evolve toward states where the compressible components tend to become predominant. The propagating MHD pulse asymptotically reaches a quasi-steady state with the final levels of density perturbation weakly depending on the degree of non-planeness of the pulse in the homogeneous transverse direction and somewhat stronger depending on plasma beta. Our study demonstrates the necessity of incorporation of compressible and 3D effects in theory of Alfven wave phase mixing. However, as far as the dynamics of weakly non-plane Alfven waves is concerned it can still be qualitatively understood in terms of the previous 2.5D models.

We investigate the interaction of magnetohydrodynamic waves with plasma density inhomogeneities. Our numerical study of the full MHD equations shows that: (A) Plasma density inhomogeneities are a source of non-linear generation of transverse compressive waves by a plane Alfven wave, and substantially enhance (by about a factor of 2) the generation of the longitudinal compressive waves. (B) Attained maximal values of the generated transverse compressive waves are insensitive to the strength of the plasma density inhomogeneity and the initial amplitude of the Alfven wave. (C) Efficiency of the generation depends weakly upon the plasma beta parameter. The maximum generated amplitude of transverse compressive wave, up to 32% of the initial Alfven wave amplitude, is reached for about beta=0.5. The results obtained demonstrate that plasma inhomogeneities enhance the efficiency of the non-linear wave coupling.

A new mechanism for the enhanced generation of compressible fluctuations by Alfven waves is presented. A strongly nonlinear regime of Alfven wave phase-mixing is numerically simulated in a one-dimensionally inhomogeneous plasma of finite temperature. It is found that the inhomogeneity of the medium determines the efficiency of nonlinear excitation of magnetoacoustic waves. The level of the compressible fluctuations is found to be higher (up to the factor of two) in inhomogeneous regions. The amplitude of the generated magnetoacoustic wave can reach up to 30% of the source Alfven wave amplitude, and this value is practically independent of the Alfven wave amplitude and the steepness of Alfven speed profile. The highest amplitudes of compressible disturbances are reached in plasmas with beta of about 0.5. The further growth of the amplitude of compressible fluctuations is depressed by saturation.

The interaction of a weakly nonlinear Alfv\'enic pulse with an Alfv\'en speed inhomogeneity in the direction perpendicular to the magnetic field is investigated. Identically to the phase mixing experienced by a harmonic Alfv\'en wave, sharp transverse gradients are generated in the pulse by the inhomogeneity. In the initial stage of the evolution of an initially plane Alfv\'enic pulse, the transverse gradients efficiently generate transversely propagating fast magnetoacoustic waves. However, high resolution full MHD numerical simulations of the developed stage of the pulse evolution show that the generation saturates due to destructive wave interference. It is shown that the weakly non-linear description of the generated fast magnetoacoustic wave is well described by the driven wave equation proposed in Nakariakov, Roberts & Murawski (1997), and a simple numerical code (2D MacCromack), which solves it with minimal CPU resources produces identical results to those obtained from the full MHD code ({\it Lare2d}, Arber et al. 2001). A parametric study of the phenomenon is undertaken, showing that, contrary to one's expectations, steeper inhomogenities of the Alfv\'en speed do not produce higher saturation levels of the fast wave generation. There is a certain optimal gradient of the inhomogeneity that ensures the maximal efficiency of the fast wave generation.

A model interpreting variations of EUV brightness upward propagating in solar coronal loops as slow magnetoacoustic waves is developed. A loop is considered to have a non-zero plane inclination angle and offset of circular loop centre from the baseline. The model also incorporates effects of dissipation and gravitational stratification. A linear evolutionary equation is derived and applied to investigation of slow wave dynamics. Both the non-zero plane inclination and the offset decrease the growth of the wave amplitude due to stratification. It is shown that wide spectrum slow magnetoacoustic waves, consistent with currently available observations in the low frequency part of the spectrum, can provide the rate of heat deposition sufficient to heat the loop. In this scenario, the heat would be deposited near the loop footpoints which agrees with the current observational data.

The evolution of quasi-isentropic magnetohydrodynamic waves of small but finite amplitude in an optically thin plasma is analyzed. The plasma is assumed to be initially homogeneous, in thermal equilibrium and with a straight and homogeneous magnetic field frozen in. Depending on the particular form of the heating/cooling function, the plasma may act as a dissipative or active medium for magnetoacoustic waves, while Alfven waves are not directly affected. An evolutionary equation for fast and slow magnetoacoustic waves in the single wave limit, has been derived and solved, allowing us to analyse the wave modification by competition of weakly nonlinear and quasi-isentropic effects. It was shown that the sign of the quasi-isentropic term determines the scenario of the evolution, either dissipative or active. In the dissipative case, when the plasma is first order isentropically stable the magnetoacoustic waves are damped and the time for shock wave formation is delayed. However, in the active case when the plasma is isentropically overstable, the wave amplitude grows, the strength of the shock increases and the breaking time decreases. The magnitude of the above effects depends upon the angle between the wave vector and the magnetic field. For hot (T > 10^4 K) atomic plasmas with solar abundances either in the interstellar medium or in the solar atmosphere, as well as for the cold (T < 10^3 K) ISM molecular gas, the range of temperature where the plasma is isentropically unstable and the corresponding time and length-scale for wave breaking have been found.