Aims. We use high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the SST to study the magnetic field in solar spicules. Methods. The equations that result from the application of the Weak Field Approximation (WFA) to the radiative transfer equations are used to infer the LOS component of the magnetic field (BLOS). Two restrictive conditions are imposed on the Stokes I and V profiles at each pixel before they can be used in a Bayesian inversion to compute its BLOS. Results. The LOS magnetic field component has been inferred in six data sets totalling 448 spectral scans in the Ca ii 8542 A line and containing both active region and quiet Sun areas, with values of hundreds of G being abundantly inferred. There seems to be no difference, from the statistical point of view, between the magnetic field strength of spicules in the quiet Sun or near an active region. On the other hand, the BLOS distributions present smaller values on the disk than off-limb, a fact that can be explained by the effect of superposition on the chromosphere of on-disk structures. We find that spicules in the vicinity of a sunspot have a magnetic field polarity (i.e. north or south) equal to that of the sunspot. This paper also contains an analysis of the effect of off-limb overlapping structures on the observed Stokes I and V parameters and the BLOS obtained from the WFA. It is found that this value is equal to or smaller than the largest LOS magnetic field components of the two structures. In addition, using random BLOS, Doppler velocities and line intensities of these two structures leads in ~ 50% of the cases to Stokes I and V parameters unsuitable to be used with the WFA. Conclusions. Our results present a scarcity of LOS magnetic field components smaller than some 50 G, which must not be taken as evidence against the existence of these magnetic field strengths in spicules.

Flare ribbons are bright manifestations of flare energy dissipation in the lower solar atmosphere. For the first time, we report on high-resolution imaging spectroscopy observations of flare ribbons situated off-limb in the H$\beta$ and Ca II 8542 {\AA} lines and make a detailed comparison with radiative hydrodynamic simulations. Observations of the X8.2-class solar flare SOL2017-09-10T16:06 UT obtained with the Swedish Solar Telescope reveal bright horizontal emission layers in H$\beta$ line wing images located near the footpoints of the flare loops. The apparent separation between the ribbon observed in the H$\beta$ wing and the nominal photospheric limb is about 300 - 500 km. The Ca II 8542 {\AA} line wing images show much fainter ribbon emissions located right on the edge of the limb, without clear separation from the limb. RADYN models are used to investigate synthetic spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. The simulations show that, towards the limb, where the line of sight is substantially oblique with respect to the vertical direction, the flaring atmosphere model reproduces the high contrast of the off-limb H$\beta$ ribbons and their significant elevation above the photosphere. The ribbons in the Ca II 8542 {\AA} line wing images are located deeper in the lower solar atmosphere with a lower contrast. A comparison of the height deposition of electron beam energy and the intensity contribution function shows that the H$\beta$ line wing intensities can be an useful tracer of flare energy deposition in the lower solar atmosphere

Aims. We aim to improve the understanding of the physical mechanisms behind the slower than free-fall motion and the two-stage evolution (an initial phase of acceleration followed by an almost constant velocity phase) detected in coronal rain events. Methods. Using the Mancha3D code, we solve the 2D ideal MHD equations. We represent the solar corona as an isothermal vertically stratified atmosphere with a uniform vertical magnetic field and the plasma condensation as a density enhancement described by a 2D Gaussian profile. We analyse the temporal evolution of the descending plasma and study its dependence on parameters such as density and magnetic field strength. Results. We confirm previous findings that the pressure gradient is the main force that opposes the action of gravity and slows down the blob descent and that larger densities require larger pressure gradients to reach the constant speed phase. We find that the shape of a condensation with a horizontal variation of density is distorted as it falls, due to the denser parts of the blob falling faster than the lighter ones. This is explained by the fact that the duration of the initial acceleration phase, and therefore the maximum falling speed attained by the plasma, increases with the ratio of blob to coronal density. We also find that the magnetic field plays a fundamental role in the evolution of the descending condensations. A strong enough magnetic field (greater than 10 G in our simulations) forces each plasma element to follow the path given by a particular field line, which allows to describe the evolution of each vertical slice of the blob in terms of 1D dynamics, without influence of the adjacent slices. In addition, under the typical conditions of the coronal rain events, the magnetic field prevents the development of the Kelvin-Helmholtz instability.

Some high-resolution observations have revealed that the active-region solar corona is filled with myriads of thin strands even in apparently uniform regions with no resolved loops. This fine structure can host collective oscillations involving a large portion of the corona due to the coupling of the motions of the neighbouring strands. We study these vibrations and the possible observational effects. Here we theoretically investigate the collective oscillations inherent to the fine structure of the corona. We have called them fundamental vibrations because they cannot exist in a uniform medium. We use the T-matrix technique to find the normal modes of random arrangements of parallel strands. We consider an increasing number of tubes to understand the vibrations of a huge number of tubes of a large portion of the corona. We additionally generate synthetic time-distance Doppler and line broadening diagrams of the vibrations of a coronal region to compare with observations. We have found that the fundamental vibrations are in the form of clusters of tubes where not all the tubes participate in the collective mode. The periods are distributed over a wide band of values. The width of the band increases with the number of strands but rapidly reaches an approximately constant value. The frequency band associated with the fine structure of the corona depends on the minimum separation between strands. We have found that the coupling between the strands is of large extent. The synthetic Dopplergrams and line-broadening maps show signatures of collective vibrations, not present in the case of purely random individual kink vibrations. We conclude that the fundamental vibrations of the corona can contribute to the energy budget of the corona and they may have an observational signature.

The purpose of this work is to develop a procedure to obtain the normal modes of a coronal loop from time-dependent numerical simulations with the aim of better understanding observed transverse loop oscillations. To achieve this goal, in this paper we present a new method and test its performance with a problem for which the normal modes can be computed analytically. In a follow-up paper, the application to the simulations of \citet{rial2013} is tackled. The method proceeds iteratively and at each step consists of (i) a time-dependent numerical simulation followed by (ii) the Complex Empirical Orthogonal Function (CEOF) analysis of the simulation results. The CEOF analysis provides an approximation to the normal mode eigenfunctions that can be used to set up the initial conditions for the numerical simulation of the following iteration, in which an improved normal mode approximation is obtained. The iterative process is stopped once the global difference between successive approximate eigenfunctions is below a prescribed threshold. The equilibrium used in this paper contains material discontinuities that result in one eigenfunction with a jump across these discontinuities and two eigenfunctions whose normal derivatives are discontinuous there. After 6 iterations, the approximation to the frequency and eigenfunctions are accurate to $\lesssim 0.7$\% except for the eigenfunction with discontinuities, which displays a much larger error at these positions.

