Source author record

V. M. Nakariakov

V. M. Nakariakov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

High Frequency Waves in Chromospheric Spicules

Using high cadence observations from the Hydrogen-alpha Rapid Dynamics camera imaging system on the Dunn Solar Telescope, we present an investigation of the statistical properties of transverse oscillations in spicules captured above the solar limb. At five equally separated atmospheric heights, spanning approximately 4900-7500 km, we have detected a total of 15 959 individual wave events, with a mean displacement amplitude of 151 +/- 124 km, a mean period of 54 +/- 45 s, and a mean projected velocity amplitude of 21 +/- 13 km s^-1. We find that both the displacement and velocity amplitudes increase with height above the solar limb, ranging from 132 +/- 111 km and 17.7 +/- 10.6 km s^-1 at 4900 km, and 168 +/- 125 km and 26.3 +/- 14.1 km s^-1 at 7500 km, respectively. Following the examination of neighboring oscillations in time and space, we find 45% of the waves to be upwardly propagating, 49% to be downwardly propagating, and 6% to be standing, with mean absolute phase velocities for the propagating waves on the order of 75-150 km s^-1. While the energy flux of the waves propagating downwards does not appear to depend on height, we find the energy flux of the upwardly propagating waves decreases with atmospheric height at a rate of -13 200 +/- 6500 W m^-2 /Mm. As a result, this decrease in energy flux as the waves propagate upwards may provide significant thermal input into the local plasma.

preprint2021arXiv

The dynamics and observability of circularly polarized kink waves

Context. Kink waves are routinely observed in coronal loops. Resonant absorption is a well-accepted mechanism that extracts energy from kink waves. Nonlinear kink waves are know to be affected by the Kelvin-Helmholtz instability. However, all previous numerical studies consider linearly polarized kink waves. Aims. We study the properties of circularly polarized kink waves on straight plasma cylinders, for both standing and propagating waves, and compare them to the properties of linearly polarized kink waves. Methods. We use the code MPI-AMRVAC to solve the full 3D Magnetohydrodynamic (MHD) equations for a straight magnetic cylinder, excited by both standing and propagating circularly polarized kink (m = 1) modes. Results. The damping due to resonant absorption is independent of the polarization state. The morphology or appearance of the induced resonant flow is different for the two polarizations, however, there are essentially no differences in the forward-modeled Doppler signals. For nonlinear oscillations, the growth rate of small scales is determined by the total energy of the oscillation rather than the perturbation amplitude. We discuss possible implications and seismological relevance.

preprint2020arXiv

Excitation of negative energy surface magnetohydrodynamic waves in an incompressible cylindrical plasma

Negative energy wave phenomena may appear in shear flows in the presence of a wave decay mechanism and external energy supply. We study the appearance of negative energy surface waves in a plasma cylinder in the incompressible limit. The cylinder is surrounded by an axial magnetic field and by a plasma of different density. Considering flow inside and viscosity outside the flux tube, we derive dispersion relations, and obtain analytical solutions for the phase speed and growth rate (increment) of the waves. It is found that the critical speed shear for the occurrence of the dissipative instability associated with negative energy waves (NEWs) and the threshold of Kelvin--Helmholtz instability (KHI) depend on the axial wavelength. The critical shear for the appearance of sausage NEW is lowest for the longest axial wavelengths, while for kink waves the minimum value of the critical shear is reached for the axial wavelength comparable to the diameter of the cylinder. The range between the critical speed of the dissipative instability and the KHI threshold is shown to depend on the difference of the Alfvén speeds inside and outside of the cylinder. For all axial wavenumbers, NEW appears for the shear flow speeds lower than the KHI threshold. It is easier to excite NEW in an underdense cylinder than in an overdense one. The negative energy surface waves can be effectively generated for azimuthal number $m=0$ with a large axial wave number and for higher modes ($m>0$) with a small axial wave number.

preprint2020arXiv

Standing kink waves in sigmoid solar coronal loops: implications for coronal seismology

