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.

Context. Transverse oscillations of coronal structures are currently intensively studied to explore the associated magnetohydrodynamic wave physics and perform seismology of the local medium. Aims. We make a first attempt to measure the thermodynamic evolution of a sample of coronal loops that undergo decaying kink oscillations in response to an eruption in the corresponding active region. Methods. Using data from the six coronal wavelengths of SDO/AIA, we performed a differential emission measure (DEM) analysis of 15 coronal loops before, during, and after the eruption and oscillation. Results. We find that the emission measure, temperature, and width of the DEM distribution undergo significant variations on time scales relevant for the study of transverse oscillations. There are no clear collective trends of increases or decreases for the parameters we analysed. The strongest variations of the parameters occur during the initial perturbation of the loops, and the influence of background structures may also account for much of this variation. Conclusions. The DEM analysis of oscillating coronal loops in erupting active regions shows evidence of evolution on time scales important for the study of the oscillations. Further work is needed to separate the various observational and physical mechanisms that may be responsible for the variations in temperature, DEM distribution width, and total emission measure.

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.

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

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