Search

Recent observational and theoretical studies indicate that the damping of solar coronal loop oscillations depends on the oscillation amplitude. We consider two mechanisms, linear resonant absorption and a nonlinear damping model. We confront theoretical predictions from these models with observed data in the plane of observables defined by the damping ratio and the oscillation amplitude. The structure of the Bayesian evidence in this plane displays a clear separation between the regions where each model is more plausible relative to the other. There is qualitative agreement between the regions of high marginal likelihood and Bayes factor for the nonlinear damping model and the arrangement of observed data. A quantitative application to 101 loop oscillation cases observed with SDO/AIA results in the marginal likelihood for the nonlinear model being larger in the majority of them. The cases with conclusive evidence for the nonlinear damping model outnumber considerably those in favor of linear resonant absorption., 15 pages, 4 figures, 1 table

The heating of the solar chromosphere and corona to the observed high temperatures, imply the presence of ongoing heating that balances the strong radiative and thermal conduction losses expected in the solar atmosphere. It has been theorized for decades that the required heating mechanisms of the chromospheric and coronal parts of the active regions, quiet-Sun, and coronal holes are associated with the solar magnetic fields. However, the exact physical process that transport and dissipate the magnetic energy which ultimately leads to the solar plasma heating are not yet fully understood. The current understanding of coronal heating relies on two main mechanism: reconnection and MHD waves that may have various degrees of importance in different coronal regions. In this review we focus on recent advances in our understanding of MHD wave heating mechanisms. First, we focus on giving an overview of observational results, where we show that different wave modes have been discovered in the corona in the last decade, many of which are associated with a significant energy flux, either generated in situ or pumped from the lower solar atmosphere. Afterwards, we summarise the recent findings of numerical modelling of waves, motivated by the observational results. Despite the advances, only 3D MHD models with Alfv\'en wave heating in an unstructured corona can explain the observed coronal temperatures compatible with the quiet Sun, while 3D MHD wave heating models including cross-field density structuring are not yet able to account for the heating of coronal loops in active regions to their observed temperature., Comment: accepted in Space Science Reviews

This paper investigates the changes in spatial properties when magnetohydrodynamic (MHD) waves undergo resonant damping in the Alfvén continuum. The analysis is carried out for a 1D cylindrical pressure-less plasma with a straight magnetic field. The effect of the damping on the spatial wave variables is determined by using complex frequencies that arise as a result of the resonant damping. Compression and vorticity are used to characterise the spatial evolution of the MHD wave. The most striking result is the huge spatial variation in the vorticity component parallel to the magnetic field. Parallel vorticity vanishes in the uniform part of the equilibrium. However, when the MHD wave moves into the non-uniform part, parallel vorticity explodes to values that are orders of magnitude higher than those attained by the transverse components in planes normal to the straight magnetic field. In the non-uniform part of the equilibrium plasma, the MHD wave is controlled by parallel vorticity and resembles an Alfvén wave, with the unfamiliar property that it has pressure variations even in the linear regime.

We use Coronal Multi-channel Polarimeter (CoMP) observations of propagating waves in the solar corona and Bayesian analysis to assess the evidence of models with resonant damping and foot-point wave power asymmetries. Two nested models are considered. The reduced model considers resonant damping as the sole cause of the measured discrepancy between outward and inward wave power. The larger model contemplates an extra source of asymmetry with origin at the foot-points. We first compute probability distributions of parameters conditional on the models and the observed data. The obtained constraints are then used to calculate the evidence for each model in view of data. We find that we need to consider the larger model to explain CoMP data and to accurately infer the damping ratio, hence, to better assess the possible contribution of the waves to coronal heating., Accepted in A&A

