Search

The multi-wavelength data from the Solar Dynamics Observatory (SDO) is extensively used in studying the physics of the Sun and its atmosphere. In this study, we estimate the formation heights of low-corona and chromospheric channels of the Atmospheric Imaging Assembly (AIA) over the atmospheres of sunspot umbrae during the quiet condition period within 20 different active regions. The upward propagating slow magnetoacoustic waves (slow MAWs) of 3-min period, which are perpetually present in sunspots, are utilized for this purpose. Employing a cross-correlation technique, the most frequent time lag between different channel pairs is measured. By combining this information with the local sound speed obtained from the characteristic formation temperatures of individual channels, we estimate the respective formation heights. The median values of formation heights obtained across all active regions in our sample are 356, 368, 858, 1180, and 1470 km, respectively, for the AIA 1600 {\AA}, 1700 {\AA}, 304 {\AA}, 131 {\AA}, and 171 {\AA} channels. The corresponding ranges in the formation heights are 247 $\--$ 453, 260 $\--$ 468, 575 $\--$ 1155, 709 $\--$ 1937, and 909 $\--$ 2585 km, respectively. These values are measured with respect to the HMI continuum. We find the formation height of UV channels is quite stable (between 250 $\--$ 500 km) and displays only a marginal difference between the AIA 1600 {\AA} and 1700 {\AA} during quiet conditions. On the other hand, the formation height of coronal channels is quite variable., Comment: accepted for publication in ApJ

One of the goals of Space Weather studies is to achieve a better understanding of impulsive phenomena, such as Coronal Mass Ejections (CMEs), in order to improve our ability to forecast them and mitigate the risk to our technologically driven society. The essential part of achieving this goal is to assess the performance of forecasting models. To this end, the quality and availability of suitable data are of paramount importance. In this work, we have merged already publicly available data of CMEs from both in-situ and remote instrumentation in order to build a database of CME properties. To evaluate the accuracy of such a database and confirm the relationship between in-situ and remote observations, we have employed the drag-based model (DBM) due to its simplicity and inexpensive cost of computational resources. In this study, we have also explored the parameter space for the drag parameter and solar wind speed using a Monte Carlo approach to evaluate how well the DBM determines the propagation of CMEs for the events in the dataset. The dataset of geoeffective CMEs constructed as a result of this work provides validation of the initial hypothesis about DBM, and solar wind speed and also yields further insight into CME features like arrival time, arrival speed, lift-off time, etc. Using a data-driven approach, this procedure allows us to present a homogeneous, reliable, and robust dataset for the investigation of CME propagation. On the other hand, possible CME events are identified where DBM approximation is not valid due to model limitations and higher uncertainties in the input parameters, those events require more thorough investigation.

We present the first quiet Sun spectropolarimetric observations obtained with the Visible SpectroPolarimeter (ViSP) at the $4-$m Daniel K. Inouye Solar Telescope (DKIST). We recorded observations in a wavelength range that includes the magnetically sensitive Fe I $6301.5/6302.5$ $\AA$ doublet. With an estimated spatial resolution of 0.08'', this represents the highest spatial resolution full-vector spectropolarimetric observations ever obtained of the quiet Sun. We identified $53$ small-scale magnetic elements, including $47$ magnetic loops and $4$ unipolar magnetic patches, with linear and circular polarisation detected in all of them. Of particular interest is a magnetic element in which the polarity of the magnetic vector appears to change three times in only $400$ km and which has linear polarisation signals throughout. We find complex Stokes $V$ profiles at the polarity inversion lines of magnetic loops and discover degenerate solutions, as we are unable to conclusively determine whether these arise due to gradients in the atmospheric parameters or smearing of opposite polarity signals. We analyse a granule which notably has linear and circular polarisation signals throughout, providing an opportunity to explore its magnetic properties. On this small scale we see the magnetic field strength range from $25$ G at the granular boundary to $2$ kG in the intergranular lane (IGL), and sanity check the values with the weak and strong field approximations. A value of $2$ kG in the IGL is among the highest measurements ever recorded for the internetwork., Comment: Accepted for publication in ApJL

Observational pre-cursors of large solar flares provide a basis for future operational systems for forecasting. Here, we study the evolution of the normalized emergence (EM), shearing (SH) and total (T) magnetic helicity flux components for 14 flaring with at least one X-class flare) and 14 non-flaring ($<$ M5-class flares) active regions (ARs) using the Spaceweather Helioseismic Magnetic Imager Active Region Patches vector magnetic field data. Each of the selected ARs contain a $\delta$-type spot. The three helicity components of these ARs were analyzed using wavelet analysis. Localised peaks of the wavelet power spectrum (WPS) were identified and statistically investigated. We find that: i) the probability density function of the identified WPS peaks for all the EM, SH and T profiles can be fitted with a set of Gaussian functions centered at distinct periods between $\sim$ 3 to 20 hours. ii) There is a noticeable difference in the distribution of periods found in the EM profiles between the flaring and non-flaring ARs, while no significant difference is found in the SH and T profiles. iii) In flaring ARs, the distributions of the shorter EM/SH/T periods ($<$ 10 hrs) split up into two groups after flares, while the longer periods ($>$ 10 hrs) do not change. iv) When the EM periodicity does not contain harmonics, the ARs do not host a large energetic flare. Finally, v) significant power at long periods ($\sim$ 20 hour) in the T and EM components may serve as pre-cursor for large energetic flares.

In review of the White Papers from the Voyage 2050 process and after the public presentation of a number of these papers in October 2019 in Madrid, we as White Paper lead authors have identified a coherent science theme that transcends the divisions around which the Topical Teams are structured. This note aims to highlight this synergistic science theme and to make the Topical Teams and the Voyage 2050 Senior Committee aware of the wide importance of these topics and the broad support that they have across the worldwide science community., Comment: 5 pages, published in Frontiers in Astronomy and Space Science

Whilst the most dynamic solar active regions (ARs) are known to flare frequently, predicting the occurrence of individual flares and their magnitude, is very much a developing field with strong potentials for machine learning applications. The present work is based on a method which is developed to define numerical measures of the mixed states of ARs with opposite polarities. The method yields compelling evidence for the assumed connection between the level of mixed states of a given AR and the level of the solar eruptive probability of this AR by employing two morphological parameters: (i) the separation parameter $S_{l-f}$ and (ii) the sum of the horizontal magnetic gradient $G_{S}$. In this work, we study the efficiency of $S_{l-f}$ and $G_{S}$ as flare predictors on a representative sample of ARs, based on the SOHO/MDI-Debrecen Data (SDD) and the SDO/HMI - Debrecen Data (HMIDD) sunspot catalogues. In particular, we investigate about 1000 ARs in order to test and validate the joint prediction capabilities of the two morphological parameters by applying the logistic regression machine learning method. Here, we confirm that the two parameters with their threshold values are, when applied together, good complementary predictors. Furthermore, the prediction probability of these predictor parameters is given at least 70\% a day before.