Using full three-dimensional magnetohydrodynamic numerical simulations, we study the effects of magnetic field sigmoidity or helicity on the properties of the fundamental kink oscillation of solar coronal loops. Our model consists of a single denser coronal loop, embedded in a plasma with dipolar force-free magnetic field with a constant alpha-parameter. For the loop with no sigmoidity, we find that the numerically determined oscillation period of the fundamental kink mode matches the theoretical period calculated using WKB theory. In contrast, with increasing sigmoidity of the loop, the actual period is increasingly smaller than the one estimated by WKB theory. Translated through coronal seismology, increasing sigmoidity results in magnetic field estimates which are increasingly shifting towards higher values, and even surpassing the average value for the highest alpha value considered. Nevertheless, the estimated range of the coronal magnetic field value lies within the mimimal/maximal limits, proving the robustness coronal seismology. We propose that the discrepancy in the estimations of the absolute value of the force-free magnetic field could be exploited seismologically to determine the free energy of coronal loops, if averages of the internal magnetic field and density can be reliably estimated by other methods.

preprint2020arXiv

The effect of magnetic field on the damping of slow waves in the solar corona

Slow magnetoacoustic waves are routinely observed in astrophysical plasma systems such as the solar corona. As a slow wave propagates through a plasma, it modifies the equilibrium quantities of density, temperature, and magnetic field. In the corona and other plasma systems, the thermal equilibrium is comprised of a balance between continuous heating and cooling processes, the magnitudes of which vary with density, temperature and magnetic field. Thus the wave may induce a misbalance between these competing processes. Its back reaction on the wave has been shown to lead to dispersion, and amplification or damping, of the wave. In this work the importance of the effect of magnetic field in the rapid damping of slow waves in the solar corona by heating/cooling misbalance is evaluated and compared to the effects of thermal conduction. The two timescales characterising the effect of misbalance are derived and calculated for plasma systems with a range of typical coronal conditions. The predicted damping times of slow waves from thermal misbalance in the solar corona are found to be of the order of 10-100 minutes, coinciding with the wave periods and damping times observed. Moreover the slow wave damping by thermal misbalance is found to be comparable to the damping by field-aligned thermal conduction. We show that in the infinite field limit, the wave dynamics is insensitive to the dependence of the heating function on the magnetic field, and this approximation is found to be valid in the corona so long as the magnetic field strength is greater than 10G for quiescent loops and plumes and 100G for hot and dense loops. In summary thermal misbalance may damp slow magnetoacoustic waves rapidly in much of the corona, and its inclusion in our understanding of slow mode damping may resolve discrepancies between observations and theory relying on compressive viscosity and thermal conduction alone.

preprint2019arXiv

Excitation of decay-less transverse oscillations of coronal loops by random motions

The relatively large-amplitude decaying regime of transverse oscillations of coronal loops has been known for two decades and interpreted in terms of MHD kink modes of cylindrical plasma waveguides. Recent observational analysis has revealed the so-called decay-less small-amplitude oscillations, with a multi-harmonic structure being detected. Several models have been proposed to explain them. In particular, decay-less oscillations have been described in terms of standing kink waves driven with continuous monoperiodic motions of loop footpoints, in terms of a simple oscillator model of forced oscillations due to harmonic external force, and as a self-oscillatory process due to the interaction of a loop with quasi-steady flows. However, an alternative mechanism is needed to explain the simultaneous excitation of several longitudinal harmonics of the oscillation. We study the mechanism of random excitation of decay-less transverse oscillations of coronal loops. With a spatially one-dimensional and time-dependent analytical model taking into account effects of the wave damping and kink speed variation along the loop, we consider transverse loop oscillations driven by random motions of footpoints. The footpoint motions are modelled by broad-band coloured noise. We have found the excitation of loop eigenmodes and analysed their frequency ratios as well as the spatial structure of the oscillations along the loop. The obtained results successfully reproduce the observed properties of decay-less oscillations. In particular, excitation of eigenmodes of a loop as a resonator can explain the observed quasi-monochromatic nature of decay-less oscillations and generation of multiple harmonics detected recently. We propose the mechanism that can interpret decay-less transverse oscillations of coronal loops in terms of kink waves randomly driven at the loop footpoints.

preprint2019arXiv

Formation of quasi-periodic slow magnetoacoustic wave trains by the heating/cooling misbalance

