Search

Funding: JR acknowledges the support of the Carnegie Trust for the Universities of Scotland. JR and AWH acknowledge the financial support of STFC through the Consolidated grant, ST/S000402/1, to the University of St Andrews. AWH acknowledges support from ERC Synergy grant ‘The Whole Sun’ (810218). CDJ acknowledges funding from the European Research Council (ERC) under grant agreement No. 647214. A ‘proof of principle’ is presented, whereby the Ohmic and viscous heating determined by a three-dimensional (3D) MHD model of a coronal avalanche are used as the coronal heating input for a series of field-aligned, one-dimensional (1D) hydrodynamic models. Three-dimensional coronal MHD models require large computational resources. For current numerical parameters, it is difficult to model both the magnetic field evolution and the energy transport along field lines for coronal temperatures much hotter than 1MK⁠, because of severe constraints on the time step from parallel thermal conduction. Using the 3D MHD heating derived from a simulation and evaluated on a single field line, the 1D models give coronal temperatures of 1MK and densities 1014--1015m−3 for a coronal loop length of 80Mm⁠. While the temperatures and densities vary smoothly along the field lines, the heating function leads to strong asymmetries in the plasma flows. The magnitudes of the velocities in the 1D model are comparable with those seen in 3D reconnection jets in our earlier work. Advantages and drawbacks of this approach for coronal modelling are discussed. Publisher PDF

Context. It is unclear why there is a rather sharp boundary in sunspots between the umbra and the penumbra. Both regions exhibit magnetoconvection, manifesting in penumbral filaments in the penumbra and in umbral dots in the umbra. Aims. Here we compare the physical properties of umbral dots and penumbral filaments. Our goal is to understand how the properties of these convective features change across the boundary between the umbra and the penumbra and how this is related to the rapid increase in brightness at the umbra-penumbra boundary. Methods. We derived ensemble averages of the physical properties of different types of convective features based on observations of two sunspots with Hinode. Results. There are strong similarities between the convective features in the outer parts of the umbra and the ones in the penumbra, with most physical parameters being smooth and continuous functions of the length of the features. Conclusions. Our results indicate that the transition in brightness from the umbra to the penumbra is solely caused by an increased effectiveness of magnetoconvection within individual convective cells. There is no significant difference in the number density of convective elements between the outer umbra and the inner penumbra. Penumbral filaments exhibit a larger area and a higher brightness compared to umbral dots. It is still unclear, how exactly the underlying magnetic field causes the increase in size and brightness of convective features in the penumbra., 8 pages, 5 figures. Accepted for publication in Astronomy & Astrophysics

We present a detailed analysis of an eruptive event that occurred on early 2019 March 8 in active region AR 12734, to which we refer as the International Women's day event. The event under study is intriguing in several aspects: 1) low-coronal eruptive signatures come in ''pairs'' (a double-peak flare, two coronal dimmings, and two EUV waves); 2) although the event is characterized by a complete chain of eruptive signatures, the corresponding coronagraphic signatures are weak; 3) although the source region of the eruption is located close to the center of the solar disc and the eruption is thus presumably Earth-directed, heliospheric signatures are very weak with little Earth-impact. We analyze a number of multi-spacecraft and multi-instrument (both remote-sensing and in situ) observations, including Soft X-ray, (extreme-) ultraviolet (E)UV), radio and white-light emission, as well as plasma, magnetic field and particle measurements. We employ 3D NLFF modeling to investigate the coronal magnetic field configuration in and around the active region, the GCS model to make a 3D reconstruction of the CME geometry and the 3D MHD numerical model EUHFORIA to model the background state of the heliosphere. Our results indicate two subsequent eruptions of two systems of sheared and twisted magnetic fields, which merge already in the upper corona and start to evolve further out as a single entity. The large-scale magnetic field significantly influences both, the early and the interplanetary evolution of the structure. During the first eruption the stability of the overlying field was disrupted which enabled the second eruption. We find that during the propagation in the interplanetary space the large-scale magnetic field, i.e. , the location of heliospheric current sheet between the AR and the Earth likely influences propagation and the evolution of the erupted structure(s)., 20 pages, 14 figures

In solar physics, a severe numerical challenge for modern simulations is properly representing a transition region between the million-degree hot corona and a much cooler plasma of about 10000 K (e.g., the upper chromosphere or a prominence). In previous 1D hydrodynamic simulations, the transition region adaptive conduction (TRAC) method has been proven to capture aspects better that are related to mass evaporation and energy exchange. We aim to extend this method to fully multidimensional magnetohydrodynamic (MHD) settings, as required for any realistic application in the solar atmosphere. Because modern MHD simulation tools efficiently exploit parallel supercomputers and can handle automated grid refinement, we design strategies for any-dimensional block grid-adaptive MHD simulations. We propose two different strategies and demonstrate their working with our open-source MPI-AMRVAC code. We benchmark both strategies on 2D prominence formation based on the evaporation--condensation scenario, where chromospheric plasma is evaporated through the transition region and then is collected and ultimately condenses in the corona. A field-line-based TRACL method and a block-based TRACB method are introduced and compared in block grid-adaptive 2D MHD simulations. Both methods yield similar results and are shown to satisfactorily correct the underestimated chromospheric evaporation, which comes from a poor spatial resolution in the transition region. Because fully resolving the transition region in multidimensional MHD settings is virtually impossible, TRACB or TRACL methods will be needed in any 2D or 3D simulations involving transition region physics., Accepted by A&A

Context. Prominence eruptions provide key observations to understand the launch of coronal mass ejections as their cold plasma traces a part of the unstable magnetic configuration. Aims. We select a well observed case to derive observational constraints for eruption models. Methods. We analyze the prominence eruption and loop expansion and contraction observed on 02 March 2015 associated with a GOES M3.7 class flare (SOL2015-03-02T15:27) using the data from Atmospheric Imaging Assembly (AIA) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We study the prominence eruption and the evolution of loops using the time-distance techniques. Results. The source region is a decaying bipolar active region where magnetic flux cancellation is present for several days before the eruption. AIA observations locate the erupting prominence within a flux rope viewed along its local axis direction. We identify and quantify the motion of loops in contraction and expansion located on the side of the erupting flux rope. Finally, RHESSI hard X-ray observations identify the loop top and two foot-point sources. Conclusions. Both AIA and RHESSI observations support the standard model of eruptive flares. The contraction occurs 19 minutes after the start of the prominence eruption indicating that this contraction is not associated with the eruption driver. Rather, this prominence eruption is compatible with an unstable flux rope where the contraction and expansion of the lateral loop is the consequence of a side vortex developing after the flux rope is launched., Comment: 16 pages, 9 figures, 1 table, and appendix

Recently, magnetic oscillations were detected in the chromosphere of a large sunspot and found to be linked to the coronal locations where a first ionization potential (FIP) effect was observed. In an attempt to shed light on the possible excitation mechanisms of these localized waves, we further investigate the same data by focusing on the relation between the spatial distribution of the magnetic wave power and the overall field geometry and plasma parameters obtained from multi-height spectropolarimetric non-local thermodynamic equilibrium (NLTE) inversions of IBIS data. We find, in correspondence with the locations where the magnetic wave energy is observed at chromospheric heights, that the magnetic fields have smaller scale heights, meaning faster expansions of the field lines, which ultimately results in stronger vertical density stratification and wave steepening. In addition, the acoustic spectrum of the oscillations at the locations where magnetic perturbations are observed is broader than that observed at other locations, which suggests an additional forcing driver to the p-modes. Analysis of the photospheric oscillations in the sunspot surroundings also reveals a broader spectrum between the two opposite polarities of the active region (the leading spot and the trailing opposite polarity plage), and on the same side where magnetic perturbations are observed in the umbra. We suggest that strong photospheric perturbations between the two polarities are responsible for this broader spectrum of oscillations, with respect to the p-mode spectrum, resulting in locally excited acoustic waves that, after crossing the equipartition layer, located close to the umbra-penumbra boundary at photopheric heights, are converted into magnetic waves and steepen due to the strong density gradient.