The evolution of magnetic helicity has a close relationship with solar eruptions and is of interest as a predictive diagnostic. In this case study, we analyse the evolution of the normalised emergence, shearing and total magnetic helicity components in the case of three flaring and three non-flaring active regions (ARs) using SHARPs (Spaceweather Helioseismic Magnetic Imager Active Region Patches) vector magnetic field data. The evolution of the three magnetic helicity components is analysed with wavelet transforms, revealing significant common periodicities of the normalised emergence, shearing and total helicity fluxes before flares in the flaring ARs. The three non-flaring ARs do not show such common periodic behaviour. This case study suggests that the presence of significant periodicities in the power spectrum of magnetic helicity components could serve as a valuable precursor for flares.

In this article, we present the application of the weighted horizontal gradient of magnetic field (WGM) flare prediction method to 3-dimensional (3D) extrapolated magnetic configurations of 13 flaring solar active regions (ARs). The main aim is to identify an optimal height range, if any, in the interface region between the photosphere and lower corona, where the flare onset time prediction capability of WGM is best exploited. The optimal height is where flare prediction, by means of the WGM method, is achieved earlier than at the photospheric level. 3D magnetic structures, based on potential and non-linear force-free field extrapolations, are constructed to study a vertical range from the photosphere up to the low corona with a 45 km step size. The WGM method is applied as a function of height to all 13 flaring AR cases that are subject to certain selection criteria. We found that applying the WGM method between 1000 and 1800 km above the solar surface would improve the prediction of the flare onset time by around 2-8 hours.Certain caveats and an outlook for future work along these lines are also discussed.

The European Solar Telescope (EST) is a research infrastructure for solar physics. It is planned to be an on-axis solar telescope with an aperture of 4 m and equipped with an innovative suite of spectro-polarimetric and imaging post-focus instrumentation. The EST project was initiated and is driven by EAST, the European Association for Solar Telescopes. EAST was founded in 2006 as an association of 14 European countries. Today, as of December 2019, EAST consists of 26 European research institutes from 18 European countries. The Preliminary Design Phase of EST was accomplished between 2008 and 2011. During this phase, in 2010, the first version of the EST Science Requirement Document (SRD) was published. After EST became a project on the ESFRI roadmap 2016, the preparatory phase started. The goal of the preparatory phase is to accomplish a final design for the telescope and the legal governance structure of EST. A major milestone on this path is to revisit and update the Science Requirement Document (SRD). The EST Science Advisory Group (SAG) has been constituted by EAST and the Board of the PRE-EST EU project in November 2017 and has been charged with the task of providing with a final statement on the science requirements for EST. Based on the conceptual design, the SRD update takes into account recent technical and scientific developments, to ensure that EST provides significant advancement beyond the current state-of-the-art. The present update of the EST SRD has been developed and discussed during a series of EST SAG meetings. The SRD develops the top-level science objectives of EST into individual science cases. Identifying critical science requirements is one of its main goals. Those requirements will define the capabilities of EST and the post-focus instrument suite. The technical requirements for the final design of EST will be derived from the SRD., Comment: 2nd edition, December 2019, 138 pages

Aims: In the present work we investigate the nature of the magnetic Rayleigh-Taylor instability at a density interface permeated by an oblique, homogeneous magnetic field in an incompressible limit. Methods: Using the system of linearised ideal incompressible magnetohydrodynamics (MHD) equations we derive the dispersion relation for perturbations of the contact discontinuity by imposing the necessary continuity conditions at the interface. The imaginary part of the frequency describes the growth rate of waves due to instability. The growth rate of waves is studied by solving numerically the dispersion relation. Results: The critical wavenumber at which waves become unstable, present for a parallel magnetic field, disappears, due to the inclination of the magnetic field and instead waves are shown to be unstable for all wavenumbers. Theoretical results are applied to diagnose the structure of the magnetic field in prominence threads. When applying our theoretical results to observed waves in prominence plumes, we obtain a wide range of field inclination angle; from 0.5 degrees up to 30 degrees. These results highlight the diagnostic possibility of our study., Comment: 9 pages, 13 figures

Numerous apparent signatures of magnetic reconnection have been reported in the solar photosphere, including inverted-Y shaped jets. The reconnection at these sites is expected to cause localised bi-directional flows and extended shock waves; however, these signatures are rarely observed as extremely high spatial-resolution data are required. Here, we use H-alpha imaging data sampled by the Swedish Solar Telescope's CRisp Imaging SpectroPolarimeter to investigate whether bi-directional flows can be detected within inverted-Y shaped jets near the solar limb. These jets are apparent in the H-alpha line wings, while no signature of either jet is observed in the H-alpha line core, implying reconnection took place below the chromospheric canopy. Asymmetries in the H-alpha line profiles along the legs of the jets indicate the presence of bi-directional flows, consistent with cartoon models of reconnection in chromospheric anemone jets. These asymmetries are present for over two minutes, longer than the lifetimes of Rapid Blue Excursions, and beyond \pm 1 \AA\ into the wings of the line indicating that flows within the inverted-Y shaped jets are responsible for the imbalance in the profiles, rather than motions in the foreground. Additionally, surges form following the occurrence of the inverted-Y shaped jets. This surge formation is consistent with models which suggest such events could be caused by the propagation of shock waves from reconnection sites in the photosphere to the upper atmosphere. Overall, our results provide evidence that magnetic reconnection in the photosphere can cause bi-directional flows within inverted-Y shaped jets and could be the driver of surges., Comment: 9 pages, 7 figures, accepted for publication in ApJ

Ever since their detection two decades ago, standing kink oscillations in coronal loops have been extensively studied both observationally and theoretically. Almost all driven coronal loop oscillations (e.g., by flares) are observed to damp through time often with Gaussian or exponential profiles. Intriguingly, however, it has been shown theoretically that the amplitudes of some oscillations could be modified from Gaussian or exponential profiles if cooling is present in the coronal loop systems. Indeed, in some cases the oscillation amplitude can even increase through time. In this article, we analyse a flare-driven coronal loop oscillation observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA) in order to investigate whether models of cooling can explain the amplitude profile of the oscillation and whether hints of cooling can be found in the intensity evolution of several SDO/AIA filters. During the oscillation of this loop system, the kink mode amplitude appears to differ from a typical Gaussian or exponential profile with some hints being present that the amplitude increases. The application of cooling coronal loop modelling allowed us to estimate the density ratio between the loop and the background plasma, with a ratio of between 2.05-2.35 being returned. Overall, our results indicate that consideration of the thermal evolution of coronal loop systems can allow us to better describe oscillations in these structures and return more accurate estimates of the physical properties of the loops (e.g., density, scale height, magnetic field strength)., Comment: 7 Figures, 14 pages, Accepted for publication in Frontiers in Astronomy and Space Science