Slow magnetoacoustic waves are omnipresent in both natural and laboratory plasma systems. The wave-induced misbalance between plasma cooling and heating processes causes the amplification or attenuation, and also dispersion, of slow magnetoacoustic waves. The wave dispersion could be attributed to the presence of characteristic time scales in the system, connected with the plasma heating or cooling due to the competition of the heating and cooling processes in the vicinity of the thermal equilibrium. We analysed linear slow magnetoacoustic waves in a plasma in a thermal equilibrium formed by a balance of optically thin radiative losses, field-align thermal conduction, and an unspecified heating. The dispersion is manifested by the dependence of the effective adiabatic index of the wave on the wave frequency, making the phase and group speeds frequency-dependent. The mutual effect of the wave amplification and dispersion is shown to result into the occurrence of an oscillatory pattern in an initially broadband slow wave, with the characteristic period determined by the thermal misbalance time scales, i.e. by the derivatives of the combined radiation loss and heating function with respect to the density and temperature, evaluated at the equilibrium. This effect is illustrated by estimating the characteristic period of the oscillatory pattern, appearing because of thermal misbalance in the plasma of the solar corona. It is found that by an order of magnitude the period is about the typical periods of slow magnetoacoustic oscillations detected in the corona.

preprint2019arXiv

Quasi-periodic pulsations of gamma-ray emissions from a solar flare on 2017 September 06

We investigate quasi-periodic pulsations (QPPs) of high-energy nonthermal emissions from an X9.3 flare (SOL2017-Sep-06T11:53), the most powerful flare since the beginning of solar cycle 24. The QPPs are identified as a series of regular and repeating peaks in the light curves in the gamma- and hard X-ray (HXR) channels recorded by the Konus-Wind, as well as the radio and microwave fluxes measured by the CALLISTO radio spectrograph during the impulsive phase. The periods are determined from the global wavelet and Fourier power spectra, as 24-30 s in the HXR and microwave channels which are associated with nonthermal electrons, and ~20 s in the gamma-ray band related to nonthermal ions. Both nonthermal electrons and ions may be accelerated by repetitive magnetic reconnection during the impulsive phase. However, we could not rule out other mechanisms such as the MHD oscillation in a sausage mode. The QPP detected in this study is useful for understanding the particle acceleration and dynamic process in solar flares and also bridging the gap between stellar and solar flares since the energy realm of the X9.3 solar flare is almost compared with a typical stellar flare.

preprint2016arXiv

Effect of a Radiation Cooling and Heating Function on Standing Longitudinal Oscillations in Coronal Loops

Standing long-period (with the periods longer than several minutes) oscillations in large hot (with the temperature higher than 3 MK) coronal loops have been observed as the quasi-periodic modulation of the EUV and microwave intensity emission and the Doppler shift of coronal emission lines, and have been interpreted as standing slow magnetoacoustic (longitudinal) oscillations. Quasi-periodic pulsations of shorter periods, detected in thermal and non-thermal emissions in solar flares could be produced by a similar mechanism. We present theoretical modelling of the standing slow magnetoacoustic mode, showing that this mode of oscillation is highly sensitive to peculiarities of the radiative cooling and heating function. We generalised the theoretical model of standing slow magnetoacoustic oscillations in a hot plasma, including the effects of the radiative losses, and accounting for plasma heating. The heating mechanism is not specified and taken empirically to compensate the cooling by radiation and thermal-conduction. It is shown that the evolution of the oscillations is described by a generalised Burgers equation. Numerical solution of an initial value problem for the evolutionary equation demonstrates that different dependences of the radiative cooling and plasma heating on the temperature lead to different regimes of the oscillations, including growing, quasi-stationary and rapidly decaying. Our findings provide a theoretical foundation for probing the coronal heating function, and may explain the observations of decayless long-period quasi-periodic pulsations in flares. The hydrodynamic approach employed in this study should be considered with caution in the modelling of non-thermal emission associated with flares, as it misses potentially important non-hydrodynamic effects.

preprint2016arXiv

Effect of a sausage oscillation on radio zebra-pattern structures in a solar flare

