Showing 1–50 of 54 results for author: De Moortel, I

MUSE observations of small-scale heating events

Authors: C. A. Breu, I. De Moortel, P. Testa

Abstract: Constraining the processes that drive coronal heating from observations is a difficult task due to the complexity of the solar atmosphere. As upcoming missions such as MUSE will provide coronal observations with unprecedented spatial and temporal resolution, numerical simulations are becoming increasingly realistic. Despite the availability of synthetic observations from numerical models, line-of-… ▽ More Constraining the processes that drive coronal heating from observations is a difficult task due to the complexity of the solar atmosphere. As upcoming missions such as MUSE will provide coronal observations with unprecedented spatial and temporal resolution, numerical simulations are becoming increasingly realistic. Despite the availability of synthetic observations from numerical models, line-of-sight effects and the complexity of the magnetic topology in a realistic setup still complicate the prediction of signatures for specific heating processes. 3D MHD simulations have shown that a significant part of the Poynting flux injected into the solar atmosphere is carried by small-scale motions, such as vortices driven by rotational flows inside intergranular lanes. MHD waves excited by these vortices have been suggested to play an important role in the energy transfer between different atmospheric layers. Using synthetic spectroscopic data generated from a coronal loop model incorporating realistic driving by magnetoconvection, we study whether signatures of energy transport by vortices and eventual dissipation can be identified with future missions such as MUSE. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 15 pages, 15 figures

Non-thermal broadening of coronal lines in a 3D MHD loop model

Authors: C. A. Breu, H. Peter, S. K. Solanki, R. Cameron, I. De Moortel

Abstract: Observed spectral profiles of emission lines from the corona are found to have widths exceeding the thermal line width. To investigate the physical mechanism, we run a 3D MHD model of a single, straightened loop in which we partially resolve turbulent motions that form in response to the driving by self-consistently evolving magneto-convection in the photosphere. The convective motions shear and t… ▽ More Observed spectral profiles of emission lines from the corona are found to have widths exceeding the thermal line width. To investigate the physical mechanism, we run a 3D MHD model of a single, straightened loop in which we partially resolve turbulent motions that form in response to the driving by self-consistently evolving magneto-convection in the photosphere. The convective motions shear and twist the magnetic field lines, leading to heating. From the model we synthesize spectral profiles of emission lines forming at temperatures around and above 1 MK. The coronal heating process generates a range of velocity amplitudes and directions structured on a scale much smaller than the resolving power of current instruments, leading to a broadening of the spectral lines. Our model includes the mass exchange between corona and chromosphere, thus we also capture flows parallel to the loop axis. We find that the spectral lines show a non-thermal line broadening roughly consistent with observations for a viewing angle perpendicular to the axis. The broadening through field-parallel flows is comparable, although slightly smaller. The line broadening is independent of the instrument resolution for a perpendicular line-of-sight. We can connect the non-thermal line broadening to heating events and flows. While small-scale velocities along the line-of-sight are mainly responsible for the broadening observed perpendicular to the loop, chromospheric evaporation is important for the line broadening observed along the loop. The model reproduces observed values for non-thermal line widths. In the model these result from continuous driving by magnetoconvection, without imposing driving motions or starting from an already braided field. △ Less

Submitted 29 March, 2024; originally announced April 2024.

Comments: 18 pages, 18 figures

The effects of driving time scales on coronal heating in a stratified atmosphere

Authors: Thomas Howson, Ineke De Moortel

Abstract: Aims: We investigate the atmospheric response to coronal heating driven by random flows with different characteristic time scales and amplitudes. Methods: We conducted a series of 3D MHD simulations of random driving imposed on a gravitationally stratified model of the solar atmosphere. In order to understand differences between alternating current (AC) and direct current (DC) heating, we consider… ▽ More Aims: We investigate the atmospheric response to coronal heating driven by random flows with different characteristic time scales and amplitudes. Methods: We conducted a series of 3D MHD simulations of random driving imposed on a gravitationally stratified model of the solar atmosphere. In order to understand differences between alternating current (AC) and direct current (DC) heating, we considered the effects of changing the characteristic time scales of the imposed velocities. We also investigated the effects of the magnitude of the velocity driving. Results: Complex foot point motions lead to a proliferation of current sheets and energy dissipation throughout the corona. For a given amplitude, DC driving typically leads to a greater rate of energy injection when compared to AC driving. This leads to the formation of larger currents, increased heating rates and higher temperatures in DC simulations. There is no difference in the spatial distribution of energy dissipation across simulations, however, energy release events in AC cases tend to be more frequent and last for less time. Higher velocity driving is associated with larger currents, higher temperatures and the corona occupying a larger fraction of the simulation volume. In all cases, most of heating is associated with small energy release events, which occur much more frequently than large events. Conclusions: When combined with observational results showing a greater abundance of power in low frequency modes, these findings suggest that energy release in the corona is likely to be driven by long time scale motions. In the corona, AC and DC driving will occur concurrently and their effects remain difficult to isolate. The distribution of field line temperatures and the asymmetry of temperature profiles may reveal the frequency and longevity of energy release events and therefore the relative importance of AC and DC heating. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Comments: 14 pages, 15 figures