We present results derived from the analysis of spectropolarimetric measurements of active region AR12546, which represents one of the largest sunspots to have emerged onto the solar surface over the last $20$ years. The region was observed with full-Stokes scans of the Fe I 617.3 nm and Ca II 854.2 nm lines with the Interferometric BIdimensional Spectrometer (IBIS) instrument at the Dunn Solar Telescope over an uncommon, extremely long time interval exceeding three hours. Clear circular polarization (CP) oscillations localized at the umbra-penumbra boundary of the observed region were detected. Furthermore, the multi-height data allowed us to detect the downward propagation of both CP and intensity disturbances at $2.5-3$~mHz, which was identified by a phase delay between these two quantities. These results are interpreted as a propagating magneto-hydrodynamic surface mode in the observed sunspot., Comment: 8 Pages, 4 Figures, ApJ Accepted

We investigate the characteristics of magneto-acoustic surface waves propagating at a single density interface, in the presence of an inclined magnetic field. For linear wave propagation, dispersion relation is obtained and analytical solutions are derived for small inclination angle. The inclination of the field renders the frequency of waves to be complex, where the imaginary part describes wave attenuation, due to lateral energy leakage., Comment: 19 pages, 5 figures

Torsional motions are ubiquitous in the solar atmosphere. In this work, we perform 3D numerical simulations which mimic a vortex-type photospheric driver with a Gaussian spatial profile. This driver is implemented to excite MHD waves in an axially symmetric, 3D magnetic flux tube embedded in a realistic solar atmosphere. The Gaussian width of the driver is varied and the resulting perturbations are compared. Velocity vectors were decomposed into parallel, perpendicular and azimuthal components with respect to pre-defined magnetic flux surfaces. These components correspond broadly to the fast, slow and Alfven modes, respectively. From these velocities the corresponding wave energy fluxes are calculated, allowing us to estimate the contribution of each mode to the energy flux. For the narrowest driver ($0.15$ Mm) the parallel component accounts for $\sim 55 - 65\%$ of the flux. This contribution increases smoothly with driver width up to nearly $90\%$ for the widest driver ($0.35$ Mm). The relative importance of the perpendicular and azimuthal components decrease at similar rates. The azimuthal energy flux varied between $\sim 35\%$ for the narrowest driver and $< 10\%$ for the widest one. Similarly, the perpendicular flux was $\sim 25 - 10\%$. We also demonstrate that the fast mode corresponds to the sausage wave in our simulations. Our results therefore show that the fast sausage wave is easily excited by this driver and that it carries the majority of the energy transported. For this vortex-type driver the Alfven wave does not contribute a significant amount of energy.

The temporal recurrence of micro-flare events is studied in a time interval before and after of major solar flares. Our sample is based on the x-ray flare observations by the Geostationary Operational Environmental Satellite (GOES) and Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The analysed data contains 1330/301 M- and X-class GOES/RHESSI energetic solar flares and 4062/4119 GOES/RHESSI micro-flares covering the period elapsed since 2002. The temporal analysis of recurrence, by Fast Fourier Transform (FFT), of the micro-flares shows multiple significant periods. Based on the GOES and RHESSI data, the temporal analysis also demonstrates that multiple periods manifest simultaneously in both statistical samples without any significant shift over time. In the GOES sample, the detected significant periods are: $11.33$, $5.61$, $3.75$, $2.80$ and $2.24$ minutes. The RHESSI data shows similar significant periods at $8.54$, $5.28$, $3.66$, $2.88$ and $2.19$ minutes. The periods are interpreted as signatures of standing oscillations, with the longest period ($P_{1}$) being the fundamental and others as higher harmonic modes. The period ratio of the fundamental and higher harmonics ($P_{1}/P_{N}$) is also analysed. The standing modes may be signatures of global oscillations of the entire solar atmosphere encompassing magnetised plasma from photosphere to corona in active regions.

Here, we test the weighted horizontal magnetic gradient ($WG_M$) as a flare precursor, introduced by Korsos et. al. (2015, ApJ, 802, L21), by applying it to a magneto-hydrodynamic (MHD) simulation of solar-like flares Chatterjee et. al. (2016, Physical Review Letters, 116, 10, 101101). The pre-flare evolution of the $WG_M$ and the behavior of the distance parameter between the area-weighted barycenters of opposite polarity sunspots at various heights is investigated in the simulated $\delta$-type sunspot. Four flares emanated from this sunspot. We found the optimum heights above the photosphere where the flare precursors of the $WG_M$ method are identifiable prior to each flare. These optimum heights agree reasonably well with the heights of the occurrence of flares identified from the analysis of their thermal and Ohmic heating signatures in the simulation. We also estimated the expected time of the flare onsets from the duration of the approaching-receding motion of the barycenters of opposite polarities before each single flare. The estimated onset time and the actual time of occurrence of each flare are in good agreement at the corresponding optimum heights. This numerical experiment further supports the use of flare precursors based on the $WG_M$ method., Comment: 12 figures, 2 tables

This study aims to provide further evidence for the potential influence of the global solar magnetic field on localised chromospheric jets, the macrospicules (MS). To find a connection between the long-term variation of properties of MS and other solar activity proxies, including e.g. the temporal variation of the frequency shift of solar global oscillations, sunspot area, etc., a database overarching seven years of observations was built up. This database contains 362 MS, based on observations at the 30.4 nm of the Atmospheric Imaging Assembly (AIA) on-board the Solar Dynamics Observatory (SDO). Three of the five investigated physical properties of MS show a clear long-term temporal variation after smoothing the raw data. Wavelet analysis of the temporal variation of maximum length, maximum area and average velocity is carried out. The results reveal a strong pattern of periodicities at around 2-year (also referred to as Quasi-Biennial Oscillations -- QBOs). Comparison to solar activity proxies, that also possess the properties of QBOs, provides some interesting features: the minima and maxima of QBOs of MS properties occur at around the same epoch as the minima and maxima of these activity proxies. For most of the time span investigated, the oscillations are out-of-phase. This out-of-phase behaviour was also corroborated by a cross-correlation analysis. These results suggest that the physical processes, that generate and drive the long-term evolution of the global solar activity proxies, may be coupled to the short-term local physical processes driving the macrospicules, and, therefore modulate the properties of local dynamics.

Coronal loops exist ubiquitously in the solar atmosphere. These loops puzzle astronomers over half a century. Solar magneto-seismology (SMS) provides a unique way to constrain the physical parameters of coronal loops. Here, we study the evolution of oscillations of a coronal loop observed by the Atmospheric Imaging Assembly (AIA). We measure geometric and physical parameters of the loop oscillations. In particular, we find that the mean period of the oscillations increased from 1048 to 1264 s during three oscillatory cycles. We employ the differential emission measure method and apply the tools of SMS. The evolution of densities inside and outside the loop is analyzed. We found that an increase of density inside the loop and decrease of the magnetic field strength along the loop are the main reasons for the increase in the period during the oscillations. Besides, we also found that the amplitude profile of the loop is different from a profile would it be a homogeneous loop. It is proposed that the distribution of magnetic strength along the loop rather than density stratification is responsible for this deviation. The variation in period and distribution of amplitude provide, in terms of SMS, a new and unprecedented insight into coronal loop diagnostics.