Context. The recent discovery of much greater magnetic flux cancellation taking place at the photosphere than previously realised has led us in our previous works to suggest magnetic reconnection driven by flux cancellation as the cause of a wide range of dynamic phenomena, including jets of various kinds and solar atmospheric heating. Aims. Previously, the theory considered energy release at a two-dimensional current sheet. Here we develop the theory further by extending it to an axisymmetric current sheet in three dimensions without resorting to complex variable theory. Methods. We analytically study reconnection and treat the current sheet as a three-dimensional structure. We apply the theory to the cancellation of two fragments of equal but opposite flux that approach each another and are located in an overlying horizontal magnetic field. Results. The energy release occurs in two phases. During Phase 1, a separator is formed and reconnection is driven at it as it rises to a maximum height and then moves back down to the photosphere, heating the plasma and accelerating a plasma jet as it does so. During Phase 2 the fluxes cancel in the photosphere and accelerate a mixture of cool and hot plasma upwards., Comment: 10 pages, 10 figures

The suitability of solar pores as magnetic wave guides has been a key topic of discussion in recent years. Here we present observational evidence of propagating magnetohydrodynamic wave activity in a group of five photospheric solar pores. Employing data obtained by the Facility Infrared Spectropolarimeter at the Dunn Solar Telescope, oscillations with periods on the order of 5 minutes were detected at varying atmospheric heights by examining Si I 10827 \r{A} line bisector velocities. Spectropolarimetric inversions, coupled with the spatially resolved root mean square bisector velocities, allowed the wave energy fluxes to be estimated as a function of atmospheric height for each pore. We find propagating magnetoacoustic sausage mode waves with energy fluxes on the order of 30 kW/m$^{2}$ at an atmospheric height of 100 km, dropping to approximately 2 kW/m$^{2}$ at an atmospheric height of around 500 km. The cross-sectional structuring of the energy fluxes reveals the presence of both body- and surface-mode sausage waves. Examination of the energy flux decay with atmospheric height provides an estimate of the damping length, found to have an average value across all 5 pores of $L_d \approx 268$km, similar to the photospheric density scale height. We find the damping lengths are longer for body mode waves, suggesting that surface mode sausage oscillations are able to more readily dissipate their embedded wave energies. This work verifies the suitability of solar pores to act as efficient conduits when guiding magnetoacoustic wave energy upwards into the outer solar atmosphere., Comment: 20 pages, 9 figures, Accepted for publication in PHILOS T R SOC A