The Kelvin-Helmholtz instability has been proposed as a mechanism to extract energy from magnetohydrodynamic (MHD) kink waves in flux tubes, and to drive dissipation of this wave energy through turbulence. It is therefore a potentially important process in heating the solar corona. However, it is unclear how the instability is influenced by the oscillatory shear flow associated with an MHD wave. We investigate the linear stability of a discontinuous oscillatory shear flow in the presence of a horizontal magnetic field within a Cartesian framework that captures the essential features of MHD oscillations in flux tubes. We derive a Mathieu equation for the Lagrangian displacement of the interface and analyse its properties, identifying two different instabilities: a Kelvin-Helmholtz instability and a parametric instability involving resonance between the oscillatory shear flow and two surface Alfv\'{e}n waves. The latter occurs when the system is Kelvin-Helmholtz stable, thus favouring modes that vary along the flux tube, and as a consequence provides an important and additional mechanism to extract energy. When applied to flows with the characteristic properties of kink waves in the solar corona, both instabilities can grow, with the parametric instability capable of generating smaller scale disturbances along the magnetic field than possible via the Kelvin-Helmholtz instability. The characteristic time-scale for these instabilities is $\sim 100$ s, for wavelengths of $200$ km. The parametric instability is more likely to occur for smaller density contrasts and larger velocity shears, making its development more likely on coronal loops than on prominence threads., Comment: 12 pages, 4 figures, accepted for publication MNRAS

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., Accepted for publication in the Astronomy and Astrophysics

This paper investigates the mixed properties of MHD waves in a non-uniform plasma. It starts with a short revision of MHD waves in a uniform plasma of infinite extent. In that case the MHD waves do not have mixed properties. They can be separated in Alfv\'{e}n waves and magneto-sonic waves. The Alfv\'{e}n waves propagate parallel vorticity and are incompressible. In addition they have no parallel displacement component. The magneto-sonic waves are compressible and in general do have a parallel component of displacement but do not propagate parallel vorticity. This clear separation has been the reason why there has been a strong inclination in the literature to use this classification in the study of MHD waves in non-uniform plasmas. The main part of this paper is concerned with MHD waves in a non-uniform plasma. It is shown that the MHD waves in that situation in general propagate both vorticity and compression and hence have mixed properties. Finally, the close connection between resonant absorption and MHD waves with mixed properties is discussed.

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.

In contrast to the situation in a laboratory, the study of the solar atmosphere has to be pursued without direct access to the physical conditions of interest. Information is therefore incomplete and uncertain and inference methods need to be employed to diagnose the physical conditions and processes. One of such methods, solar atmospheric seismology, makes use of observed and theoretically predicted properties of waves to infer plasma and magnetic field properties. A recent development in solar atmospheric seismology consists in the use of inversion and model comparison methods based on Bayesian analysis. In this paper, the philosophy and methodology of Bayesian analysis are first explained. Then, we provide an account of what has been achieved so far from the application of these techniques to solar atmospheric seismology and a prospect of possible future extensions., 19 pages, accepted in Advances in Space Research

Recent observations show cool, oscillating prominence threads fading when observed in cool spectral lines and appearing in warm spectral lines. A proposed mechanism to explain this evolution is that the threads were heated by turbulence driven by the Kelvin--Helmholtz instability that developed as a result of wave-driven shear flows on the surface of the thread. As the Kelvin--Helmholtz instability is an instability that works to mix the fluids, in the solar corona it can be expected to work by mixing the cool prominence material with that of the hot corona to form a warm boundary layer. In this paper we develop a simple phenomenological model of nonlinear Kelvin--Helmholtz mixing, using it to determine the characteristic density and temperature of the mixing layer, which for the case under study with constant pressure across the two fluids are $\rho_{\rm mixed}=\sqrt{\rho_1\rho_2}$ and $T_{\rm mixed}=\sqrt{T_1T_2}$. One result from the model is that it provides an accurate, as determined by comparison with simulation results, determination of the kinetic energy in the mean velocity field. A consequence of this is that the magnitude of turbulence, and with it the energy that can be dissipated on fast time-scales, as driven by this instability can be determined. For the prominence-corona system, the mean temperature rise possible from turbulent heating is estimated to be less than 1\% of the characteristic temperature (which is found to be $10^5$\,K). These results highlight that mixing, and not heating, are likely to be the cause of the observed transition between cool to warm material in Okamoto et. al (2015). One consequence of this result is that the mixing creates a region with higher radiative loss rates on average than either of the original fluids, meaning that this instability could contribute a net loss of thermal energy from the corona, i.e. coronal cooling., Comment: 17 pages, 9 figures, accepted for publication in ApJ