Vortex motions are frequently observed on the solar photosphere. These motions may play a key role in the transport of energy and momentum from the lower atmosphere into the upper solar atmosphere, contributing to coronal heating. The lower solar atmosphere also consists of complex networks of flux tubes that expand and merge throughout the chromosphere and upper atmosphere. We perform numerical simulations to investigate the behaviour of vortex driven waves propagating in a pair of such flux tubes in a non-force-free equilibrium with a realistically modelled solar atmosphere. The two flux tubes are independently perturbed at their footpoints by counter-rotating vortex motions. When the flux tubes merge, the vortex motions interact both linearly and nonlinearly. The linear interactions generate many small-scale transient magnetic substructures due to the magnetic stress imposed by the vortex motions. Thus, an initially monolithic tube is separated into a complex multi-threaded tube due to the photospheric vortex motions. The wave interactions also drive a superposition that increases in amplitude until it exceeds the local Mach number and produces shocks that propagate upwards with speeds of approximately $ 50$ km s$^{-1}$. The shocks act as conduits transporting momentum and energy upwards, and heating the local plasma by more than an order of magnitude, with peak temperature approximately $60,000$ K. Therefore, we present a new mechanism for the generation of magnetic waveguides from the lower solar atmosphere to the solar corona. This wave guide appears as the result of interacting perturbations in neighbouring flux tubes. Thus, the interactions of photospheric vortex motions is a potentially significant mechanism for energy transfer from the lower to upper solar atmosphere.

In this paper, we investigate the inhomogeneous spatial distribution of solar faculae. The focus is on the latitudinal and longitudinal distributions of these highly localised features covering ubiquitously the solar surface. The statistical analysis is based on white light observations of the Solar and Heliospheric Observatory (SOHO) and Solar Dynamics Observatory (SDO) between 1996 and 2014. We found that the fine structure of the latitudinal distribution of faculae displays a quasi-biennial oscillatory pattern. Furthermore, the longitudinal distribution of photospheric solar faculae does not show homogeneous behaviour either. In particular, the non-axisymmetric behaviour of these events show similar properties as that of the active longitude (AL) found in the distribution of sunspots. Our results, preliminary though, may provide a valuable observational constrain for developing the next-generation solar dynamo model.

We present a combined analysis of the applications of the weighted horizontal magnetic gradient (denoted as WG_M in Korsos et al., ApJ, 802, L21, 2015) method and the magnetic helicity tool (Berger & Field, JFM, 147, 133, 1984) employed for three active regions (ARs), namely NOAA AR11261, AR11283 and AR11429. All three active regions produced series flares and CMEs. We followed the evolution of the components of the WG_M and the magnetic helicity before the flare and CME occurrences. We found an unique and mutually shared behavior, called the U-shaped pattern, of the weighted distance component of WG_M and of the shearing component of the helicity flux before the flare and CME eruptions. This common pattern is associated with the decreasing-receding phase yet reported only known to be a necessary feature prior to solar flare eruption(s), but found now at the same time in the evolution of the shearing helicity parameter. This result leads to the conclusion that (i) the shearing motion of photospheric magnetic field may be a key driver for the solar eruption in addition to the flux emerging process, and that (ii) the found decreasing-approaching pattern in the evolution of shearing helicity may be another precursor indicator for improving the forecasting of solar eruptions.

We apply a novel pre-flare tracking of sunspot groups towards improving the estimation of flare onset time by focusing on the evolution of the 3D magnetic field construction of AR 11429. The 3D magnetic structure is based on potential field extrapolation encompassing a vertical range from the photosphere through the chromosphere and transition region into the low corona. The basis of our proxy measure of activity prediction is the so-called weighted horizontal gradient of magnetic field (WG_M) defined between spots of opposite polarities close to the polarity inversion line of an active region. The temporal variation of the distance of the barycenter of the opposite polarities is also found to possess potentially important diagnostic information about the flare onset time estimation as function of height similar to its counterpart introduced initially in an application at the photosphere only in Korsos et al. (2015). We apply the photospheric pre-flare behavioural patterns of sunspot groups to the evolution of their associated 3D-constructed AR 11429 as function of height. We found that at a certain height in the lower solar atmosphere the onset time may be estimated much earlier than at the photosphere or at any other heights. Therefore, we present a tool and recipe that may potentially identify the optimum height for flare prognostic in the solar atmosphere allowing to improve our flare prediction capability and capacity.

We propose to apply the weighted horizontal magnetic gradient ($WG_{M}$), introduced in Korsos et al (2015), for analysing the pre-flare and pre-CME behaviour and evolution of Active Regions (ARs) using the SDO/HMI-Debrecen Data catalogue. To demonstrate the power of investigative capabilities of the $WG_M$ method, in terms of flare and CME eruptions, we studied two typical ARs, namely, AR 12158 and AR 12192. The choice of ARs represent canonical cases. AR 12158 produced an X1.6 flare with fast "halo" CME ($v_{linear}$=1267 $km s^{-1}$) while in AR 12192 there occurred a range of powerful X-class eruptions, i.e. X1.1, X1.6, X3.1, X1.0, X2.0 and X2.0-class energetic flares, interestingly, none with an accompanying CME. The value itself and temporal variation of $WG_{M}$ is found to possess potentially important diagnostic information about the intensity of the expected flare class. Furthermore, we have also estimated the flare onset time from the relationship of duration of converging and diverging motions of the area-weighted barycenters of two subgroups of opposite magnetic polarities. This test turns out not only to provide information about the intensity of the expected flare-class and the flare onset time but may also indicate whether a flare will occur with/without fast CME. We have also found that, in the case when the negative polarity barycenter has moved around and the positive one "remained" at the same coordinates preceding eruption, the flare occurred with fast "halo" CME. Otherwise, when both the negative and the positive polarity barycenters have moved around, the AR produced flares without CME. If these properties found for the movement of the barycenters are generic pre-cursors of CME eruption (or lack of it), identifying them may serve as an excellent pre-condition for refining the forecast of the lift-off of CMEs.