For moderate and slow rotation, magnetic activity of solar-like stars is observed to strongly depend on rotation. These observations do not yet have a solid explanation in terms of dynamo theory. We aim to find such an explanation by numerically investigated the rotational dependency of dynamo drivers in solar-like stars. We ran semi-global convection simulations of stars with rotation rates from 0 to 30 times the solar value, corresponding to Coriolis numbers, Co, of 0 to 110. We measured the turbulent transport coefficients describing the magnetic field evolution with the help of the test-field method, and compared with the dynamo effect arising from the differential rotation. The trace of the $\alpha$ tensor increases for moderate rotation rates with Co$^{0.5}$ and levels off for rapid rotation. This behavior agrees with the kinetic $\alpha$, if one considers the decrease of the convective scale with increasing rotation. The $\alpha$ tensor becomes highly anisotropic for Co$\gtrsim 1$. Furthermore, $\alpha_{rr}$ dominates for moderate rotation (1, Comment: 16 pages with 13 figures, plus 3 pages with 4 figures in the appendix, accepted for publication in A&A

All authors: Rouillard, A. P.; Pinto, R. F.; Vourlidas, A.; De Groof, A.; Thompson, W. T.; Bemporad, A.; Dolei, S.; Indurain, M.; Buchlin, E.; Sasso, C.; Spadaro, D.; Dalmasse, K.; Hirzberger, J.; Zouganelis, I.; Strugarek, A.; Brun, A. S.; Alexandre, M.; Berghmans, D.; Raouafi, N. E.; Wiegelmann, T.; Pagano, P.; Arge, C. N.; Nieves-Chinchilla, T.; Lavarra, M.; Poirier, N.; Amari, T.; Aran, A.; Andretta, V.; Antonucci, E.; Anastasiadis, A.; Auchère, F.; Bellot Rubio, L.; Nicula, B.; Bonnin, X.; Bouchemit, M.; Budnik, E.; Caminade, S.; Cecconi, B.; Carlyle, J.; Cernuda, I.; Davila, J. M.; Etesi, L.; Espinosa Lara, F.; Fedorov, A.; Fineschi, S.; Fludra, A.; Génot, V.; Georgoulis, M. K.; Gilbert, H. R.; Giunta, A.; Gomez-Herrero, R.; Guest, S.; Haberreiter, M.; Hassler, D.; Henney, C. J.; Howard, R. A.; Horbury, T. S.; Janvier, M.; Jones, S. I.; Kozarev, K.; Kraaikamp, E.; Kouloumvakos, A.; Krucker, S.; Lagg, A.; Linker, J.; Lavraud, B.; Louarn, P.; Maksimovic, M.; Maloney, S.; Mann, G.; Masson, A.; Müller, D.; Önel, H.; Osuna, P.; Orozco Suarez, D.; Owen, C. J.; Papaioannou, A.; Pérez-Suárez, D.; Rodriguez-Pacheco, J.; Parenti, S.; Pariat, E.; Peter, H.; Plunkett, S.; Pomoell, J.; Raines, J. M.; Riethmüller, T. L.; Rich, N.; Rodriguez, L.; Romoli, M.; Sanchez, L.; Solanki, S. K.; St Cyr, O. C.; Straus, T.; Susino, R.; Teriaca, L.; del Toro Iniesta, J. C.; Ventura, R.; Verbeeck, C.; Vilmer, N.; Warmuth, A.; Walsh, A. P.; Watson, C.; Williams, D.; Wu, Y.; Zhukov, A. N.-- Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. The Solar Orbiter spacecraft will be equipped with a wide range of remote-sensing (RS) and in situ (IS) instruments to record novel and unprecedented measurements of the solar atmosphere and the inner heliosphere. To take full advantage of these new datasets, tools and techniques must be developed to ease multi-instrument and multi-spacecraft studies. In particular the currently inaccessible low solar corona below two solar radii can only be observed remotely. Furthermore techniques must be used to retrieve coronal plasma properties in time and in three dimensional (3D) space. Solar Orbiter will run complex observation campaigns that provide interesting opportunities to maximise the likelihood of linking IS data to their source region near the Sun. Several RS instruments can be directed to specific targets situated on the solar disk just days before data acquisition. To compare IS and RS, data we must improve our understanding of how heliospheric probes magnetically connect to the solar disk. Aims. The aim of the present paper is to briefly review how the current modelling of the Sun and its atmosphere can support Solar Orbiter science. We describe the results of a community-led effort by European Space Agency's Modelling and Data Analysis Working Group (MADAWG) to develop different models, tools, and techniques deemed necessary to test different theories for the physical processes that may occur in the solar plasma. The focus here is on the large scales and little is described with regards to kinetic processes. To exploit future IS and RS data fully, many techniques have been adapted to model the evolving 3D solar magneto-plasma from the solar interior to the solar wind. A particular focus in the paper is placed on techniques that can estimate how Solar Orbiter will connect magnetically through the complex coronal magnetic fields to various photospheric and coronal features in support of spacecraft operations and future scientific studies. Methods. Recent missions such as STEREO, provided great opportunities for RS, IS, and multi-spacecraft studies. We summarise the achievements and highlight the challenges faced during these investigations, many of which motivated the Solar Orbiter mission. We present the new tools and techniques developed by the MADAWG to support the science operations and the analysis of the data from the many instruments on Solar Orbiter. Results. This article reviews current modelling and tool developments that ease the comparison of model results with RS and IS data made available by current and upcoming missions. It also describes the modelling strategy to support the science operations and subsequent exploitation of Solar Orbiter data in order to maximise the scientific output of the mission. Conclusions. The on-going community effort presented in this paper has provided new models and tools necessary to support mission operations as well as the science exploitation of the Solar Orbiter data. The tools and techniques will no doubt evolve significantly as we refine our procedure and methodology during the first year of operations of this highly promising mission. © 2020 A. P. Rouillard et al., Solar Orbiter is a joint ESA and NASA mission. A. Vourlidas' Solar Orbiter effort is supported by NRL grant N00173-16-1-G029. P. Pagano would like to thank D. H. Mackay and S. L. Yardley for their valuable contributions, the European Research Council (ERC) under the European Union Horizon 2020 research and innovation program (grant agreement No. 647214), and the DiRAC Data Centric system at Durham University, operated by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (http://www.dirac.ac.uk).This equipment was funded by a BIS National E-infrastructure capital grant ST/K00042X/1, STFC capital grant ST/K00087X/1, DiRAC Operations grant ST/K003267/1 and Durham University. DiRAC is part of the National E-Infrastructure. A. Rouillard, V. Genot, M. Janvier, Elie Soubrier, F. Auchere, E. Buchlin and E. Pariat acknowledge support from the French space agency (Centre National d'Etudes Spatiales; CNES; https://cnes.fr/fr) that funds activity in plasma physics data center (Centre de Donnees de la Physique des Plasmas; CDPP; http://cdpp.eu/) and the Multi Experiment Data and Operation Center (MEDOC; https://idoc.ias.u-psud.fr/MEDOC), and the space weather team in Toulouse (Solar-Terrestrial Observations and Modelling Service; STORMS; http://storms-service.irap.omp.eu/).This includes funding for Gaia-DEM, the data mining tools AMDA (http://amda.cdpp.eu/), CLWEB (clweb.cesr. fr/) and the propagation tool (http://propagationtool.cdpp.eu).R.Pinto, M. Lavarra, Y. Wu and A. Kouloumvakos acknowledge financial support from the ANR project SLOW_ SOURCE No. ANR-17-CE31-0006-01, ANR project COROSHOCK No. ANR-18-ERC1-0006-01 and FP7 HELCATS project https://www.helcats-fp7.eu/under the FP7 EU contract number 606692. A. Warmuth acknowledges the support by DLR under grant No. 50 QL 0001. The STEREO SECCHI data are produced by a consortium of RAL (UK), NRL (USA), LMSAL (USA), GSFC (USA), MPS (Germany), CSL (Belgium), IOTA (France) and IAS (France). The ACE data were obtained from the ACE science center. The WIND data were obtained from the Space Physics Data Facility. Javier Rodriguez-Pacheco acknowledges Spanish Project: FEDER/MCIU-AEI/Project ESP2017-88436-R.

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., 14 pages, 16 figures

This paper investigates the changes in spatial properties when magnetohydrodynamic (MHD) waves undergo resonant damping in the Alfvén continuum. The analysis is carried out for a 1D cylindrical pressure-less plasma with a straight magnetic field. The effect of the damping on the spatial wave variables is determined by using complex frequencies that arise as a result of the resonant damping. Compression and vorticity are used to characterise the spatial evolution of the MHD wave. The most striking result is the huge spatial variation in the vorticity component parallel to the magnetic field. Parallel vorticity vanishes in the uniform part of the equilibrium. However, when the MHD wave moves into the non-uniform part, parallel vorticity explodes to values that are orders of magnitude higher than those attained by the transverse components in planes normal to the straight magnetic field. In the non-uniform part of the equilibrium plasma, the MHD wave is controlled by parallel vorticity and resembles an Alfvén wave, with the unfamiliar property that it has pressure variations even in the linear regime.

A single case study of a CBP in an equatorial coronal hole with an exceptionally large size is investigated to extend our understanding of the formation of mini-filaments, their destabilisation and the origin of the eruption triggering the formation of jet-like features recorded in the extreme-ultraviolet (EUV) and X-ray emission. We aim to explore the nature of the so-called micro-flares in CBPs associated with jets in coronal holes and mini coronal mass ejections in the quiet Sun. Co-observations from the Atmospheric Imaging Assembly (AIA) and Helioseismic Magnetic Imager (HMI) on board the Solar Dynamics Observatory, and GONG Halpha images are used together with a Non-Linear Force Free Field (NLFFF) relaxation approach, where the latter is based on a time series of HMI line-of-sight magnetograms. A mini-filament (MF) that formed beneath the CBP arcade around 3-4 h before the eruption is seen in the Halpha and EUV AIA images to lift up and erupt triggering the formation of an X-ray jet. No significant photospheric magnetic flux concentration displacement (convergence) is observed and neither is magnetic flux cancellation between the two main magnetic polarities forming the CBP in the time period leading to the MF liftoff. The CBP micro-flare is associated with three flare kernels that formed shortly after the MF liftoff. No observational signature is found for reconnection beneath the erupting MF. The applied NLFFF modelling successfully reproduces both the CBP loop complex as well as the magnetic flux rope that hosts the MF during the build-up to the eruption., 13 pages, 11 figures, the abstract has been reduced from its original length, accepted by A&A

The large field-of-view of the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) instrument on board the PRoject for Onboard Autonomy 2 (PROBA2) spacecraft provides a unique opportunity to study extended coronal structures observed in EUV in conjunction with global coronal magnetic field simulations. A global non-potential magnetic field model is used to simulate the evolution of the global corona from 1 September 2014 to 31 March 2015, driven by newly emerging bipolar active regions determined from Helioseismic and Magnetic Imager (HMI) magnetograms. We compare the large-scale structure of the simulated magnetic field with structures seen off-limb in SWAP EUV observations. In particular, we investigate how successful the model is in reproducing regions of closed and open structures; the scale of structures; and compare the evolution of a coronal fan observed over several rotations. The model is found to accurately reproduce observed large-scale off-limb structures. When discrepancies do arise they mainly occur off the east solar limb due to active regions emerging on the far side of the Sun, which cannot be incorporated into the model until they are observed on the Earth-facing side. When such ``late'' active region emergences are incorporated into the model, we find that the simulated corona self-corrects within a few days, so that simulated structures off the west limb more closely match what is observed. Where the model is less successful, we consider how this may be addressed, through model developments or additional observational products., Comment: 26 pages, 11 figures

The study of lithium (Li) surface abundance in the Sun and young stellar globular clusters which are seemingly anomalous in present-day scenarios, as well as the influence of rotation and magnetic braking (MB) on its depletion during pre-main sequence (PMS) and main sequence (MS). In this work, the effects of rotational mixing and of the rotational hydrostatic effects on Li abundances are studied by simulating several grids of PMS and MS rotating and non-rotating models. Those effects are combined with the additional impact of the MB (with magnetic field intensities ranging between 3.0 and 5.0 G). The data obtained from simulations are confronted by comparing different stellar parameters. The results show that the surface Li abundance for the Sun-like models at the end of the PMS and throughout the MS decreases when rotational effects are included, that is the Li depletion rate for rotating models is higher than for non-rotating ones. This effect is attenuated when the MB produced by a magnetic field is present. This physical phenomena impacts also the star effective temperature (T-eff) and its location in the HR diagram. The impact of MB in Li depletion is sensitive to the magnetic field intensity: the higher it is, the lower the Li destruction. A direct link between the magnetic fields and the convective one (CZ) size is observed: stronger magnetic fields produce shallower CZ's. This result suggests that MB effect must be taken into consideration during PMS if we aim to reproduce Li abundances in young clusters. © 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, We are pleased to acknowledge Jieun Choi her kindness in responding to the different questions about her work which were aimed to reproduce the results obtained as far as Li is concerned, and also for allowing us access to the MESA files she used. Also we very much appreciate the expert support provided by Elisa Delgado-Mena during the process review of this document. Similarly, we wish also to thank Matteo Cantiello and Bill Paxton for their very valuable support in the development of the MB routine in MESA. Finally, we wish to acknowledge the generous work and help offered by the MESA community. The authors acknowledge funding support from Spanish public funds (including FEDER funds) for research under projects ESP2017-87676-C5-2-R and ESP2017-87676-C5-5-R. JCS also acknowledges support from project RYC-2012-09913 under the 'Ramon y Cajal' program of the "Ministerio de Ciencia e Innovacion".

Context. Direct magnetic field measurements in solar prominences occur infrequently and are difficult to make and interpret. As a consequence, alternative methods are needed to derive the main properties of the magnetic field that supports the prominence mass. This is important for our understanding of solar prominences, but also for understanding how eruptive prominences may affect space weather. Aims. We present the first direct comparison of the magnetic field strength derived from spectro-polarimetric observations of a solar prominence, with corresponding results from a theoretical flux rope model constructed from on-disc normal component magnetograms. Methods. We first used spectro-polarimetric observations of a prominence obtained with the magnetograph THEMIS operating in the Canary Islands to derive the magnetic field of the observed prominence by inverting the Stokes parameters measured in the He D3 line. Next, we constructed two data-constrained non-linear force-free field (NLFFF) models of the same prominence. In one model we assumed a strongly twisted flux rope solution, and in the other a weakly twisted flux rope solution. Results. The physical extent of the prominence at the limb (height and length) is best reproduced with the strongly twisted flux rope solution. The line-of-sight average of the magnetic field for the strongly twisted solution results in a magnetic field that has a magnitude of within a factor of 1−2 of the observed magnetic field strength. For the peak field strength along the line of sight, an agreement to within 20% of the observations is obtained for the strongly twisted solution. The weakly twisted solution produces significantly lower magnetic field strengths and gives a poor agreement with the observations. Conclusions. The results of this first comparison are promising. We found that the flux rope insertion method of producing a NLFFF is able to deduce the overall properties of the magnetic field in an observed prominence.

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., 14 pages, 19 figures

Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging) and the stellar radius. We use a sample of 51 stars across a mass range of 0.1-1.34 M$_\odot$ to show how the helicity density relates to stellar mass, Rossby number, magnetic energy and age. We find that the large-scale helicity density increases with decreasing Rossby number $R_o$, peaking at $R_o \simeq 0.1$, with a saturation or decrease below that. For both fully- and partially-convective stars we find that the mean absolute helicity density scales with the mean squared toroidal magnetic flux density according to the power law: $|\langle{h\,}\rangle|$ $\propto$ $\langle{\rm{B_{tor}}^2_{}\,\rangle}^{0.86\,\pm\,0.04}$. The scatter in this relation is consistent with the variation across a solar cycle, which we compute using simulations and observations across solar cycles 23 and 24 respectively. We find a significant decrease in helicity density with age., 11 pages, 7 figures