A classic application of coronal seismology uses transverse oscillations of waveguides to obtain estimates of the magnetic field strength. The procedure requires information on the density of the structures. Often, it ignores the damping of the oscillations. We computed marginal posteriors for parameters such as the waveguide density; the density contrast; the transverse inhomogeneity length-scale; and the magnetic field strength, under the assumption that the oscillations can be modelled as standing magnetohydrodynamic (MHD) kink modes damped by resonant absorption. Our results show that the magnetic field strength can be properly inferred, even if the densities inside and outside the structure are largely unknown. Incorporating observational estimates of plasma density further constrains the obtained posteriors. The amount of information one is willing to include (a priori) for the density and the density contrast influences their corresponding posteriors, but very little the inferred magnetic field strength. The decision to include or leave out the information on the damping and the damping time-scales have a minimal impact on the obtained magnetic field strength. In contrast to the classic method which provides with numerical estimates with error bars or possible ranges of variation for the magnetic field strength, Bayesian methods offer the full distribution of plausibility over the considered range of possible values. The methods are applied to available datasets of observed transverse loop oscillations, can be extended to prominence fine structures or chromospheric spicules and implemented to propagating waves in addition to standing oscillations., Comment: 16 pages, A&A accepted

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., 11 pages, 2 figures, ApJL accepted

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.

Solar tornadoes are dark vertical filamentary structures observed in the extreme ultraviolet associated with prominence legs and filament barbs. Their true nature and relationship to prominences requires understanding their magnetic structure and dynamic properties. Recently, a controversy has arisen: is the magnetic field organized forming vertical, helical structures or is it dominantly horizontal? And concerning their dynamics, are tornadoes really rotating or is it just a visual illusion? Here, we analyze four consecutive spectropolarimetric scans of a prominence hosting tornadoes on its legs which help us shed some light on their magnetic and dynamical properties. We show that the magnetic field is very smooth in all the prominence, probably an intrinsic property of the coronal field. The prominence legs have vertical helical fields that show slow temporal variation probably related to the motion of the fibrils. Concerning the dynamics, we argue that 1) if rotation exists, it is intermittent, lasting no more than one hour, and 2) the observed velocity pattern is also consistent with an oscillatory velocity pattern (waves)., Comment: accepted for publication in ApJ

Context. The geometry and physical conditions in solar coronal loops are complicated and still not understood well. Recent high-resolution observations obtained with TRACE indicate the existence of sub-resolution transverse structuring not accessible to direct observation. This ingredient has not yet been taken into account in the previous theoretical models used for the study of transversal coronal loop oscillations and of their damping due to resonant conversion of energy. Aims. This study aims to assess the effect of the possibly unresolved internal structure of a coronal loop on the properties of its transverse oscillations and on the efficiency of resonant absorption as a damping mechanism of these oscillations. Methods. The equilibrium configuration of a single coronal loop with internal density structuring is modelled by considering the loop as composed of two very close, parallel, identical coronal slabs in Cartesian geometry. The period of the oscillation and the damping time are computed for the resonantly damped fundamental kink mode. These quantities are then compared to those obtained for two models for a single equivalent slab without internal density structuring. Results. We find that the period and the damping time of a coronal loop with internal density structuring change by less than 15%, when compared to the same oscillatory properties of a single coronal loop with either the same density contrast or a single coronal loop with the same total mass. Conclusions. Therefore the internal density structuring of a coronal loop does not affect its oscillatory properties very much. However, the sub-resolution structuring of a coronal loop with different densities in its components or with different widths could vary these results.

