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 report a quasi-periodic pulsation (QPP) event simultaneously detected from the spatial displacements of coronal loop at both EUV images and microwave emission during the preflare phase of a C1.1 flare on 2016 March 23. Using the motion magnification technique, a low-amplitude transverse oscillation with the growing period is discovered in a diffuse coronal loop in Atmospheric Imaging Assembly (AIA) image sequences at wavelength of 171 A, and the initial oscillation period is estimated to be ~397 s with a slow growth rate of 0.045. At the same time, a QPP with growing periods from roughly 300 s to nearly 500 s is discovered in the microwave flux in the same active region. Based on the imaging observations measured at EUV wavelengths by the AIA and at microwave 17 GHz by Nobeyama Radioheliograph, the diffuse coronal loop and the microwave radiation source are found to be connected through a hot loop seen in AIA images at wavelength of 94 A. The growing period of the QPP should be related to the modulation of LRC-circuit oscillating process in a current-carrying plasma loop. The existence of electric currents may imply the non-potentialities in the source region during the preflare phase.

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

At the Sun, the strongest magnetic fields are routinely detected at dark sunspots. The magnitude of the field is typically about 3000 G, with only a few exceptions that reported the magnetic field in excess of 5000 G. Given that the magnetic field decreases with height in the solar atmosphere, no coronal magnetic field above ~2000 G was ever reported. Here, we present imaging microwave observations of an anomalously strong magnetic field of about 4000 G at the base of the corona in solar active region NOAA 12673 on 06 September 2017. Combining the photospheric vector measurements of the magnetic field and the coronal probing, we created and validated a nonlinear force-free field coronal model, with which we quantify the record-breaking coronal magnetic field at various coronal heights.

We report hard X-ray and gamma-ray observations of the impulsive phase of the SOL2017-09-06T11:55 X9.3 solar flare. We focus on a high-energy part of the spectrum, >100 keV, and perform time resolved spectral analysis for a portion of the impulsive phase, recorded by the Konus-Wind experiment, that displayed prominent gamma-ray emission. Given a variety of possible emission components contributing to the gamma-ray emission, we employ a Bayesian inference to build the most probable fitting model. The analysis confidently revealed contributions from nuclear deexcitation lines, electron-positron annihilation line at 511 keV, and a neutron capture line at 2.223 MeV along with two components of the bremsstrahlung continuum. The revealed time evolution of the spectral components is particularly interesting. The low-energy bremsstrahlung continuum shows a soft-hard-soft pattern typical for impulsive flares, while the high-energy one shows a persistent hardening at the course of the flare. The neutron capture line emission shows an unusually short time delay relative to the nuclear deexcitation line component, which implies that the production of neutrons was significantly reduced soon after the event onset. This in turn may imply a prominent softening of the accelerated proton spectrum at the course of the flare, similar to the observed softening of the low-energy component of the accelerated electrons responsible for the low-energy bremsstrahlung continuum. We discuss possible physical scenarios, which might result in the obtained relationships between these gamma-ray components.

3D picture of the coronal magnetic field remains an outstanding problem in solar physics, particularly, in active regions. Nonlinear force-free field reconstructions that employ routinely available full-disk photospheric vector magnetograms represent state-of-the-art coronal magnetic field modeling. Such reconstructions, however, suffer from an inconsistency between a force-free coronal magnetic field and non-force-free photospheric boundary condition, from which the coronal reconstruction is performed. In this study we focus on integrating the additional chromospheric and / or coronal magnetic field data with the vector photospheric magnetograms with the goal of improving the reliability of the magnetic field reconstructions. We develop a corresponding modification of the available optimization codes described in Fleishman et al. (2017) and test their performance using a full-fledged MHD model obtained from the Bifrost code by performing a `voxel-by-voxel' comparison between the reconstructed and the model magnetic fields. We demonstrate that adding even an incomplete set of chromospheric magnetic field data can measurably improve the reconstruction of the coronal magnetic field, greatly improve reconstructions of the magnetic connectivity and of the coronal electric current.

Quantifying coronal magnetic field remains a central problem in solar physics. Nowadays the coronal magnetic field is often modelled using nonlinear force-free field (NLFFF) reconstructions, whose accuracy has not yet been comprehensively assessed. Here we perform a detailed casting of the NLFFF reconstruction tools, such as pi-disambiguation, photospheric field preprocessing, and volume reconstruction methods using a 3D snapshot of the publicly available full-fledged radiative MHD model. Specifically, from the MHD model we know the magnetic field vector in the entire 3D domain, which enables us to perform "voxel-by-voxel" comparison of the restored and the true magnetic field in the 3D model volume. Our tests show that the available pi-disambiguation methods often fail at the quiet sun areas dominated by small-scale magnetic elements, while they work well at the AR photosphere and (even better) chromosphere. The preprocessing of the photospheric magnetic field, although does produce a more force-free boundary condition, also results in some effective `elevation' of the magnetic field components. This `elevation' height is different for the longitudinal and transverse components, which results in a systematic error in absolute heights in the reconstructed magnetic data cube. The extrapolations performed starting from actual AR photospheric magnetogram are free from this systematic error, while have other metrics comparable with those for extrapolations from the preprocessed magnetograms. This finding favors the use of extrapolations from the original photospheric magnetogram without preprocessing. Our tests further suggest that extrapolations from a force-free chromospheric boundary produce measurably better results, than those from the photospheric boundary.