International audience; The magnetometer instrument on the Solar Orbiter mission is designed to measure the magnetic field local to the spacecraft continuously for the entire mission duration. The need to characterise not only the background magnetic field but also its variations on scales from far above to well below the proton gyroscale result in challenging requirements on stability, precision, and noise, as well as magnetic and operational limitations on both the spacecraft and other instruments. The challenging vibration and thermal environment has led to significant development of the mechanical sensor design. The overall instrument design, performance, data products, and operational strategy are described.

Context: Hot coronal jets are a basic observed feature of the solar atmosphere whose physical origin is still being actively debated. Aims: We study six recurrent jets occurring in the active region NOAA 12644 on April 04, 2017. They are observed in all the hot filters of AIA as well as cool surges in IRIS slit-jaw high spatial and temporal resolution images. Methods: The AIA filters allow us to study the temperature and the emission measure of the jets using the filter ratio method. We study the pre-jet phases by analyzing the intensity oscillations at the base of the jets with the wavelet technique. Results: A fine co-alignment of the AIA and IRIS data shows that the jets are initiated at the top of a canopy-like, double-chambered structure with cool emission on one side and hot emission in the other. The hot jets are collimated in the hot temperature filters, have high velocities (around 250 km/s) and accompanied by the cool surges and ejected kernels both moving at about 45 km/s. In the pre-phase of the jets, at their base we find quasi-periodic intensity oscillations in phase with small ejections; they have a period between 2 and 6 minutes and are reminiscent of acoustic or MHD waves. Conclusions: This series of jets and surges provides a good case-study to test the 2D and 3D magnetohydrodynamic (MHD) models that result from magnetic flux emergence. The double-chambered structure found in the observations corresponds to the cold and hot loop regions found in the models beneath the current sheet that contains the reconnection site. The cool surge with kernels is comparable with the cool ejection and plasmoids that naturally appear in the models., Comment: 12 Pages, 11 Figures

Context. Penumbral microjets (PMJs) is the name given to elongated jet-like brightenings observed in the chromosphere above sunspot penumbrae. They are transient events that last from a few seconds to several minutes, and their origin is presumed to be related to magnetic reconnection processes. Previous studies have mainly focused on their morphological and spectral characteristics, and more recently on their spectropolarimetric signals during the maximum brightness stage. Studies addressing the temporal evolution of PMJs have also been carried out, but they are based on spatial and spectral time variations only. Aims. Here we investigate, for the first time, the temporal evolution of the polarization signals produced by short-lived PMJs (lifetimes < 2 min) to infer how the magnetic field vector evolves in the upper photosphere and mid-chromosphere. Methods. We use fast-cadence spectropolarimetric observations of the Ca II 854.2 nm line taken with the CRisp Imaging Spectropolarimeter at the Swedish 1 m Solar Telescope. The weak-field approximation (WFA) is used to estimate the strength and inclination of the magnetic field vector. By separating the Ca II 854.2 nm line into two different wavelength domains to account for the chromospheric origin of the line core and the photospheric contribution to the wings, we infer the height variation of the magnetic field vector. Results. The WFA reveals larger magnetic field changes in the upper photosphere than in the chromosphere during the PMJ maximum brightness stage. In the photosphere, the magnetic field inclination and strength undergo a transient increase for most PMJs, but in 25% of the cases the field strength decreases during the brightening. In the chromosphere, the magnetic field tends to be slightly stronger during the PMJs. Conclusions. The propagation of compressive perturbation fronts followed by a rarefaction phase in the aftershock region may explain the observed behavior of the magnetic field vector. The fact that such behavior varies among the analyzed PMJs could be a consequence of the limited temporal resolution of the observations and the fast-evolving nature of the PMJs. © 2020 ESO., This work has been funded by the Spanish Ministry of Science and Innovation through project RTI2018-096886-B-C51, including a percentage from FEDER funds, and through the Centro de Excelencia Severo Ochoa grant SEV-2017-0709 awarded to the Instituto de Astrofísica de Andalucía in the period 2018–2022. The Swedish 1 m Solar Telescope is operated on the island of La Palma by the Institute for Solar Physics of Stockholm University in the Spanish Observatory del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. The Institute for Solar Physics is supported by a grant for research infrastructures of national importance from the Swedish Research Council (registration number 2017-00625).

Coronal plasma in the cores of solar active regions is impulsively heated to more than 5 MK. The nature and location of the magnetic energy source responsible for such impulsive heating is poorly understood. Using observations of seven active regions from the Solar Dynamics Observatory, we found that a majority of coronal loops hosting hot plasma have at least one footpoint rooted in regions of interacting mixed magnetic polarity at the solar surface. In cases when co-temporal observations from the Interface Region Imaging Spectrograph space mission are available, we found spectroscopic evidence for magnetic reconnection at the base of the hot coronal loops. Our analysis suggests that interactions of magnetic patches of opposite polarity at the solar surface and the associated energy release during reconnection are key to impulsive coronal heating., Comment: Published in Astronomy & Astrophysics. Online movies available at https://www.aanda.org/articles/aa/olm/2020/12/aa39099-20/aa39099-20.html

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., Comment: 15 pages and 22 figures

The study of the evolution of coronal holes (CHs) is especially important in the context of high-speed solar wind streams (HSS) emanating from them. Stream interaction regions may deliver large amount of energy into the Earths system, cause geomagnetic storms, and shape interplanetary space. By statistically analysing 16 long-living CHs observed by the SDO, we focus on coronal, morphological and underlying photospheric magnetic field characteristics as well as investigate the evolution of the associated HSSs. We use CATCH to extract and analyse CHs using observations taken by AIA and HMI. We derive changes in the CH properties and correlate them to the CH evolution. Further we analyse the properties of the HSS signatures near 1au from OMNI data by manually extracting the peak bulk velocity of the solar wind plasma. We find that the area evolution of CHs mostly shows a rough trend of growing to a maximum followed by a decay. No correlation of the area evolution to the evolution of the signed magnetic flux and signed magnetic flux density enclosed in the projected CH area was found. From this we conclude that the magnetic flux within the extracted CH boundaries is not the main cause for its area evolution. We derive CH area change rates (growth and decay) of 14.2 +/- 15.0 * 10^8 km^2/day showing a reasonable anti-correlation (cc =-0.48) to the solar activity, approximated by the sunspot number. The change rates of the signed mean magnetic flux density (27.3 +/- 32.2 mG/day) and the signed magnetic flux (30.3 +/- 31.5 * 10^18 Mx/day) were also found to be dependent on solar activity (cc =0.50 and cc =0.69 respectively) rather than on the individual CH evolutions. Further we find that the CH area-to-HSS peak velocity relation is valid for each CH over its evolution but revealing significant variations in the slopes of the regression lines., Accepted at A&A