The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement of the spatial distribution of density, however, is difficult and indirect seismology inversion methods have been suggested as an alternative. We applied Bayesian inference, model comparison, and model-averaging techniques to the inference of the cross-field density structuring in solar magnetic waveguides using information on periods and damping times for resonantly damped magnetohydrodynamic (MHD) transverse kink oscillations. Three commonly employed alternative profiles were used to model the variation of the mass density across the waveguide boundary. Parameter inference enabled us to obtain information on physical quantities such as the Alfv\'en travel time, the density contrast, and the transverse inhomogeneity length scale. The inference results from alternative density models were compared and their differences quantified. Then, the relative plausibility of the considered models was assessed by performing model comparison. Our results indicate that the evidence in favor of any of the three models is minimal, unless the oscillations are strongly damped. In such a circumstance, the application of model-averaging techniques enables the computation of an evidence-weighted inference that takes into account the plausibility of each model in the calculation of a combined inversion for the unknown physical parameters., Comment: 12 pages, 9 figures, 3 tables

Aims. We combine the magnetohydrodynamic (MHD) theory of resonantly damped quasi-mode kink oscillations with observational estimates of the period and damping of transverse coronal loop oscillations to extract information on physical parameters in oscillating loops. Methods. A numerical study of the quasi-mode period and damping, in one-dimensional fully non-uniform flux tubes, is used to obtain equilibrium models that reproduce the observed periods and damping rates. This scheme is applied to 11 loop oscillation events. Results. When only the damping rate is used. the valid equilibrium models form a one-dimensional solution curve in the two-dimensional parameter space (density contrast, transverse inhomogeneity length-scale). Lower limits to the transverse inhomogeneity are obtained in the limit of high contrast loops. When both the period and the damping rate are used, the equilibrium Alfven speed (or Alfven travel time) comes into play. The valid equilibrium models then form a one-dimensional solution curve in the three-dimensional parameter space (density contrast, transverse inhomogeneity length-scale, Alfven speed or Alfven travel time). The projection of these solutions onto the Alfven speed axis is found to be constrained to a rather limited interval. Upper limits to the internal Alfven speed are derived for 9 of the 1 1 analysed events.

In the last decade we have been overwhelmed by an avalanche of discoveries of magnetohydrodynamic (MHD) waves by the Solar and Heliospheric Observatory and Transition Region and Coronal Explorer observatories. Both standing and propagating versions of fast magnetoacoustic and slow magnetoacoustic MHD waves have been detected. Information on the damping times and damping distances of these waves is less detailed and less accurate than that on periods and amplitudes. Nevertheless, observations show the damping times and damping lengths are often short. Also, different types of MHD waves in different types of magnetic structures likely require different damping mechanisms. The phenomenon of fast damping is well documented for the standing fast magnetosonic kink waves in coronal loops. This paper concentrates on standing fast magnetosonic waves. It reports on results on periods and damping times due to resonant absorption in one-dimensional and two-dimensional models of coronal loops. Special attention is given to multiple modes.

We will discuss the observed, heavily damped transversal oscillations of coronal loops. These oscillations are often modeled as transversal kink oscillations in a cylinder. Several features are added to the classical cylindrical model. In our models we include loop curvature, longitudinal density stratification, and highly inhomogeneous radial density profiles. In this paper, we will first give an overview of recently obtained results, both analytically and numerically. After that, we shed a light on the computational aspects of the modeling process. In particular, we will focus on the parallellization of the numerical codes.

Resonantly damped fast kink quasi-modes are computed in fully resistive magnetohydrodynamics (MHD) for two-dimensional equilibrium models. The equilibrium model is a straight cylindrically symmetric flux tube with a plasma density that is non-uniform both across and along the loop. The non-uniform layer across the loop is not restricted to be thin, but its thickness can reach values up to the loop diameter. Our results indicate that the period and damping of coronal loop oscillations mainly depend on the density contrast and the inhomogeneity length-scale and are independent of the details of longitudinal stratification, depending on the weighted mean density, weighted with the wave energy. For fully non-uniform loops, quasi-modes can interact with resistive Alfven eigenmodes leading to avoided crossings and gaps in the complex frequency plane. The present study extends previous studies on coronal loop oscillations in one-dimensional equilibrium models with thick boundary layers and in equilibria with longitudinally stratified loops under the thin boundary approximation, and allow for a better comparison between observations and theory raising the prospect of coronal seismology using the time damping of coronal loop oscillations.

