Aims. We investigate the formation of small scales and the dissipation of MHD wave energy through non-linear interactions of counter-propagating, phase-mixed Alfvenic waves in a complex magnetic field. Methods. We conducted fully 3-D, non-ideal MHD simulations of transverse waves in complex magnetic fields. Continuous wave drivers were imposed on the foot points of magnetic field lines and the system was evolved for several Alfven travel times. Phase-mixed waves were allowed to reflect off the upper boundary and the interactions between the resultant counter-streaming wave packets were analysed. Results. The complex nature of the background magnetic field encourages the development of phase mixing, leading to a growth in currents and vorticities. Counter-propagating phase-mixed waves induce a cascade of energy to small scales and result in more efficient energy dissipation. This effect is enhanced in simulations with more complex magnetic fields. High-frequency drivers excite localised field line resonances and produce efficient wave heating. However, this relies on the formation of large amplitude oscillations on resonant field lines. Drivers with smaller frequencies than the fundamental frequencies of field lines are not able to excite resonances and thus do not inject sufficient Poynting flux to power coronal heating. Even in the case of high-frequency oscillations, the rate of dissipation is likely too slow to balance coronal energy losses, even within the quiet Sun. Conclusions. For the generalised phase-mixing presented here, complex background field structures enhance the rate of wave energy dissipation. However, it remains difficult for realistic wave drivers to inject sufficient Poynting flux to heat the corona. Indeed, significant heating only occurs in cases which exhibit amplitudes that are much larger than those currently observed in the solar atmosphere.

Modelling the solar Transition Region with the use of an Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, capturing the enthalpy exchange between the corona and transition region, when the corona undergoes impulsive heating. The TRAC method eliminates the need for highly resolved numerical grids in the transition region and the commensurate very short time steps that are required for numerical stability. When employed with coarse spatial resolutions, typically achieved in multi-dimensional magnetohydrodynamic codes, the errors at peak density are less than 5% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. This paper presents further examples that demonstrate the versatility and robustness of the method over a range of heating events, including impulsive and quasi-steady footpoint heating. A detailed analytical assessment of the TRAC method is also presented, showing that the approach works through all phases of an impulsive heating event because (i) the total radiative losses and (ii) the total heating when integrated over the transition region are both preserved at all temperatures under the broadening modifications of the method. The results from the numerical simulations complement this conclusion.

Phase mixing of Alfv\'en waves has been studied extensively as a possible coronal heating mechanism but without the full thermodynamic consequences considered self-consistently. It has been argued that in some cases, the thermodynamic feedback of the heating could substantially affect the transverse density gradient and even inhibit the phase mixing process. In this paper, we use MHD simulations with the appropriate thermodynamical terms included to quantify the evaporation following heating by phase mixing of Alfv\'en waves in a coronal loop and the effect of this evaporation on the transverse density profile. The numerical simulations were performed using the Lare2D code. We set up a 2D loop model consisting of a field-aligned thermodynamic equilibrium and a cross-field (background) heating profile. A continuous, sinusoidal, high-frequency Alfv\'en wave driver was implemented. As the Alfv\'en waves propagate along the field, they undergo phase mixing due to the cross-field density gradient in the coronal part of the loop. We investigated the presence of field-aligned flows, heating from the dissipation of the phase-mixed Alfv\'en waves, and the subsequent evaporation from the lower atmosphere. We find that phase mixing of Alfv\'en waves leads to modest heating in the shell regions of the loop and evaporation of chromospheric material into the corona with upflows of the order of only 5-20 m/s. Although the evaporation leads to a mass increase in the shell regions of the loop, the effect on the density gradient and, hence, on the phase mixing process, is insignificant. This paper self-consistently investigates the effect of chromospheric evaporation on the cross-field density gradient and the phase mixing process in a coronal loop. We found that the effects in our particular setup (small amplitude, high frequency waves) are too small to significantly change the density gradient.

Aims. We investigate the transfer of energy between a standing kink mode and azimuthal Alfven waves within an expanding flux tube. We consider the process of resonant absorption in a loop with a non-uniform Alfven frequency profile but in the absence of a radial density gradient. Methods. Using the MHD code, Lare3d, we model an oscillating magnetic flux tube that expands radially with height. An initially straight loop structure with a magnetic field enhancement is relaxed towards a force-free state before a standing kink mode is introduced. Results. We demonstrate that the transverse gradient in Alfven frequency required for the existence of resonant field lines can be associated with the expansion of a high field-strength flux tube from concentrated flux patches. This allows for the conversion of energy between wave modes even in the absence of the density profile typically assumed in wave heating models. As with standing modes in straight flux tubes, small scales are dominated by the vorticity at the loop apex and by currents at the foot points. The azimuthal wave exhibits the structure of the expanded flux tube and is therefore associated with smaller length scales close to the foot points of the flux tube. Conclusions. Resonant absorption can proceed throughout the corona, even in the absence of dense loop structures. The flux tube and MHD waves considered are difficult to observe and our model highlights how estimating hidden wave power within the solar atmosphere can be problematic. We highlight that, for standing modes, the global properties of field lines are important for resonant absorption and coronal conditions at a single altitude will not fully determine the nature of MHD resonances. We provide a new model in partial response to the criticism that wave heating models cannot self-consistently generate or sustain the density profile upon which they typically rely.