Context. Coronal loops are the basic building blocks of the solar corona, which are related to the mass supply and heating of solar plasmas in the corona. However, their fundamental magnetic structures are still not well understood. Most coronal loops do not expand significantly, whereas the diverging magnetic field would have an expansion factor of about 5-10 over one pressure scale height. Aims. In this study, we investigate a unique coronal loop with a roughly constant cross section, it is ultra long and quite thin. A coronal loop model with magnetic helicity is presented to explain the small expansion of the loop width. Methods. This coronal loop was predominantly detectable in the 171 A channel of the Atmospheric Imaging Assembly (AIA). Then, the local magnetic field line was extrapolated by a Potential-Field-Source-Surface model. Finally, the differential emission measure analysis made from six AIA bandpasses was applied to obtain the thermal properties of this loop. Results. This coronal loop has a projected length of roughly 130 Mm, a width of about 1.5 +(-) 0.5 Mm and a lifetime of around 90 minutes. It follows an open magnetic field line. The cross section expanded very little (i.e., 1.5-2.0) along the loop length during its whole lifetime. This loop has a nearly constant temperature at about 0.7 +(-) 0.2 MK, whereas its density exhibits the typical structure of a stratified atmosphere. Conclusions. We use a thin twisted flux tube theory to construct a model for this non-expanding loop, and find that indeed with sufficient twist a coronal loop can attain equilibrium. However, we can not rule out other possibilities such as footpoint heating by small-scale reconnection, elevated scale height by a steady flow along the loop etc., Comment: 8 pages, 5 figures, 1 tables

The Kelvin-Helmholtz instability has been proposed as a mechanism to extract energy from magnetohydrodynamic (MHD) kink waves in flux tubes, and to drive dissipation of this wave energy through turbulence. It is therefore a potentially important process in heating the solar corona. However, it is unclear how the instability is influenced by the oscillatory shear flow associated with an MHD wave. We investigate the linear stability of a discontinuous oscillatory shear flow in the presence of a horizontal magnetic field within a Cartesian framework that captures the essential features of MHD oscillations in flux tubes. We derive a Mathieu equation for the Lagrangian displacement of the interface and analyse its properties, identifying two different instabilities: a Kelvin-Helmholtz instability and a parametric instability involving resonance between the oscillatory shear flow and two surface Alfv\'{e}n waves. The latter occurs when the system is Kelvin-Helmholtz stable, thus favouring modes that vary along the flux tube, and as a consequence provides an important and additional mechanism to extract energy. When applied to flows with the characteristic properties of kink waves in the solar corona, both instabilities can grow, with the parametric instability capable of generating smaller scale disturbances along the magnetic field than possible via the Kelvin-Helmholtz instability. The characteristic time-scale for these instabilities is $\sim 100$ s, for wavelengths of $200$ km. The parametric instability is more likely to occur for smaller density contrasts and larger velocity shears, making its development more likely on coronal loops than on prominence threads., Comment: 12 pages, 4 figures, accepted for publication MNRAS

Funding: Carnegie Trust for the Universities of Scotland; Science and Technology Facilities Council (grants ST/N000609/1 and ST/P000320/1). The plasma heating associated with an avalanche involving three twisted magnetic threads within a coronal loop is investigated using three-dimensional magnetohydrodynamic simulations. The avalanche is triggered by the kink instability of one thread, with the others being engulfed as a consequence. The heating as a function of both time and location along the strands is evaluated. It is shown to be bursty at all times but to have no preferred spatial location. While there appears to be a level of "background" heating, this is shown to be comprised of individual, small heating events. A comparison between viscous and resistive (Ohmic) heating demonstrates that the strongest heating events are largely associated with the Ohmic heating that arises when the current exceeds a critical value. Viscous heating is largely (but not entirely) associated with smaller events. Ohmic heating dominates viscous heating only at the time of the initial kink instability. It is also demonstrated that a variety of viscous models lead to similar heating rates, suggesting that the system adjusts to dissipate the same amount of energy. Postprint

The magnetic network extending from the photosphere (solar radius $\simeq R_\odot$) to lower corona ($R_\odot + 10$ Mm) plays an important role in the heating mechanisms of the solar atmosphere. Here we further develop the models with realistic open magnetic flux tubes of the authors in order to model more complicated configurations. Closed magnetic loops, and combinations of closed and open magnetic flux tubes are modelled. These are embedded within a stratified atmosphere, derived from observationally motivated semi-empirical and data-driven models subject to solar gravity and capable of spanning from the photosphere up into the chromosphere and lower corona. Constructing a magnetic field comprising self-similar magnetic flux tubes, an analytic solution for the kinetic pressure and plasma density is derived. Combining flux tubes of opposite polarity it is possible to create a steady background magnetic field configuration modelling a solar atmosphere exhibiting realistic stratification. The result can be applied to SOHO/MDI and SDO/HMI and other magnetograms from the solar surface, upon which photospheric motions can be simulated to explore the mechanism of energy transport. We demonstrate this powerful and versatile method with an application to Helioseismic and Magnetic Imager data., 7 pages, 5 figures

Coronal mass ejections (CMEs) that exhibit weak or no eruption signatures in the low corona, known as stealth CMEs, are problematic as upon arrival at Earth they can lead to geomagnetic disturbances that were not predicted by space weather forecasters. We investigate the origin and eruption of a stealth event that occurred on 2015 January 3 that was responsible for a strong geomagnetic storm upon its arrival at Earth. To simulate the coronal magnetic field and plasma parameters of the eruption we use a coupled approach. This approach combines an evolutionary nonlinear force-free field model of the global corona with a MHD simulation. The combined simulation approach accurately reproduces the stealth event and suggests that sympathetic eruptions occur. In the combined simulation we found that three flux ropes form and then erupt. The first two flux ropes, which are connected to a large AR complex behind the east limb, erupt first producing two near-simultaneous CMEs. These CMEs are closely followed by a third, weaker flux rope eruption in the simulation that originated between the periphery of AR 12252 and the southern polar coronal hole. The third eruption coincides with a faint coronal dimming, which appears in the SDO/AIA 211 A observations, that is attributed as the source responsible for the stealth event and later the geomagnetic disturbance at 1 AU. The incorrect interpretation of the stealth event being linked to the occurrence of a single partial halo CME observed by LASCO/C2 is mainly due to the lack of STEREO observations being available at the time of the CMEs. The simulation also shows that the LASCO CME is not a single event but rather two near-simultaneous CMEs. These results show the significance of the coupled data-driven simulation approach in interpreting the eruption and that an operational L5 mission is crucial for space weather forecasting., Comment: Accepted in A&A, 17 pages, 12 figures

We study the properties of slow magneto-acoustic waves that are naturally excited due to turbulent convection and investigate their role in the energy balance of a plage region using three dimensional (3D) radiation-MHD simulations. We calculate the horizontally averaged (over the whole domain) frequency power spectra for both longitudinal and vertical (i.e. the component perpendicular to the surface) components of velocity. To compare our results with the observations we degrade the simulation data with Gaussian kernels having FWHM of 100 km and 200 km, and calculate horizontally averaged power spectra for the vertical component of velocity. The power spectra of the longitudinal component of velocity, averaged over field lines in the core of a kG magnetic flux concentration, reveal that the dominant period of oscillations shifts from around 6.5 minutes in the photosphere to around 4 minutes in the chromosphere. At the same time, the velocity power spectra, averaged horizontally over the whole domain, show that low frequency waves (approximately 6.5 minute period) may reach well into the chromosphere. Importantly, waves with frequencies above 5 mHz propagating along different field lines are found to be out of phase with each other even within a single magnetic concentration. The horizontally averaged power spectra of the vertical component of velocity at various effective resolutions show that the observed acoustic wave energy fluxes are underestimated, by a factor of three even if determined from observations carried out at a high spatial resolution of 200 km. Our results show that longitudinal waves carry (just) sufficient energy to heat the chromosphere in solar plage. We conjecture that current observations (with spatial resolution around 200 km) underestimate the energy flux by roughly a factor of three, or more if the observations have lower spatial resolution., Accepted for publication in A & A