Sausage modes that are axisymmetric fast magnetoacoustic oscillations of solar coronal loops are characterized by variation of the plasma density and magnetic field, and hence cause time variations of the electron plasma frequency and cyclotron frequency. The latter parameters determine the condition for the double plasma resonance (DPR), which is responsible for the appearance of zebra-pattern (ZP) structures in time spectra of solar type IV radio bursts. We perform numerical simulations of standing and propagating sausage oscillations in a coronal loop modeled as a straight, field-aligned plasma slab, and determine the time variation of the DPR layer locations. Instant values of the plasma density and magnetic field at the DPR layers allowed us to construct skeletons of the time variation of ZP stripes in radio spectra. In the presence of a sausage oscillation, the ZP structures are shown to have characteristic wiggles with the time period prescribed by the sausage oscillation. Standing and propagating sausage oscillations are found to have different signatures in ZP patterns. We conclude that ZP wiggles can be used for the detection of short-period sausage oscillations and the exploitation of their seismological potential.

preprint2016arXiv

Nonlinear waves in the terrestrial quasi-parallel foreshock

We study the applicability of the derivative nonlinear Schrödinger (DNLS) equation, for the evolution of high frequency nonlinear waves, observed at the foreshock region of the terrestrial quasi-parallel bow shock. The use of a pseudo-potential is elucidated and, in particular, the importance of canonical representation in the correct interpretation of solutions in this formulation is discussed. Numerical solutions of the DNLS equation are then compared directly with the wave forms observed by Cluster spacecraft. Non harmonic slow variations are filtered out by applying the empirical mode decomposition. We find large amplitude nonlinear wave trains at frequencies above the proton cyclotron frequency, followed in time by nearly harmonic low amplitude fluctuations. The approximate phase speed of these nonlinear waves, indicated by the parameters of numerical solutions, is of the order of the local Alfvén speed.

preprint2016arXiv

Observation of quasi-periodic solar radio bursts associated with propagating fast-mode waves

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, $δt_\mathrm{r}$=1.78$\pm$0.04 min, matches the period of the fast wave train observed at 171 $Å$, $δ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.

preprint2016arXiv

Solar science with the Atacama Large Millimeter/submillimeter Array - A new view of our Sun

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

preprint2016arXiv

Statistical Properties of Quasi-Periodic Pulsations in White-Light Flares Observed With Kepler

We embark on a study of quasi-periodic pulsations (QPPs) in the decay phase of white-light stellar flares observed by Kepler. Out of the 1439 flares on 216 different stars detected in the short-cadence data using an automated search, 56 flares are found to have pronounced QPP-like signatures in the light curve, of which 11 have stable decaying oscillations. No correlation is found between the QPP period and the stellar temperature, radius, rotation period and surface gravity, suggesting that the QPPs are independent of global stellar parameters. Hence they are likely to be the result of processes occurring in the local environment. There is also no significant correlation between the QPP period and flare energy, however there is evidence that the period scales with the QPP decay time for the Gaussian damping scenario, but not to a significant degree for the exponentially damped case. This same scaling has been observed for MHD oscillations on the Sun, suggesting that they could be the cause of the QPPs in those flares. Scaling laws of the flare energy are also investigated, supporting previous reports of a strong correlation between the flare energy and stellar temperature/radius. A negative correlation between the flare energy and stellar surface gravity is also found.

preprint2016arXiv

Statistically determined dispersion relations of magnetic field fluctuations in the terrestrial foreshock

We obtain dispersion relations of magnetic field fluctuations for two crossings of the terrestrial foreshock by Cluster spacecraft. These crossings cover plasma conditions which differ significantly in their plasma $β$ but not in the properties of the encountered ion population, both showing shell-like distribution function. Dispersion relations are reconstructed using two-point instantaneous wave number estimations from pairs of Cluster spacecraft. The accessible range of wave vectors, limited by the available spacecraft separations, extends to $\approx \!2 \times 10^4$ km. Results show multiple branches of dispersion relation, associated with different power of magnetic field fluctuations. We find that Sunward propagating fast magnetosonic waves and beam resonant modes are dominant for the high plasma $β$ interval, while the dispersions of magnetically dominated interval include Alfvén and fast magnetosonic modes propagating Sunward and unit-Sunward.

preprint2015arXiv

A Comparison Between Global Proxies of the Sun's Magnetic Activity Cycle: Inferences from Helioseismology

