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

We present a model for the intensity of optically thin EUV emission for a plasma atmosphere. We apply our model to the solar corona as observed using the six optically thin EUV channels of the SDO/AIA instrument. The emissivity of the plasma is calculated from the density and temperature using CHIANTI tables and the intensity is then determined by integration along the line of sight. We consider several different profiles for the radial density and temperature profiles, each of which are constrained by the observational data alone with no further physical assumptions. We demonstrate the method first by applying it to a quiet region of the corona, and then use it as the background component of a model including coronal holes, allowing the plasma densities and temperatures inside and outside the hole to be estimated. We compare our results with differential emission measure (DEM) inversions. More accurate estimates for the coronal density and temperature profiles have the potential to help constrain plasma properties such as the magnetic field strength when used in combination with methods such as seismology.

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

Recent observations with the Atmospheric Imaging Assembly (AIA) instrument on the SDO spacecraft have revealed the existence of decayless coronal kink oscillations. These transverse oscillations are not connected to any external phenomena like flares or coronal mass ejections, and show significantly lower amplitudes than the externally excited decaying oscillations. Numerical studies have managed to reproduce such decayless oscillations in the form of footpoint driven standing waves in coronal loops, and to treat them as a possible mechanism for wave heating of the solar corona. Our aim is to investigate the correlation between the observed amplitudes of the oscillations and input the energy flux from different drivers. We perform 3D MHD simulations in single, straight, density-enhanced coronal flux tubes for different drivers, in the presence of gravity. Synthetic images at different spectral lines are constructed with the use of the FoMo code. The development of the Kelvin-Helmholtz instability leads to mixing of plasma between the flux tube and the hot corona. Once the KHI is fully developed, the amplitudes of the decayless oscillations show only a weak correlation with the driver strength. We find that low amplitude decayless kink oscillations may correspond to significant energy fluxes of the order of the radiative losses for the Quiet Sun. A clear correlation between the input energy flux and the observed amplitudes from our synthetic imaging data cannot be established. Stronger drivers lead to higher vales of the line width estimated energy fluxes. Finally, estimations of the energy fluxes by spectroscopic data are affected by the LOS angle, favoring combined analysis of imaging and spectroscopic data for single oscillating loops.

Recent developments in the observation and modelling of kink oscillations of coronal loops have led to heightened interest over the last few years. The modification of the Transverse Density Profile (TDP) of oscillating coronal loops by non-linear effects, in particular the Kelvin-Helmholtz Instability (KHI), is investigated. How this evolution may be detected is established, in particular, when the KHI vortices may not be observed directly. A model for the loop's TDP is used which includes a finite inhomogeneous layer and homogeneous core, with a linear transition between them. The evolution of the loop's transverse intensity profile from numerical simulations of kink oscillations is analysed. Bayesian inference and forward modelling techniques are applied to infer the evolution of the TDP from the intensity profiles, in a manner which may be applied to observations. The strongest observational evidence for the development of the KHI is found to be a widening of the loop's inhomogeneous layer, which may be inferred for sufficiently well resolved loops, i.e $>$ 15 data points across the loop. The main signatures when observing the core of the loop (for this specific loop model) during the oscillation are: a widening inhomogeneous layer, decreasing intensity, an unchanged radius, and visible fine transverse structuring when the resolution is sufficient. The appearance of these signatures are delayed for loops with wider inhomogeneous layers, and quicker for loops oscillating at higher amplitudes. These cases should also result in stronger observational signatures, with visible transverse structuring appearing for wide loops observed at SDO/AIA resolution.

Kink oscillations of a coronal loop are observed and studied in detail because they provide a unique probe into the structure of coronal loops through MHD seismology and a potential test of coronal heating through the phase-mixing of Alfv\'en waves. In particular, recent observations show that standing oscillations of loops often involve also higher harmonics, beside the fundamental mode. The damping of these kink oscillations is explained by mode coupling with Alfv\'en waves. We investigate the consequences for wave-based coronal heating of higher harmonics and what coronal heating observational signatures we may use to infer the presence of higher harmonic kink oscillations. We perform a set of non-ideal MHD simulations where the damping of the kink oscillation of a flux tube via mode coupling is modelled. Our MHD simulation parameters are based on the seismological inversion of an observation for which the first three harmonics are detected. We study the phase-mixing of Alfv\'en waves that leads to the deposition of heat in the system, and we apply the seismological inversion techniques to the MHD simulation output. We find that the heating due to phase-mixing of the Alfv\'en waves triggered by the damping of the kink oscillation is relatively small, however we can illustrate i) how the heating location drifts due to the subsequent damping of lower order harmonics. We also address the role of the higher order harmonics and the width of the boundary shell in the energy deposition. We conclude that the coronal heating due to phase-mixing seems not to provide enough energy to maintain the thermal structure of the solar corona even when multi-harmonics oscillations are included, and these oscillations play an inhibiting role in the development of smaller scale structures.

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

This short review paper gives an overview of recently observed transverse, propagating velocity perturbations in coronal loops. These ubiquitous perturbations are observed to undergo strong damping as they propagate. Using 3D numerical simulations of footpoint-driven transverse waves propagating in a coronal plasma with a cylindrical density structure, in combination with analytical modelling, it is demonstrated that the observed velocity perturbations can be understood in terms of coupling of different wave modes in the inhomogeneous boundaries of the loops. Mode coupling in the inhomogeneous boundary layers of the loops leads to the coupling of the transversal (kink) mode to the azimuthal (Alfven) mode, observed as the decay of the transverse kink oscillations. Both the numerical and analytical results show the spatial profile of the damped wave has a Gaussian shape to begin with, before switching to exponential decay at large heights. In addition, recent analysis of CoMP (Coronal Multi-channel Polarimeter) Doppler shift observations of large, off- limb, trans-equatorial loops shows that Fourier power at the apex appears to be higher in the high-frequency part of the spectrum than expected from theoretical models. This excess high-frequency FFT power could be tentative evidence for the onset of a cascade of the low-to-mid frequency waves into (Alfvenic) turbulence.

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

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 different 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.

We use images of high spatial and temporal resolution, obtained with the Rapid Oscillations in the Solar Atmosphere instrument at the Dunn Solar Telescope, to reveal how the generation of transverse waves in Type I spicules is a direct result of longitudinal oscillations occurring in the photosphere. Here we show how pressure oscillations, with periodicities in the range 130 - 440 s, manifest in small-scale photospheric magnetic bright points, and generate kink waves in the Sun's outer atmosphere with transverse velocities approaching the local sound speed. Through comparison of our observations with advanced two-dimensional magneto-hydrodynamic simulations, we provide evidence for how magnetoacoustic oscillations, generated at the solar surface, funnel upwards along Type I spicule structures, before undergoing longitudinal-to-transverse mode conversion into waves at twice the initial driving frequency. The resulting kink modes are visible in chromospheric plasma, with periodicities of 65 -220 s, and amplitudes often exceeding 400 km. A sausage mode oscillation also arises as a consequence of the photospheric driver, which is visible in both simulated and observational time series. We conclude that the mode conversion and period modification is a direct consequence of the 90 degree phase shift encompassing opposite sides of the photospheric driver. The chromospheric energy flux of these waves are estimated to be approximately 300,000 W/m^2, which indicates that they are sufficiently energetic to accelerate the solar wind and heat the localized corona to its multi-million degree temperatures.