|Home | About | Journals | Submit | Contact Us | Français|
The existence of the Sun’s hot atmosphere and the solar wind acceleration continues to be an outstanding problem in solar-astrophysics. Although magnetohydrodynamic (MHD) modes and dissipation of magnetic energy contribute to heating and the mass cycle of the solar atmosphere, yet direct evidence of such processes often generates debate. Ground-based 1-m Swedish Solar Telescope (SST)/CRISP, Hα 6562.8Å observations reveal, for the first time, the ubiquitous presence of high frequency (~12–42mHz) torsional motions in thin spicular-type structures in the chromosphere. We detect numerous oscillating flux tubes on 10 June 2014 between 07:17 UT to 08:08 UT in a quiet-Sun field-of-view of 60”×60” (1”=725km). Stringent numerical model shows that these observations resemble torsional Alfvén waves associated with high frequency drivers which contain a huge amount of energy (~105Wm−2) in the chromosphere. Even after partial reflection from the transition region, a significant amount of energy (~103Wm−2) is transferred onto the overlying corona. We find that oscillating tubes serve as substantial sources of Alfvén wave generation that provide sufficient Poynting flux not only to heat the corona but also to originate the supersonic solar wind.
Continuous generation of radiation and supersonic wind from the Sun’s chromosphere and corona requires a large input of energy (~102–104Wm−2) to balance these losses1. The role of magnetohydrodynamic (MHD) waves and small-scale magnetic reconnection causing nano-flare heating have been explored as primary candidates to energize the solar atmosphere. However, direct evidence of energy sources and their dissipation are not yet fully understood2,3,4. In the era of high resolution space and ground-based observations, it is now revealed that energy and mass transport in the quiescent solar atmosphere are associated with localized static and flowing flux tubes (e.g., network & inter-network magnetic fields, spicules, vortices etc) possessing various plasma and wave processes5,6,7,8. Here, we observe directly, for the first time, the ubiquitous presence of high frequency (~12–42mHz) torsional oscillations at apparent surfaces composed of thin spicular-type structures rooted in the quiet-Sun magnetic network. These observations are described by torsional Alfvén waves associated with high frequency drivers transferring ~103Wm−2 energy into the overlying corona. These oscillating tubes serve as substantial sources of Alfvén wave generation providing sufficient Poynting flux to heat the solar corona and in originating the nascent solar wind.
Quiet-sun magnetic networks are the locations where field lines fan out into the outer atmosphere supporting waves and exotic plasma dynamics9,10. The magnetic skeleton of the bundle of fine structured small-scale flux tubes becomes visible when remnants of plasma flows (e.g., spicules, jets, surges) are confined within their boundaries. Various other similar structures are prevalent in the solar chromosphere and well resolved with modern day instruments, e.g., on-disk counterparts of type-II spicules, chromospheric counterparts of the transition region network jets, etc.11,12. A tube with its fine structures, each 120–215km wide, is observed using CRISP on the Swedish Solar Telescope (see yellow and green expanding cylinders in Fig. 1A.1, respectively). In an integrated view, the tube’s projected height and width at the top are respectively ~4Mm and ~1.5Mm. Using visualization software, we measure the structure’s length, while the width is an average of the values measured at three locations along its length. The fine structures show collective motions and behave as part of an integrated self-contained magnetic flux tube. This is evident on two clearly visible fine structures showing periodic reversal of velocity sign, i.e., blue on top & red on bottom to red on top & blue on bottom. Thereafter, the previous original condition of velocity signs is restored, i.e., blue on top & red on bottom (Fig. 1A.2). It should be noted that blue and red are the signs of line-of-sight (LOS) components of respectively outward and downward motions, which are tangled with each other. This is the observational signature of torsional oscillations of the fine structured tube over an apparent surface.
RBEs/RREs have three instantaneous chromospheric motions, i.e., up-flows, swaying, and torsional motion13. The present observations differ from RBEs /RREs in two aspects: (i) an absence of flows along their length, and (ii) simultaneous red and blue-shifted emission. The adjacent red-blue shift pattern in the chromospheric line was also observed in newly discovered small-scale twisted flux tubes, although the time evolution of the Doppler shift pattern could not be revealed14.
As mentioned above this paper invokes solar chromospheric observational data to infer the presence of high-frequency Alfvén waves in the solar atmosphere. The observations are obtained with Crisp Imaging Spectropolarimeter (CRISP) on the ground-based Swedish 1-m Solar Telescope15,16. CRISP has a field-of-view of 60”×60” (1”=725km) and a pixel scale of 0.0592”. The data is obtained with the cadence of 3.9s on 10 June 2014 between 07:17 UT to 08:08 UT. The FOV of these observations was centered at (Xc, Yc)=(403”, −211”). It contains a pore and quiet-Sun region in the southward side of NOAA AR 12080. The adaptive optics (AO) system was running on SST and used the pore as a tracking point during the acquisition of the data. Nine Hα line positions are sampled in sequence, with eight images being collected at each of −1032, −774, −516, −258, 0 +258, +516, +774, +1032mÅ w.r.t. the line core at 6562.8Å before 36 frames were collected at each of −774mÅ and the line core. After multi-object & multi-frame blind deconvolution (MOMFBD) reconstruction17, the desired data is obtained with resolution of 0.14”. All the observational data (cases) and their details are outlined in Supplementary Table 1 (see in Supplementary Material Information).
The average line profile is defined by taking the average of the FOV and then the line profile is calculated using the pixels along the length of the event. Doppler shifts (Fig. 1A.2) are estimated w.r.t. the line center at Hα 6562.8Å. Therefore, the velocities are estimated from a simple Doppler formula. We do not measure the full-width-at-half-maximum (FWHM) of the spectral-line as our event is observed in the Hα wings, which is sensitive to the temperature, opacity and velocity gradients, in addition to mass motions18. Moreover, we limit the spectral resolution by defining the line-positions to gain the required temporal resolution to observe high-frequency oscillations in the observed fast evolving chromospheric features.
The transverse motion of one of these fine-structures on the tube is presented by choosing slits P1, P2, ..P10 at various spatial positions (Fig. 1‘B’). In a 2-D projected plane, it exhibits lateral transversal motions at various heights. In the time-distance plot at various locations (Fig. 1‘C’), the position of this dark structure and its displacement in time is estimated by locating the minimum intensity at each temporal step. The spatial position of the minimum intensity at each time step is estimated with sub-pixel accuracy using an FFT cross-correlation technique to disentangle the lateral fluctuations and the bulk motion of the fine structure, and fitting linearly its true transverse motions. Careful investigation shows almost negligible time-lag between the major peak and valley of the transverse motions of this structure at P1 and P10 (Fig. 1‘C’). A wavelet analysis of the transverse oscillations shows the presence of significant power associated with ~47sec period (~21mHz)19 (see Supplementary Material Information). Overall, this indicates the direct presence of long-wavelength high-frequency torsional waves that already moved to the top of the tube20. The periodic Doppler motions of the flux tube clearly show the torsional oscillations as pronounced on an apparent surface, which is the first direct observation of high-frequency (<50sec) torsional Alfvén waves in the chromosphere. It should be noted that the first evidence of long-period (126–700s) torsional Alfvén waves in terms of the oscillation of full-width-at-half-maximum (FWHM) of Hα spectral line had previously been reported above magnetic bright points6. The present observations reveal the first evidence of high-frequency torsional waves in numerous fine structured chromospheric flux-tubes (see Supplementary Table 1 in Supplementary Material Information) at larger scales along with more complete energy estimates using a stringent numerical model.
Direct observations of such high frequency torsional oscillations of the fine structured flux tubes in the chromosphere are difficult due to opacity, ongoing heating and cooling processes, and shallowness of the atmosphere21. Moreover, this requires high spatial and temporal resolution observations, which were made in the present case by SST/CRISP16,17 at various line positions around Hα 6562.8Å. We observe twelve such cases of high frequency torsional oscillations on fine structured tubes. Their details and properties are given in Supplementary Table 1.
A physical model of such localized solar structures, subject to the torsional oscillations (Fig. 1), is made by a simple but realistic magnetic flux tube using the FLASH code (Fig. 2). The model is an axis-symmetric tube rooted at the solar photosphere (Bo=121G)22, gravitationally stratified, and determined by realistic temperature distribution23. For the sake of simplicity, we model only the dynamics of a single enveloping tube similar in size to the observed fine structured flux tube. A concise description of 3-D magnetohydrodynamic (MHD) equations, initial conditions of the plasma and realistic magnetic field as well as their equilibrium in a gravitationally stratified tube, and numerical method and setup, are outlined in the literature23,24,25,26.
The flux tube, where magnetic and gas pressure are initially balanced keeping it in non-force free equilibrium, is then subjected to a perturbation in the azimuthal component of velocity by a periodic high-frequency driver. The form of velocity perturbation is given as,
where Av is the amplitude of the pulse, y0=500km its vertical position, and w=300km its width. Av is set to 150kms−1, which results in an effective maximum velocity of about 2.4kms−1. The period of the driver (Pd) is taken as 50s with wave perturbations generated at the top of the photosphere.
In order to solve the 3-D MHD equations of the model23 numerically, we use the FLASH code25,26 with the third-order un-split Roe Riemann solver and the minmod slope limiter, as well as adaptive mesh refinement (AMR). We set the simulation box as (1.5, 1.5, 18) Mm3, and fix in time all plasma quantities to their equilibrium values at all six boundary surfaces. We use a static, non-uniform grid with a minimum (maximum) level of refinement set to 2 (5). Numerical results are displayed in Figs 2 and and33.
The horizontal convective motions present in the quiet-Sun photosphere can generate Alfvén waves carrying significant amount of energy to heat the corona27,28. Such motions at the photosphere and chromosphere correlate well at the spatial scales of super-granular cells, therefore, indicating the transfer of energy seen as velocity fluctuations even up to the chromosphere29. Exclusively observed oscillating chromospheric flux tubes are visible even in the blue-wing (−774mÅ) of Hα 6562.8Å, which may form at least 200km above the photosphere30. This indicates that the origin of transverse perturbations is somewhere at the top of the photosphere either by the transfer of gigantic velocity fluctuations or by impulsive magnetic reconnection31,32,33. Such impulses result in an azimuthal velocity perturbations (Equation 1) to excite the torsional Alfvén waves in the fine structured flux tubes.
Energy flux (W=0.5ρ VA V2θ in W m−2) is calculated using the estimated density (ρ), wave velocity amplitude (Vθ), and local Alfvén velocity (VA) with height averaged over the horizontal spatial scales in the modeled flux-tube (−0.3Mm<x<+0.3Mm) by considering the tube as a whole. Therefore, the resultant velocity amplitude and thus derived energy flux at each height signify the net contribution of Alfvén waves excited in the tube. Energy flux oscillations result from the periodic driver which varies in time with the period, and its spatial averaging also contributes to the oscillations. As the energy flux is determined by V2θ, its averaging does not lead to a zero in the energy flux. The flux tube at its equilibrium is radially inhomogeneous. Therefore, Alfvén waves and their energy estimates weakly depend on a magnitude of equilibrium magnetic field. For a stronger magnetic field, the Alfvén speed may attain a larger value but it has little effect on the energy flux which depends on the equilibrium mass density and V2θ. The high-frequency Alfvén waves carry ~104W m−2 energy in the chromosphere, and fulfill the huge requirement of energy there34. After partial reflection from the solar transition region, it still contains enough energy (~103W m−2) to compensate the inner coronal energy losses1 (Fig. 3). This very basic analysis shows that the observed high-frequency torsional Alfvén waves can act as a substantial source to channel sufficient energy to heat the corona. Velocity fluctuations related to these Alfvén waves penetrate into the corona carrying associated energy (Fig. 3), while magnetic field perturbations are mostly concentrated up to the upper transition region and inner corona (Fig. 2). These structures are ubiquitous in the quiet-Sun chromosphere (Fig. 1). For the stronger magnetic field at the tube’s foot-point, the nature of the wave is the same. However, it may pump slightly more energy to the overlying atmosphere20.
The structuring of the plasma and magnetic field along the observed flux tube may lead to an amplitude variation as is evident in the standard-deviation (SDEV) measurement of the time-distance profile of transverse motions at its different locations (see Supplementary Material Information. This amplitude variation may further introduce a variation in energy flux. This is in qualitative agreement with the calculation of energy flux at different heights of the modeled flux tube (Fig. 3).
During MHD wave evolution in realistic model flux tube coupling the solar photosphere to the corona, various other important physical processes are evident. There is a generation of vertical plasma flows due to pondermotive forces under non-linear conditions. However, these up and down flows of the plasma are mostly trapped in the solar atmosphere and do not launch the solar wind unless the escape velocity is achieved. Another effect is the evolution of multiple concentric magnetic shells within the tube where Alfvén perturbations are pronounced with time creating complex velocity fields. Essentially, the velocity perturbations over different magnetic shells (or surfaces) are not in phase, which is a fundamental property of Alfvén waves in stratified inhomogeneous tubes.
The frequency range of the observed torsional Alfvén waves in various magnetic flux tubes lies between 12–42mHz (Supplementary Table 1). Eight of these observed structures, associated with strong amplitudes, exhibit one full period of the traversal oscillations, while three show the half-period of the oscillations and disappear quickly (Fig. 1 and Supplementary Table 1). Our theoretical model depicts the wave propagation in Eulerian formalism and shows a clear evolution of the torsional oscillations. The observed chromospheric structures are the fast evolving features carrying similar high-frequency oscillations, but they fade away quickly. The major event presented in the paper (Fig. 1) exhibits almost two cycles of the oscillations vis-à-vis matching with the model. However, many of these incompressible waves, having large-amplitude, may undergo mode-conversion to compressible waves in a non-linear regime in the presence of a pondermotive force due to total magnetic pressure variations and may dissipate their energy quickly35. The fading of the oscillatory tubes may also be due to the strong velocity gradient at their edges and thereafter evolution of the Kelvin-Helmholtz instability36. It should be noted that these fast waves also seem to be associated with long wavelengths compared to the length of the tube itself. This is the reason why every tube shows torsional motion almost as a whole. These high-frequency torsional Alfvén waves are detected directly in such localized small-scale flux tubes.
Strongly homogeneous tubes in the chromosphere may be subject to coupling of the kink modes to the surface Alfvén waves. However, this may take a longer time compared to the dynamics of the observed structures here, while a discontinuity layer in the tube’s boundary is established with sufficient thickness37,38. Keeping in view the complexity (e.g., twist) in flux tube, the axis-symmetric sausage mode signature may also contribute to the Doppler velocity variations39. But, there is no evidence of an initial twist in the observed fine-structured tubes. Therefore, we rule out these possibilities in the present context. In conclusion, the SST/CRISP field-of-view of 60”×60” (1”=725km) detects at least 40 oscillating flux tubes in the chromosphere. They seem to be ubiquitous in various regions of the Sun’s chromosphere40 indicating the presence of substantial energy sources beneath the corona. Torsional waves through many of these flux tubes (‘12’) carry sufficient energy up to the TR/inner corona, which may be potentially used in heating the localized atmosphere and energizing the supersonic wind by its dissipation due to energy cascade at smaller spatial scales41. These oscillating chromospheric structures will likely be detected in abundant measure with ultra high resolution observations with the next generation solar telescopes, e.g., 2-m National Large Indian Telescope (NLST), 4-m DKIST, European Solar Telescope (EST), 8-m Chinese Giant Telescope etc, and new information will further enlighten our understanding on their dynamics and potential role in the solar atmosphere42,43,44.
How to cite this article: Srivastava, A. K. et al. High-frequency torsional Alfvén waves as an energy source for coronal heating. Sci. Rep. 7, 43147; doi: 10.1038/srep43147 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A.K.S. & B.N.D. acknowledge the RESPOND-ISRO (DOS/PAOGIA205-16/130/602) project. A.K.S. acknowledges the SERB-DST project (YSS/2015/000621) grant, the Advanced Solar Computational & Analyses Laboratory (ASCAL), and infrastructural facilities at the Department of Physics, IIT (BHU) to pursue this research. Armagh Observatory is grant-aided by the N.Ireland Department for Communities. The authors acknowledge the DJEI/DES/SFI/HEA Irish Centre for High-End Computing (ICHEC) for the provision of computing facilities and support. We also like to thank STFC for a studentship and PATT T&S and the Solarnet project which is supported by the European Commission’s FP7 Capacities Programme under Grant Agreement number 312495 for T&S. J.S. is funded by the Leverhulme Trust. M.S. aknowledges support by the “Progetti di ricerca INAF di Rilevante Interesse Nazionale” (PRIN-INAF 2014) and PRIN MIUR 2012 (prot. 2012P2HRCR) entitled “Il sole attivoed i suoi effetti sul clima dello spazio e della terra” grants, funded by the Italian National Institute for Astrophysics (INAF) and Ministry of Education, Universities and Research (MIUR), respectively. KM’s and DW’s work was done in the frame-work of the project from the National Science Centre, Poland Grant No: 2014/15/B/ST9/00106. 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 Observatorio del Roque de los Muchachos of the Instituto de Astrofsica de Canarias.
The authors declare no competing financial interests.
Author Contributions A.K.S. leads the project by defining the novel science case and writing the paper. J.S. has prepared the observational data, while K.M. & D.W. provided the numerical runs based on science case. J.G.D. and M.S. have contributed on SST data and its most plausible connection with the numerical model. E.S. & T.R. have provided the SST/CRISP observations. J.G.D. & B.N.D. have an overview on the write-up and consistency on presenting the scientific results.