We consider a series of MHD simulations in which a small amplitude, transverse velocity perturbation is introduced into a complex magnetic field. We analysed the deformation of the wave fronts as the perturbation propagates through the braided magnetic structures and explore the nature of Alfv\'enic wave phase mixing in this regime. Spatial gradients in the local Alfv\'en speed and variations in the length of magnetic field lines ensure that small scales form throughout the propagating wave front due to phase mixing. Additionally, the presence of complex, intricate current sheets associated with the background field locally modifies the polarisation of the wave front. The combination of these two effects enhances the rate of viscous dissipation, particularly in more complex field configurations. Unlike in classical phase mixing configurations, the greater spatial extent of Alfv\'en speed gradients ensures that wave energy is deposited over a larger cross-section of the magnetic structure. Further, the complexity of the background magnetic field ensures that small gradients in a wave driver can map to large gradients within the coronal plasma. The phase mixing of MHD waves in a complex magnetic field will progress throughout the braided volume. As a result, in a non-ideal regime wave energy will be dissipated over a greater cross-section than in classical phase mixing models. The formation rate of small spatial scales in a propagating wave front is a function of the complexity of the background magnetic field. As such, if the coronal field is sufficiently complex it remains plausible that phase mixing induced wave heating can contribute to maintaining observed temperatures. Furthermore, the weak compressibility of the transverse wave and the observed phase mixing pattern may provide seismological information about the nature of the background plasma.

MHD waves are ubiquitous in the solar corona where the highly structured magnetic fields provide efficient wave guides for their propagation. While MHD waves have been observed originating from lower layers of the solar atmosphere, recent studies have shown that some can be generated in-situ by the collision of dense counter propagating flows. In this theoretical study, we analyse the mechanism that triggers the propagation of kink and sausage modes in the solar corona, following the collision of counter-propagating flows and how the properties of the flows affect the properties of the generated waves. In order to study in detail such a mechanism we run a series of ideal 2D and 3D MHD simulations where we vary the properties of the counter-propagating flows and by means of a simple technique to estimate the amplitudes of the kink and sausage modes, we investigate their role on the generation and propagation of the MHD waves. We find that the amplitude of the waves is largely dependent on the kinetic energy of the flows and that the onset of kink or sausage modes depends on the asymmetries between the colliding blobs. Moreover, the initial wavelength of the MHD waves is associated with the magnetic configuration resulting from the collision of the flows. We also find that genuine 3D systems respond with smaller wave amplitudes. In this study, we present a parameter space description of the mechanism that leads to the generation of MHD waves from the collision of flows in the corona. Future observations of these waves can be used to understand the properties of the plasma and magnetic field of the solar corona.

Thermal non-equilibrium (TNE) is believed to be a potentially important process in understanding some properties of the magnetically closed solar corona. Through one-dimensional hydrodynamic models, this paper addresses the importance of the numerical spatial resolution, footpoint heating timescales and background heating on TNE. Inadequate transition region (TR) resolution can lead to significant discrepancies in TNE cycle behaviour, with TNE being suppressed in under-resolved loops. A convergence on the periodicity and plasma properties associated with TNE required spatial resolutions of less than 2 km for a loop of length 180 Mm. These numerical problems can be resolved using an approximate method that models the TR as a discontinuity using a jump condition, as proposed by Johnston et al. (2017a,b). The resolution requirements (and so computational cost) are greatly reduced while retaining good agreement with fully resolved results. Using this approximate method we (i) identify different regimes for the response of coronal loops to time-dependent footpoint heating including one where TNE does not arise and (ii) demonstrate that TNE in a loop with footpoint heating is suppressed unless the background heating is sufficiently small. The implications for the generality of TNE are discussed.

Aims: This paper presents 2.5D numerical experiments of Alfv\'en wave phase mixing and aims to assess the effects of nonlinearities on wave behaviour and dissipation. In addition, this paper aims to quantify how effective the model presented in this work is at providing energy to the coronal volume. Methods: The model is presented and explored through the use of several numerical experiments which were carried out using the Lare2D code. The experiments study footpoint driven Alfv\'en waves in the neighbourhood of a two-dimensional x-type null point with initially uniform density and plasma pressure. A continuous sinusoidal driver with a constant frequency is used. Each experiment uses different driver amplitudes to compare weakly nonlinear experiments with linear experiments. Results: We find that the wave trains phase-mix owing to variations in the length of each field line and variations in the field strength. The nonlinearities reduce the amount of energy entering the domain, as they reduce the effectiveness of the driver, but they have relatively little effect on the damping rate (for the range of amplitudes studied). The nonlinearities produce density structures which change the natural frequencies of the field lines and hence cause the resonant locations to move. The shifting of the resonant location causes the Poynting flux associated with the driver to decrease. Reducing the magnetic diffusivity increases the energy build-up on the resonant field lines, however, it has little effect on the total amount of energy entering the system. From an order of magnitude estimate, we show that the Poynting flux in our experiments is comparable to the energy requirements of the quiet Sun corona. However a (possibly unphysically) large amount of magnetic diffusion was used however and it remains unclear if the model is able to provide enough energy under actual coronal conditions.