Coronal mass ejections (CMEs), which are among the most magnificent solar eruptions, are a major driver of space weather and can thus affect diverse human technologies. Different processes have been proposed to explain the initiation and release of CMEs from solar active regions (ARs), without reaching consensus on which is the predominant scenario, and thus rendering impossible to accurately predict when a CME is going to erupt from a given AR. To investigate AR magnetic properties that favor CMEs production, we employ multi-spacecraft data to analyze a long duration AR (NOAA 11089, 11100, 11106, 11112 and 11121) throughout its complete lifetime, spanning five Carrington rotations from July to November 2010. We use data from the Solar Dynamics Observatory to study the evolution of the AR magnetic properties during the five near-side passages, and a proxy to follow the magnetic flux changes when no magnetograms are available, i.e. during far-side transits. We characterized the angular widths, speeds and masses of 108 CMEs that we associated to the AR, when examining a 124-day period. Such an ejectivity tracking was possible thanks to the multiviewpoint images provided by the STEREO and SOHO in a quasiquadrature configuration. We also inspected the X-ray flares registered by the GOES satellite and found 162 to be associated to the AR under study. Given the substantial number of ejections studied, we use a statistical approach instead of a single-event analysis. We found three well defined periods of very high CMEs activity and two periods with no mass ejections that are preceded or accompanied by characteristic changes in the AR magnetic flux, free magnetic energy and/or presence of electric currents. Our large sample of CMEs and long term study of a single AR, provide further evidence relating AR magnetic activity to CME and Flare production., Accepted for publication in Advances in Space Research

Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $\delta$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted flux tube from the solar interior into a non-magnetized stratified atmosphere. Results. We emerge two $\Omega$-shaped loops forming a quadrupolar region. The emerging flux forms two initially separated bipoles, that later come in contact creating a $\delta$-spot central region. Above the two bipoles, two magnetic lobes expand and interact through a series of current sheets at the interface between them. Two recurrent confined eruptions are produced. In both cases, reconnection between sheared low-lying field lines forms a flux rope. Reconnection between the two lobes high in the atmosphere forms retracting field lines that push against the flux rope, creating a current sheet between them. It also forms a third magnetic lobe between the two emerged ones, that later acts as a strapping field. The flux rope eruptions are triggered when the reconnection between the flux ropes and the field above them becomes efficient enough to remove the tension of the overlying field. These reconnection events occur internally in the system. The first erupting flux rope almost fully reconnects with the overlying field. The second eruption is confined by the overlying field. During the confined eruption, the flux rope is enhanced in size, flux and twist, similar to confined-flare-to-flux-rope observations. Proxies of the emission reveal the two erupting filaments channels. A flare arcade is formed only in the second eruption due to the longer-lasting and more efficient reconnection below the flux rope., Comment: 15 pages, 19 figures, to appear in Astronomy and Astrophysics

K.A.M. and D.H.M. gratefully acknowledge the support of the Leverhulme Trust, the STFC and the Carnegie Trust for the Universities of Scotland. K.A.M., A.S.S., and E.E.D. thank ISSI (Bern) for their support of the team "Understanding Solar Jets and their Role in Atmospheric Structure and Dynamics." D.H.M. thanks the EU for their support under FP7. Coronal jets occur frequently on the Sun, and may contribute significantly to the solar wind. With the suite of instruments available now, we can observe these phenomena in greater detail than ever before. Modeling and simulations can assist further with understanding the dynamic processes involved, but previous studies tended to consider only one mechanism (e.g., emergence or rotation) for the origin of the jet. In this study we model a series of idealized archetypal jet configurations and follow the evolution of the coronal magnetic field. This is a step toward understanding these idealized situations before considering their observational counterparts. Several simple situations are set up for the evolution of the photospheric magnetic field: a single parasitic polarity rotating or moving in a circular path; as well as opposite polarity pairs involved in flyby (shearing), cancellation or emergence; all in the presence of a uniform, open background magnetic field. The coronal magnetic field is evolved in time using a magnetofrictional relaxation method. While magnetofriction cannot accurately reproduce the dynamics of an eruptive phase, the structure of the coronal magnetic field, as well as the buildup of electric currents and free magnetic energy are instructive. Certain configurations and motions produce a flux rope and allow the significant buildup of free energy, reminiscent of the progenitors of so-called blowout jets, whereas other, simpler configurations are more comparable to the standard jet model. The next stage is a comparison with observed coronal jet structures and their corresponding photospheric evolution. Postprint Publisher PDF

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., 15 pages, 10 figures, 1 table, accepted for publication in A&A

This paper investigates the mixed properties of MHD waves in a non-uniform plasma. It starts with a short revision of MHD waves in a uniform plasma of infinite extent. In that case the MHD waves do not have mixed properties. They can be separated in Alfv\'{e}n waves and magneto-sonic waves. The Alfv\'{e}n waves propagate parallel vorticity and are incompressible. In addition they have no parallel displacement component. The magneto-sonic waves are compressible and in general do have a parallel component of displacement but do not propagate parallel vorticity. This clear separation has been the reason why there has been a strong inclination in the literature to use this classification in the study of MHD waves in non-uniform plasmas. The main part of this paper is concerned with MHD waves in a non-uniform plasma. It is shown that the MHD waves in that situation in general propagate both vorticity and compression and hence have mixed properties. Finally, the close connection between resonant absorption and MHD waves with mixed properties is discussed.

AH acknowledges funding from the Science and Technology Facilities Council (UK) through the consolidated grant ST/N000609/1. DP and TVD were supported by the GOA-2015-014 (KU Leuven) and the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 724326). The transverse structure of coronal loops plays a key role in the physics but the small transverse scales can be difficult to observe directly. For wider loops the density profile may be estimated by forward modelling of the transverse intensity profile. The transverse density profile may also be estimated seismologically using kink oscillations in coronal loops. The strong damping of kink oscillations is attributed to resonant absorption and the damping profile contains information about the transverse structure of the loop. However, the analytical descriptions for damping by resonant absorption presently only describe the behaviour for thin inhomogeneous layers. Previous numerical studies have demonstrated that this thin boundary approximation produces poor estimates of the damping behaviour in loops with wider inhomogeneous layers. Both the seismological and forward modelling approaches suggest loops have a range of layer widths and so there is a need for a description of the damping behaviour that accurately describes such loops. We perform a parametric study of the damping of standing kink oscillations by resonant absorption for a wide range of inhomogeneous layer widths and density contrast ratios, with a focus on the values most relevant to observational cases. We describe the damping profile produced by our numerical simulations without prior assumption of its shape and compile our results into a lookup table which may be used to produce accurate seismological estimates for kink oscillation observations. Publisher PDF

This research has received funding from the Science and Technology Facilities Council (UK) through the consolidated grant ST/N000609/1 and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 647214). Aims . This paper presents 2.5D numerical experiments of Alfvén wave phase mixing and aims to assess the effects of nonlinearities on the wave behaviour and dissipation. In addition, this paper aims to quantify how effective the model presented here 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 outusing the Lare2D code (Arber et al. 2001). 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 . It was found that the wave trains phase-mix due to variations in the length of each field line as well asvariations 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 increased 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, it was shown that the Poynting flux in our experiments was comparable to the energy requirements of the Quiet Sun corona, although a (possibly unphysically) large amount of magnetic diffusion was used, and it remains unclear if the model is able to provide enough energy under actual coronal conditions. Postprint