This present study deals with the dissipative instability that appears in a compressible partially ionised plasma slab embedded in a uniform magnetic field, modelling the state of the plasma in solar prominences. In the partially ionised plasma, the dominant dissipative effect is the Cowling resistivity. The regions outside the slab (modelling the solar corona) are fully ionised, and the dominant mechanism of dissipation is viscosity. Analytical solutions to the extended magnetohydrodynamic (MHD) equations are found inside and outside of the slab and solutions are matched at the boundaries of the slab. The dispersion relation is derived and solutions are found analytically in the slender slab limit, while the conditions necessary for the appearance of the instability is investigated numerically for the entire parameter space. Our study is focussed on the effect of the compressibility on the generation and evolution of instabilities. We find that compressibility reduces the threshold of the equilibrium flow, where waves can be unstable, to a level that is comparable to the internal cusp speed, which is of the same order of flow speeds that are currently observed in solar prominences. Our study addresses only the slow waves, as these are the most likely perturbations to become unstable, however the time-scales of the instability are found to be rather large ranging from $10^5$-$10^7$ seconds. It is determined that the instability threshold is further influenced by the concentration of neutrals and the strength of the viscosity of the corona. Interestingly, these two latter aspects have opposite effects. Our numerical analysis shows that the interplay between the equilibrium flow, neutrals and dispersion can change considerably the nature of waves.

The longitudinal distribution of solar active regions shows non-homogeneous spatial behaviour, which is often referred to as Active Longitude (AL). Evidence for a significant statistical relationships between the AL and the longitudinal distribution of flare and coronal mass ejections (CME) occurrences is found in Gyenge et al, 2017 (ApJ, 838, 18). The present work forecasts the spatial position of AL, hence the most flare/CME capable active regions are also predictable. Our forecast method applies Autoregressive Integrated Moving Average model for the next 2 years time period. We estimated the dates when the solar flare/CME capable longitudinal belts face towards Earth.