Journal ref: A&A 661, A144 (2022)

Forward modelling of heating within a coronal arcade

Authors: Lianne Fyfe, Thomas Howson, Ineke De Moortel

Abstract: Aims. We investigate the synthetic observational signatures from numerical models of coronal heating in an arcade to determine what features are associated with such heating, and what tools can be used to identify them. Methods. We consider two simulations of coronal arcades driven by footpoint motions with different characteristic timescales. Forward modelling is then used, and the synthetic em… ▽ More Aims. We investigate the synthetic observational signatures from numerical models of coronal heating in an arcade to determine what features are associated with such heating, and what tools can be used to identify them. Methods. We consider two simulations of coronal arcades driven by footpoint motions with different characteristic timescales. Forward modelling is then used, and the synthetic emission data is analysed. Results. The total intensity and Doppler velocities clearly show the magnetic structure of the coronal arcade. Contrasts in the local Doppler shift also highlight the locations of separatrix surfaces. The distinguishing feature of the AC and DC models is that of the frequencies. Through FFT analysis of the Doppler velocities, when short timescale footpoint motions are present, higher frequencies are observed. For longer timescale motions, the dominant signal is that of lower frequencies; however, higher frequencies were also detected, which matched the natural Alfvén frequency of the background magnetic field. Alfvénic wave signatures were identified in both models with fast wave signatures observable in the AC model. Finally, the estimates of the kinetic energy using the Doppler shifts were found to significantly underestimate within these models. Conclusions. Observables identified within this article were from features such as Alfvén waves, fast waves, the arcade structure and separatrix surfaces. The two models were differentiated by examining the frequencies present. The Doppler velocities cannot provide accurate estimates of the total kinetic energy or the component parallel to the LOS. This is due to some plasma outside the formation temperature range of the ion, the multi-directional driver, and cancellation of the velocity along the LOS. The impact each factor has on the estimation is dependent on the set up of the model and the chosen emission line. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 15 pages and 18 figures

Journal ref: A&A 656, A120 (2021)

Investigating coronal wave energy estimates using synthetic non-thermal line widths

Authors: Lianne Fyfe, Thomas Howson, Ineke De Moortel, Vaibhav Pant, Tom Van Doorsselaere

Abstract: Aims. Estimates of coronal wave energy remain uncertain as a large fraction of the energy is likely hidden in the non-thermal line widths of emission lines. In order to estimate these wave energies, many previous studies have considered the root mean squared wave amplitudes to be a factor of $\sqrt{2}$ greater than the non-thermal line widths. However, other studies have used different factors. To… ▽ More Aims. Estimates of coronal wave energy remain uncertain as a large fraction of the energy is likely hidden in the non-thermal line widths of emission lines. In order to estimate these wave energies, many previous studies have considered the root mean squared wave amplitudes to be a factor of $\sqrt{2}$ greater than the non-thermal line widths. However, other studies have used different factors. To investigate this problem, we consider the relation between wave amplitudes and the non-thermal line widths within a variety of 3D magnetohydrodynamic (MHD) simulations. Methods. We consider the following 3D numerical models: Alfvén waves in a uniform magnetic field, transverse waves in a complex braided magnetic field, and two simulations of coronal heating in an arcade. We applied the forward modelling code FoMo to generate the synthetic emission data required to analyse the non-thermal line widths. Results. Determining a single value for the ratio between the non-thermal line widths and the root mean squared wave amplitudes is not possible across multiple simulations. It was found to depend on a variety of factors, including line-of-sight angles, velocity magnitudes, wave interference, and exposure time. Indeed, some of our models achieved the values claimed in recent articles while other more complex models deviated from these ratios. Conclusions. To estimate wave energies, an appropriate relation between the non-thermal line widths and root mean squared wave amplitudes is required. However, evaluating this ratio to be a singular value, or even providing a lower or upper bound on it, is not realistically possible given its sensitivity to various MHD models and factors. As the ratio between wave amplitudes and non-thermal line widths is not constant across our models, we suggest that this widely used method for estimating wave energy is not robust. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: 11 pages and 11 figures

Journal ref: A&A 656, A56 (2021)

Magnetic reconnection and the Kelvin-Helmholtz instability in the solar corona

Authors: Thomas Howson, Ineke De Moortel, David Pontin