Whilst there are observational indications that transverse MHD waves carry enough energy to maintain the thermal structure of the solar corona, it is not clear whether such energy can be efficiently and effectively converted into heating. Phase-mixing of Alfv\'en waves is considered a candidate mechanism, as it can develop transverse gradient where magnetic energy can be converted into thermal energy. However, phase-mixing is a process that crucially depends on the amplitude and period of the transverse oscillations, and only recently have we obtained a complete measurement of the power spectrum for transverse oscillations in the corona. We aim to investigate the heating generated by phase-mixing of transverse oscillations triggered by buffeting of a coronal loop that follows from the observed coronal power spectrum as well as the impact of these persistent oscillations on the structure of coronal loops. We consider a 3D MHD model of an active region coronal loop and we perturb its footpoints with a 2D horizontal driver that represents a random buffeting motion of the loop footpoints. Our driver is composed of 1000 pulses superimposed to generate the observed power spectrum. We find that the heating supply from the observed power spectrum in the solar corona through phase mixing is not sufficient to maintain the million degree active region solar corona. We also find that the development of Kelvin-Helmholtz instabilities could be a common phenomenon in coronal loops that could affect their apparent life time. This study concludes that is unlikely that phase-mixing of Alfv\'en waves resulting from an observed power spectrum of transverse coronal loop oscillations can heat the active region solar corona. However, transverse waves could play an important role in the development of small scale structures.

Recent numerical studies revealed that transverse motions of coronal loops can induce the Kelvin-Helmholtz Instability (KHI). This process could be important in coronal heating because it leads to dissipation of energy at small spatial-scale plasma interactions. Meanwhile, small amplitude decayless oscillations in coronal loops have been discovered recently in observations of SDO/AIA. We model such oscillations in coronal loops and study wave heating effects, considering a kink and Alfv{\'e}n driver separately and a mixed driver at the bottom of flux tubes. Both the transverse and Alfv{\'e}n oscillations can lead to the KHI. Meanwhile, the Alfv{\'e}n oscillations established in loops will experience phase mixing. Both processes will generate small spatial-scale structures, which can help the dissipation of wave energy. Indeed, we observe the increase of internal energy and temperature in loop regions. The heating is more pronounced for the simulation containing the mixed kink and Alfv{\'e}n driver. This means that the mixed wave modes can lead to a more efficient energy dissipation in the turbulent state of the plasma and that the KHI eddies act as an agent to dissipate energy in other wave modes. Furthermore, we also obtained forward modelling results using the FoMo code. We obtained forward models which are very similar to the observations of decayless oscillations. Due to the limited resolution of instruments, neither Alfv{\'e}n modes nor the fine structures are observable. Therefore, this numerical study shows that Alfv{\'e}n modes probably can co-exist with kink modes, leading to enhanced heating.

Transverse MHD waves permeate the solar atmosphere and are a candidate for coronal heating. However, the origin of these waves is still unclear. In this work, we analyse coordinated observations from \textit{Hinode}/SOT and \textit{IRIS} of a prominence/coronal rain loop-like structure at the limb of the Sun. Cool and dense downflows and upflows are observed along the structure. A collision between a downward and an upward flow with an estimated energy flux of $10^{7}-10^{8}$~erg~cm$^{-2}$~s$^{-1}$ is observed to generate oscillatory transverse perturbations of the strands with an estimated $\approx40~$km~s$^{-1}$ total amplitude, and a short-lived brightening event with the plasma temperature increasing to at least $10^{5}~$K. We interpret this response as sausage and kink transverse MHD waves based on 2D MHD simulations of plasma flow collision. The lengths, density and velocity differences between the colliding clumps and the strength of the magnetic field are major parameters defining the response to the collision. The presence of asymmetry between the clumps (angle of impact surface and/or offset of flowing axis) is crucial to generate a kink mode. Using the observed values we successfully reproduce the observed transverse perturbations and brightening, and show adiabatic heating to coronal temperatures. The numerical modelling indicates that the plasma $\beta$ in this loop-like structure is confined between $0.09$ and $0.36$. These results suggest that such collisions from counter-streaming flows can be a source of in-situ transverse MHD waves, and that for cool and dense prominence conditions such waves could have significant amplitudes.

The role of type II spicules in the corona has been a much debated topic in recent years. This paper aims to shed light on the impact of type II spicules in the corona using novel 2.5D radiative MHD simulations including ion-neutral interaction effects with the Bifrost code. We find that the formation of simulated type II spicules, driven by the release of magnetic tension, impacts the corona in various manners. Associated with the formation of spicules, the corona exhibits 1) magneto-acoustic shocks and flows which supply mass to coronal loops, and 2) transversal magnetic waves and electric currents that propagate at Alfv\'en speeds. The transversal waves and electric currents, generated by the spicule's driver and lasting for many minutes, are dissipated and heat the associated loop. These complex interactions in the corona can be connected with blue shifted secondary components in coronal spectral lines (Red-Blue asymmetries) observed with Hinode/EIS and SOHO/SUMER, as well as the EUV counterpart of type II spicules and propagating coronal disturbances (PCDs) observed with the 171~\AA\ and 193~\AA\ SDO/AIA channels.

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

In addition to their jet-like dynamic behaviour, spicules usually exhibit strong transverse speeds, multi-stranded structure and heating from chromospheric to transition region temperatures. In this work we first analyse \textit{Hinode} \& \textit{IRIS} observations of spicules and find different behaviours in terms of their Doppler velocity evolution and collective motion of their sub-structure. Some have a Doppler shift sign change that is rather fixed along the spicule axis, and lack coherence in the oscillatory motion of strand-like structure, matching rotation models or long wavelength torsional Alfv\'en waves. Others exhibit a Doppler shift sign change at maximum displacement and coherent motion of their strands, suggesting a collective MHD wave. By comparing with an idealised 3-D MHD simulation combined with radiative transfer modelling, we analyse the role of transverse MHD waves and associated instabilities in spicule-like features. We find that Transverse Wave Induced Kelvin-Helmholtz (TWIKH) rolls lead to coherence of strand-like structure in imaging and spectral maps, as seen in some observations. The rapid transverse dynamics and the density and temperature gradients at the spicule boundary lead to ring-shaped \ion{Mg}{2} k and \ion{Ca}{2} H source functions in the transverse cross-section, potentially allowing IRIS to capture the KHI dynamics. Twists and currents propagate along the spicule at Alfv\'enic speeds, and the temperature variations within TWIKH rolls produce sudden appearance / disappearance of strands seen in Doppler velocity and in \ion{Ca}{2} H intensity. However, only a mild intensity increase in higher temperature lines is obtained, suggesting there is an additional heating mechanism at work in spicules.