This paper introduces the Sheffield Magnetohydrodynamics Algorithm Using GPUs (SMAUG+), an advanced numerical code for solving magnetohydrodynamic (MHD) problems, using multi-GPU systems. Multi-GPU systems facilitate the development of accelerated codes and enable us to investigate larger model sizes and/or more detailed computational domain resolutions. This is a significant advancement over the parent single-GPU MHD code, SMAUG (Griffiths, M., Fedun, V., and Erd\'elyi, R. (2015). A fast MHD code for gravitationally stratified media using graphical processing units: SMAUG. Journal of Astrophysics and Astronomy,36(1):197-223). Here, we demonstrate the validity of the SMAUG+ code, describe the parallelisation techniques and investigate performance benchmarks. The initial configuration of the Orszag-Tang vortex simulations are distributed among 4, 16, 64 and 100 GPUs. Furthermore, different simulation box resolutions are applied: $1000 \times 1000$, $2044 \times 2044$, $4000 \times 4000$ and $8000 \times 8000$. We also tested the code with the Brio-Wu shock tube simulations with model size of 800 employing up to 10 GPUs. Based on the test results, we observed speed ups and slow downs, depending on the granularity and the communication overhead of certain parallel tasks. The main aim of the code development is to provide massively parallel code without the memory limitation of a single GPU. By using our code, the applied model size could be significantly increased. We demonstrate that we are able to successfully compute numerically valid and large 2D MHD problems.

High spatial and temporal resolution images of a sunspot, obtained simultaneously in multiple optical and UV wavelengths, are employed to study the propagation and damping characteristics of slow magnetoacoustic waves up to transition region heights. Power spectra are generated from intensity oscillations in sunspot umbra, across multiple atmospheric heights, for frequencies up to a few hundred mHz. It is observed that the power spectra display a power-law dependence over the entire frequency range, with a significant enhancement around 5.5 mHz found for the chromospheric channels. The phase-difference spectra reveal a cutoff frequency near 3 mHz, up to which the oscillations are evanescent, while those with higher frequencies propagate upwards. The power-law index appears to increase with atmospheric height. Also, shorter damping lengths are observed for oscillations with higher frequencies suggesting frequency-dependent damping. Using the relative amplitudes of the 5.5 mHz (3 minute) oscillations, we estimate the energy flux at different heights, which seems to decay gradually from the photosphere, in agreement with recent numerical simulations. Furthermore, a comparison of power spectra across the umbral radius highlights an enhancement of high-frequency waves near the umbral center, which does not seem to be related to magnetic field inclination angle effects., Comment: Accepted for publication in ApJ

Ellerman bombs (EBs) have been widely studied over the past two decades; however, only recently have counterparts of these events been observed in the quiet-Sun. The aim of this article is to further understand small-scale quiet-Sun Ellerman-like brightenings (QSEBs) through research into their spectral signatures, including investigating whether the hot signatures associated with some EBs are also visible co-spatial to any QSEBs. We combine H$\alpha$ and Ca II $8542$ \AA\ line scans at the solar limb with spectral and imaging data sampled by the Interface Region Imaging Spectrograph (IRIS). Twenty one QSEBs were identified with average lifetimes, lengths, and widths measured to be around $120$ s, $0.63$", and $0.35$", respectively. Three of these QSEBs displayed clear repetitive flaring through their lifetimes, comparable to the behaviour of EBs in Active Regions (ARs). Two QSEBs in this sample occurred co-spatial with increased emission in SDO/AIA $1600$ \AA\ and IRIS slit-jaw imager $1400$ \AA\ data, however, these intensity increases were smaller than reported co-spatial to EBs. One QSEB was also sampled by the IRIS slit during its lifetime, displaying increases in intensity in the Si IV $1393$ \AA\ and Si IV $1403$ \AA\ cores as well as the C II and Mg II line wings, analogous to IRIS bursts (IBs). Using RADYN simulations, we are unable to reproduce the observed QSEB H$\alpha$ and Ca II $8542$ \AA\ line profiles leaving the question of the temperature stratification of QSEBs open. Our results imply that some QSEBs could be heated to Transition Region temperatures, suggesting that IB profiles should be observed throughout the quiet-Sun., Comment: 12 pages, 10 figures, accepted by ApJ

Jets, whatever small (e.g. spicules) or large (e.g. macrospicules) their size, may play a key role in momentum and energy transport from photosphere to chromosphere and at least to the low corona. Here, we investigate the properties of abundant, large-scale dynamic jets observable in the solar atmosphere: the macrospicules (MS). These jets are observationally more distinct phenomena than their little, and perhaps more ubiquitous, cousins, the spicules. Investigation of long-term variation of the properties of macrospicules may help to a better understanding of their underlying physics of generation and role in coronal heating. Taking advantage of the high temporal and spatial resolution of the Solar Dynamics Observatory, a new dataset, with several hundreds of macrospicules, was constructed encompassing a period of observations over six years. Here, we analyse the measured properties and relations between these properties of macrospicules as function of time during the observed time interval. We found that cross-correlations of several of these macrospicule properties display a strong oscillatory pattern. Next, wavelet analysis is used to provide more detailed information about the temporal behaviour of the various properties of MS. For coronal hole macrospicules, a significant peak is found at around 2-year period. This peak also exists partially or is shifted to longer period, in the case of quiet Sun macrospicules. These observed findings may be rooted in the underlying mechanism generating the solar magnetic field, i.e. the global solar dynamo.

The spatial inhomogeneity of the distribution of coronal mass ejection (CME) occurrences in the solar atmosphere could provide a tool to estimate the longitudinal position of the most probable CME-capable active regions in the Sun. The anomaly in the longitudinal distribution of active regions themselves is often referred to as active longitude (AL). In order to reveal the connection between the AL and CME spatial occurrences, here we investigate the morphological properties of active regions. The first morphological property studied is the separateness parameter, which is able to characterize the probability of the occurrence of an energetic event, such as a solar flare or CME. The second morphological property is the sunspot tilt angle. The tilt angle of sunspot groups allows us to estimate the helicity of active regions. The increased helicity leads to a more complex buildup of the magnetic structure and also can cause CME eruption. We found that the most complex active regions appear near the AL and that the AL itself is associated with the most tilted active regions. Therefore, the number of CME occurrences is higher within the AL. The origin of the fast CMEs is also found to be associated with this region. We concluded that the source of the most probably CME-capable active regions is at the AL. By applying this method, we can potentially forecast a flare and/or CME source several Carrington rotations in advance. This finding also provides new information for solar dynamo modeling.

Macrospicules (MS) are localised small-scale jet-like phenomena in the solar atmosphere, which have the potential to transport considerable amount of momentum and energy from the lower solar atmospheric regions to the Transition Region and the low corona. A detailed statistical analysis of their temporal behaviour and spatial properties is carried out in this work. By means of state-of-the-art spatial and temporal resolution observations, yielded by the Atmospheric Imaging Assembly (AIA) of Solar Dynamics Observatory (SDO), we constructed a database covering a 5.5-year long period, containing 301 macrospicules that occurred between June 2010 and December 2015 detected at 30.4 nm wavelength. Here, we report the long-term variation of the height, length, average speed and width of MS in Coronal Holes and Quiet Sun areas both in the northern and southern hemisphere of the Sun. This new database helps to refine our knowledge about the physical properties of MS. Cross-correlation of these properties show a relatively strong correlation, but not always a dominant one. However, a more detailed analysis indicates a wave-like signature in the behaviour of MS properties in time. The period of these long-term oscillatory behaviours are just under two years. Also, in terms of solar north/south hemispheres, a strong asymmetry was found in the spatial distribution of MS properties, which may be accounted for the solar dynamo. This latter feature may then indicate a strong and rather intrinsic link between global internal and local atmospheric phenomena in the Sun.

Burst-like events with signatures in the UV are often observed co-spatial to strong line-of-sight photospheric magnetic fields. Several authors, for example, have noted the spatial relationship between Ellerman bombs (EBs) and Moving Magnetic Features (MMFs), regions of flux which disconnect from a sunspot or pore before propagating away in the moat flow and often displaying evidence of cancellation. In this article, data collected by the Solar Dynamics Observatory's Helioseismic and Magnetic Imager and Atmospheric Imaging Assembly are analysed in an attempt to understand the potential links between such cancellation and UV burst formation. Two MMFs from AR 11579, three bi-poles from AR 11765, and six bi-poles (four of which were co-spatial to IRIS bursts) in AR 11850 were identified for analysis. All of these cancellation features were found to have lifetimes of the order hours and cancellation rates of the order 10^14-10^15 Mx s^-1. H-alpha line wing data from the Dunn Solar Telescope's Interferometric BIdimensional Spectrometer were also available for AR 11579 facilitating a discussion of links between MMFs and EBs. Using an algebraic model of photospheric magnetic reconnection, the measured cancellation rates are then used to ascertain estimates of certain quantities (such as up-flow speeds, jet extents, and potential energy releases) which compared reasonably to the properties of EBs reported within the literature. Our results suggest that cancellation rates of the order measured here are capable of supplying enough energy to drive certain UV bursts (including EBs), however, they are not a guaranteeing condition for burst formation., Comment: 13 pages, 6 figures, 2 tables in Monthly Notices of the Royal Astronomical Society 2016

In this paper the longitudinal and latitudinal spatial distribution of macrospicules is examined. We found a statistical relationship between the active longitude determined by sunspot groups and the longitudinal distribution of macrospicules. This distribution of macrospicules shows an inhomogeneity and non-axysimmetrical behaviour in the time interval from June $2010$ until December $2012$ covered by observations of the Solar Dynamic Observatory (SDO) satellite. The enhanced positions of the activity and its time variation has been calculated. The migration of the longitudinal distribution of macrospicules shows a similar behaviour as that of the sunspot groups.

We present observational evidence of compressible magnetohydrodynamic wave modes propagating from the solar photosphere through to the base of the transition region in a solar magnetic pore. High cadence images were obtained simultaneously across four wavelength bands using the Dunn Solar Telescope. Employing Fourier and wavelet techniques, sausage-mode oscillations displaying significant power were detected in both intensity and area fluctuations. The intensity and area fluctuations exhibit a range of periods from 181-412s, with an average period ~290s, consistent with the global p-mode spectrum. Intensity and area oscillations present in adjacent bandpasses were found to be out-of-phase with one another, displaying phase angles of 6.12 degrees, 5.82 degrees and 15.97 degrees between 4170 Angstrom continuum - G-band, G-band - Na I D1 and Na I D1 - Ca II K heights, respectively, reiterating the presence of upwardly-propagating sausage-mode waves. A phase relationship of ~0 degrees between same-bandpass emission and area perturbations of the pore best categorises the waves as belonging to the `slow' regime of a dispersion diagram. Theoretical calculations reveal that the waves are surface modes, with initial photospheric energies in excess of 35000 W/m^2. The wave energetics indicate a substantial decrease in energy with atmospheric height, confirming that magnetic pores are able to transport waves that exhibit appreciable energy damping, which may release considerable energy into the local chromospheric plasma., Comment: 14 pages, 5 figures, accepted into ApJ

We analyse high temporal and spatial resolution time-series of spectral scans of the Halpha line obtained with the CRisp Imaging SpectroPolarimeter (CRISP) instrument mounted on the Swedish Solar Telescope. The data reveal highly dynamic, dark, short-lived structures known as Rapid Redshifted and Blueshifted Excursions (RREs, RBEs) that are on-disk absorption features observed in the red and blue wings of spectral lines formed in the chromosphere. We study the dynamics of RREs and RBEs by tracking their evolution in space and time, measuring the speed of the apparent motion, line-of-sight Doppler velocity, and transverse velocity of individual structures. A statistical study of their measured properties shows that RREs and RBEs have similar occurrence rates, lifetimes, lengths, and widths. They also display non-periodic, non-linear transverse motions perpendicular to their axes at speeds of 4 - 31 km/s. Furthermore, both types of structures either appear as high speed jets and blobs that are directed outwardly from a magnetic bright point with speeds of 50 - 150 km/s, or emerge within a few seconds. A study of the different velocity components suggests that the transverse motions along the line-of-sight of the chromospheric flux tubes are responsible for the formation and appearance of these redshifted/blueshifted structures. The short lifetime and fast disappearance of the RREs/RBEs suggests that, similar to type II spicules, they are rapidly heated to transition region or even coronal temperatures. We speculate that the Kelvin-Helmholtz instability triggered by observed transverse motions of these structures may be a viable mechanism for their heating., Comment: 8 pages, 6 figures, accepted in ApJ

High-resolution observations of the solar photosphere have identified a wide variety of spiralling motions in the solar plasma. These spirals vary in properties, but are observed to be abundant at the solar surface. In this work these spirals are studied for their potential as magnetohydrodynamic (MHD) wave generation mechanisms. The inter-granular lanes, where these spirals are commonly observed, are also regions where the magnetic field strength is higher than average. This combination of magnetic field and spiralling plasma is a recipe for the generation of Alfv\'en waves and other MHD waves. This work employs numerical simulations of a self-similar magnetic flux tube embedded in a realistic, gravitationally stratified, solar atmosphere to study the effects of a single magnetic flux tube perturbed by a logarithmic velocity spiral driver. The expansion factor of the logarithmic spiral driver is varied and multiple simulations are run for a range of values of the expansion factor centred around observational constraints. The simulations are analysed using `flux surfaces' constructed from the magnetic field lines so that the vectors perpendicular, parallel and azimuthal to the local magnetic field vector can be calculated. The results of this analysis show that the Alfv\'en wave is the dominant wave for lower values of the expansion factor, whereas, for the higher values the parallel component is dominant. This transition occurs within the range of the observational constraints, meaning that spiral drivers, as observed in the solar photosphere, have the potential to generate a variety of MHD wave modes.

We report the observation of the first two harmonics of the horizontally polarized kink waves excited in a coronal loop system lying at south-east of AR 11719 on 2013 April 11. The detected periods of the fundamental mode ($P_1$), its first overtone ($P_2$) in the northern half, and that in the southern one are $530.2 \pm 13.3$, $300.4 \pm 27.7$, and $334.7 \pm 22.1$ s, respectively. The periods of the first overtone in the two halves are the same considering uncertainties in the measurement. We estimate the average electron density, temperature, and length of the loop system as $(5.1 \pm 0.8) \times 10^8$ cm$^{-3}$, $0.65 \pm 0.06$ MK, and $203.8 \pm 13.8$ Mm, respectively. As a zeroth order estimation, the magnetic field strength, $B = 8.2 \pm 1.0$ G, derived by the coronal seismology using the fundamental kink mode matches with that derived by a potential field model. The extrapolation model also shows the asymmetric and nonuniform distribution of the magnetic field along the coronal loop. Using the amplitude profile distributions of both the fundamental mode and its first overtone, we observe that the antinode positions of both the fundamental mode and its first overtone shift towards the weak field region along the coronal loop. The results indicate that the density stratification and the temperature difference effects are larger than the magnetic field variation effect on the period ratio. On the other hand, the magnetic field variation has a greater effect on the eigenfunction of the first overtone than the density stratification does for this case., Comment: 24 pages, 6 figures, 1 table, accepted for publication in ApJ

Ellerman bombs (EBs) have been widely studied in recent years due to their dynamic, explosive nature and apparent links to the underlying photospheric magnetic field implying that they may be formed by magnetic reconnection in the photosphere. Despite a plethora of researches discussing the morphologies of EBs, there has been a limited investigation of how these events appear at the limb, specifically, whether they manifest as vertical extensions away from the disc. In this article, we make use of high-resolution, high-cadence observations of an AR at the solar limb, collected by the CRISP instrument, to identify EBs and infer their physical properties. The upper atmosphere is also probed using the SDO/AIA. We analyse 22 EB events evident within these data, finding that 20 appear to follow a parabolic path away from the solar surface at an average speed of 9 km s^(-1), extending away from their source by 580 km, before retreating back at a similar speed. These results show strong evidence of vertical motions associated with EBs, possibly explaining the dynamical `flaring' (changing in area and intensity) observed in on-disc events. Two in-depth case studies are also presented which highlight the unique dynamical nature of EBs within the lower solar atmosphere. The viewing angle of these observations allows for a direct linkage between these EBs and other small-scale events in the H-alpha line wings, including a potential flux emergence scenario. The findings presented here suggest that EBs could have a wider-reaching influence on the solar atmosphere than previously thought, as we reveal a direct linkage between EBs and an emerging small-scale loop, and other near-by small-scale explosive events. However, as previous research found, these extensions do not appear to impact upon the H-alpha line core, and are not observed by the SDO/AIA EUV filters., Comment: 10 pages, 7 figures, accepted to ApJ

There have been ubiquitous observations of wave-like motions in the solar atmosphere for decades. Recent improvements to space- and ground-based observatories have allowed the focus to shift to smaller magnetic structures on the solar surface. In this paper, high-resolution ground-based data taken using the Swedish 1 m Solar Telescope is combined with co-spatial and co-temporal data from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) satellite to analyze running penumbral waves (RPWs). RPWs have always been thought to be radial wave propagation that occurs within sunspots. Recent research has suggested that they are in fact upwardly propagating field-aligned waves (UPWs). Here, RPWs within a solar pore are observed for the first time and are interpreted as UPWs due to the lack of a penumbra that is required to support RPWs. These UPWs are also observed co-spatially and co-temporally within several SDO/AIA elemental lines that sample the transition region and low corona. The observed UPWs are traveling at a horizontal velocity of around 17 +- 0.5 km s-1 and a minimum vertical velocity of 42 +- 21 km s-1. The estimated energy of the waves is around 150 W m-2, which is on the lower bound required to heat the quiet-Sun corona. This is a new, yet unconsidered source of wave energy within the solar chromosphere and low corona., Comment: 7 pages, 4 figures

In this paper, we consider the complex modified Korteweg-de Vries (mKdV) equation as a model of few-cycle optical pulses. Using the Lax pair, we construct a generalized Darboux transformation and systematically generate the first-, second- and third-order rogue wave solutions and analyze the nature of evolution of higher-order rogue waves in detail. Based on detailed numerical and analytical investigations, we classify the higher-order rogue waves with respect to their intrinsic structure, namely, fundamental pattern, triangular pattern, and ring pattern. We also present several new patterns of the rogue wave according to the standard and non-standard decomposition. The results of this paper explain the generalization of higher-order rogue waves in terms of rational solutions. We apply the contour line method to obtain the analytical formulas of the length and width of the first-order RW of the complex mKdV and the NLS equations. In nonlinear optics, the higher-order rogue wave solutions found here will be very useful to generate high-power few-cycle optical pulses which will be applicable in the area of ultra-short pulse technology., Comment: 31 pages, 22 figures, accepted by Phys.Rev.E

The launch of the Hinode satellite has allowed for seeing-free observations at high-resolution and high-cadence making it well suited to study the dynamics of quiescent prominences. In recent years it has become clear that quiescent prominences support small-amplitude transverse oscillations, however, sample sizes are usually too small for general conclusions to be drawn. We remedy this by providing a statistical study of transverse oscillations in vertical prominence threads. Over a three-hour period of observations it was possible to measure the properties of 3436 waves, finding periods from 50 to 6000 s with typical velocity amplitudes ranging between 0.2 to 23 km s$^{-1}$. The large number of observed waves allows the determination of the frequency dependence of the wave properties and derivation of the velocity power spectrum for the transverse waves. For frequencies less than 7 mHz, the frequency-dependence of the velocity power is consistent with the velocity power spectra generated from observations of the horizontal motions of magnetic elements in the photosphere, suggesting that the prominence transverse waves are driven by photospheric motions. However, at higher frequencies the two distributions significantly diverge, with relatively more power found at higher frequencies in the prominence oscillations. These results highlight that waves over a large frequency range are ubiquitous in prominences, and that a significant amount of the wave energy is found at higher frequency., Comment: 6 pages, 5 figures, accepted for publication in ApJ letters

The presence of photospheric magnetic reconnection has long been thought to give rise to short and impulsive events, such as Ellerman bombs (EBs) and Type II spicules. In this article, we combine high-resolution, high-cadence observations from the Interferometric BIdimensional Spectrometer (IBIS) and Rapid Oscillations in the Solar Atmosphere (ROSA) instruments at the Dunn Solar Telescope, National Solar Observatory, New Mexico with co-aligned Atmospheric Imaging Assembly (SDO/AIA) and Solar Optical Telescope (Hinode/SOT) data to observe small-scale events situated within an active region. These data are then compared with state-of-the-art numerical simulations of the lower atmosphere made using the MURaM code. It is found that brightenings, in both the observations and the simulations, of the wings of the H alpha line profile, interpreted as EBs, are often spatially correlated with increases in the intensity of the FeI 6302.5A line core. Bi-polar regions inferred from Hinode/SOT magnetic field data show evidence of flux cancellation associated, co-spatially, with these EBs, suggesting magnetic reconnection could be a driver of these high-energy events. Through the analysis of similar events in the simulated lower atmosphere, we are able to infer that line profiles analogous to the observations occur co-spatially with regions of strong opposite polarity magnetic flux. These observed events and their simulated counterparts are interpreted as evidence of photospheric magnetic reconnection at scales observable using current observational instrumentation., Comment: 10 pages, 7 figures, accepted to ApJ

The source of the non-thermal energy required for the heating of the upper solar atmosphere to temperatures in excess of a million degrees and the acceleration of the solar wind to hundreds of kilometres per second is still unclear. One such mechanism for providing the required energy flux is incompressible torsional Alfv\'en and kink magnetohydrodynamic (MHD) waves, which are magnetically dominated waves supported by the Sun's pervasive and complex magnetic field. In particular, propagating MHD kink waves have recently been observed to be ubiquitous throughout the solar atmosphere, but, until now, critical details of the transport of the kink wave energy throughout the Sun's atmosphere were unclear. Here, the ubiquity of the waves is exploited for statistical studies in the highly dynamic solar chromosphere. This large-scale investigation allows for the determination of the chromospheric kink wave velocity power spectra, a missing link necessary for determining the energy transport between the photosphere and corona. Crucially, the power spectra contains evidence for horizontal photospheric motions being the main mechanism for kink wave generation in the quiescent Sun. In addition, a comparison to measured coronal power spectra is provided, revealing frequency-dependent transmission profiles suggesting there is enhanced damping of kink waves in the lower corona., Comment: In Press ApJ - March 2014

Using data obtained by the high temporal and spatial resolution Rapid Oscillations in the Solar Atmosphere (ROSA) instrument on the Dunn Solar Telescope, we investigate at an unprecedented level of detail transverse oscillations in chromospheric fine structures near the solar disk center. The oscillations are interpreted in terms of propagating and standing magnetohydrodynamic kink waves. Wave characteristics including the maximum transverse velocity amplitude and the phase speed are measured as a function of distance along the structure's length. Solar magneto-seismology is applied to these measured parameters to obtain diagnostic information on key plasma parameters (e.g., magnetic field, density, temperature, flow speed) of these localised waveguides. The magnetic field strength of the mottle along the $\sim$2 Mm length is found to decrease by a factor of 12, while the local plasma density scale height is $\sim280\pm$80 km., Comment: 8 pages, 7 figures, accepted in ApJ

The details of the mechanism(s) responsible for the observed heating and dynamics of the solar atmosphere still remain a mystery. Magnetohydrodynamic (MHD) waves are thought to play a vital role in this process. Although it has been shown that incompressible waves are ubiquitous in off-limb solar atmospheric observations their energy cannot be readily dissipated. We provide here, for the first time, on-disk observation and identification of concurrent MHD wave modes, both compressible and incompressible, in the solar chromosphere. The observed ubiquity and estimated energy flux associated with the detected MHD waves suggest the chromosphere is a vast reservoir of wave energy with the potential to meet chromospheric and coronal heating requirements. We are also able to propose an upper bound on the flux of the observed wave energy that is able to reach the corona based on observational constraints, which has important implications for the suggested mechanism(s) for quiescent coronal heating., Comment: Published in Nature Communications, 3, 1315 (2012)

The solar atmosphere is a dynamic environment, constantly evolving to form a wide range of magnetically dominated structures (coronal loops, spicules, prominences, etc.) which cover a significant percentage of the surface at any one time. Oscillations and waves in many of these structures are now widely observed and have led to the new analytic technique of solar magneto-seismology, where inferences of the background conditions of the plasma can be deduced by studying magneto-hydrodynamic (MHD) waves. Here, we generalise a novel magneto-seismological method designed to infer the density distribution of a bounded plasma structure from the relationship of its fundamental and first harmonics. Observations of the solar atmosphere have emphatically shown that stratification, leading to complex density profiles within plasma structures, is common thereby rendering this work instantly accessible to solar physics. We show, in a dynamic waveguide, how the period ratio differs from the idealised harmonic ratios prevalent in homogeneous structures. These ratios show strong agreement with recent observational work. Next, anti-node shifts are also analysed. Using typical scaling parameters for bulk flows within atmospheric waveguides, e.g., coronal loops, it is found that significant anti-node shifts can be predicted, even to the order of 10 Mm. It would be highly encouraged to design specific observations to confirm the predicted anti-node shifts and apply the developed theory of solar magneto-seismology to gain more accurate waveguide diagnostics of the solar atmosphere., Comment: 18 pages, 6 figures