Annual oscillations have been detected in many indices of solar activity during many cycles. Recent multi spacecraft observations of coronal bright points revealed slow retrograde toroidal phase drift (with the speed of 3 m/s of 1 yr oscillations, which naturally suggested their connection with Rossby-type waves in the interior. We have studied from a theoretical point of view the dynamics of global magneto-Kelvin and magneto-Rossby waves in the solar tachocline with toroidal magnetic field. Using spherical coordinates, the dispersion relations of the waves and latitudinal structure of solutions were obtained analytically. We have also obtained the spectrum of unstable magneto-Rossby wave harmonics in the presence of the latitudinal differential rotation. Estimated periods and phase speeds show that the magneto-Rossby waves rather than the Kelvin waves match with the observations of 1 yr oscillations. On the other hand, Morlet wavelet analysis of Greenwich Royal Observatory sunspot areas for the solar cycle 23 has revealed multiple periodicities with periods of 450-460 days, 370-380 days, 310-320 days, 240-270 days, and 150-175 days in hemispheric and full disk data. Comparison of theoretical results with the observations allow us to conclude that the global magneto-Kelvin waves in the upper overshoot tachocline may be responsible for the periodicity of 450-460 days (1.3 yrs), while the remaining periods can be connected with different harmonics of global fast magneto-Rossby waves.

Here we report on the unique observation of flaring coronal loops at the solar limb using high resolution imaging spectropolarimetry from the Swedish 1-meter Solar Telescope. The vantage position, orientation and nature of the chromospheric material that filled the flare loops allowed us to determine their magnetic field with unprecedented accuracy using the weak-field approximation method. Our analysis reveals coronal magnetic field strengths as high as 350 Gauss at heights up to 25 Mm above the solar limb. These measurements are substantially higher than a number of previous estimates and may have considerable implications for our current understanding of the extended solar atmosphere.

In the solar atmosphere, Alfv\'en waves are believed to play an important role in the transfer of energy from the photosphere to the corona and solar wind, and in the heating of the chromosphere. We perform numerical computations to investigate energy transport and dissipation associated with torsional Alfv\'en waves propagating in magnetic flux tubes that expand from the photosphere to the corona in quiet-Sun conditions. We place a broadband driver at the photosphere that injects a wave energy flux of $10^7$ erg cm$^{-2}$ s$^{-1}$ and consider Ohm's magnetic diffusion and ion-neutral collisions as dissipation mechanisms. We find that only a small fraction of the driven flux, $\sim 10^5$ erg cm$^{-2}$ s$^{-1}$, is able to reach coronal heights, but it may be sufficient to partly compensate the total coronal energy loss. The frequency of maximal transmittance is $\sim 5$ mHz for a photospheric field strength of 1 kG and is shifted to smaller/larger frequencies for weaker/stronger fields. Lower frequencies are reflected at the transition region, while higher frequencies are dissipated producing enough heat to balance chromospheric radiative losses. Heating in the low and middle chromosphere is due to Ohmic dissipation, while ion-neutral friction dominates in the high chromosphere. Ohmic diffusion is enhanced by phase mixing because of the expansion of the magnetic field. This effect has the important consequence of increasing the chromospheric dissipation and, therefore, reducing the energy flux that reaches the corona. We provide empirical fits of the transmission coefficient that could be used as input for coronal models.

We report on the discovery of periodic coronal rain in an off-limb sequence of {\it Solar Dynamics Observatory}/Atmospheric Imaging Assembly images. The showers are co-spatial and in phase with periodic (6.6~hr) intensity pulsations of coronal loops of the sort described by Auchere et al. (2014) and Froment et al. (2015, 2017). These new observations make possible a unified description of both phenomena. Coronal rain and periodic intensity pulsations of loops are two manifestations of the same physical process: evaporation / condensation cycles resulting from a state of thermal nonequilibrium (TNE). The fluctuations around coronal temperatures produce the intensity pulsations of loops, and rain falls along their legs if thermal runaway cools the periodic condensations down and below transition-region (TR) temperatures. This scenario is in line with the predictions of numerical models of quasi-steadily and footpoint heated loops. The presence of coronal rain -- albeit non-periodic -- in several other structures within the studied field of view implies that this type of heating is at play on a large scale.

Rieger-type periodicity has been detected in different activity indices over many solar cycles. It was recently shown that the periodicity correlates with solar activity having a shorter period during stronger cycles. Solar activity level is generally asymmetric between northern and southern hemispheres, which could suggest the presence of a similar behavior in the Rieger-type periodicity. We analyse the sunspot area/number and the total magnetic flux data for northern and southern hemispheres during solar cycles 19-23 which had remarkable north-south asymmetry. Using wavelet analysis of sunspot area and number during the north-dominated cycles (19-20) we obtained the periodicity of 160-165 days in the stronger northern hemisphere and 180-190 days in the weaker southern hemisphere. On the other hand, south-dominated cycles (21-23) display the periodicity of 155-160 days in the stronger southern hemisphere and 175-188 days in the weaker northern hemisphere. Therefore, the Rieger-type periodicity has the north-south asymmetry in sunspot area/number data during solar cycles with strong hemispheric asymmetry. We suggest that the periodicity is caused by magnetic Rossby waves in the internal dynamo layer. Using the dispersion relation of magnetic Rossby waves and observed Rieger periodicity we estimated the magnetic field strength in the layer as 45-49 kG in more active hemispheres (north during the cycles 19-20 and south during the cycles 21-23) and 33-40 kG in weaker hemispheres. The estimated difference in the hemispheric field strength is around 10 kG, which provides a challenge for dynamo models. Total magnetic flux data during the cycle 20-23 reveals no clear north-south asymmetry which needs to be explained in the future.

Ellerman bombs (EBs) have been widely studied over the past two decades; however, only recently have counterparts of these events been observed in the quiet-Sun. The aim of this article is to further understand small-scale quiet-Sun Ellerman-like brightenings (QSEBs) through research into their spectral signatures, including investigating whether the hot signatures associated with some EBs are also visible co-spatial to any QSEBs. We combine H$\alpha$ and Ca II $8542$ \AA\ line scans at the solar limb with spectral and imaging data sampled by the Interface Region Imaging Spectrograph (IRIS). Twenty one QSEBs were identified with average lifetimes, lengths, and widths measured to be around $120$ s, $0.63$", and $0.35$", respectively. Three of these QSEBs displayed clear repetitive flaring through their lifetimes, comparable to the behaviour of EBs in Active Regions (ARs). Two QSEBs in this sample occurred co-spatial with increased emission in SDO/AIA $1600$ \AA\ and IRIS slit-jaw imager $1400$ \AA\ data, however, these intensity increases were smaller than reported co-spatial to EBs. One QSEB was also sampled by the IRIS slit during its lifetime, displaying increases in intensity in the Si IV $1393$ \AA\ and Si IV $1403$ \AA\ cores as well as the C II and Mg II line wings, analogous to IRIS bursts (IBs). Using RADYN simulations, we are unable to reproduce the observed QSEB H$\alpha$ and Ca II $8542$ \AA\ line profiles leaving the question of the temperature stratification of QSEBs open. Our results imply that some QSEBs could be heated to Transition Region temperatures, suggesting that IB profiles should be observed throughout the quiet-Sun.

We investigate the nature of dissipative instability appearing in a prominence planar thread filled with partially ionised plasma in the incompressible limit. The importance of partial ionisation is investigated in terms of the ionisation factor and wavelength of waves propagating in the slab. To highlight the role of partial ionisation, we have constructed models describing various situations we can meet in solar prominence fine structure. Matching the solutions for the transversal component of the velocity and total pressure at the interfaces between the prominence slab and surrounding plasmas, we derived a dispersion relation whose imaginary part describes the evolution of the instability. Results are obtained in the limit of weak dissipation. We have investigated the appearance of instabilities in prominence dark plumes using single and two-fluid approximations. We show that dissipative instabilities appear for flow speeds that are less than the Kelvin-Helmholtz instability threshold. The onset of instability is determined by the equilibrium flow strength, the ionisation factor of the plasma, the wavelength of waves and the ion-neutral collisional rate. For a given wavelength and for ionisation degrees closer to a neutral gas, the propagating waves become unstable for a narrow band of flow speeds, meaning that neutrals have a stabilising effect. Our results show that the partially ionised plasma describing prominence dark plumes becomes unstable only in a two-fluid (charged particles-neutrals) model, that is for periods that are smaller than the ion-neutral collision time. The present study improves our understanding of stability of solar prominences and the role of partial ionisation in destabilising the plasma. We show the necessity of two-fluid approximation when discussing the nature of instabilities: waves in a single fluid approximation show a great deal of stability.

Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the gravity projected along the dipped field lines of the magnetic structure, best explains these oscillations. However, several factors can influence the longitudinal oscillations, potentially invalidating the pendulum model. The aim of this work is to study the influence of large-scale variations in the magnetic field strength along the field lines, i.e., variations of the cross-sectional area along the flux tubes supporting prominence threads. We studied the normal modes of several flux tube configurations, using linear perturbation analysis, to assess the influence of different geometrical parameters on the oscillation properties. We found that the influence of the symmetric and asymmetric expansion factors on longitudinal oscillations is small.}{We conclude that the longitudinal oscillations are not significantly influenced by variations of the cross-section of the flux tubes, validating the pendulum model in this context.

Observations have shown that magnetohydrodynamic waves over a large frequency range are ubiquitous in solar prominences. The waves are probably driven by photospheric motions and may transport energy up to prominences suspended in the corona. Dissipation of wave energy can lead to heating of the cool prominence plasma, so contributing to the local energy balance within the prominence. Here we discuss the role of Alfv\'en wave dissipation as a heating mechanism for the prominence plasma. We consider a slab-like quiescent prominence model with a transverse magnetic field embedded in the solar corona. The prominence medium is modelled as a partially ionized plasma composed of a charged ion-electron single fluid and two separate neutral fluids corresponding to neutral hydrogen and neutral helium. Friction between the three fluids acts as a dissipative mechanism for the waves. The heating caused by externally-driven Alfv\'en waves incident on the prominence slab is analytically explored. We find that the dense prominence slab acts as a resonant cavity for the waves. The fraction of incident wave energy that is channelled into the slab strongly depends upon the wave period, $P$. Using typical prominence conditions, we obtain that wave energy trapping and associated heating are negligible when $P \gtrsim 100$ s, so that it is unlikely that those waves have a relevant influence on prominence energetics. When $1$ s $\lesssim P \lesssim 100$ s the energy absorption into the slab shows several sharp and narrow peaks, that can reach up to 100%, when the incident wave frequency matches a cavity resonance of the slab. Wave heating is enhanced at those resonant frequencies. Conversely, when $P \lesssim 1$ s cavity resonances are absent, but the waves are heavily damped by the strong dissipation. We estimate that wave heating may compensate for about 10% of radiative losses of the prominence plasma.

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

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.

The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov Demoulin (1999) under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, characterized by the excitation of standing vertical oscillations. These oscillations show a strong attenuation with time produced by the mechanism of continuum damping due to the inhomogeneous transition between the prominence and solar corona. The characteristic period of the vertical oscillations does not depend strongly on the twist of the flux rope. Nonlinearity is the responsible for triggering the Kelvin-Helmholtz instability associated to the vertical oscillations and that eventually produces horizontal structures. Contrary to other configurations in which the longitudinal axis of the prominence is permeated by a perpendicular magnetic field, like in unsheared arcades, the orientation of the prominence along the flux rope axis prevents the development of Rayleigh-Taylor instabilities and therefore the appearance of vertical structuring along this axis.

Coronal rain clumps and prominence knots are dense condensations with chromospheric to transition region temperatures that fall down in the much hotter corona. Their typical speeds are in the range 30--150~km~s$^{-1}$ and of the order of 10--30~km~s$^{-1}$, respectively, i.e., they are considerably smaller than free fall velocities. These cold blobs contain a mixture of ionized and neutral material that must be dynamically coupled in order to fall together, as observed. We investigate this coupling by means of hydrodynamic simulations in which the coupling arises from the friction between ions and neutrals. The numerical simulations presented here are an extension of those of \citet{oliver2014} to the partially ionized case. We find that, although the relative drift speed between the two species is smaller than 1~m~s$^{-1}$ at the blob center, it is sufficient to produce the forces required to strongly couple charged particles and neutrals. The ionization degree has no discernible effect on the main results of our previous work for a fully ionized plasma: the condensation has an initial acceleration phase followed by a period with roughly constant velocity and, in addition, the maximum descending speed is clearly correlated with the ratio of initial blob to environment density.

Long-term records of sunspot number and concentrations of cosmogenic radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's magnetic activity over hundreds and thousands of years. We identify several clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100 years. We found that the periods of the first five spherical harmonics of the slow magnetic Rossby mode in the presence of a steady toroidal magnetic field of 1200-1300 G in the lower tachocline are in perfect agreement with the time scales of observed variations. The steady toroidal magnetic field can be generated in the lower tachocline either due to the steady dynamo magnetic field for low magnetic diffusivity or due to the action of the latitudinal differential rotation on the weak poloidal primordial magnetic field, which penetrates from the radiative interior. The slow magnetic Rossby waves lead to variations of the steady toroidal magnetic field in the lower tachocline, which modulate the dynamo magnetic field and consequently the solar cycle strength. This result constitutes a key point for long-term prediction of the cycle strength. According to our model, the next deep minimum in solar activity is expected during the first half of this century.

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

We investigate the nature of dissipative instability at the boundary (seen here as tangential discontinuity) between the viscous corona and the partially ionised prominence plasma in the incompressible limit. The importance of the partial ionisation is investigated in terms of the ionisation fraction. Matching the solutions for the transversal component of the velocity and total pressure at the interface between the prominence and coronal plasmas, we derive a dispersion relation whose imaginary part describes the evolution of the instability. Results are obtained in the limit of weak dissipation. Using simple analytical methods, we show that dissipative instabilities appear for flow speeds that are lower than the Kelvin-Helmholtz instability threshold. While viscosity tends to destabilise the plasma, the effect of partial ionisation (through the Cowling resistivity) will act towards stabilising the interface. For ionisation degrees closer to a neutral gas the interface will be unstable for larger values of equilibrium flow. The same principle is assumed when studying the appearance of instability at the interface between prominences and dark plumes. The unstable mode appearing in this case has a very small growth rate and dissipative instability cannot explain the appearance of flows in plumes. The present study improves our understanding of the complexity of dynamical processes at the interface of solar prominences and solar corona, and the role partial ionisation can have on the stability of the plasma. Our results clearly show that the problem of partial ionisation introduces new aspects of plasma stability with consequences on the evolution of solar prominences.

Cool and dense prominences found in the solar atmosphere are known to be partially ionized because of their relative low temperature. In this Letter, we address the long-standing problem of how the neutral component of the plasma in prominences is supported against gravity. Using the multiple fluid approach we solve the time-dependent equations in two dimensions considering the frictional coupling between the neutral and ionized components of the magnetized plasma representative of a solar prominence embedded in a hot coronal environment. We demonstrate that given an initial density enhancement in the two fluids, representing the body of the prominence, the system is able to relax in the vicinity of magnetic dips to a stationary state in which both neutrals and ionized species are dynamically suspended above the photosphere. Two different coupling processes are considered in this study, collisions between ions and neutrals and charge exchange interactions. We find that for realistic conditions ions are essentially static while neutrals have a very small downflow velocity. The coupling between ions and neutrals is so strong at the prominence body that the behavior is similar to that of a single fluid with an effective density equal to the sum of the ion and neutral species. We also find that the charge exchange mechanism is about three times more efficient sustaining neutrals than elastic scattering of ions with neutrals.

A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we considered a transverse initial perturbation, whereas the temporal evolution of an axisymmetric one is examined here. In both papers we use a method based on Fourier integrals to solve the initial value problem. Previous studies on wave propagation in magnetic wave guides have emphasized that the wave train dispersion is influenced by the particular dependence of the group velocity on the longitudinal wavenumber. Here we also find that long initial perturbations result in low amplitude wave packets and that large values of the magnetic tube to environment density ratio yield longer wave trains. To test the detectability of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling of the perturbations must be carried out. This is left for a future work.

In this paper we present a numerical study of the time evolution of solar prominences embedded in sheared magnetic arcades. The prominence is represented by a density enhancement in a background stratified atmosphere and is connected to the photosphere through the magnetic field. By solving the ideal magnetohydrodynamic (MHD) equations in three dimensions we study the dynamics for a range of parameters representative of real prominences. Depending on the parameters considered, we find prominences that are suspended above the photosphere, i.e., detached prominences, but also configurations resembling curtain or hedgerow prominences whose material continuously connects to the photosphere. The plasma$-\beta$ is an important parameter that determines the shape of the structure. In many cases magnetic Rayleigh-Taylor (MRT) instabilities and oscillatory phenomena develop. Fingers and plumes are generated, affecting the whole prominence body and producing vertical structures in an essentially horizontal magnetic field. However, magnetic shear is able to reduce or even to suppress this instability.

The dispersion of small amplitude, impulsively excited wave trains propagating along a magnetic flux tube is investigated. The initial disturbance is a localized transverse displacement of the tube that excites a fast kink wave packet. The spatial and temporal evolution of the perturbed variables (density, plasma displacement, velocity, ...) is given by an analytical expression containing an integral that is computed numerically. We find that the dispersion of fast kink wave trains is more important for shorter initial disturbances (i.e. more concentrated in the longitudinal direction) and for larger density ratios (i.e. for larger contrasts of the tube density with respect to the environment density). This type of excitation generates a wave train whose signature at a fixed position along a coronal loop is a short event (duration ~ 20 s) in which the velocity and density oscillate very rapidly with typical periods of the order of a few seconds. The oscillatory period is not constant but gradually declines during the course of this event. Peak values of the velocity are of the order of 10 km/s and are accompanied by maximum density variations of the order of 10-15% the unperturbed loop density.

Context. Prominence oscillations have been mostly detected using Doppler velocity, although there are also claimed detections by means of periodic variations in half-width or line intensity. However, scarce observational evidence exists about simultaneous detection of oscillations in several spectral indicators. Aims. Our main aim here is to explore the relationship between spectral indicators, such as Doppler shift, line intensity, and line half-width, and the linear perturbations excited in a simple prominence model. Methods. Our equilibrium background model consists of a bounded, homogeneous slab, which is permeated by a transverse magnetic field, having prominence-like physical properties. Assuming linear perturbations, the dispersion relation for fast and slow modes has been derived, as well as the perturbations for the different physical quantities. These perturbations have been used as the input variables in a one-dimensional radiative transfer code, which calculates the full spectral profile of the hydrogen H-alpha and H-beta lines. Results. We have found that different oscillatory modes produce spectral indicator variations in different magnitudes. Detectable variations in the Doppler velocity were found for the fundamental slow mode only. Substantial variations in the H-beta line intensity were found for specific modes. Other modes lead to lower and even undetectable parameter variations. Conclusions. To perform prominence seismology, analysis of the H-alpha and H-beta spectral line parameters could be a good tool to detect and identify oscillatory modes.

Seismology of coronal loops using observations of damped transverse oscillations in combination with results from theoretical models is a tool to indirectly infer physical parameters in the solar atmospheric plasma. Existing seismology schemes based on approximations to the period and damping time of kink oscillations are often used beyond their theoretical range of applicability. These approximations assume that the variation of density across the loop is confined to a nonuniform layer much thinner than the radius of the loop, but the results of the inversion problem often do not satisfy this preliminary hypothesis. Here, we determine the accuracy of the analytic approximations to the period and damping time, and its impact on seismology estimates, when largely nonuniform loops are considered. We find that the accuracy of the approximations when used beyond their range of applicability is strongly affected by the form of the density profile across the loop, that is observationally unknown and so must be arbitrarily imposed as part of the theoretical model. The error associated with the analytic approximations can be larger than 50% even for relatively thin nonuniform layers. This error directly affects the accuracy of approximate seismology estimates compared to actual numerical inversions. In addition, assuming different density profiles can produce noncoincident intervals of the seismic variables in inversions of the same event. The ignorance about the true shape of density variation across the loop is an important source of error that may dispute the reliability of parameters seismically inferred assuming an ad hoc density profile.

Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that results in the formation of a large pressure gradient that opposes gravity. Eventually this pressure gradient becomes so large that the blob acceleration vanishes or even points upwards. Then, the blob descent is characterised by an initial acceleration phase followed by an essentially constant velocity phase. These two stages can be identified in published time-distance diagrams of coronal rain events. Both the duration of the first stage and the velocity attained by the blob increase for larger values of the ratio of blob to coronal density, for larger blob mass, and for smaller coronal temperature. Dense blobs are characterised by a detectable density growth (up to 60% in our calculations) and by a steepening of the density in their lower part, that could lead to the formation of a shock. They also emit sound waves that could be detected as small intensity changes with periods of the order of 100 s and lasting between a few and about ten periods. Finally, the curvature of the falling path is only relevant when a very dense blob falls along inclined magnetic field lines.

Solar prominence models used so far in the analysis of MHD waves in such structures are quite elementary. In this work, we calculate numerically magnetohydrostatic models in two-dimensional configurations under the presence of gravity. Our interest is in models that connect the magnetic field to the photosphere and include an overlying arcade. The method used here is based on a relaxation process and requires solving the time-dependent nonlinear ideal MHD equations. Once a prominence model is obtained, we investigate the properties of MHD waves superimposed on the structure. We concentrate on motions purely two-dimensional neglecting propagation in the ignorable direction. We demonstrate how by using different numerical tools we can determine the period of oscillation of stable waves. We find that vertical oscillations, linked to fast MHD waves, are always stable and have periods in the 4-10 min range. Longitudinal oscillations, related to slow magnetoacoustic-gravity waves, have longer periods in the range of 28-40 min. These longitudinal oscillations are strongly influenced by the gravity force and become unstable for short magnetic arcades.

Magnetohydrodynamic (MHD) waves are ubiquitously observed in the solar atmosphere. Kink waves are a type of transverse MHD waves in magnetic flux tubes that are damped due to resonant absorption. The theoretical study of kink MHD waves in solar flux tubes is usually based on the simplification that the transverse variation of density is confined to a nonuniform layer much thinner than the radius of the tube, i.e., the so-called thin boundary approximation. Here, we develop a general analytic method to compute the dispersion relation and the eigenfunctions of ideal MHD waves in pressureless flux tubes with transversely nonuniform layers of arbitrary thickness. Results for kink waves are produced and are compared with fully numerical resistive MHD eigenvalue computations in the limit of small resistivity. We find that the frequency and resonant damping rate are the same in both ideal and resistive cases. The actual results for thick nonuniform layers deviate from the behavior predicted in the thin boundary approximation and strongly depend on the shape of the nonuniform layer. The eigenfunctions in ideal MHD are very different from those in resistive MHD. The ideal eigenfunctions display a global character regardless of the thickness of the nonuniform layer, while the resistive eigenfunctions are localized around the resonance and are indistinguishable from those of ordinary resistive Alfv\'en modes. Consequently, the spatial distribution of wave energy in the ideal and resistive cases is dramatically different. This poses a fundamental theoretical problem with clear observational consequences.

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

Many solar filaments and prominences show short-lived horizontal threads lying parallel to the photosphere. In this work the possible link between Rayleigh-Taylor instabilities and thread lifetimes is investigated. This is done by calculating the eigenmodes of a thread modelled as a Cartesian slab under the presence of gravity. An analytical dispersion relation is derived using the incompressible assumption for the magnetohydrodynamic (MHD) perturbations. The system allows a mode that is always stable, independently of the value of the Alfv\'en speed in the thread. The character of this mode varies from being localised at the upper interface of the slab when the magnetic field is weak, to having a global nature and resembling the transverse kink mode when the magnetic field is strong. On the contrary, the slab model permits another mode that is unstable and localised at the lower interface when the magnetic field is weak. The growth rates of this mode can be very short, of the order of minutes for typical thread conditions. This Rayleigh-Taylor unstable mode becomes stable when the magnetic field is increased, and in the limit of strong magnetic field it is essentially a sausage magnetic mode. The gravity force might have a strong effect on the modes of oscillation of threads, depending on the value of the Alfv\'en speed. In the case of threads in quiescent filaments, where the Alfv\'en speed is presumably low, very short lifetimes are expected according to the slab model. In active region prominences, the stabilising effect of the magnetic tension might be enough to suppress the Rayleigh-Taylor instability for a wide range of wavelengths.

Given the difficulty in directly determining prominence physical parameters from observations, prominence seismology stands as an alternative method to probe the nature of these structures. We show recent examples of the application of magnetohydrodynamic (MHD) seismology techniques to infer physical parameters in prominence plasmas. They are based on the application of inversion techniques using observed periods, damping times, and plasma flow speeds of prominence thread oscillations. The contribution of Hinode to the subject has been of central importance. We show an example based on data obtained with Hinode's Solar Optical Telescope. Observations show an active region limb prominence, composed by a myriad of thin horizontal threads that flow following a path parallel to the photosphere and display synchronous vertical oscillations. The coexistence of waves and flows can be firmly established. By making use of an interpretation based on transverse MHD kink oscillations, a seismological analysis of this event is performed. It is shown that the combination of high quality Hinode observations and proper theoretical models allows flows and waves to become two useful characteristics for our understanding of the nature of solar prominences.

We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordinates and the time evolution of the initial perturbations has been studied in the zero-beta approximation. For propagation constrained to the plane of the arcade, the considered initial perturbations do not excite trapped modes of the system. This weakness of the model is overcome by the inclusion of wave propagation in the ignorable direction. The inclusion of perpendicular propagation produces two main results. First, damping by wave leakage is less efficient because the loop is able to act as a wave trap of vertical oscillations. Second, the consideration of an inhomogeneous corona enables the resonant damping of vertical oscillations and the energy transfer from the interior of the loop to the external coronal medium.

Rapidly rotating stars show short-period oscillations in magnetic activity and polar appearance of starspots. The aim of this paper is to study large-scale shallow water waves in the tachoclines of rapidly rotating stars and their connection to the periodicity and the formation of starspots at high latitudes. Shallow-water magnetohydrodynamic equations were used to study the dynamics of large-scale waves at the rapidly rotating stellar tachoclines in the presence of toroidal magnetic field. Dispersion relations and latitudinal distribution of wave modes were derived. We found that low-frequency magnetic Rossby waves tend to be located at poles, but high-frequency magnetic Poincare waves are concentrated near the equator in rapidly rotating stars. These results have important implications for the evolution of the stellar wind in young Sun-like stars. Unstable magnetic Rossby waves may lead to the local enhancement of magnetic flux at high latitudes of tachoclines in rapidly rotating stars. The enhanced magnetic flux may rise upwards owing to the magnetic buoyancy in the form of tubes and appear as starspots at polar regions. Magnetic Rossby waves may also cause observed short-term periodicity in the stellar magnetic activity.

Quasi-biennial oscillations (QBO) are frequently observed in the solar activity indices. However, no clear physical mechanism for the observed variations has been suggested so far. Here we study the stability of magnetic Rossby waves in the solar tachocline using the shallow water magnetohydrodynamic approximation. Our analysis shows that the combination of typical differential rotation and a toroidal magnetic field with a strength > 10^5 G triggers the instability of the m=1 magnetic Rossby wave harmonic with a period of 2 years. This harmonic is antisymmetric with respect to the equator and its period (and growth rate) depends on the differential rotation parameters and the magnetic field strength. The oscillations may cause a periodic magnetic flux emergence at the solar surface and consequently may lead to the observed QBO in the solar activity features. The period of QBO may change throughout the cycle, and from cycle to cycle, due to variations of the mean magnetic field and differential rotation in the tachocline.

Transverse oscillations and propagating waves are frequently observed in threads of solar prominences/filaments and have been interpreted as kink magnetohydrodynamic (MHD) modes. We investigate the spatial damping of propagating kink MHD waves in transversely nonuniform and partially ionized prominence threads. Resonant absorption and ion-neutral collisions (Cowling's diffusion) are the damping mechanisms taken into account. The dispersion relation of resonant kink waves in a partially ionized magnetic flux tube is numerically solved by considering prominence conditions. Analytical expressions of the wavelength and damping length as functions of the kink mode frequency are obtained in the Thin Tube and Thin Boundary approximations. For typically reported periods of thread oscillations, resonant absorption is an efficient mechanism for the kink mode spatial damping, while ion-neutral collisions have a minor role. Cowling's diffusion dominates both the propagation and damping for periods much shorter than those observed. Resonant absorption may explain the observed spatial damping of kink waves in prominence threads. The transverse inhomogeneity length scale of the threads can be estimated by comparing the observed wavelengths and damping lengths with the theoretically predicted values. However, the ignorance of the form of the density profile in the transversely nonuniform layer introduces inaccuracies in the determination of the inhomogeneity length scale.

We investigate standing kink magnetohydrodynamic (MHD) oscillations in a prominence fine structure modeled as a straight and cylindrical magnetic tube only partially filled with the prominence material, and with its ends fixed at two rigid walls representing the solar photosphere. The prominence plasma is partially ionized and a transverse inhomogeneous transitional layer is included between the prominence thread and the coronal medium. Thus, ion-neutral collisions and resonant absorption are the considered damping mechanisms. Approximate analytical expressions of the period, the damping time, and their ratio are derived for the fundamental mode in the thin tube and thin boundary approximations. We find that the dominant damping mechanism is resonant absorption, which provides damping ratios in agreement with the observations, whereas ion-neutral collisions are irrelevant for the damping. The values of the damping ratio are independent of both the prominence thread length and its position within the magnetic tube, and coincide with the values for a tube fully filled with the prominence plasma. The implications of our results in the context of the MHD seismology technique are discussed, pointing out that the reported short-period (2 - 10 min) and short-wavelength (700 - 8,000 km) thread oscillations may not be consistent with a standing mode interpretation and could be related to propagating waves. Finally, we show that the inversion of some prominence physical parameters, e.g., Alfv\'en speed, magnetic field strength, transverse inhomogeneity length-scale, etc., is possible using observationally determined values of the period and damping time of the oscillations along with the analytical approximations of these quantities.

We investigate the transverse oscillations of a line-tied multi-stranded coronal loop composed of several parallel cylindrical strands. First, the collective fast normal modes of the loop are found with the T-matrix theory. There is a huge quantity of normal modes with very different frequencies and a complex structure of the associated magnetic pressure perturbation and velocity field. The modes can be classified as bottom, middle, and top according to their frequencies and spatial structure. Second, the temporal evolution of the velocity and magnetic pressure perturbation after an initial disturbance are analyzed. We find complex motions of the strands. The frequency analysis reveals that these motions are a combination of low and high frequency modes. The complexity of the strand motions produces a strong modulation of the whole tube movement. We conclude that the presumed internal fine structure of a loop influences its transverse oscillations and so its transverse dynamics cannot be properly described by those of an equivalent monolithic loop.

We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent and therefore only coupled MHD fast and Alfven modes remain. Two types of numerical experiments are performed. First, the resonant wave energy transfer between a fast normal mode of the system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, the fast wave with global character transfers its energy to Alfvenic oscillations localized around a particular magnetic surface within the arcade, thus producing the damping of the initial fast MHD mode. Second, the time evolution of a localized impulsive excitation, trying to mimic a nearby coronal disturbance, is considered. In this case, the generated fast wavefront leaves its energy on several magnetic surfaces within the arcade. The system is therefore able to trap energy in the form of Alfvenic oscillations, even in the absence of a density enhancement such as that of a coronal loop. These local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of wave energy trapped by the system via wave energy conversion strongly depends on the wavelength of perturbations in the perpendicular direction, but is almost independent from the ratio of the magnetic to density scale heights.

Solar prominences are partially ionised plasmas displaying flows and oscillations. These oscillations show time and spatial damping and, commonly, have been explained in terms of magnetohydrodynamic (MHD) waves. We study the spatial damping of linear non-adiabatic MHD waves in a flowing partially ionised plasma, having prominence-like physical properties. We consider single fluid equations for a partially ionised hydrogen plasma including in the energy equation optically thin radiation, thermal conduction by electrons and neutrals, and heating. Keeping the frequency real and fixed, we have solved the obtained dispersion relations for the complex wavenumber, k, and have analysed the behaviour of the damping length, wavelength and the ratio of the damping length to the wavelength, versus period, for Alfven, fast, slow and thermal waves.

Transverse oscillations of coronal loops are often observed and have been theoretically interpreted as kink magnetohydrodynamic (MHD) modes. Numerical simulations by Terradas et al. (2008, ApJ 687, L115) suggest that shear flows generated at the loop boundary during kink oscillations could give rise to a Kelvin-Helmholtz instability (KHI). Here, we investigate the linear stage of the KHI in a cylindrical magnetic flux tube in the presence of azimuthal shear motions. We consider the basic, linearized MHD equations in the beta = 0 approximation, and apply them to a straight and homogeneous cylindrical flux tube model embedded in a coronal environment. Azimuthal shear flows with a sharp jump of the velocity at the cylinder boundary are included in the model. We obtain an analytical expression for the dispersion relation of the unstable MHD modes supported by the configuration, and compute analytical approximations of the critical velocity shear and the KHI growth rate in the thin tube limit. A parametric study of the KHI growth rates is performed by numerically solving the full dispersion relation. We find that fluting-like modes can develop a KHI in time-scales comparable to the period of kink oscillations of the flux tube. The KHI growth rates increase with the value of the azimuthal wavenumber and decrease with the longitudinal wavenumber. However, the presence of a small azimuthal component of the magnetic field can suppress the KHI. Azimuthal motions related to kink oscillations of untwisted coronal loops may trigger a KHI, but this phenomenon has not been observed to date. We propose that the azimuthal component of the magnetic field is responsible for suppressing the KHI in a stable coronal loop. The required twist is small enough to prevent the development of the pinch instability.

Apart from the 11-year solar cycle, another periodicity around 155-160 days was discovered during solar cycle 21 in high energy solar flares, and its presence in sunspot areas and strong magnetic flux has been also reported. This periodicity has an elusive and enigmatic character, since it usually appears only near the maxima of solar cycles, and seems to be related with a periodic emergence of strong magnetic flux at the solar surface. Therefore, it is probably connected with the tachocline, a thin layer located near the base of the solar convection zone, where strong dynamo magnetic field is stored. We study the dynamics of Rossby waves in the tachocline in the presence of a toroidal magnetic field and latitudinal differential rotation. Our analysis shows that the magnetic Rossby waves are generally unstable and that the growth rates are sensitive to the magnetic field strength and to the latitudinal differential rotation parameters. Variation of the differential rotation and the magnetic field strength throughout the solar cycle enhance the growth rate of a particular harmonic in the upper part of the tachocline around the maximum of the solar cycle. This harmonic is symmetric with respect to the equator and has a period of 155-160 days. A rapid increase of the wave amplitude could give place to a magnetic flux emergence leading to observed periodicities in solar activity indicators related with magnetic flux.

Prominences are partially ionized, magnetized plasmas embedded in the solar corona. Damped oscillations and propagating waves are commonly observed. These oscillations have been interpreted in terms of magnetohydrodynamic (MHD) waves. Ion-neutral collisions and non-adiabatic effects (radiation losses and thermal conduction) have been proposed as damping mechanisms. We study the effect of the presence of helium on the time damping of non-adiabatic MHD waves in a plasma composed by electrons, protons, neutral hydrogen, neutral helium (He I), and singly ionized helium (He II) in the single-fluid approximation. The dispersion relation of linear non-adiabatic MHD waves in a homogeneous, unbounded, and partially ionized prominence medium is derived. The period and the damping time of Alfven, slow, fast, and thermal waves are computed. A parametric study of the ratio of the damping time to the period with respect to the helium abundance is performed. The efficiency of ion-neutral collisions as well as thermal conduction is increased by the presence of helium. However, if realistic abundances of helium in prominences (~10%) are considered, this effect has a minor influence on the wave damping. The presence of helium can be safely neglected in studies of MHD waves in partially ionized prominence plasmas.

Transverse oscillations of solar filament and prominence threads have been frequently reported. These oscillations have the common features of being of short period (2-10 min) and being damped after a few periods. Kink magnetohydrodynamic (MHD) wave modes have been proposed as responsible for the observed oscillations, whereas resonant absorption in the Alfven continuum and ion-neutral collisions are the best candidates to be the damping mechanisms. Here, we study both analytically and numerically the time damping of kink MHD waves in a cylindrical, partially ionized filament thread embedded in a coronal environment. The thread model is composed of a straight and thin, homogeneous filament plasma, with a transverse inhomogeneous transitional layer where the plasma physical properties vary continuously from filament to coronal conditions. The magnetic field is homogeneous and parallel to the thread axis. We find that the kink mode is efficiently damped by resonant absorption for typical wavelengths of filament oscillations, the damping times being compatible with the observations. Partial ionization does not affect the process of resonant absorption, and the filament plasma ionization degree is only important for the damping for wavelengths much shorter than those observed. To our knowledge, this is the first time that the phenomenon of resonant absorption is studied in a partially ionized plasma.

From recent high resolution observations obtained with the Swedish 1 m Solar Telescope in La Palma, we detect swaying motions of individual filament threads in the plane of the sky. The oscillatory character of these motions are comparable with oscillatory Doppler signals obtained from corresponding filament threads. Simultaneous recordings of motions in the line of sight and in the plane of the sky give information about the orientation of the oscillatory plane. These oscillations are interpreted in the context of the magnetohydrodynamic theory. Kink magnetohydrodynamic waves supported by the thread body are proposed as an explanation of the observed thread oscillations. On the basis of this interpretation and by means of seismological arguments, we give an estimation of the thread Alfv\'en speed and magnetic field strength by means of seismological arguments.

Oscillations and propagating waves are commonly seen in high-resolution observations of filament threads, i.e., the fine-structures of solar filaments/prominences. Since the temperature of prominences is typically of the order of 10^4 K, the prominence plasma is only partially ionized. In this paper, we study the effect of neutrals on the wave propagation in a filament thread modeled as a partially ionized homogeneous magnetic flux tube embedded in an homogeneous and fully ionized coronal plasma. Ohmic and ambipolar magnetic diffusion are considered in the basic resistive MHD equations. We numerically compute the eigenfrequencies of kink, slow, and Alfven linear MHD modes, and obtain analytical approximations in some cases. We find that the existence of propagating modes is constrained by the presence of critical values of the longitudinal wavenumber. In particular, the lower and upper frequency cut-offs of kink and Alfven waves owe their existence to magnetic diffusion parallel and perpendicular to magnetic field lines, respectively. The slow mode only has a lower frequency cut-off, which is caused by perpendicular magnetic diffusion and is significantly affected by the ionization degree. In addition, ion-neutral collisions is the most efficient damping mechanism for short wavelengths while ohmic diffusion dominates in the long-wavelength regime.

Transverse oscillations of small amplitude are commonly seen in high-resolution observations of filament threads, i.e. the fine-structures of solar filaments/prominences, and are typically damped in a few periods. Kink wave modes supported by the thread body offer a consistent explanation of these observed oscillations. Among the proposed mechanisms to explain the kink mode damping, resonant absorption in the Alfven continuum seems to be the most efficient as it produces damping times of about 3 periods. However, for a nonzero-beta plasma and typical prominence conditions, the kink mode is also resonantly coupled to slow (or cusp) continuum modes, which could further reduce the damping time. In this Letter, we explore for the first time both analytically and numerically the effect of the slow continuum on the damping of transverse thread oscillations. The thread model is composed of a homogeneous and straight cylindrical plasma, an inhomogeneous transitional layer, and the homogeneous coronal plasma. We find that the damping of the kink mode due to the slow resonance is much less efficient than that due to the Alfven resonance.

Quiescent prominences are thin slabs of cold, dense plasma embedded in the much hotter and rarer solar corona. Although their global shape is rather irregular, they are often characterised by an internal structure consisting of a large number of thin, parallel threads piled together. Prominences often display periodic disturbances mostly observed in the Doppler displacement of spectral lines and with an amplitude typically of the order of or smaller than 2--3 km s$^{-1}$, a value which seems to be much smaller than the characteristic speeds of the prominence plasma (namely the Alfv\'en and sound velocities). Two particular features of these small amplitude prominence oscillations is that they seem to damp in a few periods and that they seem not to affect the whole prominence structure. In addition, in high spatial resolution observations, in which threads can be discerned, small amplitude oscillations appear to be clearly associated to these fine structure constituents. Prominence seismology tries to bring together the results from these observations (e.g. periods, wavelengths, damping times) and their theoretical modeling (by means of the magnetohydrodynamic theory) to gain insight into physical properties of prominences that cannot be derived from direct observation. In this paper we discuss works that have not been described in previous reviews, namely the first seismological application to solar prominences and theoretical advances on the attenuation of prominence oscillations.

We derive analytical solutions and dispersion relations of global magnetic Poincar\'e (magneto-gravity) and magnetic Rossby waves in the approximation of shallow water magnetohydrodynamics. The solutions are obtained in a rotating spherical coordinate system for strongly and weakly stable stratification separately in the presence of toroidal magnetic field. In both cases magnetic Rossby waves split into fast and slow magnetic Rossby modes. In the case of strongly stable stratification (valid in the radiative part of the tachocline) all waves are slightly affected by the layer thickness and the toroidal magnetic field, while in the case of weakly stable stratification (valid in the upper overshoot layer of the tachocline) magnetic Poincar\'e and fast magnetic Rossby waves are found to be concentrated near the solar equator, leading to equatorially trapped waves. However, slow magnetic Rossby waves tend to concentrate near the poles, leading to polar trapped waves. The frequencies of all waves are smaller in the upper weakly stable stratification region than in the lower strongly stable stratification one.

High resolution observations of solar filaments suggest the presence of groups of prominence threads, i.e. the fine-structures of prominences, which oscillate coherently (in phase). In addition, mass flows along threads have been often observed. Here, we investigate the effect of mass flows on the collective fast and slow nonadiabatic magnetoacoustic wave modes supported by systems of prominence threads. Prominence fine-structures are modeled as parallel, homogeneous and infinite cylinders embedded in a coronal environment. The magnetic field is uniform and parallel to the axis of threads. Configurations of identical and nonidentical threads are both explored. We apply the T-matrix theory of acoustic scattering to obtain the oscillatory frequency and the eigenfunctions of linear magnetosonic disturbances. We find that the existence of wave modes with a collective dynamics, i.e. those that produce significant perturbations in all threads, is only possible when the Doppler-shifted individual frequencies of threads are very similar. This can be only achieved for very particular values of the plasma physical conditions and flow velocities within threads.

We study the collective kinklike normal modes of a system of several cylindrical loops using the T-matrix theory. Loops that have similar kink frequencies oscillate collectively with a frequency which is slightly different from that of the individual kink mode. On the other hand, if the kink frequency of a loop is different from that of the others, it oscillates individually with its own frequency. Since the individual kink frequency depends on the loop density but not on its radius for typical 1 MK coronal loops, a coupling between kink oscillations of neighboring loops take place when they have similar densities. The relevance of these results in the interpretation of the oscillations studied by \citet{schrijver2000} and \citet{verwichte2004}, in which transverse collective loop oscillations seem to be detected, is discussed. In the first case, two loops oscillating in antiphase are observed; interpreting this motion as a collective kink mode suggests that their densities are roughly equal. In the second case, there are almost three groups of tubes that oscillate with similar periods and therefore their dynamics can be collective, which again seems to indicate that the loops of each group share a similar density. All the other loops seem to oscillate individually and their densities can be different from the rest.

First results from a high-resolution three-dimensional nonlinear numerical study of the kink oscillation are presented. We show in detail the development of a shear instability in an untwisted line-tied magnetic flux tube. The instability produces significant deformations of the tube boundary. An extended transition layer may naturally evolve as a result of the shear instability at a sharp transition between the flux tube and the external medium. We also discuss the possible effects of the instability on the process of resonant absorption when an inhomogeneous layer is included in the model. One of the implications of these results is that the azimuthal component of the magnetic field of a stable flux tube in the solar corona, needed to prevent the shear instability, is probably constrained to be in a very specific range.

Small-amplitude prominence oscillations are usually damped after a few periods. We study the attenuation of non-adiabatic magnetoacoustic waves in a slab prominence embedded in the coronal medium. We assume an equilibrium configuration with a transverse magnetic field to the slab axis and investigate wave damping by thermal conduction and radiative losses. The differential MHD equations that govern linear slow and fast modes are numerically solved to obtain the complex oscillatory frequency and the corresponding eigenfunctions. We find that coronal thermal conduction and radiative losses from the prominence plasma reveal as the most relevant damping mechanisms. Both mechanisms govern together the attenuation of hybrid modes, whereas prominence radiation is responsible for the damping of internal modes and coronal conduction essentially dominates the attenuation of external modes. In addition, the energy transfer between the prominence and the corona caused by thermal conduction has a noticeable effect on the wave stability, radiative losses from the prominence plasma being of paramount importance for the thermal stability of fast modes. We conclude that slow modes are efficiently damped, with damping times compatible with observations. On the contrary, fast modes are less attenuated by non-adiabatic effects and their damping times are several orders of magnitude larger than those observed. The presence of the corona causes a decrease of the damping times with respect to those of an isolated prominence slab, but its effect is still insufficient to obtain damping times of the order of the period in the case of fast modes.

High-resolution observations provide evidence about the existence of small-amplitude transverse oscillations in solar filament fine structures. These oscillations are believed to represent fast magnetohydrodynamic (MHD) waves and the disturbances are seen to be damped in short timescales of the order of 1 to 4 periods. In this Letter we propose that, due to the highly inhomogeneous nature of the filament plasma at the fine structure spatial scale, the phenomenon of resonant absorption is likely to operate in the temporal attenuation of fast MHD oscillations. By considering transverse inhomogeneity in a straight flux tube model we find that, for density inhomogeneities typical of filament threads, the decay times are of a few oscillatory periods only.

High-resolution observations show that oscillations and waves in prominence threads are common and that they are attenuated in a few periods. In addition, observers have also reported the presence of material flows in such prominence fine-structures. Here we investigate the time damping of non-leaky oscillations supported by a homogeneous cylindrical prominence thread embedded in an unbounded corona and with a steady mass flow. Thermal conduction and radiative losses are taken into account as damping mechanisms, and the effect of these non-ideal effects and the steady flow on the attenuation of oscillations is assessed. We solve the general dispersion relation for linear, non-adiabatic magnetoacoustic and thermal waves supported by the model, and find that slow and thermal modes are efficiently attenuated by non-adiabatic mechanisms. On the contrary, fast kink modes are much less affected and their damping times are much larger than those observed. The presence of flow has no effect on the damping of slow and thermal waves, whereas fast kink waves are more (less) attenuated when they propagate parallel (anti-parallel) to the flow direction. Although the presence of steady mass flows improves the efficiency of non-adiabatic mechanisms on the attenuation of transverse, kink oscillations for parallel propagation to the flow, its effect is still not enough to obtain damping times compatible with observations.

Recent observations with the Hinode Solar Optical Telescope display an active region prominence whose fine threads oscillate in the vertical direction as they move along a path parallel to the photosphere. A seismological analysis of this event is carried out by taking advantage of the small radius of these structures compared to the total length of magnetic field lines, i.e. by using the thin tube approximation. This analysis reveals that the oscillatory period is only slightly modified by the existence of the flow and that the difference between the period of a flowing thread and a static one is below the error bars of these observations. Moreover, although it is not possible to obtain values of the physical parameters, a lower bound for the Alfv\'en speed (ranging between 120 km s$^{-1}$ and 350 km s$^{-1}$) is obtained for each of the threads. Such Alfv\'en speeds agree with the intense magnetic fields and large densities usually found in active region prominences.

We study the excitation and damping of transverse oscillations in a multi-stranded model of a straight line-tied coronal loop. The transverse geometry of our equilibrium configuration is quite irregular and more realistic than the usual cylindrical loop model. By numerically solving the time-dependent ideal magnetohydrodynamic equations in two dimensions we show how the global motion of the whole bundle of strands, excited by an external disturbance, is converted into localized Alfv\'enic motions due to the process of resonant absorption. This process produces the attenuation of the transverse oscillations. At any location in the structure two dominant frequencies are found, the frequency of the global mode, or quasi-mode, and the local Alfv\'en frequency. We find that the mechanism of mode conversion, due to the coupling between fast and Alfv\'en waves, is not compromised by the complicated geometry of the model. We also show that it is possible to have energy conversion not only at the external edge of the composite loop but also inside the structure. The implications of these results and their relationship with the observations are discussed.

We study transverse fast magnetohydrodynamic waves in a system of two coronal loops modeled as smoothed, dense plasma cylinders in a uniform magnetic field. The collective oscillatory properties of the system due to the interaction between the individual loops are investigated from two points of view. Firstly, the frequency and spatial structure of the normal modes are studied. The system supports four trapped normal modes in which the loops move rigidly in the transverse direction. The direction of the motions is either parallel or perpendicular to the plane containing the axes of the loops. Two of these modes correspond to oscillations of the loops in phase, while in the other two they move in antiphase. Thus, these solutions are the generalization of the kink mode of a single cylinder to the double cylinder case. Secondly, we analyze the time-dependent problem of the excitation of the pair of tubes. We find that depending on the shape and location of the initial disturbance, different normal modes can be excited. The frequencies of normal modes are accurately recovered from the numerical simulations. In some cases, because of the simultaneous excitation of several eigenmodes, the system shows beating and the phase lag between the loops is $\pi/2$.

Oscillations of magnetic structures in the solar corona have often been interpreted in terms of magnetohydrodynamic waves. We study the adiabatic magnetoacoustic modes of a prominence plasma slab with a uniform longitudinal magnetic field, surrounded by a prominence-corona transition region (PCTR) and a coronal medium. Considering linear small-amplitude oscillations, the dispersion relation for the magnetoacoustic slow and fast modes is deduced assuming evanescent-like perturbations in the coronal medium. In the system without PCTR, a classification of the oscillatory modes according to the polarisation of their eigenfunctions is made in order to distinguish modes with fast-like or slow-like properties. Internal and external slow modes are governed by the prominence and coronal properties respectively, and fast modes are mostly dominated by prominence conditions for the observed wavelengths. In addition, the inclusion of an isothermal PCTR does not substantially influence the mode frequencies, but new solutions (PCTR slow modes) are present.

The main aim of this study is to point out the difficulties found when trying to assess the statistical significance of the North-South asymmetry (hereafter SSNSA) of the most usually considered time series of solar activity. First of all, we distinguish between solar activity time series composed by integer or non-integer and dimensionless data, or composed by non-integer and dimensional data. For each of these cases, we discuss the most suitable statistical tests which can be applied and highlight the difficulties to obtain valid information about the statistical significance of solar activity time series. Our results suggest that, apart from the need to apply the suitable statistical tests, other effects such as the data binning, the considered units and the need, in some tests, to consider groups of data, affect substantially the determination of the statistical significance of the asymmetry. Our main conclusion is that the assessment of the statistical significance of the N-S asymmetry of solar activity is a difficult matter and that an absolute answer cannot be given, since many different effects influence the results given by the statistical tests. In summary, the quantitative results about the statistical significance of the N-S asymmetry of solar activity provided by different authors, as well as the studies about its behaviour, must be considered with care because they depend from the chosen values of different parameters or from the considered units.

The theory of magnetohydrodynamic (MHD) waves in solar coronal slabs in a zero-$\beta$ configuration and for parallel propagation of waves does not allow the existence of surface waves. When oblique propagation of perturbations is considered both surface and body waves are able to propagate. When the perpendicular wave number is larger than a certain value, the body kink mode becomes a surface wave. In addition, a sausage surface mode is found below the internal cut-off frequency. When non-uniformity in the equilibrium is included, surface and body modes are damped due to resonant absorption. In this paper, first, a normal-mode analysis is performed and the period, the damping rate, and the spatial structure of eigenfunctions are obtained. Then, the time-dependent problem is solved, and the conditions under which one or the other type of mode is excited are investigated.

Observations of transversal coronal loop oscillations very often show the excitation and damping of oscillations in groups of coronal loops rather than in individual and isolated structures. We present results on the oscillatory properties (periods, damping rates, and spatial distribution of perturbations) for resonantly damped oscillations in a system of two inhomogeneous coronal slabs and compare them to the properties found in single slab loop models. A system of two identical coronal loops is modeled, in Cartesian geometry, as being composed by two density enhancements. The linear magnetohydrodynamic (MHD) wave equations for oblique propagation of waves are solved and the damping of the different solutions, due to the transversal inhomogeneity of the density profile, is computed. The physics of the obtained results is analyzed by an examination of the perturbed physical variables. We find that, due to the interaction between the loops, the normal modes of oscillation present in a single slab split into symmetric and antisymmetric oscillations when a system of two identical slabs is considered. The frequencies of these solutions may differ from the single slab results when the distance between the loops is of the order of a few slab widths. Oblique propagation of waves weakens this interaction, since solutions become more confined to the edges of the slabs. The damping is strong for surface-like oscillations, while sausage body-like solutions are unaffected. For some solutions, and small slab separations, the damping in a system of two loops differs substantially from the damping of a single loop.

Context. One of the typical features shown by observations of solar prominence oscillations is that they are damped in time and that the values of the damping times are usually between one and three times the corresponding oscillatory period. However, the mechanism responsible for the attenuation is still not well-known. Aims. Thermal conduction, optically thin or thick radiation and heating are taken into account in the energy equation, and their role on the attenuation of prominence oscillations is evaluated. Methods. The dispersion relation for linear non-adiabatic magnetoacoustic waves is derived considering an equilibrium made of a prominence plasma slab embedded in an unbounded corona. The magnetic field is orientated along the direction parallel to the slab axis and has the same strength in all regions. By solving the dispersion relation for a fixed wavenumber, a complex oscillatory frequency is obtained, and the period and the damping time are computed. Results. The effect of conduction and radiation losses is different for each magnetoacoustic mode and depends on the wavenumber. In the observed range of wavelengths the internal slow mode is attenuated by radiation from the prominence plasma, the fast mode by the combination of prominence radiation and coronal conduction and the external slow mode by coronal conduction. The consideration of the external corona is of paramount importance in the case of the fast and external slow modes, whereas it does not affect the internal slow modes at all. Conclusions. Non-adiabatic effects are efficient damping mechanisms for magnetoacoustic modes, and the values of the obtained damping times are compatible with those observed.

The influence of a toroidal magnetic field on the dynamics of Rossby waves in a thin layer of ideal conductive fluid on a rotating sphere is studied in the "shallow water" magnetohydrodynamic approximation for the first time. Dispersion relations for magnetic Rossby waves are derived analytically in Cartesian and spherical coordinates. It is shown that the magnetic field causes the splitting of low order (long wavelength) Rossby waves into two different modes, here denoted fast and slow {\em magnetic Rossby waves}. The high frequency mode (the fast magnetic Rossby mode) corresponds to an ordinary hydrodynamic Rossby wave slightly modified by the magnetic field, while the low frequency mode (the slow magnetic Rossby mode) has new and interesting properties since its frequency is significantly smaller than that of the same harmonics of pure Rossby and Alfv{\'e}n waves.

Observations reveal an increase with height of the line width of several coronal spectral lines probably caused by outwardly propagating Alfv{\'e}n waves. However, the spectral line width sometimes shows a sudden decrease at a height 0.1-0.2 R, where the ratio of sound to Alfven speeds may approach unity. Qualitative analysis shows that the resonant energy conversion from Alfven to acoustic waves near the region of the corona where the plasma $\beta$ approaches unity may explain the observed spectral line width reduction.