The recent detection of multiple harmonic standing transverse oscillations in coronal loops by Verwichte et al. is of special importance, as it allows one to obtain information on the longitudinal density variation in loops. Verwichte et al. detected the simultaneous presence of both the fundamental and the first-overtone mode in two coronal loops. Here we point out that the ratio of the period of the fundamental mode to the period of the overtone mode differs from 2 in loops with longitudinal density stratification. Conversely, the difference between this ratio and 2 can be used as a seismological tool to obtain information about the density scale height in loops.

The observed coronal loop oscillations and their damping are often theoretically described by the use of a very simple coronal loop model, viz. a straight, longitudinally invariant, axi-symmetric, and pressureless flux tube with a different density inside and outside of the loop. In this paper we generalize the model by including longitudinal density stratification and we examine how the longitudinal density stratification alters the linear eigenmodes of the system, their oscillation frequencies, and the damping rates by resonant absorption.

We present a method to assess the plausibility of alternative mechanisms to explain the damping of magnetohydrodynamic (MHD) transverse waves in solar coronal loops. The considered mechanisms are resonant absorption of kink waves in the Alfv\'en continuum, phase-mixing of Alfv\'en waves, and wave leakage. Our methods make use of Bayesian inference and model comparison techniques. We first infer the values for the physical parameters that control the wave damping, under the assumption of a particular mechanism, for typically observed damping time-scales. Then, the computation of marginal likelihoods and Bayes factors enable us to quantify the relative plausibility between the alternative mechanisms. We find that, in general, the evidence is not large enough to support a single particular damping mechanism as the most plausible one. Resonant absorption and wave leakage offer the most probable explanations in strong damping regimes, while phase mixing is the best candidate for weak/moderate damping. When applied to a selection of 89 observed transverse loop oscillations, with their corresponding measurements of damping times scales and taking into account data uncertainties, we find that only in a few cases positive evidence for a given damping mechanism is available.

Time and spatial damping of transverse magnetohydrodynamic (MHD) kink oscillations is a source of information on the cross-field variation of the plasma density in coronal waveguides. We show that a probabilistic approach to the problem of determining the density structuring from the observed damping of transverse oscillations enables us to obtain information on the two parameters that characterise the cross-field density profile. The inference is performed by computing the marginal posterior distributions for density contrast and transverse inhomo- geneity length-scale using Bayesian analysis and damping ratios for transverse oscillations under the assumption that damping is produced by resonant absorption. The obtained distributions show that, for damping times of a few oscillatory periods, low density contrasts and short inho- mogeneity length scales are more plausible in explaining observations. This means that valuable information on the cross-field density profile can be obtained even if the inversion problem, with two unknowns and one observable, is a mathematically ill-posed problem., 5 pages, 3 figures, accepted

Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. Alfven waves and magneto-sonic waves are particular classes of MHD waves. These wave modes are clearly different and have pure properties in uniform plasmas of infinite extent only. Due to plasma non-uniformity MHD waves have mixed properties and cannot be classified as pure Alfven or magneto-sonic waves. However, vorticity is a quantity unequivocally related to Alfven waves as compression is for magneto-sonic waves. Here, we investigate MHD waves superimposed on a one-dimensional non-uniform straight cylinder with constant magnetic field. For a piecewise constant density profile we find that the fundamental radial modes of the non-axisymmetric waves have the same properties as surface Alfven waves at a true discontinuity in density. Contrary to the classic Alfveen waves in a uniform plasma of infinite extent, vorticity is zero everywhere except at the cylinder boundary. If the discontinuity in density is replaced with a continuous variation of density, vorticity is spread out over the whole interval with non-uniform density. The fundamental radial modes of the non-axisymmetric waves do not need compression to exist unlike the radial overtones. In thin magnetic cylinders the fundamental radial modes of the non-axisymmetric waves with phase velocities between the internal and the external Alfven velocities can be considered as surface Alfven waves. On the contrary, the radial overtones can be related to fast-like magneto-sonic modes., Accepted for publication in ApJ

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.

