ORIGINAL RESEARCH article
Significance of Cooling Effect on Comprehension of Kink Oscillations of Coronal Loops
- 1Solar Physics and Space Plasma Research Centre (SP2RC), School of Mathematics and Statistics, University of Sheffield, Sheffield, United Kingdom
- 2School of Build Environment, Engineering and Computing, Leeds Beckett University, Leeds, United Kingdom
- 3Astrophysics Research Centre (ARC), School of Mathematics and Physics, Queen’s University, Belfast, United Kingdom
- 4Laboratory of Interplanetary Medium, Space Research Institute (IKI) Russian Academy of Sciences, Moscow, Russia
- 5Department of Astronomy, Eötvös Loránd University, Budapest, Hungary
Kink oscillations of coronal loops have been widely studied, both observationally and theoretically, over the past few decades. It has been shown that the majority of observed driven coronal loop oscillations appear to damp with either exponential or Gaussian profiles and a range of mechanisms have been proposed to account for this. However, some driven oscillations seem to evolve in manners which cannot be modeled with purely Gaussian or exponential profiles, with amplification of oscillations even being observed on occasions. Recent research has shown that incorporating the combined effects of coronal loop expansion, resonant absorption, and cooling can cause significant deviations from Gaussian and exponential profiles in damping profiles, potentially explaining increases in oscillation amplitude through time in some cases. In this article, we analyze 10 driven kink oscillations in coronal loops to further investigate the ability of expansion and cooling to explain complex damping profiles. Our results do not rely on fitting a periodicity to the oscillations meaning complexities in both temporal (period changes) and spatial (amplitude changes) can be accounted for in an elegant and simple way. Furthermore, this approach could also allow us to infer some important diagnostic information (such as, for example, the density ratio at the loop foot-points) from the oscillation profile alone, without detailed measurements of the loop and without complex numerical methods. Our results imply the existence of correlations between the density ratio at the loop foot-points and the amplitudes and periods of the oscillations. Finally, we compare our results to previous models, namely purely Gaussian and purely exponential damping profiles, through the calculation of
Standing, driven kink-mode oscillations within coronal loops have been extensively studied by the community ever since they were detected in the solar atmosphere at the end of the 20th Century (Aschwanden et al., 1999; Nakariakov et al., 1999) and for reviews see e.g. Andries et al. (2009); Ruderman and Erdélyi (2009). It is now known that, typically, flare-driven kink-mode oscillations are observed to damp within only a few periods, faster than would routinely be expected given the large amplitudes present (see, for example, Aschwanden and Terradas 2008), implying the presence of complex physics during these relatively routine events. The damping profiles of the majority of driven kink-mode oscillations can be well fitted by either Gaussian or exponential profiles (as shown in the statistical study of Goddard et al., 2016) and numerous physical mechanisms have been proposed to account for such decay. Resonant absorption is one such mechanism and has been widely analyzed in numerous analytical works with the aim of better understanding damping within kink-mode oscillations. (e.g., Goossens et al., 2002; Ruderman and Roberts 2002; Dymova and Ruderman 2006; Shukhobodskiy and Ruderman 2018). Interestingly, resonant absorption was proposed well before the first direct detection of kink-modes in coronal loops in the solar atmosphere (Hollweg and Yang 1988; Goossens et al., 1992; Erdelyi and Goossens 1994, Erdelyi and Goossens 1995). The potential transition from Gaussian to exponential damping during the course of a single coronal loop oscillation for propagating kink waves was derived by Hood et al. (2013) in the absence of cooling and for standing kink waves by Ruderman and Terradas (2013). For a review on resonant absorption without complex profiles in a solar context see Goossens et al. (2011).
The improvement in both the spatial and temporal resolutions of coronal imaging which resulted from the launch of the Solar Dynamics Observatory’s Atmospheric Imaging Assembly (SDO/AIA; Lemen et al., 2012) in 2010 has facilitated major advances in understanding of coronal loop oscillations over the past decade. For example, it is now known that multiple harmonics can be present in kink-mode oscillations (Pascoe et al., 2016a), that decay-less observations are present in some coronal loops (Anfinogentov et al., 2013), and that coronal loops can be multi-thermal in nature (Krishna Prasad et al., 2017). One of the main benefits of the continuous, full-disk observations provided by the SDO/AIA instrument is that large statistical studies of coronal loop oscillations can now be conducted (see, e.g., Zimovets and Nakariakov 2015; Goddard et al., 2016). Goddard et al. (2016), for example, studied 120 kink-mode oscillations finding parameters such as amplitudes, damping times, and periods. Ground-based facilities also enable such statistical investigations in chromospheric wave-guides, as were conducted by Kuridze et al. (2012), Morton et al. (2012) who studied the properties of transverse oscillations in mottles. Furthermore, for the example of undamped kink oscillations presented by Aschwanden and Schrijver (2011), Kumar et al. (2013) showed that such phenomena could be triggered by the combination of fast and slow MHD waves present within the system. Moreover, Wang et al. (2012) suggested that amplification could occur due to additional energy input, potentially associated with flaring during the lifetime of the initial oscillation (Pascoe et al., 2020).
Such studies are highly valuable and provide important constraints which can help with the verification of numerical and analytical theories. Interestingly, although the evolution of the majority of coronal loop oscillations analyzed by Goddard et al. (2016) could be modeled using Gaussian, exponential, or Gaussian followed by exponential fits, 21 decaying kink-mode oscillations were identified which had combinations of both exponential and non-exponential damping profiles co-temporally. The presence of these complex or non-standard damping profiles, which cannot be approximated by a few Gaussian profiles and exponential profiles attached to each other in a definite order, sometimes including amplification of the oscillations through time, implies further work is required to fully understand what effects are important in defining the evolution of kink-modes in some coronal loops.
The first analytical models of kink oscillations of coronal loops were simple and considered homogeneous magnetic flux tubes (e.g. Ryutov and Ryutova 1976; Edwin and Roberts 1983). More recently, the idea of analytically studying cooling effect on MHD waves (see e.g. Morton et al., 2010) and in particular on the seismological properties for kink oscillations of coronal loops (see Morton and Erdélyi 2009; Morton and Erdélyi 2010) was proposed. Morton and Erdélyi (2010); Ruderman (2011a) concluded that in the absence of damping due to resonant absorption cooling can cause amplification of coronal loops. Furthermore Ruderman (2011b) concluded that combining cooling and resonant damping can result in the amplitude not varying in time (i.e., the oscillation being decay-less). Ruderman et al. (2017) studied the combination of cooling and expansion for non damped kink oscillation, it was shown that expansion of coronal loops acts in favor of amplification. Shukhobodskiy et al. (2018) considered the similar problem in the presence of damping due to resonant absorption and found that the combination of cooling and expansion can lead to an increase in the oscillation amplitude even in presence of resonant damping. The theory was later tested observationally by Nelson et al. (2019) on one kink-mode oscillation. It was found that increases in the amplitude of the oscillation through time could be explained if cooling was considered, without the need of external forces (e.g., additional flaring) affecting the oscillating system in order to sustain oscillation.
Cooling has been observed in many coronal loops during their lifetimes and has numerous observational signatures, both direct and indirect. The most direct signature is the visible transition of the loops from channels sampling hotter plasma to channels sampling cooler plasma through time (as was shown by, for example, Winebarger and Warren 2005; López Fuentes et al., 2007; Aschwanden and Terradas 2008). More indirect mechanisms include coronal rain (e.g., Antolin et al., 2015) and associated process such as downflows in transition region spectra at the foot-points of coronal loops (Kleint et al., 2014; Ishikawa et al., 2020). Both coronal rain and transition region downflows have been shown to be stable (Antolin et al., 2015; Straus et al., 2015) and variable (Antolin et al., 2015; Nelson et al., 2020a) over short time-scales meaning the thermal evolution of coronal loops may be an on-going and variable process throughout their lifetimes. It is known that transition region downflows occur above almost all sunspots (Samanta et al., 2018; Nelson et al., 2020b) implying that a high number of coronal loops do not maintain constant temperatures through time. The effects of thermal evolution (both heating and cooling) within loops should, therefore, be considered when analyzing the oscillations of coronal loops (see e.g., Morton et al., 2010; Erdélyi et al., 2011; Al-Ghafri and Erdélyi 2013; Al-Ghafri et al., 2014 to name a few studies).
In this article, we analyze 10 examples of kink mode oscillations in coronal loops which damp with non-Gaussian and non-exponential profiles, as identified in the statistical study of Goddard et al. (2016). Our aim is to understand whether a range of damping profiles can be explained using the theories develop in Shukhobodskiy et al. (2018) which accounts for cooling within the loops during the oscillations. We are not claiming that cooling occurs in the studied loops (this would need to be studied in more detail in a separate study) but are, instead, displaying the flexibility of the theory and offering a discussion of when it may be useful in understanding the solar atmosphere. Our work is laid out as follows: In Section 2 we detail the coronal loops studied here and the data used to analyze them; In Section 3 we present our results; Before in Section 4 we draw our conclusions.
2 Observations—Event Selection and Model Fitting
The aim of this article is to understand whether a range of damping profiles of kink mode oscillations can be explained by the effects of cooling and to compare this model to several other previously studied models. To tackle this aim, we selected 10 events randomly from the sample identified to damp with both exponential and non-exponential components in the statistical study of Goddard et al. (2016) for analysis. The combination of exponential and non-exponential components implies a level of complexity in the amplitude profiles through time which could be an observational signature of cooling. Each of these 10 examples of kink-mode oscillations were sampled by the SDO/AIA instrument at discrete times between November 3, 2010 and July 18, 2013, with the periods of these oscillations ranging from slightly over 2 min up to nearly 16 min. Here, we specifically study data from the 171 Å filter which samples the coronal plasma at temperatures of around 0.6 MK with a cadence of 12 s and a pixel scale of 0.6″. A 300″×300″ field-of-view (FOV) around each oscillation was downloaded for a 1-h time-series beginning 15 min before the oscillation was deemed to have begun in Goddard et al. (2016) using the ssw_cutout_service.pro routine. Frames from the SDO/AIA instrument can have reduced exposure times during flares resulting in frames being dropped during download. Therefore, we filled any gaps in the time-series with synthetic images generated by averaging the intensities in the previous and following frames at each pixel. Since no consecutive frames were dropped, this should have no effect on the identification and analysis of oscillations with periods of the order minutes. Time-distance diagrams were then constructed to mimic the slits studied in Goddard et al. (2016), with a width of five pixels. In Table 1, we detail the basic information of the 10 oscillations studied in this article. We should also note that Pascoe et al. (2017) have extensively studied density contrast for the event 40 02. Furthermore, comprehensive temporal analysis was performed by Goddard and Nisticò (2020) for the event 03 01.
In Figure 1, we plot the time-distance diagrams constructed for each of the events studied here (using the routine detailed in Krishna Prasad et al., 2012) over the course of 1 h. Each slit has a different length, however, they have all been scaled to the same y-axis to emphasize the oscillations. The blue dots over-laid on each panel indicate the edges of the loop, identified using the Canny edge-detection method in Wolfram Mathematica 12. The red dots in each panel indicate the mid-points between the detected edges that we infer to be the central axis of the oscillating loop here. The background trend was then accounted for by fitting a function, using spline fitting, through the averages between local extremum points, similar to the method described in Pascoe et al. (2016b). Once the background trend had been subtracted, the amplitudes of the oscillations were normalized and fitted with a sum of Gaussian profiles of the form:
FIGURE 1. Time-space diagrams for each of the selected events studied here. The x-axis of each panel corresponds to 1 h. The y-axes are all different in lengths but are scaled independently to improve the visibility of the oscillations. Exact timing and distance information is provided in Table 1. The results of the Canny edge-detection algorithm are over-laid as blue crosses. The central position between these points, assumed to be the center of the loop axis through time, is signified by the red crosses. The vertical white dashed lines denote the time-period plotted in Figure 2. We note that Event 40 04 could correspond to two individual loops over-lying each other, however, we do not remove it as the complex amplitude profile allows us to show-case both the strengths and limitations of the method.
3.1 Theoretical Model
In this article, we expand on the work of Nelson et al. (2019) by applying the model proposed by Shukhobodskiy et al. (2018) to analyze 10 complex kink-mode oscillations. Here, let us provide a brief summary of the theoretical model for completeness. We consider an expanding and cooling loop of semi-circular shape surrounded by an annulus layer which is fixed in the dense photosphere. We note that the effects of curvature are only applied to the density distribution meaning, essentially, we consider the loop to be a straight magnetic flux tube with radius
The variation of the loop cross-section,
where z is the height of the tube (set to be equal to 0 at the apex),
With the aid of Eqs 2–4, Shukhobodskiy et al. (2018) derived the relation for the dimensionless amplitude,
We note that to obtain results for comparison with observed amplitude profiles, the arbitrary constant,
3.2 Results of the Model Fitting
In order to fit these observations with the model, we first need to obtain the damping profiles for each of the 10 events studied here. In Figure 2, we plot each of the oscillations identified in Figure 1 with red dots. Each of these oscillations appears to be qualitatively different with different periods and amplitudes, with some having large amplitudes and some only appearing to have small amplitudes (only several pixels). This variety in the studied oscillations allows us to test the model in a more dynamic way than if all oscillations displayed similar behavior. With this in mind, we keep the fitting of Event 40 04 despite the apparent overlap between two independent loops in the slit and label this as a “low-confidence” [LC] event in the remainder of this article. The green and blue lines in Figure 2 mark the spline fits (following the method of Pascoe et al., 2016b) and the results of applying Eq. 1, respectively, for each of the events. The dashed black vertical lines indicate the temporal regions of interest within which the oscillations were studied in the remainder of this article. The next step in our analysis was to remove the effects of the background trends and normalize the amplitudes such that the damping profiles could be calculated. The results of this procedure are plotted in Figure 3. It is immediately evident that the damping profiles of Events 03 01, 04 01, 40 02, 40 04 (LC), 40 07, 40 08 and 43 05 clearly deviate from both typical exponential and Gaussian damping. Some of these events display an increase in the amplitude through time whilst some others have stagnation periods where the amplitude remains constant for a while before further amplitude reduction takes place. Such damping profiles are difficult to explain using standard methods.
FIGURE 2. Displacement profile of each coronal loop studied in this article. The red points correspond to the red points over-laid on Figure 1 and track the measured amplitudes of the oscillation. The green lines display the background trend estimated using spline fitting. The blue lines are the results of fitting Eq. 1 to each oscillation. Note the ability of this method to deal with variable periodicities, best seen in Events 40 02 and 40 07. The dashed black vertical lines on each panel outline the time-period plotted in Figure 3.
FIGURE 3. The absolute value of the displacement of the Gaussian fitted position profiles (blue line from Figure 2) from the spline fitted background trend (green line from Figure 2) through time for each of the events studied here.
Next, we applied the model proposed by Shukhobodskiy et al. (2018) to the damping profiles, varying the four parameters discussed in the previous section to minimize the departure from the measured amplitude profile (as was previously done by Nelson et al., 2019). In Figure 4 we plot the amplitudes of each peak of each oscillation as green dots, with spline fits over-laid as blue lines. The optimized results of the model fitting are over-laid as red lines. In order to compare with other single models, we also plot fits calculated using both exponential (dashed) and Gaussian (dot-dashed) profiles, as discussed in Pascoe et al. (2016b), which study the damping of a fundamental mode. The parameters used for the cooling model provide sufficiently good fits in a qualitative sense in the majority of cases, as can be seen in Figure 4, with only Event 40 04 (LC) being visually poorly fitted. Essentially, including the effects of cooling in the fitting of kink-mode damping profiles can produce good fits for a wide variety of profiles but even this method can struggle when low amplitude, highly variable damping profiles are considered. Deviations from standard Gaussian or exponential damping profiles are most prominent in Events 03 01, 03 02, 04 01, 26 01, 40 07 and 43 05 where several saddle points are evident in the fits. It is worth noting that this model is also capable of capturing features of both the Gaussian and exponential damping profiles by varying its parameters, with near-Gaussian damping profiles being returned for Events 40 08 and 48 01 and near-exponential damping being returned for Event 40 02. We note that Gaussian damping profiles can only be obtained for low values of ζ (
FIGURE 4. The fitting of the model, introduced by Shukhobodskiy et al. (2018), to the coronal loop oscillations studied here. The blue lines plot spline fits through the maxima of each peak plotted in Figure 3. The red lines plot the approximation of the oscillations obtained through fitting the theoretical model. Reasonable fits are achieved for all events except Event 40 04 (LC). The dot-dashed lines represents Gaussian fit and the dashed lines correspond to exponential fits. Finally, the red crosses indicate the peaks from Figure 3.
Highly complex amplitude profiles, for example Event 43 05, display both decreases and increases in the amplitude through time. Fortunately, this complexity can be explained by theory once cooling is considered. For lower values of ζ less than 1.2 the theoretical profile reduces to Gaussian decay. Whereas, for higher values of density ratios, the initial part would start with exponential decrease first and then after some time the local increase may appear after which the Gaussian decay profile appears. This behavior provides higher degrees of freedom in fitting the observations meaning more amplitude profiles are now theoretically permitted. This fact makes it ultimately both convenient and useful to consider the evolution of temperature during oscillations for accurate approximation of seismological parameters. Goddard and Nakariakov (2016) discussed the importance of the projection angle in order to obtain more accurate information about coronal loops, which could cause a cascade of associated errors and require the need to perform multiple angle analysis of the same phenomena. Despite this fact, the model used in this article is not dependent of the absolute values of the damping profile, rather it depends on the normalized shape of oscillation, making it resistant to the above errors in obtaining the dimensionless parameters. As such, that could lead to faster analysis of kink oscillations of coronal loops. However, such an error does limit the estimation of parameters with dimensions in the model used here.
We should note that ζ is responsive to the position where the analytical solution changes it shape (e.g., where amplitude increases are detected). Any increases in the assumed cooling time results in a stretching of the analytical solution and lowers the value of ζ, however, such effects becomes slower or negligible as ζ becomes closer to 1 as seen in Figure 5. Additionally, increases in the cooling time lead to increases of κ in order to fit the observed shape. Nevertheless, in all occasions the model is able to capture the observed profile for various lengths of cooling. Essentially, this introduces some degeneracy in our solutions which will need to be further studied through future research. Moreover, oscillations with low values of ζ are barely affected by the cooling time, where kink oscillations with high values of ζ could lead to overestimation of this parameter and require additional measurements to determine
FIGURE 5. Examination of the effect of cooling time on the fits of the analytical model to observations. The κ values for each of the cooling times (from quickest to slowest) are
The inferred, best fit, parameters (defined in the previous section) for each of the oscillations are presented in Table 2. To obtain accurate approximations for ζ, care should be taken in accurately positioning the dashed black lines in Figure 2 due to its high sensitivity to temporal variations. In Figure 6, we plot the relationship between the density ratio at the loop foot-points, ζ, and the oscillation amplitudes (left panel) and periods (right panel) as measured by Goddard et al. (2016). It is immediately evident that the higher amplitude and longer period oscillations appear to have lower density ratios between the loop and the external environment at the foot-points. Although we are hesitant to draw strong conclusions about this due to the small statistical size and lack of direct density diagnostics, this result does suggest further research should be done to investigate such links in the future. Some important steps forward could come through analyzing the ratios between the O IV lines sampled by the Interface Region Imaging Spectrograph (IRIS; De Pontieu et al., 2014) at the foot-points of oscillating loops. Some limitations of this model are evident in Figure 4 where deviations between the fitted theoretical model and the measured observed values are evident. It is unclear at this stage whether these limitations are caused by differences in the observed and modeled structures or due to measurement errors. It should be noted that the recent work of Goddard et al. (2017) shows that some of the loops might not have thin boundaries, with Pascoe et al. (2019) discussing the possible limitations of the thin boundary approximation used within theoretical approach. This could possibly lead to underestimation of ζ, with the current model supporting
TABLE 2. Comparative analysis for goodness of fit from the observed data for the theoretical model presented by Shukhobodskiy et al. (2018) and for Gaussian and exponential fits presented by Pascoe et al. (2016b).
FIGURE 6. Scatter plots displaying the apparent relationships between the amplitude and period of the coronal loop oscillations (as measured by Goddard et al., 2016) and the density ratio inferred here.
Table 3 provides the comparison with previous models, we could see that inclusion of cooling to Events 03 02, 26 01, 40 02, 40 07 with oscillation time been equal to the oscillation time could lead the cooling theoretical model to have better approximations to the damping profiles than exponential and Gaussian fits, in case only extremum points (i.e. amplitude maxima) are taken as the point of reference. If instead we analyze the overall damping profile (i.e. calculate the
TABLE 3. Best-fit parameters fitted by the theoretical model for each of the loop oscillations studied here.
In this article, we have further demonstrated the promise of the theory developed by Shukhobodskiy et al. (2018) by applying it to 10 driven kink-mode oscillations. These events were sampled by the SDO/AIA instrument and were randomly selected from the population identified to have complex damping profiles in the statistical study of Goddard et al. (2016). Space-time diagrams were constructed for each oscillation before the edges of the loops through time were determined by employing the Canny edged-detection method (see Figure 1). The mid-point between the edges was defined as the central axis of the oscillating loop for each time-step. The background trends were then subtracted from each oscillation using spline fitting before the detrended data were then fitted with a series of Gaussian functions. The results of the application of these steps to the oscillations can be seen in Figure 2.
Next, we plotted the absolute values of the amplitude variation of the kink-mode oscillations through time in Figure 3. Some oscillations displayed clear deviations from the usual Gaussian or exponential damping profiles even at this stage. Specifically, five (Events 04 01, 40 02, 40 04 (LC), 40 07, and 43 05) of the 10 oscillations studied here displayed evidence of increases in amplitude through time. Such amplification was shown to be possible during oscillations of cooling loops by Ruderman (2011b), Shukhobodskiy and Ruderman (2018) and has been reported before in observations by Nelson et al. (2019). Using an iterative approach, we then constructed best fits for the parameters defined in Eq. 5 such that the difference between the observed and theoretical damping profiles was minimized. In Figure 4 it was shown that this theoretical model provides qualitatively reasonable approximations for the damping of kink oscillations of coronal loops for nine of the 10 cases (all except Event 40 04 [LC]). This result indicates that cooling should be studied in more detail in future research. Additionally, we have also shown that the model is capable of capturing the properties of fitting both simple Gaussian and exponential profiles, depending on the ratio of external and internal densities, ζ, at the loop foot-points. Furthermore, for sufficiently large loops with
It was shown by Shukhobodskiy et al. (2018) that higher expansion factors of coronal loops could result in overall amplification of the first part of the oscillation eventually resulting in Gaussian decay of the amplitude profile. Thus studying complex damping profiles could allow us to extract interesting seismological information, such as the ratios of densities at foot-points. Furthermore, the information presented in Figure 6 suggests, that lower ratios of internal and external densities have higher amplitude and period. One of the reasons why this could happen is that the coronal loops investigated are lighter and less disturbance created by the flare is needed to disturb the system more prominently. A larger statistical study and more targeted observations with spectral instruments suitable for providing density diagnostics (such as IRIS) should be used to study this potenatial result in more detail in the future.
Overall, we presented an alternative approach to model more complex damping profiles that deviate from the widely studied Gaussian and exponential decays. A number of examples of such deviation were addressed in this study, thus strengthening the claim made by Nelson et al. (2019) that the effect of cooling should be further considered for analyzing oscillations of coronal loops. This is in agreement with some earlier studies, e.g. by Morton and Erdélyi (2009), Morton and Erdélyi (2010). Furthermore, in the case of strong damping the theoretical model proposed by Shukhobodskiy et al. (2018) is able to provide seismological information, just from the damping profile alone. For some events, the cooling model demonstrated significant improvement over standard Gaussian on exponential fitting when
Although little evidence is presented about the actual cooling for the events studied here, it is possible that the gradual fading of the intensity of several coronal loops presented in Figure 1 (specifically for Events 03 01, 04 01, and 43 05) may imply that some thermal evolution (either heating or cooling) is present. In a future work, it would be worthwhile to analyze oscillations within loops which are undergoing specific processes linked to cooling, such as coronal rain (Antolin et al., 2015) or the associated high velocity transition region downflows (see, for example, Nelson et al., 2020b, Nelson et al., 2020a and literature within), to understand whether abnormal amplitude profiles through time (such as those studied here) are present. Additionally, it would be worthwhile to develop further the analytical methods, which then could help to extract additional information about the kink oscillation of coronal loops, in order to optimize this theory for the application of solar magneto-seismology to kink-mode oscillations.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://jsoc.stanford.edu.
AS and CN conceived the study. DS and AS conducted the data analysis, with DS having a major contribution to the analysis. All authors consulted on the interpretation of results. DS, AS, and CN drafted the manuscript. All authors read, discussed and commented on the draft. RE lead the research overall.
DS is grateful to the University of Sheffield for the support received. AS is thankful to Interreg Northwest Europe for the support received to conduct this research through grant number: NWE847. CN is thankful to the Science and Technology Facilities Council (STFC) for the support received to conduct this research through grant number: ST/P000304/1. RE is grateful to STFC (UK, grant number ST/M000826/1). RE also acknowledges support from the Chinese Academy of Sciences President’s International Fellowship Initiative (PIFI, grant number 2019VMA0052) and The Royal Society (grant nr IE161153).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
SDO/AIA data are courtesy of NASA/SDO and the AIA science team. We thank S. Krishna Prasad for making their slit creation code available to us.
Al-Ghafri, K. S., Ruderman, M. S., Williamson, A., and Erdélyi, R. (2014). Longitudinal Magnetohydrodynamics oscillations in dissipative, cooling coronal loops. Astrophys. J. 786, 36. doi:10.1088/0004-637X/786/1/36
Andries, J., van Doorsselaere, T., Roberts, B., Verth, G., Verwichte, E., and Erdélyi, R. (2009). Coronal seismology by means of kink oscillation overtones. Space Sci. Rev. 149, 3–29. doi:10.1007/s11214-009-9561-2
Antolin, P., Vissers, G., Pereira, T. M. D., Rouppe van der Voort, L., and Scullion, E. (2015). The multithermal and multi-stranded nature of coronal rain. Astrophys. J. 806, 81. doi:10.1088/0004-637X/806/1/81
Aschwanden, M. J., Fletcher, L., Schrijver, C. J., and Alexander, D. (1999). Coronal loop oscillations observed with the transition region and coronal explorer. Astrophys. J. 520, 880–894. doi:10.1086/307502
Aschwanden, M. J., and Schrijver, C. J. (2011). Coronal loop oscillations observed with atmospheric imaging assembly—kink mode with cross-sectional and density oscillations. Astrophys. J. 736, 102. doi:10.1088/0004-637X/736/2/102
De Pontieu, B., Title, A. M., Lemen, J. R., Kushner, G. D., Akin, D. J., Allard, B., et al. (2014). The Interface region imaging Spectrograph (IRIS). Sol. Phys. 289, 2733–2779. doi:10.1007/s11207-014-0485-y
Goddard, C. R., Nisticò, G., Nakariakov, V. M., and Zimovets, I. V. (2016). A statistical study of decaying kink oscillations detected using SDO/AIA. Astron. Astrophys. 585, A137. doi:10.1051/0004-6361/201527341
Goddard, C. R., Pascoe, D. J., Anfinogentov, S., and Nakariakov, V. M. (2017). A statistical study of the inferred transverse density profile of coronal loop threads observed with SDO/AIA. Astron. Astrophys. 605, A65. doi:10.1051/0004-6361/201731023
Goossens, M., Andries, J., and Aschwanden, M. J. (2002). Coronal loop oscillations. An interpretation in terms of resonant absorption of quasi-mode kink oscillations. Astron. Astrophys. 394, L39. doi:10.1051/0004-6361:20021378
Hood, A. W., Ruderman, M., Pascoe, D. J., De Moortel, I., Terradas, J., and Wright, A. N. (2013). Damping of kink waves by mode coupling. I. Analytical treatment. Astron. Astrophys. 551, A39. doi:10.1051/0004-6361/201220617
Ishikawa, R. T., Katsukawa, Y., Antolin, P., and Toriumi, S. (2020). Temporal and spatial scales in coronal rain revealed by UV imaging and spectroscopic observations. Sol. Phys. 295, 53. doi:10.1007/s11207-020-01617-z
Kleint, L., Antolin, P., Tian, H., Judge, P., Testa, P., De Pontieu, B., et al. (2014). Detection of supersonic downflows and associated heating events in the transition region above sunspots. Astrophys. J. Lett. 789, L42. doi:10.1088/2041-8205/789/2/L42
Krishna Prasad, S., Jess, D. B., Klimchuk, J. A., and Banerjee, D. (2017). Unravelling the components of a multi-thermal coronal loop using magnetohydrodynamic seismology. Astrophys. J. 834, 103. doi:10.3847/1538-4357/834/2/103
Kumar, P., Cho, K. S., Chen, P. F., Bong, S. C., and Park, S. H. (2013). Multiwavelength study of a solar eruption from AR NOAA 11112: II. Large-scale coronal wave and loop oscillation. Sol. Phys. 282, 523–541. doi:10.1007/s11207-012-0158-7
Kuridze, D., Morton, R. J., Erdélyi, R., Dorrian, G. D., Mathioudakis, M., Jess, D. B., et al. (2012). Transverse oscillations in chromospheric mottles. Astrophys. J. 750, 51. doi:10.1088/0004-637X/750/1/51
Lemen, J. R., Title, A. M., Akin, D. J., Boerner, P. F., Chou, C., Drake, J. F., et al. (2012). The atmospheric imaging assembly (AIA) on the solar dynamics observatory (SDO). Sol. Phys. 275, 17–40. doi:10.1007/s11207-011-9776-8
Morton, R. J., and Erdélyi, R. (2010). Application of the theory of damping of kink oscillations by radiative cooling of coronal loop plasma. Astron. Astrophys. 519, A43. doi:10.1051/0004-6361/201014504
Morton, R. J., Verth, G., Jess, D. B., Kuridze, D., Ruderman, M. S., Mathioudakis, M., et al. (2012). Observations of ubiquitous compressive waves in the Sun’s chromosphere. Nat. Commun. 3, 1315. doi:10.1038/ncomms2324
Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., and Davila, J. M. (1999). TRACE observation of damped coronal loop oscillations: implications for coronal heating. Science. 285, 862–864. doi:10.1126/science.285.5429.862
Nelson, C. J., Krishna Prasad, S., and Mathioudakis, M. (2020a). Evolution of downflows in the transition region above a sunspot over short time-scales. Astron. Astrophys. 640, 13. doi:10.1051/0004-6361/202038155
Nelson, C. J., Shukhobodskiy, A. A., Erdélyi, R., and Mathioudakis, M. (2019). The effect of cooling on driven kink oscillations of coronal loops. Frontiers in Astronomy and Space Sciences. 6, 45. doi:10.3389/fspas.2019.00045
Pascoe, D. J., Goddard, C. R., and Nakariakov, V. M. (2016a). Spatially resolved observation of the fundamental and second harmonic standing kink modes using SDO/AIA. Astron. Astrophys. 593, A53. doi:10.1051/0004-6361/201628784
Pascoe, D. J., Goddard, C. R., Nisticò, G., Anfinogentov, S., and Nakariakov, V. M. (2016b). Damping profile of standing kink oscillations observed by SDO/AIA. Astron. Astrophys. 585, L6. doi:10.1051/0004-6361/201527835
Pascoe, D. J., Hood, A. W., and Van Doorsselaere, T. (2019). Coronal loop seismology using standing kink oscillations with a lookup table. Frontiers in Astronomy and Space Sciences. 6, 22. doi:10.3389/fspas.2019.00022
Pascoe, D. J., Russell, A. J. B., Anfinogentov, S. A., Simões, P. J. A., Goddard, C. R., Nakariakov, V. M., et al. (2017). Seismology of contracting and expanding coronal loops using damping of kink oscillations by mode coupling. Astron. Astrophys. 607, A8. doi:10.1051/0004-6361/201730915
Ruderman, M. S., Shukhobodskaya, D., and Shukhobodskiy, A. A. (2019). Resonant damping of propagating kink waves in non-stationary, longitudinally stratified, and expanding solar waveguides. Front. Astron. Space Sci. 6, 10. doi:10.3389/fspas.2019.00010
Shukhobodskiy, A. A., Ruderman, M. S., and Erdélyi, R. (2018). Resonant damping of kink oscillations of thin cooling and expanding coronal magnetic loops. Astron. Astrophys. 619, A173. doi:10.1051/0004-6361/201833714
Keywords: solar corona, coronal loop oscillations, magnetohydrodynamics, kink oscillations, waves
Citation: Shukhobodskaia D, Shukhobodskiy AA, Nelson CJ, Ruderman MS and Erdélyi R (2021) Significance of Cooling Effect on Comprehension of Kink Oscillations of Coronal Loops. Front. Astron. Space Sci. 7:579585. doi: 10.3389/fspas.2020.579585
Received: 02 July 2020; Accepted: 24 November 2020;
Published: 14 January 2021.
Edited by:Valery M. Nakariakov, University of Warwick, United Kingdom
Reviewed by:Pankaj Kumar, Goddard Space Flight Center, National Aeronautics and Space Administration, United States
David James Pascoe, KU Leuven, Belgium
Copyright © 2021 Shukhobodskaia, Shukhobodskiy, Nelson, Ruderman and Erdélyi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Robert Erdélyi, firstname.lastname@example.org