The last solar minimum was, by recent standards, unusually deep and long. We are now close to the maximum of the subsequent solar cycle, which is relatively weak. In this article we make comparisons between different global (unresolved) measures of the Sun's magnetic activity, to investigate how they are responding to this weak-activity epoch. We focus on helioseismic data, which are sensitive to conditions, including the characteristics of the magnetic field, in the solar interior. Also considered are measures of the magnetic field in the photosphere (sunspot number and sunspot area), the chromosphere and corona (10.7cm radio flux and 530.3nm green coronal index), and two measures of the Sun's magnetic activity closer to Earth (the interplanetary magnetic field and the galactic cosmic-ray intensity). Scaled versions of the activity proxies diverge from the helioseismic data around 2000, indicating a change in relationship between the proxies. The degree of divergence varies from proxy to proxy with sunspot area and 10.7cm flux showing only small deviations, while sunspot number, coronal index, and the two interplanetary proxies show much larger departures. In Cycle 24 the deviations in the solar proxies and the helioseismic data decrease, raising the possibility that the deviations observed in Cycle 23 are just symptomatic of a 22-year Hale cycle. However, the deviations in the helioseismic data and the interplanetary proxies increase in Cycle 24. Interestingly the divergence in the solar proxies and the helioseismic data are not reflected in the shorter-term variations (often referred to as quasi-biennial oscillations) observed on top of the dominant 11-year solar cycle. However, despite being highly correlated in Cycle 22, the short-term variations in the interplanetary proxies show very little correlation with the helioseismic data during Cycles 23 and 24.

preprint2015arXiv

A statistical study of decaying kink oscillations detected using SDO/AIA

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.

preprint2015arXiv

Decayless low-amplitude kink oscillations: a common phenomenon in the solar corona?

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.

preprint2015arXiv

Fast magnetoacoustic wave trains of sausage symmetry in cylindrical waveguides of the solar corona

Fast magnetoacoustic waves guided along the magnetic field by plasma non-uniformities, in particular coronal loops, fibrils and plumes, are known to be highly dispersive, which leads to the formation of quasi-periodic wave trains excited by a broadband impulsive driver, e.g. a solar flare. We investigated effects of cylindrical geometry on the fast sausage wave train formation. We performed magnetohydrodynamic numerical simulations of fast magnetoacoustic perturbations of a sausage symmetry, propagating from a localised impulsive source along a field-aligned plasma cylinder with a smooth radial profile of the fast speed. The wave trains are found to have pronounced period modulation, with the longer instant period seen in the beginning of the wave train. The wave trains have also a pronounced amplitude modulation. Wavelet spectra of the wave trains have characteristic tadpole features, with the broadband large-amplitude heads preceding low-amplitude quasi-monochromatic tails. The mean period of the wave train is about the transverse fast magnetoacoustic transit time across the cylinder. The mean parallel wavelength is about the diameter of the waveguiding plasma cylinder. Instant periods are longer than the sausage wave cutoff period. The wave train characteristics depend on the fast magnetoacoustic speed in both the internal and external media, and the smoothness of the transverse profile of the equilibrium quantities, and also the spatial size of the initial perturbation. If the initial perturbation is localised at the axis of the cylinder, the wave trains contain higher radial harmonics that have shorter periods.

preprint2015arXiv

Oscillations in stellar superflares