Abstract: Context. The magnetic Kelvin-Helmholtz instability (KHI) has been proposed as a means of generating magnetohy- drodynamic turbulence and encouraging wave energy dissipation in the solar corona, particularly within transversely oscillating loops. Aims. Our goal is to determine whether the KHI encourages magnetic reconnection in oscillating flux tubes in the solar corona. This will establish whether… ▽ More Context. The magnetic Kelvin-Helmholtz instability (KHI) has been proposed as a means of generating magnetohy- drodynamic turbulence and encouraging wave energy dissipation in the solar corona, particularly within transversely oscillating loops. Aims. Our goal is to determine whether the KHI encourages magnetic reconnection in oscillating flux tubes in the solar corona. This will establish whether the instability enhances the dissipation rate of energy stored in the magnetic field. Methods. We conducted a series of three-dimensional magnetohydrodynamic simulations of the KHI excited by an oscillating velocity shear. We investigated the effects of numerical resolution, field line length, and background currents on the growth rate of the KHI and on the subsequent rate of magnetic reconnection. Results. The KHI is able to trigger magnetic reconnection in all cases, with the highest rates occurring during the initial growth phase. Reconnection is found to occur preferentially along the boundaries of Kelvin-Helmholtz vortices, where the shear in the velocity and magnetic fields is greatest. The estimated rate of reconnection is found to be lowest in simulations where the KHI growth rate is reduced. For example, this is the case for shorter field lines or due to shear in the background field. Conclusions. In non-ideal regimes, the onset of the instability causes the local reconnection of magnetic field lines and enhances the rate of coronal wave heating. However, we found that if the equilibrium magnetic field is sheared across the Kelvin-Helmholtz mixing layer, the instability does not significantly enhance the rate of reconnection of the background field, despite the free energy associated with the non-potential field. △ Less

Submitted 30 September, 2021; originally announced September 2021.

Comments: 17 pages, 19 figures

Probing the physics of the solar atmosphere with the Multi-slit Solar Explorer (MUSE): I. Coronal Heating

Authors: Bart De Pontieu, Paola Testa, Juan Martinez-Sykora, Patrick Antolin, Konstantinos Karampelas, Viggo Hansteen, Matthias Rempel, Mark C. M. Cheung, Fabio Reale, Sanja Danilovic, Paolo Pagano, Vanessa Polito, Ineke De Moortel, Daniel Nobrega-Siverio, Tom Van Doorsselaere, Antonino Petralia, Mahboubeh Asgari-Targhi, Paul Boerner, Mats Carlsson, Georgios Chintzoglou, Adrian Daw, Ed DeLuca, Leon Golub, Takuma Matsumoto, Ignacio Ugarte-Urra , et al. (2 additional authors not shown)