We report on our three-dimensional (3D) magnetohydrodynamic (MHD) simulations of cylindrical weakly twisted flux tubes emerging from 18 Mm below the photosphere. We perform a parametric study, by varying the initial magnetic field strength ($B_0$), radius ($R$), twist ($\alpha)$ and length of the emerging part of the flux tube ($\lambda$) to investigate how these parameters affect the transfer of the magnetic field from the convection zone to the photosphere. We show that the efficiency of emergence at the photosphere (i.e. how strong the photospheric field will be in comparison to $B_0$) depends not only on the $B_0$ but also the morphology of the emerging field and the twist. We show that parameters such as $B_0$ and magnetic flux cannot alone determine whether a flux tube will emerge to the solar surface. For instance, high-$B_0$ (weak-$B_0$) fields may fail (succeed) to emerge at the photosphere, depending on their geometrical properties. We also show that the photospheric magnetic field strength can vary greatly for flux tubes with the same $B_0$ but different geometric properties. Moreover, in some cases we have found scaling laws, whereby the magnetic field strength scales with the local density as $B\propto \rho^\kappa$, where $\kappa \approx 1$ deeper in the convection zone and $\kappa, Comment: 22 pages, 10 pages

Quiescent solar prominences are observed to exist within the solar atmosphere for up to several solar rotations. Their eruption is commonly preceded by a slow increase in height that can last from hours to days. This increase in the prominence height is believed to be due to their host magnetic flux rope transitioning through a series of neighbouring quasi-equilibria before the main loss-of-equilibrium that drives the eruption. Recent work suggests that the removal of prominence mass from a stable, quiescent flux rope is one possible cause for this change in height. However, these conclusions are drawn from observations and are subject to interpretation. Here we present a simple model to quantify the effect of "mass-draining" during the pre-eruptive height-evolution of a solar flux rope. The flux rope is modeled as a line current suspended within a background potential magnetic field. We first show that the inclusion of mass, up to $10^{12}$~kg, can modify the height at which the line current experiences loss-of-equilibrium by up to 14\%. Next, we show that the rapid removal of mass prior to the loss-of-equilibrium can allow the height of the flux rope to increase sharply and without upper bound as it approaches its loss-of-equilibrium point. This indicates that the critical height for the loss-of-equilibrium can occur at a range of heights depending explicitly on the amount and evolution of mass within the flux rope. Finally, we demonstrate that for the same amount of drained mass, the effect on the height of the flux rope is up to two order of magnitude larger for quiescent than for active region prominences., 13 pages, 6 figures, Accepted for publication in the Astrophysical Journal

© 2019. The American Astronomical Society. All rights reserved. Recent numerical simulations have demonstrated that transverse coronal loop oscillations are susceptible to the Kelvin-Helmholtz (KH) instability due to the counterstreaming motions at the loop boundary. We present the first analytical model of this phenomenon. The region at the loop boundary where the shearing motions are greatest is treated as a straight interface separating time-periodic counterstreaming flows. In order to consider a twisted tube, the magnetic field at one side of the interface is inclined. We show that the evolution of the displacement at the interface is governed by Mathieu's equation, and we use this equation to study the stability of the interface. We prove that the interface is always unstable and that, under certain conditions, the magnetic shear may reduce the instability growth rate. The result, that the magnetic shear cannot stabilize the interface, explains the numerically found fact that the magnetic twist does not prevent the onset of the KH instability at the boundary of an oscillating magnetic tube. We also introduce the notion of the loop σ-stability. We say that a transversally oscillating loop is σ-stable if the KH instability growth time is larger than the damping time of the kink oscillation. We show that even relatively weakly twisted loops are σ-stable. ispartof: ASTROPHYSICAL JOURNAL vol:870 issue:2 status: published

The combination of magnetohydrodynamic (MHD) simulation and multi-wavelength observations is an effective way to study mechanisms of magnetic flux rope eruption. We develop a data-driven MHD model using the zero-$\beta$ approximation. The initial condition is provided by nonlinear force-free field derived by the magneto-frictional method based on vector magnetic field observed by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO). The bottom boundary uses observed time series of the vector magnetic field and the vector velocity derived by the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM). We apply the data-driven model to active region 11123 observed from 06:00 UT on 2011 November 11 to about 2 hours later. The evolution of the magnetic field topology coincides with the flare ribbons observed in the 304 and 1600 \AA $ $ wavebands by the Atmospheric Imaging Assembly. The morphology, propagation path, and propagation range of the flux rope are comparable with the observations in 304 \AA . We also find that a data-constrained boundary condition, where the bottom boundary is fixed to the initial values, reproduces a similar simulation result. This model can reproduce the evolution of a magnetic flux rope in its dynamic eruptive phase., Comment: 17 pages, 5 figures, accepted by ApJL

Recent observations at high spatial resolution have shown that magnetic flux cancellation occurs on the solar surface much more frequently than previously thought, and so this led Priest et al 2018 to propose magnetic reconnection driven by photospheric flux cancellation as a mechanism for chromospheric and coronal heating. In particular, they estimated analytically the amount of energy released as heat and the height of the energy release during flux cancellation. In the present work, we take the next step in the theory by setting up a two-dimensional resistive MHD simulation of two canceling polarities in the presence of a horizontal external field and a stratified atmosphere in order to check and improve upon the analytical estimates. Computational evaluation of the energy release during reconnection is found to be in good qualitative agreement with the analytical estimates. In addition, we go further and undertake an initial study of the atmospheric response to reconnection. We find that, during the cancellation, either hot ejections or cool ones or a combination of both hot and cool ejections can be formed, depending on the height of the reconnection location. The hot structures can have the density and temperature of coronal loops, while the cooler structures are suggestive of surges and large spicules., Comment: 19 pages, 13 figures

Context. The magnetic polarities of bipolar active regions (ARs) exhibit elongations in line-of-sight magnetograms during their emergence. These elongations are referred to as magnetic tongues and attributed to the presence of twist in the emerging magnetic flux-ropes (FRs) that form ARs. Aims. The presence of magnetic tongues affects the measurement of any AR characteristic that depends on its magnetic flux distribution. The AR tilt-angle is one of them. We aim to develop a method to isolate and remove the flux associated with the tongues to determine the AR tilt-angle with as much precision as possible. Methods. As a first approach, we used a simple emergence model of a FR. This allowed us to develop and test our aim based on a method to remove the effects of magnetic tongues. Then, using the experience gained from the analysis of the model, we applied our method to photospheric observations of bipolar ARs that show clear magnetic tongues. Results. Using the developed procedure on the FR model, we can reduce the deviation in the tilt estimation by more than 60%. Next we illustrate the performance of the method with four examples of bipolar ARs selected for their large magnetic tongues. The new method efficiently removes the spurious rotation of the bipole. This correction is mostly independent of the method input parameters and significant since it is larger than all the estimated tilt errors. Conclusions. We have developed a method to isolate the magnetic flux associated with the FR core during the emergence of bipolar ARs. This allows us to compute the AR tilt-angle and its evolution as precisely as possible. We suggest that the high dispersion observed in the determination of AR tilt-angles in studies that massively compute them from line-of sight magnetograms can be partly due to the existence of magnetic tongues whose presence is not sufficiently acknowledged. Fil: Poisson, Mariano. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; Argentina Fil: López Fuentes, M. C.. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; Argentina Fil: Mandrini, Cristina Hemilse. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; Argentina Fil: Démoulin, Pascal. Centre National de la Recherche Scientifique. Observatoire de Paris; Francia Fil: MacCormack, C.. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; Argentina

Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Open Access funding provided by Max Planck Society., Context. Recently, there have been some reports of unusually strong photospheric magnetic fields (which can reach values of over 7 kG) inferred from Hinode SOT/SP sunspot observations within penumbral regions. These superstrong penumbral fields are even larger than the strongest umbral fields on record and appear to be associated with supersonic downflows. The finding of such fields has been controversial since they seem to show up only when spatially coupled inversions are performed. Aims. Here, we investigate and discuss the reliability of those findings by studying in detail observed spectra associated with particularly strong magnetic fields at the inner edge of the penumbra of active region 10930. Methods. We applied classical diagnostic methods and various inversions with different model atmospheres to the observed Stokes profiles in two selected pixels with superstrong magnetic fields, and compared the results with a magnetohydrodynamic simulation of a sunspot whose penumbra contains localized regions with strong fields (nearly 5 kG at τ = 1) associated with supersonic downflows. Results. The different inversions provide different results: while the SPINOR 2D inversions consider a height-dependent single-component model and return B > 7 kG and supersonic positive vLOS (corresponding to a counter-Evershed flow), height-dependent two-component inversions suggest the presence of an umbral component (almost at rest) with field strengths ∼4 - 4.2 kG and a penumbral component with vLOS ∼ 16 - 18 km s-1 and field strengths up to ∼5.8 kG. Likewise, height-independent two-component inversions find a solution for an umbral component and a strongly redshifted (vLOS ∼ 15 - 17 km s-1) penumbral component with B ∼ 4 kG. According to a Bayesian information criterion, the inversions providing a better balance between the quality of the fits and the number of free parameters considered by the models are the height-independent two-component inversions, but they lie only slightly above the SPINOR 2D inversions. Since it is expected that the physical parameters all display considerable gradients with height, as supported by magnetohydrodynamic (MHD) sunspot simulations, the SPINOR 2D inversions are the preferred ones. Conclusions. According to the MHD sunspot simulation analyzed here, the presence of counter-Evershed flows in the photospheric penumbra can lead to the necessary conditions for the observation of ∼5 kG fields at the inner penumbra. Although a definite conclusion about the potential existence of fields in excess of 7 kG cannot be given, their nature could be explained (based on the simulation results) as the consequence of the extreme dynamical effects introduced by highly supersonic counter-Evershed flows (vLOS > 10 km s-1 and up to ∼30 km s-1 according to SPINOR 2D). The latter are much faster and more compressive downflows than those found in the MHD simulations and therefore could lead to field intensification up to considerably stronger fields. Also, a lower gas density would lead to a deeper depression of the τ = 1 surface, making possible the observation of deeper-lying stronger fields. The superstrong magnetic fields are expected to be nearly force-free, meaning that they can attain much larger strengths than expected when considering only balance between magnetic pressure and the local gas pressure. © A. Siu-Tapia et al.., This work was partly carried out in the framework of the International Max Planck Research School (IMPRS) for Solar System Science at the Max Planck Institute for Solar System Research (MPS) supported by the Max Planck Society. It was also supported by the Spanish Science Ministry “Centro de Excelencia Severo Ochoa” Program under grant SEV-2017-0709 and project RTI2018-096886-B-C51. The National Center for Atmospheric Research is sponsored by the National Science Foundation. Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as a domestic partner and NASA and STFC (UK) as international partners. It is operated by these agencies in co-operation with ESA and NSC (Norway).