Magnetohydrodynamic (MHD) wave activity is ubiquitous in the solar atmosphere. MHD seismology aims to determine difficult to measure physical parameters in solar atmospheric magnetic and plasma structures by a combination of observed and theoretical properties of MHD waves and oscillations. This technique, similar to seismology or helio-seismology, demands the solution of two problems. The direct problem involves the computation of wave properties of given theoretical models. The inverse problem implies the calculation of unknown physical parameters, by means of a comparison of observed and theoretical wave properties. Solar atmospheric seismology has been successfully applied to different structures such as coronal loops, prominence plasmas, spicules, or jets. However, it is still in its infancy. Far more is there to come. We present an overview of recent results, with particular emphasis in the inversion procedure.

There are ubiquitous flows observed in the solar atmosphere of sub-Alfv\'{e}nic speeds, however after flaring and coronal mass ejection events flows can become Alfv\'{e}nic. In this Letter, we derive an expression for the standing kink mode frequency due to siphon flow in coronal loops, valid for both low and high speed regimes. It is found that siphon flow introduces a linear spatially dependent phase shift along coronal loops and asymmetric eigenfunctions. We demonstrate how this theory can be used to determine the kink and flow speed of oscillating coronal loops with reference to an observational case study. It is shown that the presence of siphon flow can cause the underestimation of magnetic field strength in coronal loops using the traditional seismological methods.

We analyse the oscillatory properties of resonantly damped transverse kink oscillations in two-dimensional prominence threads. The fine structures are modelled as cylindrically symmetric magnetic flux tubes with a dense central part with prominence plasma properties and an evacuated part, both surrounded by coronal plasma. The equilibrium density is allowed to vary non-uniformly in both the transverse and the longitudinal directions.We examine the influence of longitudinal density structuring on periods, damping times, and damping rates for transverse kink modes computed by numerically solving the linear resistive magnetohydrodynamic (MHD) equations. The relevant parameters are the length of the thread and the density in the evacuated part of the tube, two quantities that are difficult to directly estimate from observations. We find that both of them strongly influence the oscillatory periods and damping times, and to a lesser extent the damping ratios. The analysis of the spatial distribution of perturbations and of the energy flux into the resonances allows us to explain the obtained damping times. Implications for prominence seismology, the physics of resonantly damped kink modes in two-dimensional magnetic flux tubes, and the heating of prominence plasmas are discussed., 12 pages, 9 figures, A&A accepted

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., Comment: Accepted for publication in ApJ

Small amplitude oscillations are a commonly observed feature in prominences/filaments. These oscillations appear to be of local nature, are associated to the fine structure of prominence plasmas, and simultaneous flows and counterflows are also present. The existing observational evidence reveals that small amplitude oscillations, after excited, are damped in short spatial and temporal scales by some as yet not well determined physical mechanism(s). Commonly, these oscillations have been interpreted in terms of linear magnetohydrodynamic (MHD) waves, and this paper reviews the theoretical damping mechanisms that have been recently put forward in order to explain the observed attenuation scales. These mechanisms include thermal effects, through non-adiabatic processes, mass flows, resonant damping in non-uniform media, and partial ionization effects. The relevance of each mechanism is assessed by comparing the spatial and time scales produced by each of them with those obtained from observations. Also, the application of the latest theoretical results to perform prominence seismology is discussed, aiming to determine physical parameters in prominence plasmas that are difficult to measure by direct means., Comment: 36 pages, 16 figures, Space Science Reviews (accepted)

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., Comment: 27 pages, 11 figures

There is not as yet full agreement on the mechanism that causes the rapid damping of the oscillations observed by TRACE in coronal loops. It has been suggested that the variation of the observed values of the damping time as function of the corresponding observed values of the period contains information on the possible damping mechanism. The aim of this Letter is to show that, for resonant absorption, this is definitely not the case unless detailed a priori information on the individual loops is available.

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.

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.