Spicules have been proposed as significant contributors to the mass and energy balance of the corona. While previous observations have provided a glimpse of short-lived transient brightenings in the corona that are associated with spicules, these observations have been contested and are the subject of a vigorous debate both on the modeling and the observational side. Therefore, it remains unclear whether plasma is heated to coronal temperatures in association with spicules. We use high-resolution observations of the chromosphere and transition region with the Interface Region Imaging Spectrograph (IRIS) and of the corona with the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) to show evidence of the formation of coronal structures associated with spicular mass ejections and heating of plasma to transition region and coronal temperatures. Our observations suggest that a significant fraction of the highly dynamic loop fan environment associated with plage regions may be the result of the formation of such new coronal strands, a process that previously had been interpreted as the propagation of transient propagating coronal disturbances (PCD)s. Our observations are supported by 2.5D radiative MHD simulations that show heating to coronal temperatures in association with spicules. Our results suggest that heating and strong flows play an important role in maintaining the substructure of loop fans, in addition to the waves that permeate this low coronal environment.

Aims: We quantify the effects of twisted magnetic fields on the development of the magnetic Kelvin-Helmholtz instability (KHI) in transversely oscillating coronal loops. Methods: We modelled a fundamental standing kink mode in a straight, density-enhanced magnetic flux tube using the magnetohydrodynamics code, Lare3d. In order to evaluate the impact of an azimuthal component of the magnetic field, various degrees of twist were included within the flux tube's magnetic field. Results: The process of resonant absorption is only weakly affected by the presence of a twisted magnetic field. However, the subsequent evolution of the KHI is sensitive to the strength of the azimuthal component of the field. Increased twist values inhibit the deformation of the loop's density profile, which is associated with the growth of the instability. Despite this, much smaller scales in the magnetic field are generated when there is a non-zero azimuthal component present. Hence, the instability is more energetic in cases with (even weakly) twisted fields. Field aligned flows at the loop apex are established in a twisted regime once the instability has formed. Conclusions: The KHI may have implications for wave heating in the solar atmosphere due to the creation of small length scales and the generation of a turbulent regime. Whilst magnetic twist does suppress the development of the vortices associated with the instability, the formation of the KHI in a twisted regime will be accompanied by greater Ohmic dissipation due to the larger currents that are produced, even if only weak twist is present. The presence of magnetic twist will likely make the instability more difficult to detect in the corona, but will enhance its contribution to heating the solar atmosphere.

We proposed that the use of an approximate "jump condition" at the solar transition region permits fast and accurate numerical solutions of the one dimensional hydrodynamic equations when the corona undergoes impulsive heating. In particular, it eliminates the need for the very short timesteps imposed by a highly resolved numerical grid. This paper presents further examples of the applicability of the method for cases of non-uniform heating, in particular, nanoflare trains (uniform in space but non-uniform in time) and spatially localised impulsive heating, including at the loop apex and base of the transition region. In all cases the overall behaviour of the coronal density and temperature shows good agreement with a fully resolved one dimensional model and is significantly better than the equivalent results from a 1D code run without using the jump condition but with the same coarse grid. A detailed assessment of the errors introduced by the jump condition is presented showing that the causes of discrepancy with the fully resolved code are (i) the neglect of the terms corresponding to the rate of change of total energy in the unresolved atmsophere, (ii) mass motions at the base of the transition region and (iii) for some cases with footpoint heating, an over-estimation of the radiative losses in the transition region.

We present a new computational approach that addresses the difficulty of obtaining the correct interaction between the solar corona and the transition region in response to rapid heating events. In the coupled corona, transition region and chromosphere system, an enhanced downward conductive flux results in an upflow (chromospheric evaporation). However, obtaining the correct upflow generally requires high spatial resolution in order to resolve the transition region. With an unresolved transition region, artificially low coronal densities are obtained because the downward heat flux jumps across the unresolved region to the chromosphere, underestimating the upflows. Here, we treat the lower transition region as a discontinuity that responds to changing coronal conditions through the imposition of a jump condition that is derived from an integrated form of energy conservation. To illustrate and benchmark this approach against a fully resolved one-dimensional model, we present field-aligned simulations of coronal loops in response to a range of impulsive (spatially uniform) heating events. We show that our approach leads to a significant improvement in the coronal density evolution than just when using coarse spatial resolutions insufficient to resolve the lower transition region. Our approach compensates for the jumping of the heat flux by imposing a velocity correction that ensures that the energy from the heat flux goes into driving the transition region dynamics, rather than being lost through radiation. Hence, it is possible to obtain improved coronal densities. The advantages of using this approach in both one-dimensional hydrodynamic and three-dimensional magnetohydrodynamic simulations are discussed.