Two different mechanisms may act to induce quasi-periodic pulsations (QPP) in whole-disk observations of stellar flares. One mechanism may be magneto-hydromagnetic (MHD) forces and other processes acting on flare loops as seen in the Sun. The other mechanism may be forced local acoustic oscillations due to the high-energy particle impulse generated by the flare (known as `sunquakes' in the Sun). We analyze short-cadence Kepler data of 257 flares in 75 stars to search for QPP in the flare decay branch or post-flare oscillations which may be attributed to either of these two mechanisms. About 18 percent of stellar flares show a distinct bump in the flare decay branch of unknown origin. The bump does not seem to be a highly-damped global oscillation because the periods of the bumps derived from wavelet analysis do not correlate with any stellar parameter. We detected damped oscillations covering several cycles (QPP), in seven flares on five stars. The periods of these oscillations also do not correlate with any stellar parameter, suggesting that these may be a due to flare loop oscillations. We searched for forced global oscillations which might result after a strong flare. To this end, we investigated the behaviour of the amplitudes of solar-like oscillations in eight stars before and after a flare. However, no clear amplitude change could be detected. We also analyzed the amplitudes of the self-excited pulsations in two delta Scuti stars and one gamma Doradus star before and after a flare. Again, no clear amplitude changes were found. Our conclusions are that a new process needs to be found to explain the high incidence of bumps in stellar flare light curves, that flare loop oscillations may have been detected in a few stars and that no conclusive evidence exists as yet for flare induced global acoustic oscillations (starquakes).

preprint2015arXiv

SSALMON - The Solar Simulations for the Atacama Large Millimeter Observatory Network

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

preprint2014arXiv

Evolution of fast magnetoacoustic pulses in randomly structured coronal plasmas

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-$β$ 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.

preprint2014arXiv

Multi-height observations of magnetoacoustic cut-off frequency in a sunspot atmosphere

We used the Pixelised Wavelet Filtering (PWF) method to compute narrow-band power maps of SDO/AIA imaging datasets in the 1700 Å, 1600 Å and 304 Å bandpasses that correspond to different heights. The cut-off frequency was defined as contours where the spectral power droped to the median level. It was measured as a function of the spatial location. We inferred the magnetic field inclination according to the MAG wave theory in the low-$β$ limit and compared it with the potential field extrapolation.{We analysed intensity oscillations in a symmetric sunspot AR11131(08 Dec 2010) and an asymmetric sunspot AR11330 (27 Oct 2011). We reconstructed the magnetic field inclination in the radial direction for the symmetric sunspot and in both radial and azimuthal directions for the asymmetric sunspot. We observed 3D variation of the main oscillation periods in sunspots. We found that shorter-period oscillations were mostly constrained in sunspot umbrae, while longer-period oscillations formed an annular shape enclosing the umbra. Longer periods are found to be distributed further away from the sunspot centre. Our results indicate that 3-min oscillation are generated in the chromosphere, possibly by the acoustic resonator model, while 5-min and longer-period oscillations seemed to originate in a level under the photosphere. The reconstructed field inclinations gives the values of the field inclination that are systematically larger than the values obtained by the potential field extrapolation. The inclined magnetic field line can account for 60-80% of cut-offfrequency lowering only.

preprint2014arXiv

Wave dynamics in a sunspot umbra

The high spatial and time resolution data obtained with SDO/AIA for the sunspot in active region NOAA 11131 on 08 December 2010 were analysed with the time-distance plot technique and the pixelised wavelet filtering method. Oscillations in the 3 min band dominate in the umbra. The integrated spectrum of umbral oscillations contains distinct narrowband peaks at 1.9 min, 2.3 min, and 2.8 min. The power significantly varies in time, forming distinct oscillation trains. The oscillation power distribution over the sunspot in the horizontal plane reveals that the enhancements of the oscillation amplitude, or wave fronts, have a distinct structure consisting of an evolving two-armed spiral and a stationary circular patch at the spiral origin, situated near the umbra centre. This structure is seen from the temperature minimum to the corona. In time, the spiral rotates anti-clockwise. The wave front spirality is most pronounced during the maximum amplitude phases of the oscillations. In the low-amplitude phases the spiral breaks into arc-shaped patches. The 2D cross-correlation function shows that the oscillations at higher atmospheric levels occur later than at lower layers. The phase speed is estimated to be about 100 km/s. The fine spectral analysis shows that the central patch corresponds to the high-frequency oscillations, while the spiral arms highlight the lower-frequency oscillations in the 3-min band. The vertical and horizontal radial structure of the oscillations is consistent with the model that interprets umbral oscillations as slow magnetoacoustic waves filtered by the atmospheric temperature non-uniformity in the presence of the magnetic field inclination from the vertical. The mechanism for the polar-angle structure of the oscillations, in particular the spirality of the wave fronts, needs to be revealed.

preprint2013arXiv

Characteristics of magnetoacoustic sausage modes

Aims: We perform an advanced study of the fast magnetoacoustic sausage oscillations of coronal loops in the context of MHD coronal seismology to establish the dependence of the sausage mode period and cut-off wavenumber on the plasma-beta of the loop-filling plasma. A parametric study of the ratios for different harmonics of the mode is also carried out. Methods: Full magnetohydrodynamic numerical simulations were performed using Lare2d, simulating hot, dense loops in a magnetic slab environment. The symmetric Epstein profile and a simple step-function profile were both used to model the density structure of the simulated loops. Analytical expressions for the cut-off wavenumber and the harmonic ratio between the second longitudinal harmonic and the fundamental were also examined. Results: It was established that the period of the global sausage mode is only very weakly dependent on the value of the plasma-beta inside a coronal loop, which justifies the application of this model to hot flaring loops. The cut-off wavenumber kc for the global mode was found to be dependent on both internal and external values of the plasma-beta, again only weakly. By far the most important factor in this case was the value of the density contrast ratio between the loop and the surroundings. Finally, the deviation of the harmonic ratio P1/2P2 from the ideal non-dispersive case was shown to be considerable at low k, again strongly dependent on plasma density. Quantifying the behaviour of the cut-off wavenumber and the harmonic ratio has significant applications to the field of coronal seismology.

preprint2013arXiv

Quasi-periodic wiggles of microwave zebra structures in a solar flare

Quasi-periodic wiggles of microwave zebra pattern structures with period range from about 0.5 s to 1.5 s are found in a X-class solar flare on 2006 December 13 at the 2.6-3.8 GHz with the Chinese Solar Broadband Radio Spectrometer (SBRS/Huairou). Periodogram and correlation analysis show that the wiggles have two-three significant periodicities and almost in phase between stripes at different frequency. The Alfven speed estimated from the zebra pattern structures is about 700 Km/s. We obtain the spatial size of the waveguiding plasma structure to be about 1 Mm with the detected period of about 1 s. It suggests the ZP wiggles can be associated with the fast mag- netoacoustic oscillations in the flaring active region. The lack of a significant phase shift between wiggles of different stripes suggests that the ZP wiggles are caused by a standing sausage oscillation.

preprint2013arXiv

Slow Magnetoacoustic Oscillations in the Microwave Emission of Solar Flares

Analysis of the microwave data, obtained in the 17 GHz channel of the Nobeyama Radioheliograph during the M1.6 flare on 4th Nov 2010, revealed the presence of 11.8-min oscillations of the emitting plasma density. The oscilla- tions decayed with the characteristic time of about 25-min. These oscillations are also well-seen in the variation of EUV emission intensity measured in the 335 A channel of SDO/AIA. The observed properties of the oscillations are consistent with the properties of so-called SUMER oscillations, observed in the EUV and soft X-ray bands usually as a periodic Doppler shift. The accepted interpretation of SUMER oscillations is a standing slow magnetoacoustic wave. Our analysis presents the first direct observation of the slow magnetoacoustic oscillations in the microwave emission of a solar flare.

preprint2011arXiv

Frequency drifts of 3-min oscillations in microwave and EUV emission above sunspots

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.

preprint2011arXiv

Height distribution of the power of 3-min oscillations over sunspots

Context. The height structure of 3-min oscillations over sunspots is studied in the context of the recently discovered effect of height inversion: over the umbra, the spatial location of the maximum of chromospheric 3-min oscillation power corresponds to the relative decrease in the power of photospheric oscillations. Aims. We investigate whether the height inversion of the power of 3-min oscillations is a common feature of the spatial structure of the oscillations for the majority of sunspots. Methods. Spectrogram sequences of H-alpha 6563Å and FeI 6569Å over sunspots, acquired with very high cadency (about 2 s or better) are obtained. The distribution of the oscillation power of the line-of-sight velocity signal is studied by using methods of wavelet frequency filtration and Fourier analysis. Results. The effect of the height inversion is found in 9 of 11 analyzed active regions. The interpretation of this effect is possibly connected to both the decrease in the level of photosphere in sunspot umbrae and the magnetic field topology.

preprint2011arXiv

Instrumental oscillations in RHESSI count rates during solar flares

Aims: We seek to illustrate the analysis problems posed by RHESSI spacecraft motion by studying persistent instrumental oscillations found in the lightcurves measured by RHESSI's X-ray detectors in the 6-12 keV and 12-25 keV energy range during the decay phase of the flares of 2004 November 4 and 6. Methods: The various motions of the RHESSI spacecraft which may contribute to the manifestation of oscillations are studied. The response of each detector in turn is also investigated. Results: We find that on 2004 November 6 the observed oscillations correspond to the nutation period of the RHESSI instrument. These oscillations are also of greatest amplitude for detector 5, while in the lightcurves of many other detectors the oscillations are small or undetectable. We also find that the variation in detector pointing is much larger during this flare than the counterexample of 2004 November 4. Conclusions: Sufficiently large nutation motions of the RHESSI spacecraft lead to clearly observable oscillations in count rates, posing a significant hazard for data analysis. This issue is particularly problematic for detector 5 due to its design characteristics. Dynamic correction of the RHESSI counts, accounting for the livetime, data gaps, and the transmission of the bi-grid collimator of each detector, is required to overcome this issue. These corrections should be applied to all future oscillation studies.

preprint2011arXiv

Three-minute oscillations above sunspot umbra observed with SDO/AIA and NoRH

Three-minute oscillations over sunspot's umbra in AR 11131 were observed simultaneously in UV/EUV emission by SDO/AIA and in radio emission by Nobeyama Radioheliograph (NoRH). We use 24-hours series of SDO and 8-hours series of NoRH observations to study spectral, spatial and temporal variations of pulsations in the 5-9 mHz frequency range at different layers of the solar atmosphere. High spatial and temporal resolution of SDO/AIA in combination with long-duration observations allowed us to trace the variations of the cut-off frequency and spectrum of oscillations across the umbra. We found that higher frequency oscillations are more pronounced closer to the umbra's center, while the lower frequencies concentrate to the peripheral parts. We interpreted this discovery as a manifestation of variation of the magnetic field inclination across the umbra at the level of temperature-minimum. Possible implications of this interpretation for the diagnostics of sunspot atmospheres is discussed.

preprint2011arXiv

Transient induced MHD oscillations : A tool to probe the solar active regions

Solar transients and eruptive phenomena which are ubiquitous in the solar atmosphere, can shed new light to the understanding of the outstanding problems like coronal heating and the solar wind acceleration. Observations in the entire electromagnetic spectrum of such dynamical processes of large and small-scale transient/eruptive events, with highly dynamic magnetic field configuration, and energetic particles, provide crucial information about the plasma processes at mega-Kelvin temperature embedded in a complex magnetic field, and also energy build-up/energy-release processes, taking place in such events. One of the most important phenomenological aspects of solar eruptive phenomena is the induced magnetohydrodynamic (MHD) waves generated during these energetic processes, which carry a potential signature to probing the solar active regions. In this paper, we briefly review the recent trends of the transient (e.g., flares) induced quasi-periodic oscillations in the solar atmosphere and discuss their implications in diagnosing the solar active regions, providing the clue to understanding local plasma dynamics and heating.

preprint2010arXiv

Oscillatory processes in solar flares

Electromagnetic (radio, visible-light, UV, EUV, X-ray and gamma-ray) emission generated by solar and stellar flares often contains pronounced quasi-periodic pulsations (QPP). Physical mechanisms responsible for the generation of long-period QPP (with the periods longer than one second) are likely to be associated with MHD processes. The observed modulation depths, periods and anharmonicity of QPP suggest that they can be linked with some kind of MHD auto-oscillations, e.g. an oscillatory regime of magnetic reconnection. Such regimes, of both spontaneous and induced nature, have been observed in resistive-MHD numerical simulations. The oscillations are essentially nonlinear and non-stationary. We demonstrate that a promising novel method for their analysis is the Empirical Mode Decomposition technique.

preprint2010arXiv

Relationship between wave processes in sunspots and quasi-periodic pulsations in active region flares

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.

preprint2010arXiv

Web-Based Data Processing System for Automated Detection of Oscillations with Applications to the Solar Atmosphere

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.

preprint2002arXiv

Line-of-sight Effects on Observability of Kink and Sausage Modes in Coronal Structures with Imaging Telescopes

Kink modes of solar coronal structures, perturbing the loop in the direction along the line-of-sight (LOS), can be observed as emission intensity disturbances propagating along the loop provided the angle between the LOS and the structure is not ninety degrees. The effect is based upon the change of the column depth of the loop (along the LOS) by the wave. The observed amplitude of the emission intensity variations can be larger than the actual amplitude of the wave by a factor of two and there is an optimal angle maximizing the observed amplitude. For other angles this effect can also attenuate the observed wave amplitude. The observed amplitude depends upon the ratio of the wave length of kink perturbations to the width of the structure and on the angle between the LOS and the axis of the structure. Sausage modes are always affected negatively from the observational point of view, as the observed amplitude is always less than the actual one. This effect should be taken into account in the interpretation of wave phenomena observed in the corona with space-borne and ground-based imaging telescopes.