We present an analytic approximate seismic inversion scheme for damped transverse coronal loop oscillations based on the thin tube and thin boundary approximation for computing the period and the damping time. Asymptotic expressions for the period and damping rate are used to illustrate the process of seismological inversion in a simple and easy to follow manner. The inversion procedure is formulated in terms of two simple functions, which are given by simple closed expressions. The analytic seismic inversion shows that an infinite amount of 1-dimensional equilibrium models can reproduce the observed periods and damping times. It predicts a specific range of allowable values for the Alfven travel time and lower bounds for the density contrast and the inhomogeneity length scale. When the results of the present analytic seismic inversion are compared with those of a previous numerical inversion, excellent agreement is found up to the point that the analytic seismic inversion emerges as a tool for validating results of numerical inversions. Actually it helped us to identify and correct inaccuracies in a previous numerical investigation., Comment: 7 pages, 1 figure, A&A, accepted

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., Comment: 19 pages, 9 figures, accepted for publication in Solar Physics

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., Comment: 25 pages, 9 figures

Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? The finding and explanation of dislocations may shed light on the nature and characteristics of the propagating waves, their interaction in the corona and in general on the plasma dynamics. We positively identify dislocations in coronal waves observed by the Coronal Multi-channel Polarimeter (CoMP) as singularities in the Doppler shifts of emission coronal lines. We study the possible singularities that can be expected in coronal waves and try to reproduce the observed dislocations in terms of localization and frequency of appearance. The observed dislocations can only be explained by the interference of a kink and a sausage wave modes propagating with different frequencies along the coronal magnetic field. In the plane transverse to the propagation, the cross-section of the oscillating plasma must be smaller than the spatial resolution, and the two waves result in net longitudinal and transverse velocity components that are mixed through projection onto the line of sight. Alfv\'en waves can be responsible of the kink mode, but a magnetoacoustic sausage mode is necessary in all cases. Higher (flute) modes are excluded. The kink mode has a pressure amplitude that is smaller than the pressure amplitude of the sausage mode, though its observed velocity is larger. This concentrates dislocations on the top of the loop. To explain dislocations, any model of coronal waves must include the simultaneous propagation and interference of kink and sausage wave modes of comparable but different frequencies, with a sausage wave amplitude much smaller than the kink one., Comment: 11 pages. 5 figures. Accepted for publication in A&A

Context. Quiescent solar prominence fine structures are typically modelled as density enhancements, called threads, which occupy a fraction of a longer magnetic flux tube. This is justified from the spatial distribution of the imaged plasma emission or absorption of prominences at small spatial scales. The profile of the mass density along the magnetic field is unknown, however, and several arbitrary alternatives are employed in prominence wave studies. The identification and measurement of period ratios from multiple harmonics in standing transverse thread oscillations offer a remote diagnostics method to probe the density variation of these structures. Aims. We present a comparison of theoretical models for the field-aligned density along prominence fine structures. They aim to imitate density distributions in which the plasma is more or less concentrated around the centre of the magnetic flux tube. We consider Lorentzian, Gaussian, and parabolic profiles. We compare theoretical predictions based on these profiles for the period ratio between the fundamental transverse kink mode and the first overtone to obtain estimates for the density ratios between the central part of the tube and its foot-points and to assess which one would better explain observed period ratio data. Methods. Bayesian parameter inference and model comparison techniques were developed and applied. To infer the parameters, we computed the posterior distribution for the density gradient parameter that depends on the observable period ratio. The model comparison involved computing the marginal likelihood as a function of the period ratio to obtain the plausibility of each density model as a function of the observable. We also computed the Bayes factors to quantify the relative evidence for each model, given a period ratio observation. Results. A Lorentzian density profile, with plasma density concentrated around the centre of the tube, seems to offer the most plausible inversion result. A Gaussian profile would require unrealistically high values of the density gradient parameter, and a parabolic density distribution does not enable us to obtain well-constrained posterior probability distributions. However, our model comparison results indicate that the evidence points to the Gaussian and parabolic profiles for period ratios in between 2 and 3, while the Lorentzian profile is preferred for higher period ratio values. The method we present can be used to obtain information on the plasma structure along threads, provided period ratio measurements become widely available.