We present a seismological method to probe the solar atmosphere above sunspot umbrae with three-minute oscillations. Our technique allows us to estimate both the vertical distance between atmospheric layers and the wave-propagation speed without specifying any extra parameters, in particular, the phase speed of the wave or the emission formation heights. Our method uses the projected wave paths of slow MHD waves propagating through the atmospheric layers of different heights and guided by the magnetic field. The length of the projected wave path depends upon the distance between the layers and the inclination angle of the magnetic field with respect to the line of sight, allowing us to estimate the distance between the layers from measured projected wave paths and the local magnetic-field vector. In turn, the wave-propagation delay registered at different heights allows for the calculation of the phase speed. We estimated the vertical distance between the emission layers at the temperature minimum (1600 {\AA}) and transition region (304 {\AA}), as well as the average phase speed above the sunspot umbrae, for three active regions. We found that the distance between the 1600 {\AA} emission layer and the transition region above the sunspot umbrae lies in the range of 500 -- 800 km. The average phase speed between these layers was found to be about 30 km/s giving the sound speed of 6 km/s. The temperature between the layers has been roughly estimated as 3000 K and corresponds to the region of the temperature minimum. The results obtained are consistent with the semiempirical model of the sunspot umbrae atmosphere by Fontenla, 2009. The manuscript has been accepted for publication in Solar Physics journal.

Calculation of the potential field inside a three-dimensional box with the normal magnetic field component given on all boundaries is needed for estimation of important quantities related to the magnetic field such as free energy and relative helicity. In this work we present an analysis of three methods for calculating potential field inside a three-dimensional box. The accuracy and performance of the methods are tested on artificial models with a priori known solutions.

We study the writhe, twist and magnetic helicity of different magnetic flux ropes, based on models of the solar coronal magnetic field structure. These include an analytical force-free Titov--D\'emoulin equilibrium solution, non force-free magnetohydrodynamic simulations, and nonlinear force-free magnetic field models. The geometrical boundary of the magnetic flux rope is determined by the quasi-separatrix layer and the bottom surface, and the axis curve of the flux rope is determined by its overall orientation. The twist is computed by the Berger--Prior formula that is suitable for arbitrary geometry and both force-free and non-force-free models. The magnetic helicity is estimated by the twist multiplied by the square of the axial magnetic flux. We compare the obtained values with those derived by a finite volume helicity estimation method. We find that the magnetic helicity obtained with the twist method agrees with the helicity carried by the purely current-carrying part of the field within uncertainties for most test cases. It is also found that the current-carrying part of the model field is relatively significant at the very location of the magnetic flux rope. This qualitatively explains the agreement between the magnetic helicity computed by the twist method and the helicity contributed purely by the current-carrying magnetic field.

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

Magnetic helicity is a conserved quantity of ideal magneto-hydrodynamics characterized by an inverse turbulent cascade. Accordingly, it is often invoked as one of the basic physical quantities driving the generation and structuring of magnetic fields in a variety of astrophysical and laboratory plasmas. We provide here the first systematic comparison of six existing methods for the estimation of the helicity of magnetic fields known in a finite volume. All such methods are reviewed, benchmarked, and compared with each other, and specifically tested for accuracy and sensitivity to errors. To that purpose, we consider four groups of numerical tests, ranging from solutions of the three-dimensional, force-free equilibrium, to magneto-hydrodynamical numerical simulations. Almost all methods are found to produce the same value of magnetic helicity within few percent in all tests. In the more solar-relevant and realistic of the tests employed here, the simulation of an eruptive flux rope, the spread in the computed values obtained by all but one method is only 3%, indicating the reliability and mutual consistency of such methods in appropriate parameter ranges. However, methods show differences in the sensitivity to numerical resolution and to errors in the solenoidal property of the input fields. In addition to finite volume methods, we also briefly discuss a method that estimates helicity from the field lines' twist, and one that exploits the field's value at one boundary and a coronal minimal connectivity instead of a pre-defined three-dimensional magnetic-field solution.

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