It has been found that the Kelvin-Helmholtz instability (KHI) induced by both transverse and torsional oscillations in coronal loops can reinforce the effects of wave heating. In this study, we model a coronal loop as a system of individual strands, and we study wave heating effects by considering a combined transverse and torsional driver at the loop footpoint. We deposit the same energy into the multi-stranded loop and an equivalent monolithic loop, and then observe a faster increase in the internal energy and temperature in the multi-stranded model. Therefore, the multi-stranded model is more efficient in starting the heating process. Moreover, higher temperature is observed near the footpoint in the multi-stranded loop and near the apex in the monolithic loop. The apparent heating location in the multi-stranded loop agrees with the previous predictions and observations. Given the differences in the results from our multi-stranded loop and monolithic loop simulations, and given that coronal loops are suggested to be multi-stranded on both theoretical and observational grounds, our results suggest that the multi-strandedness of coronal loops needs to be incorporated in future wave-based heating mechanisms., Comment: 8 figures, accepted for publication in ApJ

We use the synoptic maps of the photospheric magnetic field observed at Wilcox Solar Observatory, Mount Wilson Observatory, Kitt Peak, SOHO/MDI, SOLIS/VSM, and SDO/HMI to study the distribution of weak photospheric magnetic field values in 1974–2018. We fit the histogram distribution of weak field values for each synoptic map of the six data sets separately with a parameterized Gaussian function in order to calculate the possible shift (to be called here the weak-field asymmetry) of the maximum of the Gaussian distribution from zero. We estimate the statistical significance of the weak-field asymmetry for each rotation. We also calculate several versions of lower-resolution synoptic maps from the high-resolution maps and calculate their rotational weak-field asymmetries. We find that the weak-field asymmetries increase with decreasing map resolution. A very large fraction of weak-field asymmetries are statistically significant, with the fraction of significant weak-field asymmetries increasing with decreasing resolution. Significant weak-field asymmetries of high- and low-resolution maps mainly occur at the same times and have the same sign. Weak-field asymmetries for the different data sets and resolutions vary quite similarly in time, and their mutual correlations are very high, especially for low-resolution maps. These results give strong evidence for weak-field asymmetries reflecting a real feature of weak field values, which is best seen in medium- and low-resolution synoptic maps and is most likely related to the supergranulation scale of the photospheric field.

Context. The Sun features on its surface typical flow patterns called the granulation, mesogranulation, and supergranulation. These patterns arise due to convective flows transporting energy from the interior of the Sun to its surface. The other well known elements structuring the solar photosphere are magnetic fields arranged from single, isolated, small-scale flux tubes to large and extended regions visible as sunspots and active regions. Aims. In this paper we will shed light on the interaction between the convective flows in large-scale cells as well as the large-scale magnetic fields in active regions, and investigate in detail the statistical distribution of flow velocities during the evolution and formation of National Oceanic and Atmospheric Administration active region 11190. Methods. To do so, we employed local correlation tracking methods on data obtained by the Solar Dynamics Observatory in the continuum as well as on processed line-of-sight magnetograms. Results. We find that the flow fields in an active region can be modelled by a two-component distribution. One component is very stable, follows a Rayleigh distribution, and can be assigned to the background flows, whilst the other component is variable in strength and velocity range and can be attributed to the flux emergence visible both in the continuum maps as well as magnetograms. Generally, the plasma flows, as seen by the distribution of the magnitude of the velocity, follow a Rayleigh distribution even through the time of formation of active regions. However, at certain moments of large-scale fast flux emergence, a second component featuring higher velocities is formed in the velocity magnitudes distribution. Conclusions. The plasma flows are generally highly correlated to the motion of magnetic elements and vice versa except during the times of fast magnetic flux emergence as observed by rising magnetic elements. At these times, the magnetic fields are found to move faster than the corresponding plasma.© ESO 2019., This research received support by the Austrian Science Fund (FWF) P27800. Jose Ivan Campos Rozo and Santiago Vargas Dominguez acknowledge funding from Universidad Nacional de Colombia research project code 36127: Magnetic field in the solar atmosphere. Additional funding was possible through an Odysseus grant of the Fund for Scientific Research-Flanders (FWO Vlaanderen), the IAP P7/08 CHARM (Belspo), and GOA-2015-014 (KU Leuven). This work has also received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 724326).

The formation of sunspots and starspots is not yet fully understood and is therefore one of the major open problems in solar and stellar physics. Magnetic flux concentrations can be produced by the negative effective magnetic pressure instability (NEMPI). This instability is strongly suppressed by rotation. However, the presence of an outer coronal envelope was previously found to strengthen the flux concentrations and make them more prominent. It also allows for the formation of bipolar regions (BRs). We want to know whether the presence of an outer coronal envelope also changes the excitation conditions and the rotational dependence of NEMPI. We use direct numerical simulations and mean-field simulations. We adopt a simple two-layer model of turbulence that mimics the jump between the convective turbulent and coronal layers below and above the surface of a star, respectively. The computational domain is Cartesian and located at a certain latitude of a rotating sphere. We investigate the effects of rotation on NEMPI by changing the Coriolis number, the latitude, and the box resolution. Rotation has a strong impact on the process of BR formation. Even rather slow rotation is found to suppress their formation. However, increasing the imposed magnetic field strength also makes the structures stronger and alleviates the rotational suppression somewhat. The presence of a coronal layer itself does not significantly alleviate the effects of rotational suppression., 18 pages, 24 figures, submitted to A&A

Context. Whilst there are observational indications that transverse magnetohydrodynamic (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. Aims. 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. Methods. We considered a 3D MHD model of an active region coronal loop and we perturbed its footpoints with a 2D horizontal driver that represents a random buffeting motion of the loop footpoints. Our driver was composed of 1000 pulses superimposed to generate the observed power spectrum. Results. 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. Conclusions. 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.

Funding: This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No.724326 and No.647214) (IDM). P.A. acknowledges funding from his STFC Ernest Rutherford Fellowship (No. ST/R004285/1). 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 wavemodes 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. Postprint