Abstract: The Multi-slit Solar Explorer (MUSE) is a proposed NASA MIDEX mission, currently in Phase A, composed of a multi-slit EUV spectrograph (in three narrow spectral bands centered around 171A, 284A, and 108A) and an EUV context imager (in two narrow passbands around 195A and 304A). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (<0.5 arcsec), and t… ▽ More The Multi-slit Solar Explorer (MUSE) is a proposed NASA MIDEX mission, currently in Phase A, composed of a multi-slit EUV spectrograph (in three narrow spectral bands centered around 171A, 284A, and 108A) and an EUV context imager (in two narrow passbands around 195A and 304A). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (<0.5 arcsec), and temporal resolution (down to ~0.5s) thanks to its innovative multi-slit design. By obtaining spectra in 4 bright EUV lines (Fe IX 171A , Fe XV 284A, Fe XIX-Fe XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will for the first time be able to "freeze" (at a cadence as short as 10 seconds) with a spectroscopic raster the evolution of the dynamic coronal plasma over a wide range of scales: from the spatial scales on which energy is released (~0.5 arcsec) to the large-scale often active-region size (170 arcsec x 170 arcsec) atmospheric response. We use advanced numerical modeling to showcase how MUSE will constrain the properties of the solar atmosphere on the spatio-temporal scales (~0.5 arcsec, ~20 seconds) and large field-of-view on which various state-of-the-art models of the physical processes that drive coronal heating, solar flares and coronal mass ejections (CMEs) make distinguishing and testable predictions. We describe how the synergy between MUSE, the single-slit, high-resolution Solar-C EUVST spectrograph, and ground-based observatories (DKIST and others) can address how the solar atmosphere is energized, and the critical role MUSE plays because of the multi-scale nature of the physical processes involved. In this first paper, we focus on how comparisons between MUSE observations and theoretical models will significantly further our understanding of coronal heating mechanisms. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 46 pages, 22 figures, submitted to ApJ

A fast multi-dimensional magnetohydrodynamic formulation of the transition region adaptive conduction (TRAC) method

Authors: C. D. Johnston, A. W. Hood, I. De Moortel, P. Pagano, T. A. Howson

Abstract: We have demonstrated that the Transition Region Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, successfully removing the influence of numerical resolution on the coronal density response to impulsive heating. This is achieved by adjusting the parallel thermal conductivity, radiative loss, and heating rates to broaden the… ▽ More We have demonstrated that the Transition Region Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, successfully removing the influence of numerical resolution on the coronal density response to impulsive heating. This is achieved by adjusting the parallel thermal conductivity, radiative loss, and heating rates to broaden the transition region (TR), below a global cutoff temperature, so that the steep gradients are spatially resolved even when using coarse numerical grids. Implementing the original 1D formulation of TRAC in multi-dimensional magnetohydrodynamic (MHD) models would require tracing a large number of magnetic field lines at every time step in order to prescribe a global cutoff temperature to each field line. In this paper, we present a highly efficient formulation of the TRAC method for use in multi-dimensional MHD simulations, which does not rely on field line tracing. In the TR, adaptive local cutoff temperatures are used instead of global cutoff temperatures to broaden any unresolved parts of the atmosphere. These local cutoff temperatures are calculated using only local grid cell quantities, enabling the MHD extension of TRAC to efficiently account for the magnetic field evolution, without tracing field lines. Consistent with analytical predictions, we show that this approach successfully preserves the properties of the original TRAC method. In particular, the total radiative losses and heating remain conserved under the MHD formulation. Results from 2D MHD simulations of impulsive heating in unsheared and sheared arcades of coronal loops are also presented. These simulations benchmark the MHD TRAC method against a series of 1D models and demonstrate the versatility and robustness of the method in multi-dimensional magnetic fields. We show, for the first time, that pressure differences, ... △ Less

Submitted 21 June, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 18 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 654, A2 (2021)

MHD Waves in open coronal structures

Authors: D. Banerjee, S. Krishna Prasad, V. Pant, J. A. McLaughlin, P. Antolin, N. Magyar, L. Ofman, H. Tian, T. Van Doorsselaere, I. De Moortel, T. Wang

Abstract: Modern observatories have revealed the ubiquitous presence of magnetohydrodynamic waves in the solar corona. The propagating waves (in contrast to the standing waves) are usually originated in the lower solar atmosphere which makes them particularly relevant to coronal heating. Furthermore, open coronal structures are believed to be the source regions of solar wind, therefore, the detection of MHD… ▽ More Modern observatories have revealed the ubiquitous presence of magnetohydrodynamic waves in the solar corona. The propagating waves (in contrast to the standing waves) are usually originated in the lower solar atmosphere which makes them particularly relevant to coronal heating. Furthermore, open coronal structures are believed to be the source regions of solar wind, therefore, the detection of MHD waves in these structures is also pertinent to the acceleration of solar wind. Besides, the advanced capabilities of the current generation telescopes have allowed us to extract important coronal properties through MHD seismology. The recent progress made in the detection, origin, and damping of both propagating slow mangetoacoustic waves and kink (Alfvénic) waves is presented in this review article especially in the context of open coronal structures. Where appropriate, we give an overview on associated theoretical modelling studies. A few of the important seismological applications of these waves are discussed. The possible role of Aflvénic waves in the acceleration of solar wind is also touched upon. △ Less

Submitted 21 September, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

Comments: Accepted for publication in the Space Science Reviews. 48 pages and 11 figures

Coronal heating by MHD waves

Authors: Tom Van Doorsselaere, Abhishek K. Srivastava, Patrick Antolin, Norbert Magyar, Soheil Vasheghani Farahani, Hui Tian, Dmitrii Y. Kolotkov, Leon Ofman, Mingzhe Guo, Iñigo Arregui, Ineke De Moortel, David Pascoe

Abstract: The heating of the solar chromosphere and corona to the observed high temperatures, imply the presence of ongoing heating that balances the strong radiative and thermal conduction losses expected in the solar atmosphere. It has been theorized for decades that the required heating mechanisms of the chromospheric and coronal parts of the active regions, quiet-Sun, and coronal holes are associated wi… ▽ More The heating of the solar chromosphere and corona to the observed high temperatures, imply the presence of ongoing heating that balances the strong radiative and thermal conduction losses expected in the solar atmosphere. It has been theorized for decades that the required heating mechanisms of the chromospheric and coronal parts of the active regions, quiet-Sun, and coronal holes are associated with the solar magnetic fields. However, the exact physical process that transport and dissipate the magnetic energy which ultimately leads to the solar plasma heating are not yet fully understood. The current understanding of coronal heating relies on two main mechanism: reconnection and MHD waves that may have various degrees of importance in different coronal regions. In this review we focus on recent advances in our understanding of MHD wave heating mechanisms. First, we focus on giving an overview of observational results, where we show that different wave modes have been discovered in the corona in the last decade, many of which are associated with a significant energy flux, either generated in situ or pumped from the lower solar atmosphere. Afterwards, we summarise the recent findings of numerical modelling of waves, motivated by the observational results. Despite the advances, only 3D MHD models with Alfvén wave heating in an unstructured corona can explain the observed coronal temperatures compatible with the quiet Sun, while 3D MHD wave heating models including cross-field density structuring are not yet able to account for the heating of coronal loops in active regions to their observed temperature. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: accepted in Space Science Reviews

Forward Modelling of MHD Waves in Braided Magnetic Fields

Authors: Lianne Fyfe, Thomas Howson, Ineke De Moortel

Abstract: Aim. We investigate synthetic observational signatures generated from numerical models of transverse waves propagating in braided magnetic fields. Methods. We examine two simulations with different levels of magnetic field braiding and impose a periodic, transverse wave driver at the lower boundary. After waves reflect off the top boundary, a complex pattern of wave interference forms. We analyse… ▽ More Aim. We investigate synthetic observational signatures generated from numerical models of transverse waves propagating in braided magnetic fields. Methods. We examine two simulations with different levels of magnetic field braiding and impose a periodic, transverse wave driver at the lower boundary. After waves reflect off the top boundary, a complex pattern of wave interference forms. We analyse the synthetic emission produced by the forward modelling code, FoMo. We examine line intensity, Doppler shifts and kinetic energy along several line-of-sight (LOS) angles. Results. The Doppler shifts showed the presence of transverse waves. However, when analysing the intensity, waves are less easily observed for more complex magnetic fields and may be mistaken for background noise. Depending on the LOS angle, the observable signatures of waves reflect some of the magnetic field braiding, particularly when multiple emission lines are available. In the more braided simulation, signatures of phase mixing can be identified. We identify possible ambiguities in the interpretation of wave modes based on the synthetic emission signatures. Conclusions. Most of the observables within this article behave as expected, given knowledge of the evolution of the parameters in 3D space. However, some intriguing observational signatures are present. Detecting regions of magnetic field complexity is somewhat possible when waves are present. However, simultaneous spectroscopic imaging from different lines is important to identify these locations. Care needs to be taken when interpreting intensity and Doppler velocity signatures as torsional motions as, in our setup, such signatures were a consequence of the magnetic field complexity and not torsional waves. Finally, the kinetic energy, estimated with Doppler velocities, is dependent on the polarisation of the wave, complexity of the background field and the LOS. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 15 pages and 22 figures

Journal ref: A&A 643, A86 (2020)

The effects of driving time scales on heating in a coronal arcade

Authors: Thomas Howson, Ineke De Moortel, Lianne Fyfe

Abstract: Context. The relative importance of AC and DC heating in maintaining the temperature of the corona is not well constrained. Aims. Investigate the effects of the characteristic time scales of photospheric driving on the injection and dissipation of energy within a coronal arcade. Methods. We have conducted three dimensional MHD simulations of foot point driving imposed on an arcade. We modified the… ▽ More Context. The relative importance of AC and DC heating in maintaining the temperature of the corona is not well constrained. Aims. Investigate the effects of the characteristic time scales of photospheric driving on the injection and dissipation of energy within a coronal arcade. Methods. We have conducted three dimensional MHD simulations of foot point driving imposed on an arcade. We modified the typical driving time scales to understand the efficiency of heating obtained using AC and DC drivers. We considered the implications for the injected Poynting flux and the nature of the energy release in dissipative regimes. Results. For the same driver amplitude and complexity, long time scale motions are able to inject a much greater Poynting flux into the corona. Consequently, in non-ideal regimes, slow stressing motions result in a greater increase in plasma temperature than for wave-like driving. In dissipative simulations, Ohmic heating is found to be much more significant than viscous heating. For all drivers in our parameter space, energy dissipation is greatest close to the base of the arcade where the magnetic field strength is strongest and at separatrix surfaces, where the field connectivity changes. Across all simulations, the background field is stressed with random foot point motions (in a manner typical of DC heating studies) and even for short time scale driving, the injected Poynting flux is large given the small amplitude flows considered. For long time scale driving, the rate of energy injection was comparable to the expected requirements in active regions. The heating rates were found to scale with the perturbed magnetic field strength and not the total field strength. Conclusions. Alongside recent studies which show power within the corona is dominated by low frequency motions, our results suggest that in the closed corona, DC heating is more significant than AC heating. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 14 pages, 16 figures

Journal ref: A&A 643, A85 (2020)

Effect of Coronal Loop Structure on Wave Heating by Phase Mixing

Authors: P. Pagano, I. De Moortel, R. J. Morton

Abstract: The mechanism behind coronal heating still elude direct observation and modelling of viable theoretical processes and the subsequent effect on coronal structures is one of the key tools available to assess possible heating mechanisms. Wave-heating via phase-mixing of Magnetohydrodynamics (MHD) transverse waves has been proposed as a possible way to convert magnetic energy into thermal energy but i… ▽ More The mechanism behind coronal heating still elude direct observation and modelling of viable theoretical processes and the subsequent effect on coronal structures is one of the key tools available to assess possible heating mechanisms. Wave-heating via phase-mixing of Magnetohydrodynamics (MHD) transverse waves has been proposed as a possible way to convert magnetic energy into thermal energy but increasingly, MHD models suggest this is not a sufficiently efficient mechanism. We model heating by phase-mixing of transverse MHD waves in various configurations, to investigate whether certain circumstances can enhance the heating sufficiently to sustain the million degree solar corona and to assess the impact of the propagation and phase-mixing of transverse MHD waves on the structure of the boundary shell of coronal loops. We use 3D MHD simulations of a pre-existing density enhancement in magnetised medium and a boundary driver to trigger the propagation of transverse waves with the same power spectrum as measured by the COmP. We consider different density structures, boundary conditions at the non-drive footpoint, characteristics of the driver, and different forms of magnetic resistivity. We find that different initial density structures affect the evolution of the boundary shell and some driver configurations enhance the heating generated from the dissipation of the MHD waves. In particular, drivers coherent on a larger spatial scale and higher dissipation coefficients generate significant heating, although it is still insufficient. We conclude that while phase-mixing of transverse MHD waves is unlikely to sustain the thermal structure of the corona, there are configurations that allow for an enhanced efficiency of this mechanism. We provide possible signatures to identify the presence of such configurations, such as the location of where the heating is deposited along the coronal loop △ Less

Submitted 9 September, 2020; originally announced September 2020.

Journal ref: A&A 643, A73 (2020)

Phase mixing and wave heating in a complex coronal plasma

Authors: Thomas Howson, Ineke De Moortel, Jack Reid

Abstract: 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 sys… ▽ More 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. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 14 pages, 19 figures

Journal ref: A&A 636, A40 (2020)

Chromospheric evaporation and phase mixing of Alfvén waves in coronal loops

Authors: H. J. Van Damme, I. De Moortel, P. Pagano, C. D. Johnston

Abstract: Phase mixing of Alfvén 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 wit… ▽ More Phase mixing of Alfvén 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én 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én wave driver was implemented. As the Alfvén 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én waves, and the subsequent evaporation from the lower atmosphere. We find that phase mixing of Alfvén 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. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: 11 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 635, A174 (2020)

Modelling the solar transition region using an adaptive conduction method

Authors: C. D. Johnston, P. J. Cargill, A. W. Hood, I. De Moortel, S. J. Bradshaw, A. C. Vaseekar

Abstract: 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 an… ▽ More 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. △ Less

Submitted 28 February, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: 20 pages, 9 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 635, A168 (2020)

Resonant absorption in expanding coronal magnetic flux tubes with uniform density

Authors: Thomas Howson, Ineke De Moortel, Patrick Antolin, Tom Van Doorsselaere, Andrew Wright

Abstract: 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 i… ▽ More 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. △ Less

Submitted 24 September, 2019; originally announced September 2019.

Journal ref: A&A 631, A105 (2019)

Magnetohydrodynamic waves in braided magnetic fields

Authors: Thomas Howson, Ineke De Moortel, Jack Reid, Alan Hood

Abstract: 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énic wave phase mixing in this regime. Spatial gradients in the local Alfvén speed and variations in the l… ▽ More 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énic wave phase mixing in this regime. Spatial gradients in the local Alfvén 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én 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. △ Less

Submitted 8 August, 2019; originally announced August 2019.

MHD simulations of the in-situ generation of kink and sausage waves in the solar corona by collision of dense plasma clumps

Authors: P. Pagano, H. J. Van Damme, P. Antolin, I. De Moortel

Abstract: 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 th… ▽ More 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. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Journal ref: A&A 626, A53 (2019)

The effects of numerical resolution, heating timescales and background heating on thermal non-equilibrium in coronal loops

Authors: C. D. Johnston, P. J. Cargill, P. Antolin, A. W. Hood, I. De Moortel, S. J. Bradshaw

Abstract: 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… ▽ More 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. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 15 pages, 10 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 625, A149 (2019)

Phase mixing of nonlinear Alfven waves

Authors: Alexander Prokopyszyn, Alan Hood, Ineke De Moortel

Abstract: Aims: This paper presents 2.5D numerical experiments of Alfvén 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… ▽ More Aims: This paper presents 2.5D numerical experiments of Alfvén 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én 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. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Journal ref: A&A 624, A90 (2019)

Contribution of observed multi frequency spectrum of Alfvén waves to coronal heating

Authors: Paolo Pagano, Ineke De Moortel

Abstract: 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én waves is considered a candidate mechanism, as it can develop transverse gradient where magnetic energy can be converted into thermal energ… ▽ More 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én 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én 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. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Journal ref: A&A 623, A37 (2019)

Heating effects from driven transverse and Alfvén waves in coronal loops

Authors: Mingzhe Guo, Tom Van Doorsselaere, Kostas Karampelas, Bo Li, Patrick Antolin, Ineke De Moortel

Abstract: 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 suc… ▽ More 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{é}n driver separately and a mixed driver at the bottom of flux tubes. Both the transverse and Alfv{é}n oscillations can lead to the KHI. Meanwhile, the Alfv{é}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{é}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{é}n modes nor the fine structures are observable. Therefore, this numerical study shows that Alfv{é}n modes probably can co-exist with kink modes, leading to enhanced heating. △ Less

Submitted 19 November, 2018; originally announced November 2018.

In-situ generation of transverse MHD waves from colliding flows in the solar corona

Authors: Patrick Antolin, Paolo Pagano, Ineke De Moortel, Valery M. Nakariakov

Abstract: 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 betwee… ▽ More 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 $β$ 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. △ Less

Submitted 1 July, 2018; originally announced July 2018.

Comments: Accepted for publication in The Astrophysical Journal Letters

Impact of Type II Spicules in the Corona: Simulations and Synthetic Observables

Authors: Juan Martínez-Sykora, Bart De Pontieu, Ineke De Moortel, Viggo H. Hansteen, Mats Carlsson

Abstract: 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 v… ▽ More 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én 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~Å and 193~Å SDO/AIA channels. △ Less

Submitted 16 May, 2018; originally announced May 2018.

Comments: 9 figures, 13 pages, accepted to be published in ApJ

Contribution of phase-mixing of Alfvén waves to coronal heating in multi-harmonic loop oscillations

Authors: P. Pagano, D. J. Pascoe, I. De Moortel

Abstract: 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én waves. In particular, recent observations show that standing oscillations of loops often involve also higher harmonics, beside the fundamental mode. The dampin… ▽ More 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én 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én 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én 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én 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. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Journal ref: A&A 616, A125 (2018)

Transverse Wave Induced Kelvin-Helmholtz Rolls in Spicules

Authors: Patrick Antolin, Don Schmit, Tiago M. D. Pereira, Bart De Pontieu, Ineke De Moortel

Abstract: 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.… ▽ More 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én 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énic 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. △ Less

Submitted 2 March, 2018; originally announced March 2018.

Comments: Accepted for publication in The Astrophysical Journal

Observations and numerical models of solar coronal heating associated with spicules

Authors: Bart De Pontieu, Ineke De Moortel, Juan Martinez-Sykora, Scott McIntosh

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2017; originally announced October 2017.

Comments: 9 pages, 5 figures

Journal ref: The Astrophysical Journal Letters, Volume 845, Issue 2, article id. L18, 7 pp. (2017)

Energetics of the Kelvin-Helmholtz instability induced by transverse waves in twisted coronal loops

Authors: Thomas Howson, Ineke De Moortel, Patrick Antolin

Abstract: 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,… ▽ More 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. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Journal ref: A&A 607, A77 (2017)

A new approach for modelling chromospheric evaporation in response to enhanced coronal heating: II Non-uniform heating

Authors: C. D. Johnston, A. W. Hood, P. J. Cargill, I. De Moortel

Abstract: 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… ▽ More 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. △ Less

Submitted 16 May, 2017; v1 submitted 11 May, 2017; originally announced May 2017.

Comments: 15 pages, 8 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 605, A8 (2017)

Contribution of mode coupling and phase-mixing of Alfvén waves to coronal heating

Authors: Paolo Pagano, Ineke De Moortel

Abstract: Phase-mixing of Alfvén 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 maint… ▽ More Phase-mixing of Alfvén 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énic 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én speed varies through a boundary shell and a footpoint driver is set to trigger kink modes which mode couple to torsional Alfvén modes in the boundary shell. These Alfvén 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. △ Less

Submitted 16 March, 2017; originally announced March 2017.

Journal ref: A&A 601, A107 (2017)

The effects of resistivity and viscosity on the Kelvin-Helmholtz instability in oscillating coronal loops

Authors: Thomas Howson, Ineke De Moortel, Patrick Antolin

Abstract: 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 e… ▽ More 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. △ Less

Submitted 7 March, 2017; originally announced March 2017.

Comments: 13 pages, 12 figures

Journal ref: A&A 602, A74 (2017)

JPEG2000 Image Compression on Solar EUV Images

Authors: C. E. Fischer, D. Müller, I. De Moortel

Abstract: 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. Miss… ▽ More 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. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: 25 pages, published in Solar Physics

Journal ref: Sol Phys (2017) 292: 16

Observational signatures of transverse MHD waves and associated dynamic instabilities

Authors: Patrick Antolin, Ineke De Moortel, Tom Van Doorsselaere, Takaaki Yokoyama

Abstract: 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 trans… ▽ More 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. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 24 pages, 19 figures. Accepted for publication in ApJ. In this ArXiv version, some of the figures have reduced quality due to size limits

Modelling observed decay-less oscillations as resonantly enhanced Kelvin-Helmholtz vortices from transverse MHD waves and their seismological application

Authors: Patrick Antolin, Ineke De Moortel, Tom Van Doorsselaere, Takaaki Yokoyama

Abstract: In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfvén 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… ▽ More In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfvén 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én 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. △ Less

Submitted 30 September, 2016; originally announced September 2016.

Comments: 8 pages, 4 figures; Accepted for publication in The Astrophysical Journal Letters

A new approach for modelling chromospheric evaporation in response to enhanced coronal heating: I The method

Authors: C. D. Johnston, A. W. Hood, P. J. Cargill, I. De Moortel

Abstract: 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 requ… ▽ More 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. △ Less

Submitted 12 May, 2017; v1 submitted 16 September, 2016; originally announced September 2016.

Comments: 15 pages, 11 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 597, A81 (2017)

Observational Signatures of Waves and Flows in the Solar Corona

Authors: Ineke De Moortel, Patrick Antolin, Tom Van Doorsselaere

Abstract: 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 w… ▽ More 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 Å Fe ix and the 193 Å 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 Å Fe ix and the 193 Å Fe xii spectral lines. However, these potential observational signatures need to be treated cautiously because they are likely to be model-dependent. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Journal ref: Solar Physics, Volume 290, Issue 2, pp.399-421 (2015)

Recent advances in coronal heating

Authors: Ineke De Moortel, Philippa Browning

Abstract: 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 sma… ▽ More 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. △ Less

Submitted 4 October, 2015; originally announced October 2015.

Journal ref: Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 373, issue 2042, pp. 20140269-20140269

Transverse, Propagating Velocity Perturbations in Solar Coronal Loops

Authors: I. De Moortel, D. J. Pascoe, A. N. Wright, A. W. Hood

Abstract: 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… ▽ More 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. △ Less

Submitted 4 October, 2015; originally announced October 2015.

On the Parallel and Perpendicular Propagating Motions Visible in Polar Plumes: An Incubator For (Fast) Solar Wind Acceleration?

Authors: Jiajia Liu, Scott W. McIntosh, Ineke De Moortel, Yuming Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 1 July, 2015; originally announced July 2015.

Comments: 15 pages, 4 figures and 1 table

Journal ref: The Astrophysical Journal, 806:273 (7pp), 2015 June 20

Statistical Evidence for the Existence of Alfvénic Turbulence in Solar Coronal Loops

Authors: Jiajia Liu, Scott W. McIntosh, Ineke De Moortel, James Threlfall, Christian Bethge

Abstract: 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énic turbulence may play a key role… ▽ More 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énic 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énic 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énic 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énic turbulence in coronal loops. △ Less

Submitted 18 November, 2014; originally announced November 2014.

Comments: 10 pages, 6 figures, 1 table

Journal ref: ApJ, 797, 7, 2014

The Evolving Magnetic Scales of the Outer Solar Atmosphere and Their Potential Impact on Heliospheric Turbulence

Authors: Scott W. McIntosh, Christian Bethge, James Threlfall, Ineke De Moortel, Robert J. Leamon, Hui Tian

Abstract: 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… ▽ More 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. △ Less

Submitted 11 November, 2013; originally announced November 2013.

Comments: 12 pages, 9 color figures - Submitted to Sp. Sci. Rev

First comparison of wave observations from CoMP and AIA/SDO

Authors: J. Threlfall, I. De Moortel, S. W. McIntosh, C. Bethge

Abstract: 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… ▽ More 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. △ Less

Submitted 20 June, 2013; v1 submitted 14 June, 2013; originally announced June 2013.

Comments: 14 pages, 14 figures, 1 table

The Source of Three-minute Magneto-acoustic Oscillations in Coronal Fans

Authors: D. B. Jess, I. De Moortel, M. Mathioudakis, D. J. Christian, K. P. Reardon, P. H. Keys, F. P. Keenan

Abstract: 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 surroundi… ▽ More 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. △ Less

Submitted 15 August, 2012; originally announced August 2012.

Comments: Accepted into ApJ (13 pages and 10 figures)

A Contemporary View of Coronal Heating

Authors: Clare E. Parnell, Ineke De Moortel

Abstract: 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 he… ▽ More 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. △ Less

Submitted 26 June, 2012; originally announced June 2012.

Journal ref: Phil. Trans. Roy. Soc. A, vol 370, no. 1970, 3217-3240 (2012)

Propagating Disturbances in Coronal Loops: A Detailed Analysis of Propagation Speeds

Authors: G. Kiddie, I. De Moortel, G. Del Zanna, S. W. McIntosh, I. Whittaker

Abstract: 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 a… ▽ More 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Å. △ Less

Submitted 16 May, 2012; v1 submitted 4 May, 2012; originally announced May 2012.

Comments: 26 Papges, 15 Figures

Nonlinear wave propagation and reconnection at magnetic X-points in the Hall MHD regime

Authors: J. W. Threlfall, C. E. Parnell, I. De Moortel, K. G. McClements, T. D. Arber

Abstract: 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 magnet… ▽ More 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. △ Less

Submitted 16 May, 2012; v1 submitted 16 February, 2012; originally announced February 2012.

Comments: 13 pages, 10 figures

MHD Waves and Coronal Seismology: an overview of recent results

Authors: Ineke De Moortel, Valery M. Nakariakov

Abstract: 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 perturbatio… ▽ More 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. △ Less

Submitted 9 February, 2012; originally announced February 2012.

Observed damping of the slow magnetoacoustic mode

Authors: M. S. Marsh, I. De Moortel, R. W. Walsh

Abstract: 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 t… ▽ More 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. △ Less

Submitted 6 April, 2011; originally announced April 2011.

Comments: Accepted for publication in ApJ

Phase Mixing of Nonlinear Visco-resistive Alfvén Waves

Authors: J. A. McLaughlin, I. De Moortel, A. W. Hood

Abstract: We investigate the behaviour of nonlinear, nonideal Alfvén 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én… ▽ More We investigate the behaviour of nonlinear, nonideal Alfvén 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én 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én 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. △ Less

Submitted 31 January, 2011; originally announced January 2011.

Comments: 16 pages, 17 figures