Magnetic waves are a relevant component in the dynamics of the solar atmosphere. Their significance has increased because of their potential as a remote diagnostic tool and their presumed contribution to plasma heating processes. We discuss our current understanding on coronal heating by magnetic waves, based on recent observational evidence and theoretical advances. The discussion starts with a selection of observational discoveries that have brought magnetic waves to the forefront of the coronal heating discussion. Then, our theoretical understanding on the nature and properties of the observed waves and the physical processes that have been proposed to explain observations are described. Particular attention is given to the sequence of processes that link observed wave characteristics with concealed energy transport, dissipation, and heat conversion. We conclude with a commentary on how the combination of theory and observations should help us understanding and quantifying magnetic wave heating of the solar atmosphere., Comment: Accepted for publication in Royal Society of London Philosophical Transactions Series A

Resonantly damped fast kink quasi-modes are computed in fully resistive magnetohydrodynamics for fully non-uniform two-dimensional equilibrium models. The equilibrium model is a straight cylindrically symmetric flux tube with a plasma density that is non-uniform both across and along the loop. The variation of density across the loop can cover the whole loop. Our results indicate that the period and damping of coronal loop oscillations mainly depend on the density contrast and the inhomogeneity length-scale and are independent of the details of longitudinal stratification. This study extends previous studies on coronal loop oscillations, and allows for a better comparison between observations and theory.

The magnetohydrodynamic modes of oscillation of a force-free, sheared coronal arcade are investigated. The approach is to consider the general equations for this type of equilibrium and then to introduce small perturbations. The linearised momentum equation then results in two coupled partial differential equations for the velocity components in the directions normal and perpendicular to the equilibrium magnetic field. A numerical code is used to solve the resulting system of coupled partial differential equations and its accuracy has been checked by numerically solving a case for which analytical or simple numerical solutions exist.

The analysis of quickly damped transverse oscillations of solar coronal loops using magneto-hydrodynamic seismology allow us to infer physical parameters that are difficult to measure otherwise. Under the assumption that such damped oscillations are due to the resonant conversion of global modes into Alfven oscillations of the tube surface, we carry out a global seismological analysis of a large set of coronal loops. A Bayesian hierarchical method is used to obtain distributions for coronal loop physical parameters by means of a global analysis of a large number of observations. The resulting distributions summarise global information and constitute data-favoured information that can be used for the inversion of individual events. The results strongly suggest that internal Alfven travel times along the loop are larger than 100 s and smaller than 540 s with 95% probability. Likewise, the density contrast between the loop interior and the surrounding is larger than 2.3 and below 6.9 with 95% probability., Comment: 9 pages, 6 figures, accepted for publication in A&A

We present a Bayesian seismology inversion technique for propagating magnetohydrodynamic (MHD) transverse waves observed in coronal waveguides. The technique uses theoretical predictions for the spatial damping of propagating kink waves in transversely inhomogeneous coronal waveguides. It combines wave amplitude damping length scales along the waveguide with theoretical results for resonantly damped propagating kink waves to infer the plasma density variation across the oscillating structures. Provided the spatial dependence of the velocity amplitude along the propagation direction is measured and the existence of two di erent damping regimes is identified, the technique would enable us to fully constrain the transverse density structuring, providing estimates for the density contrast and its transverse inhomogeneity length scale.

The detection of multiple mode harmonic kink oscillations in coronal loops enables to obtain information on coronal density stratification and magnetic field expansion using seismology inversion techniques. The inference is based on the measurement of the period ratio between the fundamental mode and the first overtone and theoretical results for the period ratio under the hypotheses of coronal density stratification and magnetic field expansion of the wave guide. We present a Bayesian analysis of multiple mode harmonic oscillations for the inversion of the density scale height and magnetic flux tube expansion, under each of the hypotheses. Then, the two models are compared using a Bayesian model comparison scheme to assess how plausible each one is, given our current state of knowledge., 5 pages