In this study we combine the multiwavelength ultraviolet -- optical (Solar Dynamics Observatory, SDO) and radio (Nobeyama Radioheliograph, NoRH) observations to get further insight into space-frequency distribution of oscillations at different atmospheric levels of the Sun. We processed the observational data on NOAA 11711 active region and found oscillations propagating from the photospheric level through the transition region upward into the corona. The power maps of low-frequency (1--2 mHz) oscillations reproduce well the fan-like coronal structures visible in the Fe ix 171A line. High frequency oscillations (5--7 mHz) propagate along the vertical magnetic field lines and concentrate inside small-scale elements in the umbra and at the umbra-penumbra boundary. We investigated the dependence of the dominant oscillation frequency upon the distance from the sunspot barycentre to estimate inclination of magnetic tubes in higher levels of sunspots where it cannot be measured directly, and found that this angle is close to 40 degrees above the umbra boundaries in the transition region.

The Siberian Solar Radio Telescope (SSRT) is a solar-dedicated directly-imaging interferometer observing the Sun at 5.7 GHz. The SSRT operates in the two-dimensional mode since 1996. The imaging principle of the SSRT restricts its opportunities in observations of very bright flare sources, while it is possible to use `dirty' images in studies of low brightness features, which do not overlap with side lobes from bright sources. The interactive CLEAN technique routinely used for the SSRT data provides imaging of active regions but consumes much time and efforts and does not reveal low-brightness features below the CLEAN threshold. The newly developed technique combines the CLEAN routine with the directly imaging capability of the SSRT and provides clean images with an enhanced dynamic range automatically. These elaborations considerably extend the range of tasks, which can be solved with the SSRT. We show some examples of the present opportunities of the SSRT and compare its data with the images produced by the Nobeyama Radioheliograph at 17 GHz as well as observations in different spectral ranges.

We use daily full-disk vector magnetograms from Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) system to synthesize the first Carrington maps of the photospheric vector magnetic field. We describe these maps and make a comparison of observed radial field with the radial field estimate from LOS magnetograms. Further, we employ these maps to study the hemispheric pattern of current helicity density, Hc, during the rising phase of the solar cycle 24. Longitudinal average over the 23 consecutive solar rotations shows a clear signature of the hemispheric helicity rule, i.e. Hc is predominantly negative in the North and positive in South. Although our data include the early phase of cycle 24, there appears no evidence for a possible (systematic) reversal of the hemispheric helicity rule at the beginning of cycle as predicted by some dynamo models. Further, we compute the hemispheric pattern in active region latitudes (-30 deg \le \theta \le 30 deg) separately for weak (100< |B_r| <500 G)and strong (|B_r|>1000 G) radial magnetic fields. We find that while the current helicity of strong fields follows the well-known hemispheric rule (i.e., \theta . Hc < 0), H_c of weak fields exhibits an inverse hemispheric behavior (i.e., \theta . Hc > 0) albeit with large statistical scatter. We discuss two plausible scenarios to explain the opposite hemispheric trend of helicity in weak and strong field region.

The paper presents the possibility of fast and quality azimuth disambiguation of vector magnetogram data regardless of location on the solar disc. The new Super Fast and Quality (SFQ) code of disambiguation is tried out on well-known models of Metcalf et al. (2006), Leka et al. (2009) and artificial model of fixed configuration AR 10930 (Rudenko et al., 2010). We make comparison of Hinode SOT SP vector magnetograms of AR 10930 disambiguated with three codes: SFQ, NPFC (Georgoulis, 2005), and SME (Rudenko et al., 2010). We exemplify the SFQ disambiguation of SDO/HMI measurements of the full disc. The preliminary examination indicates that the SFQ algorithm provides better quality than NPFC and is comparable to SME. In contrast to other codes, SFQ supports relatively high quality of results regardless of the magnetogram proximity to the limb (when being very close to the limb, it is efficient unlike all other algorithms).

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

Many phenomena seen in solar atmosphere are connected with underphotospheric processes. MHD oscillation and other dynamical processes observed in quit sun objects are, for example, these phenomena. In our research we use Stereo Behind EUV(171\AA and 304\AA ) observations to study dynamical processes in quiet Sun regions. We found an eruption event in a coronal bright point on the border of a small coronal hole close to the center of the Solar disk. Several oscillating loops became visible inside of the hole after the eruption. We suppose that these loops oscillation were induced by the eruption process in coronal bright point. The nature of interaction transfer agent is unclear because any propagating disturbance was not detected. For data processing we used Pixelize wavelet filtration (PWF method) and Time-distance plots. We measured the supposed interaction transfer speed and found it about 2-3 km/s. We tried to find out the nature of the interaction transfer agent.

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

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