Phase-mixing of Alfv\'en waves in the solar corona has been identified as one possible candidate to explain coronal heating. While this scenario is supported by observations of ubiquitous oscillations in the corona carrying sufficient wave energy and by theoretical models that have described the concentration of energy in small scale structures, it is still unclear whether this wave energy can maintain the million degree solar corona. The aim of this work is to assess how much energy can be converted by a phase-mixing process triggered by the propagation of Alfv\'enic waves in a cylindric coronal structure, such as a coronal loop, and to estimate the impact on the coronal heating. To this end, we run 3D MHD simulations of a magnetised cylinder where the Alfv\'en speed varies through a boundary shell and a footpoint driver is set to trigger kink modes which mode couple to torsional Alfv\'en modes in the boundary shell. These Alfv\'en waves are expected to phase-mix and the system allows us to study the subsequent thermal energy deposition. We run a reference simulation to explain the main process and then we vary simulation parameters. When we take into consideration high values of magnetic resistivity and strong footpoint drivers, we find i) that phase-mixing leads to a temperature increase of the order of $10^5$ K or less, depending on the structure of the boundary shell, ii) that this energy is able to balance the radiative losses only in the localised region involved in the heating, iii) and how the boundary layer and the persistence of the driver influence the thermal structure of the system. Our conclusion is that due to the extreme physical parameters we adopted and the moderate impact on the heating of the system, it is unlikely that phase-mixing can contribute on a global scale to the heating of the solar corona.

Aims - Investigate the effects of resistivity and viscosity on the onset and growth of the Kelvin-Helmholtz instability (KHI) in an oscillating coronal loop. Methods - We modelled a standing kink wave in a density-enhanced loop with the three dimensional (3-D), resistive magnetohydrodynamics code, Lare3d. We conducted a parameter study on the viscosity and resistivity coefficients to examine the effects of dissipation on the KHI. Results - Enhancing the viscosity and resistivity acts to suppress the KHI. Larger values delay the formation of the instability and, in some cases, prevent the onset completely. This leads to the earlier onset of heating for smaller values of the transport coefficients. We note that viscosity has a greater effect on the development of the KHI than resistivity. Furthermore, when using anomalous resistivity, the Ohmic heating rate associated with the KHI may be greater than that associated with the phase mixing that occurs in an instability-suppressed regime (using uniform resistivity). Conclusions - From our study, it is clear that the heating rate crucially depends on the formation of small length scales (influenced by the numerical resolution) as well as the values of resistivity and viscosity. As larger values of the transport coefficients suppress the KHI, the onset of heating is delayed but the heating rate is larger. As increased numerical resolution allows smaller length scales to develop, the heating rate will be higher even for the same transport coefficients.

For future solar missions as well as ground-based telescopes, efficient ways to return and process data have become increasingly important. Solar Orbiter, e.g., which is the next ESA/NASA mission to explore the Sun and the heliosphere, is a deep-space mission, which implies a limited telemetry rate that makes efficient onboard data compression a necessity to achieve the mission science goals. Missions like the Solar Dynamics Observatory (SDO) and future ground-based telescopes such as the Daniel K. Inouye Solar Telescope, on the other hand, face the challenge of making petabyte-sized solar data archives accessible to the solar community. New image compression standards address these challenges by implementing efficient and flexible compression algorithms that can be tailored to user requirements. We analyse solar images from the Atmospheric Imaging Assembly (AIA) instrument onboard SDO to study the effect of lossy JPEG2000 (from the Joint Photographic Experts Group 2000) image compression at different bit rates. To assess the quality of compressed images, we use the mean structural similarity (MSSIM) index as well as the widely used peak signal-to-noise ratio (PSNR) as metrics and compare the two in the context of solar EUV images. In addition, we perform tests to validate the scientific use of the lossily compressed images by analysing examples of an on-disk and off-limb coronal-loop oscillation time-series observed by AIA/SDO.

MHD waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints, but also by wave processes that localise the wave power in undetectable spatial scales. In this study we conduct 3D MHD simulations and forward modelling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin-Helmholtz instability (KHI), resonant absorption and phase mixing. In the presence of a cross-loop temperature gradient we find that emission lines sensitive to the loop core catch different signatures than those more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity modulation produced by the KHI mixing. Common signatures to all considered models include an intensity and loop width modulation at half the kink period, fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of $0.33\arcsec$ and spectral resolution of 25~km~s$^{-1}$, although severe over-estimation of the line width is obtained. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and the KHI motions. We estimate this hidden wave energy to be a factor of $5-10$ of the observed value.

In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfv\'en waves. Due to its localised nature, a direct detection of this mechanism is extremely difficult. Yet, it is the leading theory explaining the observed fast damping of the global transverse waves. However, at odds with this theoretical prediction, recent observations indicate that in the low amplitude regime such transverse MHD waves can also appear decay-less, a yet unsolved phenomenon. Recent numerical work has shown that Kelvin-Helmholtz instabilities (KHI) often accompany transverse MHD waves. In this work, we combine 3D MHD simulations and forward modelling to show that for currently achieved spatial resolution and observed small amplitudes, an apparent decay-less oscillation is obtained. This effect results from the combination of periodic brightenings produced by the KHI and the coherent motion of the KHI vortices amplified by resonant absorption. Such effect is especially clear in emission lines forming at temperatures that capture the boundary dynamics rather than the core, and reflects the low damping character of the local azimuthal Alfv\'en waves resonantly coupled to the kink mode. Due to phase mixing, the detected period can vary depending on the emission line, with those sensitive to the boundary having shorter periods than those sensitive to the loop core. This allows to estimate the density contrast at the boundary.

Propagating perturbations have been observed in extended coronal loop structures for a number of years, but the interpretation in terms of slow (propagating) magneto-acoustic waves and/or as quasi-periodic upflows remains unresolved. We used forward-modelling to construct observational signatures associated with a simple slow magneto-acoustic wave or periodic flow model. Observational signatures were computed for the 171 {\AA} Fe ix and the 193 {\AA} Fe xii spectral lines. Although there are many differences between the flow and wave models, we did not find any clear, robust observational characteristics that can be used in isolation ( i.e. that do not rely on a comparison between the models). For the waves model, a relatively rapid change of the average line widths as a function of (shallow) line-of-sight angles was found, whereas the ratio of the line width amplitudes to the Doppler velocity amplitudes is relatively high for the flow model. The most robust observational signature found is that the ratio of the mean to the amplitudes of the Doppler velocity is always higher than one for the flow model. This ratio is substantially higher for flows than for waves, and for the flows model used in the study is exactly the same in the 171 {\AA} Fe ix and the 193 {\AA} Fe xii spectral lines. However, these potential observational signatures need to be treated cautiously because they are likely to be model-dependent.

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

The solar corona, the tenuous outer atmosphere of the Sun, is orders of magnitude hotter than the solar surface. This 'coronal heating problem' requires the identification of a heat source to balance losses due to thermal conduction, radiation and (in some locations) convection. The review papers in this Theo Murphy meeting issue present an overview of recent observational findings, large- and small-scale numerical modelling of physical processes occurring in the solar atmosphere and other aspects which may affect our understanding of the proposed heating mechanisms. At the same time, they also set out the directions and challenges which must be tackled by future research. In this brief introduction, we summarize some of the issues and themes which reoccur throughout this issue.

We combine observations of the Coronal Multi-channel Polarimeter (CoMP) and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) to study the characteristic properties of (propagating) Alfvenic motions and quasi-periodic intensity disturbances in polar plumes. This unique combination of instruments highlights the physical richness of the processes taking place at the base of the (fast) solar wind. The (parallel) intensity perturbations with intensity enhancements around 1% have an apparent speed of 120 km/s (in both the 171A and 193A passbands) and a periodicity of 15 minutes, while the (perpendicular) Alfvenic wave motions have a velocity amplitude of 0.5 km/s, a phase speed of 830 km/s, and a shorter period of 5 minutes on the same structures. These observations illustrate a scenario where the excited Alfvenic motions are propagating along an inhomogeneously loaded magnetic field structure such that the combination could be a potential progenitor of the magnetohydrodynamic turbulence required to accelerate the fast solar wind.

Recent observations have demonstrated that waves which are capable of carrying large amounts of energy are ubiquitous throughout the solar corona. However, the question of how this wave energy is dissipated (on which time and length scales) and released into the plasma remains largely unanswered. Both analytic and numerical models have previously shown that Alfv\'enic turbulence may play a key role not only in the generation of the fast solar wind, but in the heating of coronal loops. In an effort to bridge the gap between theory and observations, we expand on a recent study [De Moortel et al., ApJL, 782:L34, 2014] by analyzing thirty-seven clearly isolated coronal loops using data from the Coronal Multi-channel Polarimeter (CoMP) instrument. We observe Alfv\'enic perturbations with phase speeds which range from 250-750 km/s and periods from 140-270 s for the chosen loops. While excesses of high frequency wave-power are observed near the apex of some loops (tentatively supporting the onset of Alfv\'enic turbulence), we show that this excess depends on loop length and the wavelength of the observed oscillations. In deriving a proportional relationship between the loop length/wavelength ratio and the enhanced wave power at the loop apex, and from the analysis of the line-widths associated with these loops, our findings are supportive of the existence of Alfv\'enic turbulence in coronal loops.

The presence of turbulent phenomena in the outer solar atmosphere is a given. However, because we are reduced to remotely sensing the atmosphere of a star with instruments of limited spatial and/or spectral resolution, we can only infer the physical progression from macroscopic to microscopic phenomena. Even so, we know that many, if not all, of the turbulent phenomena that pervade interplanetary space have physical origins at the Sun and so in this brief article we consider some recent measurements which point to sustained potential source(s) of heliospheric turbulence in the magnetic and thermal domains. In particular, we look at the scales of magnetism that are imprinted on the outer solar atmosphere by the relentless magneto-convection of the solar interior and combine state-of-the-art observations from the Solar Dynamics Observatory (SDO) and the Coronal Multi-channel Polarimeter (CoMP) which are beginning to hint at the origins of the wave/plasma interplay prevalent closer to the Earth. While linking these disparate scales of observation and understanding of their connection is near to impossible, it is clear that the constant evolution of subsurface magnetism on a host of scales guides and governs the flow of mass and energy at the smallest scales. In the near future significant progress in this area will be made by linking observations from high resolution platforms like the Interface Region Imaging Spectrograph (IRIS) and Advanced Technology Solar Telescope (ATST) with full-disk synoptic observations such as those presented herein.

Waves have long been thought to contribute to the heating of the solar corona and the generation of the solar wind. Recent observations have demonstrated evidence of quasi-periodic longitudinal disturbances and ubiquitous transverse wave propagation in many different coronal environments. This paper investigates signatures of different types of oscillatory behaviour, both above the solar limb and on-disk, by comparing findings from the Coronal Multi-channel Polarimeter (CoMP) and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) for the same active region. We study both transverse and longitudinal motion by comparing and contrasting time-distance images of parallel and perpendicular cuts along/across active region fan loops. Comparisons between parallel space-time features in CoMP Doppler velocity and transverse oscillations in AIA images are made, together with space-time analysis of propagating quasi-periodic intensity features seen near the base of loops in AIA. Signatures of transverse motions are observed along the same magnetic structure using CoMP Doppler velocity (Vphase=600-750km/s, P=3-6mins) and in AIA/SDO above the limb (P=3-8mins). Quasi-periodic intensity features (Vphase=100-200km/s, P=6-11mins) also travel along the base of the same structure. On the disk, signatures of both transverse and longitudinal intensity features were observed by AIA; both show similar properties to signatures found along structures anchored in the same active region three days earlier above the limb. Correlated features are recovered by space-time analysis of neighbouring tracks over perpendicular distances of <2.6Mm.

We use images of high spatial, spectral and temporal resolution, obtained using both ground- and space-based instrumentation, to investigate the coupling between wave phenomena observed at numerous heights in the solar atmosphere. Intensity oscillations of 3 minutes are observed to encompass photospheric umbral dot structures, with power at least three orders-of-magnitude higher than the surrounding umbra. Simultaneous chromospheric velocity and intensity time series reveal an 87 \pm 8 degree out-of-phase behavior, implying the presence of standing modes created as a result of partial wave reflection at the transition region boundary. An average blue-shifted Doppler velocity of ~1.5 km/s, in addition to a time lag between photospheric and chromospheric oscillatory phenomena, confirms the presence of upwardly-propagating slow-mode waves in the lower solar atmosphere. Propagating oscillations in EUV intensity are detected in simultaneous coronal fan structures, with a periodicity of 172 \pm 17 s and a propagation velocity of 45 \pm 7 km/s. Numerical simulations reveal that the damping of the magneto-acoustic wave trains is dominated by thermal conduction. The coronal fans are seen to anchor into the photosphere in locations where large-amplitude umbral dot oscillations manifest. Derived kinetic temperature and emission measure time-series display prominent out-of-phase characteristics, and when combined with the previously established sub-sonic wave speeds, we conclude that the observed EUV waves are the coronal counterparts of the upwardly-propagating magneto-acoustic slow-modes detected in the lower solar atmosphere. Thus, for the first time, we reveal how the propagation of 3 minute magneto-acoustic waves in solar coronal structures is a direct result of amplitude enhancements occurring in photospheric umbral dots.

Determining the heating mechanism (or mechanisms) that causes the outer atmosphere of the Sun, and many other stars, to reach temperatures orders of magnitude higher than their surface temperatures has long been a key problem. For decades the problem has been known as the coronal heating problem, but it is now clear that `coronal heating' cannot be treated or explained in isolation and that the heating of the whole solar atmosphere must be studied as a highly coupled system. The magnetic field of the star is known to play a key role, but, despite significant advancements in solar telescopes, computing power and much greater understanding of theoretical mechanisms, the question of which mechanism or mechanisms are the dominant supplier of energy to the chromosphere and corona is still open. Following substantial recent progress, we consider the most likely contenders and discuss the key factors that have made, and still make, determining the actual (coronal) heating mechanism (or mechanisms) so difficult.

The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnetic energy release in this environment. This paper investigates the role of the Hall term in the propagation and dissipation of waves, their interaction with 2D magnetic X-points and the nature of the resulting reconnection. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Depending on the ratio of ion skin depth to system size, our model demonstrates that Hall effects can play a key role in the wave-null interaction. In particular, the initial fast-wave pulse now consists of whistler and ion-cyclotron components; the dispersive nature of the whistler wave leads to (i) earlier interaction with the null, (ii) the creation of multiple additional, transient nulls and, hence, an increased number of energy release sites. In the Hall regime, the relevant timescales (such as the onset of reconnection and the period of the oscillatory relaxation) of the system are reduced significantly, and the reconnection rate is enhanced.

Quasi-periodic disturbances have been observed in the outer solar atmosphere for many years now. Although first interpreted as upflows (Schrijver et al. (1999)), they have been widely regarded as slow magnetoacoustic waves, due to observed velocities and periods. However, recent observations have questioned this interpretation, as periodic disturbances in Doppler velocity, line width and profile asymmetry were found to be in phase with the intensity oscillations (De Pontieu et al. (2010),Tian1 et al. (2011))}, suggesting the disturbances could be quasi-periodic upflows. Here we conduct a detailed analysis of the velocities of these disturbances across several wavelengths using the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We analysed 41 examples, including both sunspot and non sunspot regions of the Sun. We found that the velocities of propagating disturbances (PDs) located at sunspots are more likely to be temperature dependent, whereas the velocities of PDs at non sunspot locations do not show a clear temperature dependence. We also considered on what scale the underlying driver is affecting the properties of the PDs. Finally, we found that removing the contribution due to the cooler ions in the 193 A wavelength suggests that a substantial part of the 193 emission of sunspot PDs can be contributed to the cool component of 193\AA.

We present a comprehensive review of MHD wave behaviour in the neighbourhood of coronal null points: locations where the magnetic field, and hence the local Alfven speed, is zero. The behaviour of all three MHD wave modes, i.e. the Alfven wave and the fast and slow magnetoacoustic waves, has been investigated in the neighbourhood of 2D, 2.5D and (to a certain extent) 3D magnetic null points, for a variety of assumptions, configurations and geometries. In general, it is found that the fast magnetoacoustic wave behaviour is dictated by the Alfven-speed profile. In a $\beta=0$ plasma, the fast wave is focused towards the null point by a refraction effect and all the wave energy, and thus current density, accumulates close to the null point. Thus, null points will be locations for preferential heating by fast waves. Independently, the Alfven wave is found to propagate along magnetic fieldlines and is confined to the fieldlines it is generated on. As the wave approaches the null point, it spreads out due to the diverging fieldlines. Eventually, the Alfven wave accumulates along the separatrices (in 2D) or along the spine or fan-plane (in 3D). Hence, Alfven wave energy will be preferentially dissipated at these locations. It is clear that the magnetic field plays a fundamental role in the propagation and properties of MHD waves in the neighbourhood of coronal null points. This topic is a fundamental plasma process and results so far have also lead to critical insights into reconnection, mode-coupling, quasi-periodic pulsations and phase-mixing.

Recent observations have revealed that MHD waves and oscillations are ubiquitous in the solar atmosphere, with a wide range of periods. We give a brief review of some aspects of MHD waves and coronal seismology which have recently been the focus of intense debate or are newly emerging. In particular, we focus on four topics: (i) the current controversy surrounding propagating intensity perturbations along coronal loops, (ii) the interpretation of propagating transverse loop oscillations, (iii) the ongoing search for coronal (torsional) Alfven waves and (iv) the rapidly developing topic of quasi-periodic pulsations (QPP) in solar flares.

Spectroscopic and stereoscopic imaging observations of slow magnetoacoustic wave propagation within a coronal loop are investigated to determine the decay length scale of the slow magnetoacoustic mode in three dimensions and the density profile within the loop system. The slow wave is found to have an e-folding decay length scale of $20,000^{+4000}_{-3000}$km with a uniform density profile along the loop base. These observations place quantitive constraints on the modelling of wave propagation within coronal loops. Theoretical forward modelling suggests that magnetic field line divergence is the dominant damping factor and thermal conduction is insufficient, given the observed parameters of the coronal loop temperature, density and wave mode period.

We investigate the behaviour of nonlinear, nonideal Alfv\'en wave propagation within an inhomogeneous magnetic environment. The governing MHD equations are solved in 1D and 2D using both analytical techniques and numerical simulations. We find clear evidence for the ponderomotive effect and visco-resistive heating. The ponderomotive effect generates a longitudinal component to the transverse Alfv\'en wave, with a frequency twice that of the driving frequency. Analytical work shows the addition of resistive heating. This leads to a substantial increase in the local temperature and thus gas pressure of the plasma, resulting in material being pushed along the magnetic field. In 2D, our system exhibits phase mixing and we observe an evolution in the location of the maximum heating, i.e. we find a drifting of the heating layer. Considering Alfv\'en wave propagation in 2D with an inhomogeneous density gradient, we find that the equilibrium density profile is significantly modified by both the flow of density due to visco-resistive heating and the nonlinear response to the localised heating through phase mixing.

Aims. To determine the effect of the Hall term in the generalised Ohm's law on the damping and phase mixing of Alfven waves in the ion cyclotron range of frequencies in uniform and non-uniform equilibrium plasmas. Methods. Wave damping in a uniform plasma is treated analytically, whilst a Lagrangian remap code (Lare2d) is used to study Hall effects on damping and phase mixing in the presence of an equilibrium density gradient. Results. The magnetic energy associated with an initially Gaussian field perturbation in a uniform resistive plasma is shown to decay algebraically at a rate that is unaffected by the Hall term to leading order in k^2di^2 where k is wavenumber and di is ion skin depth. A similar algebraic decay law applies to whistler perturbations in the limit k^2di^2>>1. In a non-uniform plasma it is found that the spatially-integrated damping rate due to phase mixing is lower in Hall MHD than it is in MHD, but the reduction in the damping rate, which can be attributed to the effects of wave dispersion, tends to zero in both the weak and strong phase mixing limits.

This paper extends the models of Craig & McClymont (1991) and McLaughlin & Hood (2004) to include finite $\beta$ and nonlinear effects. We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001) and consider an initial condition in $ {\bf{v}}\times{\bf{B}} \cdot {\hat{\bf{z}}}$ (a natural variable of the system). We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a 'cusp-like' point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime.

We use a combination of analytical theory, numerical simulation, and data analysis to study the propagation of acoustic waves along coronal loops. We show that the intensity perturbation of a wave depends on a number of factors, including dissipation of the wave energy, pressure and temperature gradients in the loop atmosphere, work action between the wave and a flow, and the sensitivity properties of the observing instrument. In particular, the scale length of the intensity perturbation varies directly with the dissipation scale length (i.e., damping length) and the scale lengths of pressure, temperature, and velocity. We simulate wave propagation in three different equilibrium loop models and find that dissipation and pressure and temperature stratification are the most important effects in the low corona where the waves are most easily detected. Velocity effects are small, and cross-sectional area variations play no direct role for lines-of-sight that are normal to the loop axis. The intensity perturbation scale lengths in our simulations agree very well with the scale lengths we measure in a sample of loops observed by TRACE. The median observed value is 4.35x10^9 cm. In some cases the intensity perturbation increases with height, which is likely an indication of a temperature inversion in the loop (i.e., temperature that decreases with height). Our most important conclusion is that thermal conduction, the primary damping mechanism, is accurately described by classical transport theory. There is no need to invoke anomalous processes to explain the observations.