Issue 
A&A
Volume 585, January 2016



Article Number  A137  
Number of page(s)  9  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201527341  
Published online  11 January 2016 
A statistical study of decaying kink oscillations detected using SDO/AIA
^{1} Centre for Fusion, Space and Astrophysics, Department of Physics, University of Warwick, CV4 7AL, UK
email: c.r.goddard@warwick.ac.uk
^{2} Astronomical Observatory at Pulkovo of the Russian Academy of Sciences, 196140 St. Petersburg, Russia
^{3} School of Space Research, Kyung Hee University, 446701 Yongin, Gyeonggi, Korea
^{4} Space Research Institute (IKI) of Russian Academy of Sciences, Profsoyuznaya St. 84/32, 117997 Moscow, Russia
Received: 10 September 2015
Accepted: 4 November 2015
Context. Despite intensive studies of kink oscillations of coronal loops in the last decade, a largescale statistically significant investigation of the oscillation parameters has not been made using data from the Solar Dynamics Observatory (SDO).
Aims. We carry out a statistical study of kink oscillations using extreme ultraviolet imaging data from a previously compiled catalogue.
Methods. We analysed 58 kink oscillation events observed by the Atmospheric Imaging Assembly (AIA) on board SDO during its first four years of operation (2010–2014). Parameters of the oscillations, including the initial apparent amplitude, period, length of the oscillating loop, and damping are studied for 120 individual loop oscillations.
Results. Analysis of the initial loop displacement and oscillation amplitude leads to the conclusion that the initial loop displacement prescribes the initial amplitude of oscillation in general. The period is found to scale with the loop length, and a linear fit of the data cloud gives a kink speed of C_{k} = (1330 ± 50) km s^{1}. The main body of the data corresponds to kink speeds in the range C_{k} = (800−3300) km s^{1}. Measurements of 52 exponential damping times were made, and it was noted that at least 21 of the damping profiles may be better approximated by a combination of nonexponential and exponential profiles rather than a purely exponential damping envelope. There are nine additional cases where the profile appears to be purely nonexponential and no damping time was measured. A scaling of the exponential damping time with the period is found, following the previously established linear scaling between these two parameters.
Key words: Sun: corona / Sun: oscillations / methods: observational
© ESO, 2016
1. Introduction
Waves and oscillations in the solar corona have been observed and studied for decades. Detections at extreme ultraviolet (EUV) and Xray wavelengths have confirmed theories and models related to magnetohydrodynamic (MHD) wave modes in magnetically and gravitationally structured plasmas (De Moortel & Nakariakov 2012). The importance of these waves and oscillations stems from their relation to the local plasma parameters of the medium allowing coronal seismology to be attempted (Stepanov et al. 2012a; Liu & Ofman 2014).
Kink oscillations of coronal loops have been intensively studied since their detection with the Transition Region And Coronal Explorer (TRACE; Handy et al. 1999) in 1999 (Aschwanden et al. 1999; Nakariakov et al. 1999). Prior to their detection they were the subject of a range of theoretical and numerical studies (e.g. Zaitsev & Stepanov 1982; Edwin & Roberts 1983; Roberts et al. 1984; Murawski & Roberts 1994). This phenomenon has been clearly observed with the spatial and temporal resolution of recent EUV imagers such as TRACE and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO; Lemen et al. 2012). Standing global modes induced by flaring activity are the most commonly detected form of kink oscillation (e.g. Nakariakov et al. 1999). Other detections have included higher spatial harmonics in coronal loops (e.g. Van Doorsselaere et al. 2007), their propagating form (Tomczyk et al. 2007), oscillations of polar plumes (Thurgood et al. 2014), propagating kink waves in streamer stalks (Chen et al. 2010, 2011) and kink waves in coronal jets (Vasheghani Farahani et al. 2009).
The period of a fundamental kink oscillation is given by the expression P = 2L/C_{k}, where L is the loop length and C_{k} is the kink speed, which is the phase speed of the kink wave (Zaitsev & Stepanov 1982; Edwin & Roberts 1983). The kink speed is dependant on the internal Alfvén speed (C_{A0}) and the density contrast between the ambient plasma and the loop (ρ_{e}/ρ_{0}); in the lowβ plasma limit this can be approximated as C_{k} = (2/(1 + ρ_{e}/ρ_{0}))^{1/2}C_{A0}. The Alfvén speed is given by , meaning that if a kink speed is measured and the density contrast can be estimated then the magnetic field in the loop, B_{0}, can be estimated (e.g. Nakariakov & Ofman 2001).
Standing kink oscillations can be induced by a flaregenerated blast wave (e.g. McLaughlin & Ofman 2008), by periodic driving from a coronal wave (Selwa et al. 2010), or loop contraction due to the reconfiguration of the active region magnetic field after an eruption (Russell et al. 2015). Alternative mechanisms such as Alfvénic vortex shedding (Nakariakov et al. 2009) have also been proposed. However, it has recently been shown that the vast majority of kink oscillations of coronal loops are excited by low coronal eruptions (LCE; Zimovets & Nakariakov 2015), discussed in more detail below.
The commonly applied explanation for the observed rapid damping of these oscillations is resonant absorption. The wave is considered to be a kink mode which evolves, because of a resonance at a certain thin loop layer, from a global kink mode to a local azimuthal motion (e.g. Ruderman & Roberts 2002; Goossens et al. 2002). The local motion cannot be observed directly and it undergoes its own damping. Other proposed mechanisms include phase mixing (e.g. Ofman & Aschwanden 2002), and viscous or dissipative damping (e.g. Nakariakov et al. 1999).
For standing kink modes, it has been shown that an exponential damping profile is expected (Ruderman & Roberts 2002). Numerical simulations from Pascoe et al. (2012) revealed a change in the characteristic damping profile of propagating transverse velocity displacements, which is directly applicable to standing kink modes. At the beginning of the oscillation a Gaussian profile is present, followed by an exponential damping profile; it has been suggested that the time at which the transition between the two damping profiles occurs depends on the transverse density structuring. Purely exponential, purely Gaussian, and a Gaussian followed by exponential damping are possible depending on the parameters of the loop and background plasma. More recently Pascoe et al. (2013) performed an analytical and parametric study of kink wave damping via mode coupling; this confirmed the existence of a Gaussian phase in the damping envelope, which was found to be more suitable for seismological application than the asymptotic exponential stage.
Ofman & Aschwanden (2002) studied scaling relations for different damping mechanisms, and determined that the observations considered were consistent with the phase mixing mechanism (Heyvaerts & Priest 1983). However, as has been discussed by Arregui et al. (2008), a comparison between the observed scaling and the linear scaling typical for resonant absorption requires that all loops have the same crosssectional structuring, and this means that the scaling from resonant absorption can depart from a linear dependence between damping time and period.
There have been various applications of kink oscillation observations to perform remote diagnostics of the parameters of the surrounding plasma (e.g. Stepanov et al. 2012b). By comparing the results of observations and modelling the absolute value of the loop’s internal magnetic field has been estimated (e.g. Nakariakov & Ofman 2001), the density scale height within the loop can be inferred (e.g. Andries et al. 2005), as can the characterisation of the height variation of the loop’s minor radius (e.g. Verth & Erdélyi 2008).
Recent observational studies of kink oscillations in the corona include Verwichte et al. (2013), where results were combined to study the scaling of the damping time with period, and put constraints on the loops density contrast and inhomogeneity layer thickness. In addition, a decayless regime of kink oscillations has been detected (Nisticò et al. 2013; Anfinogentov et al. 2013) and the statistical properties of this regime have been determined (Anfinogentov et al. 2015). A linear scaling of the period with loop length was obtained, confirming their interpretation of the oscillations as standing kink waves. Recently Zimovets & Nakariakov (2015) catalogued 169 individual decaying kink oscillations using SDO/AIA data, and utilised these observations to shed light on the excitation mechanism for the oscillations. They found that 95% of oscillations corresponded to displacement of the loop by a nearby lower coronal eruption/ejection, indicating this is the most common excitation mechanism. An LCE is the emergence of an expanding plasma structure (such as a flux rope, filament, or loop arcade) or plasma ejection (mainly along magnetic field lines) into the lower corona (roughly R< 1.4 R_{⊙}), usually – but not necessarily – accompanied by a flare and a subsequent CME.
The aim of our study is to use this catalogue to investigate the distributions of the parameters and to describe the initial excitation of each oscillation, the oscillation itself, and its damping. The determination of statistical relationships between the different parameters of kink oscillations of coronal loops is also performed. In Sect. 2 the data and analysis are described, in Sect. 3 the results are presented, and a discussion and conclusion are given in Sect. 4.
2. Observations and analysis
Decaying kink oscillations of coronal loops from the catalogue in Zimovets & Nakariakov (2015) are analysed using 171 Å data from SDO/AIA. A series of images was obtained for each of the 58 kink oscillation events listed using the provided date, oscillation time, and location to define a field of view and time interval, which were subsequently submitted as a SolarSoft (SSW) cutout request. The data cubes obtained had a time span of 30 or 45 min (which was extended for longperiod oscillations as required), the standard pixel size of 0.6 arcsec, and a temporal cadence of 12 s, or 24 s in some cases when the intermediate image was not returned owing to a short exposure time.
Movies created from the data cubes were initially inspected by eye, and loop oscillations with sufficient quality for timedistance (TD) analysis were noted. Timedistance maps were created by taking linear slits 5 pixels in width perpendicular to the oscillating loop and stacking the intensities along the slit (averaged over the width to increase the signaltonoise ratio) in time. An example of an active region is shown in Fig. 1; the slits used to create TD maps of different loops are overplotted. This process resulted in 127 TD maps, which were interpolated to an equispaced temporal grid of 12 s. The ends of the slits used are listed in Table A.1, along with the event number from Zimovets & Nakariakov (2015) and a loop ID. The first endpoint of each slit is overplotted on the disk in Fig. 2 in blue, and the average slit position for each event is plotted in red, the observed clustering is due to different events occurring within the same active region.
The projected loop length was estimated for each oscillating loop. The major radius (via the apparent loop height) or diameter (via the distance between footpoints) was measured by eye, depending on the orientation of the loop with respect to our line of sight (LoS), and a semicircular loop approximation was used, L = πR. The loop lengths are listed in Table A.1. In a few cases the loop length could not be estimated because the footpoint positions or height could not be determined.
For each TD map the amplitude of the initial displacement and initial oscillation amplitude were recorded. The initial displacement is defined as the difference between the initial loop position and the first maximum, and the initial amplitude is defined between the first maximum and minimum (see Fig. 3c). In this case the displacement of the loop’s upper edge was estimated. In addition to the number of oscillation cycles observed, the start time of the oscillation was also recorded, all listed in Table A.1.
Fig. 1 Active region from event number 40 from Table A.1. The three blue lines show some of the slits used to create timedistance maps for analysing the oscillations of the corresponding loops. 
Fig. 2 Slit positions (x1, y1 from Table A.1) used to produce the sample of timedistance maps to analyse kink oscillations of coronal loops, plotted as blue asterisks. The overplotted red circles show the average slit position for each event. 
Gaussian tracking of the loop centre to record the oscillating displacement was not appropriate in many cases owing to the overlap of multiple loops, or because only the edge of the oscillating loop was clearly defined, and so the oscillations were mapped out by taking a series of points along the centre or edge of the loops by hand, and an error for the points was defined based on the clarity of the time distance map, ±1 pixel in most cases. Three examples of TD maps are shown in Fig. 3, the points used to define the oscillation are overplotted in panel 3c, where the loop edge was mapped out and points were not taken when its position could not be reliably determined.
There is a certain degree of subjectiveness and error associated with the loop length measurement, the displacement measurements, and the points taken to map out the oscillation; however, the sample size is large enough that this will not affect the overall results, and only affects individual measurements.
The data points for each oscillation were detrended by fitting with a polynomial function of the form y = c_{0} + c_{1}t + c_{2}t^{2}, and subtracting this from the data. All fitting was performed with user defined functions and the IDL rountine mpfitexpr.pro. Fitting with a sinusoidal function of the form y = Asin(2πt/P + φ) was performed for each detrended oscillation with the period (P) as one of the free parameters. The best fitting period and the corresponding error were recorded for each welldefined oscillation.
Fig. 3 Three typical time distance maps, corresponding to loop 1 from event 32, loop 4 from event 40, and loop 2 from event 48. The fits correspond to the detrending polynomial added to the sinusoidal fit and multiplied by the exponential damping profile for panels a) and b). The red points in panel c) were taken by eye to map out the oscillation and used for the subsequent fitting. The vertical blue bars represent the measurement of the initial displacement (a_{0}) and the initial oscillation amplitude (a_{osc}) of the upper loop edge. 
To analyse the damping behaviour of the oscillations the absolute value of the detrended oscillatory signal was taken and scrutinised by eye. For over 50% of the TD maps clear damping could not be seen, or the number of oscillation cycles was not sufficient to perform fitting of the damping envelope. However, these oscillations are still clearly part of the decaying, rather than decayless, regime. For oscillations with a clear exponential decaying trend a weighted fit of the form A(t) = A_{0}e^{− t/τ} was performed on the maxima of the absolute value of the detrended signal. The damping time (τ) and the corresponding error were recorded. For cases where the damping was not observed for the whole duration of the signal, or where there were clearly nonexponential regions of the damping envelope, the fit was only made for the region which was approximated well by an exponential decay. By eye it was determined whether each damping profile was best described by a purely exponential profile (see Fig. 3 panel a), or a combination of both nonexponential and exponential profiles (see Fig. 3 panel b).
The result of the whole fitting process is overplotted on the TD maps in Fig. 3, where the detrending polynomial has been added to the sinusoidal fit multiplied by the damping profile for panels a and b. For panel a the damping profile was measured for the whole signal, for panel b it was only measured for the last 2 cycles.
3. Results
Fig. 4 Distribution of the measured initial displacement of 120 kink oscillations of coronal loops. The bin size is 2 Mm. 
The periods of 120 individual kink oscillations were determined from the 58 events analysed, and for 118 of these the corresponding loop lengths were estimated. From the oscillatory signals 52 exponential damping times were obtained. In addition to these measurements, other details of the oscillations were recorded and are listed in Table A.1, including whether damping was observed and if so whether it is best described by a purely exponential, or a combination of nonexponential and exponential damping profiles. The number of oscillation cycles was also recorded. In many cases this number was limited by the oscillatory signal becoming unclear, and not the damping reducing the oscillation to an undetectable amplitude.
3.1. Oscillation parameter histograms
Analysis of the amplitudes of the initial displacement and subsequent oscillation may allow inferences to be made about the excitation mechanism and may give details of the typical scales involved. The measured initial displacements ranged from 0.6 to 31.8 Mm. In Fig. 4 a histogram of the measured initial loop displacements is shown with a bin size of 2 Mm. The distribution peaks strongly at 1–3 Mm, and 59% of the measurements are covered by the range 1–5 Mm. The distribution after this range is more uniform, but the number of cases decreases towards the upper limit.
The initial oscillation amplitudes ranged from 0.5 to 27.6 Mm. In Fig. 5 a histogram of the measured initial oscillation amplitudes is shown, with a bin size of 2 Mm. The distribution peaks at 3–5 Mm, with 57% of the measurements lying in the range 1–5 Mm. The distribution again flattens and decreases towards the upper limit of the measurements.
In Fig. 6 the initial oscillation amplitude is plotted against the initial loop displacement. A rough correlation between these two parameters is observed. A linear fit of the data cloud that passes through the origin gives a gradient of 0.86 ± 0.01.
The measured oscillation periods ranged from 1.5 to 28 min. In Fig. 7 a histogram of the measured oscillation periods is shown. The distribution peaks at 4–7 min, and drops quickly to the maximum detected period of 28 min (not shown in the histogram). No periods below 1.5 min were recorded, and there is a decrease in occurrence approaching the lower periods. In Fig. 8 a histogram of the measured loop lengths is shown. The most common length is in the range 220–260 Mm, but there is a roughly even distribution between 140 and 460 Mm, which decreases above and below this range with minimum and maximum values of 77 and 596 Mm, respectively.
3.2. Dependence of the period on length
Fig. 5 Distribution of the measured initial oscillation amplitude of 120 kink oscillations of coronal loops, recorded from the first cycle of oscillations after the initial displacement. The bin size is 2 Mm. 
Fig. 6 Initial oscillation amplitude of 120 kink oscillations of coronal loops, plotted against the initial displacement of the loop position, both measured in Mm. A linear fit of the data passing through the origin is shown by the solid black line with a gradient of 0.86 ± 0.01. 
In Fig. 9 the period is plotted against the loop length, and the period clearly increases with the length. The period errors correspond to the scaled covariance from the period fitting. An unweighed linear fit was made because the errors alone do not reflect the distribution of the data. The variation of the density contrast and Alfvén speed between different loops and active regions also introduce an intrinsic spread to the data.
The black line corresponds to an unweighted linear fit of the data. The best fitting linear function is P [min] = (0.025 ± 0.001) L[Mm], where L is the loop length and P is the period, giving a kink speed of C_{k} = (1330 ± 50) km s^{1} from the empirically determined gradient and the equation P = 2L/C_{k}. The gradient can be varied to give upper and lower bounds to the data cloud, giving a kink speed range of C_{k} = (800−3300) km s^{1}. Calculation of the kink speed for each individual data point gives the distribution shown in the inset histogram. This has a most common value of 900–1100 km s^{1}, a peak value of 1340 ± 60 km s^{1}, and a Gaussian width of 620 ± 60 km s^{1}; the last two values were obtained by fitting the observed distribution with a Gaussian model.
3.3. Relationship between the damping time and period
In Fig. 10 the damping time is plotted against the period, and a correlation between the two parameters is obtained. The damping time errors correspond to the scaled covariance from the exponential damping fits. The statistics are limited to the cases where the damping time could be measured (see Table A.1), so the figure is less populated than Fig. 9. A weighted linear fit was made, corresponding to the solid black line. The best fitting linear function is τ = (1.53 ± 0.03)P, where τ is the damping time and P is the period.
In Fig. 10 the red circles correspond to damping profiles which were determined to be exponential by eye, and the blue squares correspond to damping profiles best described by a combination of a nonexponential and exponential profile. For the latter case, which corresponds to 21 of the measurements, the damping time is determined from the exponential part of the profile. This corresponds to the values “E” and “E,NE” in the column “Damping Profile” in Table A.1. No difference between the two cases is observed. There are nine additional cases where the profile appears to be purely nonexponential, and no damping time was measured, denoted in Table A.1 as “NE”.
Fig. 7 Distribution of the measured periods of 120 individual kink oscillations of coronal loops. The bin size is 1 min. 
Fig. 8 Distribution of the measured loop lengths for 118 individual coronal loops that undergo kink oscillations. The bin size is 20 Mm. 
Fig. 9 Period plotted against loop length for 118 kink oscillations of coronal loops. The solid black line corresponds to an unweighted linear fit of the data. The best fitting linear function is P [min] = (0.025 ± 0.001) L[Mm], giving a kink speed of C_{k} = (1330 ± 50) km s^{1} from the gradient. The dashed lines correspond to kink speeds of 800 and 3300 km s^{1} for the upper and lower lines, respectively. 
Fig. 10 Damping time plotted against period for 54 kink oscillations of coronal loops. The solid black line corresponds to a weighted linear fit of the data. The best fitting linear function is τ = (1.53 ± 0.03)P. The red circles correspond to damping envelopes best described by an exponential profile, and the blue squares correspond to those best described by a combination of nonexponential and exponential. 
4. Discussion and conclusions
In this paper we provide a comprehensive statistical analysis of decaying kink oscillations of coronal loops excited by flaring events observed with SDO/AIA at 171 Å.
4.1. Oscillation parameter histograms
Details of the initial loop displacement and the subsequent oscillation have been analysed. Figure 4 indicates that initiation of kink oscillations by lower amplitude initial displacements is more common, and Fig. 5 shows that relatively low initial oscillation amplitudes are also more common. A comparison between the two histograms mentioned leads to the conclusion that the initial displacement prescribes the initial amplitude of oscillation in the majority of cases. However, there are cases where a large initial displacement leads to a new equilibrium position for the loop where it oscillates with an amplitude much smaller than the initial displacement, as has been recently discussed by Russell et al. (2015), but the exact cause of this effect is still unknown. Cases where a small initial displacement results in a larger amplitude oscillation are also seen. A simple explanation is that the passage of the LCE has left the external plasma and magnetic field more rarefied, and the perturbed loop passes through its initial equilibrium with an amplitude that is greater than the initial displacement. A linear correlation between these two parameters is seen in 6, with spreading due to measurement errors, but also reflecting the different regimes discussed above.
The above results are subject to a LoS effect. The histograms in Figs. 3 and 4 will include an effect in the distribution relating to how the observed kink oscillations are distributed over the LoS angles. If all the initial displacements and oscillation amplitudes were equal, distributions due to the different LoS angles varying the measured values would still be obtained. It can be seen in Fig. 2 that loop positions are mostly offlimb near the equator, so northsouth oscillation polarisations are better detected. Further work is required to remove the effect of the varying LoS angles from the measurements and to obtain the true displacements and oscillation amplitudes.
Some inferences can be made from the period distribution shown in Fig. 7. The dropoff in the occurrence of higher periods likely reflects a physical dropoff, as larger length loops may be less likely to be formed and are also more difficult to detect. The decrease in the distribution for lower periods may include an observational bias, as oscillations of shorter loops are more difficult to observe, but may also reflect the excitation mechanism. If, as proposed by Zimovets & Nakariakov (2015), excitation due to LCE perturbations is the dominant mechanism, this should preferentially excite loops above a certain height that have correspondingly longer lengths. The distributions of the periods and loop lengths differ, and this is likely due to the loop length estimations becoming more inaccurate for longer loops, as well as variations in the additional parameters that prescribe the period of oscillation.
If general statistics of the loop lengths were available they could be used to normalise the loop length distribution in Fig. 8. This would allow us to determine whether kink oscillations occur in all loops with the same probability or if loops of certain lengths are more likely to undergo a kink oscillation. Line of sight effects should also be included in such a study, and the varying ellipticity of the loops themselves.
4.2. Dependence of the period on length
The period scales with the loop length as expected, and the best fitting kink speed of C_{k} = (1330 ± 50) km s^{1} is in agreement with previous results. The range of kink speeds from the main data cloud (C_{k} = (800−3300) km s^{1}) should correspond to the spread of loop density contrasts, transverse density profiles, and Alfvén speeds. There are some points in Fig. 9, corresponding to short loops with larger periods than expected from the main body of the data, that give a lower kink speed value. These data points may correspond to loops and active regions with significantly lower Alfvén speeds or greater density contrasts. More statistics and analysis of the active region parameters are required to determine whether these points have a physical explanation or are due to measurement errors and/or random spread.
Comparisons can be made with the results of a similar study which focussed on the decayless regime of kink oscillations (Anfinogentov et al. 2015). They obtained a similar scaling of period with loop length, but with a lower gradient and, therefore, a higher best fitting kink speed of C_{k} = (1850 ± 70) km s^{1}. The period distribution they present is similar to the distribution we obtain, but it peaks at a lower value and their loop length distribution is significantly less uniform than ours. These differences may be due to a selection effect from their study because the data spanned one month, whereas the data presented here span four years. It may also reflect the different driving mechanisms; in particular, the excitation mechanism for the decayless regime remains unknown. The discrepancy between the statistics of decaying and decayless kink oscillations should be further investigated when a larger set of events becomes available.
4.3. Relationship between the damping time and period
Our results on the linear scaling of the damping time with the oscillation period are mainly consistent with the previously obtained results (Aschwanden et al. 2002; Verwichte et al. 2013). The spread of the data makes it unreasonable to make inferences about the damping mechanism from the observed dependence. In particular it is not possible to discriminate between a linear or power law dependence, and such an approach requires consideration of the influence of the varying crosssectional loop structuring and other parameters.
Distinguishing between the exponential damping time measured from oscillations that showed a purely exponential profile and those with a combination of a nonexponential and exponential profiles did not reveal any systematic differences. Only the exponential region was fit if there was also a nonexponential region present. This indicates that the exponential stage of the damping is the same whether a nonexponential stage is present or not. Our findings indicate that in some cases the exponential fitting of the kink oscillation damping is not sufficient to capture the whole damping profile. The amplitude of some oscillations remains approximately constant for a significant period of time, it has not been determined whether this corresponds to a slowly decaying Gaussian profile or to a periodic driver sustaining a constant amplitude oscillation.
The detected clearly nonexponential sections of the damping envelopes may be better approximated by a Gaussian profile. If this is confirmed to be the case then this is evidence for the Gaussian damping regime discussed by Pascoe et al. (2012). This will be the subject of further study as the detection of Gaussian damping envelopes allows new seismology to be performed and comparisons with the theoretical predictions to be made.
4.4. Conclusion
We summarise our main findings as follows:

The initial loop displacement prescribes the initial oscillationamplitude in the majority of cases.

The period scales linearly with the loop length, as expected, and a kink speed of C_{k} = (1330 ± 50) km s^{1} is obtained; the majority of the data points are in the range (800–3300) km s^{1}, following a Gaussian distribution.

A linear scaling of the damping time with period is observed and nonexponential damping profiles have been detected.
In conclusion, a statistically significant number of individual kink oscillations have been analysed and histograms of the measured parameters have allowed insightful inferences to be made. Details of the distribution of amplitudes, periods, and loop lengths may be useful when considering the observational capabilities of future instruments. In addition, the scaling between different parameters has been studied and the damping behaviour has been characterised, both of which – after further work – may allow seismological inferences and measurements to be made.
Acknowledgments
The work was supported by the European Research Council under the SeismoSun Research Project No. 321141 (CRG, VMN), and the STFC consolidated grant ST/L000733/1 (GN, VMN). IVZ was partially supported by the RFBR (research project No. 153221078). The data is used courtesy of the SDO/AIA team.
References
 Andries, J., Arregui, I., & Goossens, M. 2005, ApJ, 624, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Anfinogentov, S., Nisticò, G., & Nakariakov, V. M. 2013, A&A, 560, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Anfinogentov, S. A., Nakariakov, V. M., & Nisticò, G. 2015, A&A, 583, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Arregui, I., Ballester, J. L., & Goossens, M. 2008, ApJ, 676, L77 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., Fletcher, L., Schrijver, C. J., & Alexander, D. 1999, ApJ, 520, 880 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., de Pontieu, B., Schrijver, C. J., & Title, A. M. 2002, Sol. Phys., 206, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, Y., Song, H. Q., Li, B., et al. 2010, ApJ, 714, 644 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, Y., Feng, S. W., Li, B., et al. 2011, ApJ, 728, 147 [NASA ADS] [CrossRef] [Google Scholar]
 De Moortel, I., & Nakariakov, V. M. 2012, Roy. Soc. London Philosoph. Trans. Ser. A, 370, 3193 [Google Scholar]
 Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Goossens, M., Andries, J., & Aschwanden, M. J. 2002, A&A, 394, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Handy, B. N., Acton, L. W., Kankelborg, C. C., et al. 1999, Sol. Phys., 187, 229 [NASA ADS] [CrossRef] [Google Scholar]
 Heyvaerts, J., & Priest, E. R. 1983, A&A, 117, 220 [NASA ADS] [Google Scholar]
 Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, W., & Ofman, L. 2014, Sol. Phys., 289, 3233 [NASA ADS] [CrossRef] [Google Scholar]
 McLaughlin, J. A., & Ofman, L. 2008, ApJ, 682, 1338 [NASA ADS] [CrossRef] [Google Scholar]
 Murawski, K., & Roberts, B. 1994, Sol. Phys., 151, 305 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., & Ofman, L. 2001, A&A, 372, L53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., & Davila, J. M. 1999, Science, 285, 862 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Nakariakov, V. M., Aschwanden, M. J., & van Doorsselaere, T. 2009, A&A, 502, 661 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nisticò, G., Nakariakov, V. M., & Verwichte, E. 2013, A&A, 552, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ofman, L., & Aschwanden, M. J. 2002, ApJ, 576, L153 [NASA ADS] [CrossRef] [Google Scholar]
 Pascoe, D. J., Hood, A. W., de Moortel, I., & Wright, A. N. 2012, A&A, 539, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pascoe, D. J., Hood, A. W., De Moortel, I., & Wright, A. N. 2013, A&A, 551, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roberts, B., Edwin, P. M., & Benz, A. O. 1984, ApJ, 279, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S., & Roberts, B. 2002, ApJ, 577, 475 [NASA ADS] [CrossRef] [Google Scholar]
 Russell, A. J. B., Simões, P. J. A., & Fletcher, L. 2015, A&A, 581, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Selwa, M., Murawski, K., Solanki, S. K., & Ofman, L. 2010, A&A, 512, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stepanov, A. V., Zaitsev, V. V., & Nakariakov, V. M. 2012a, Stellar Coronal Seismology as a Diagnostic Tool for Flare Plasma (WileyVCH Verlag GmbH & Co. KGaA) [Google Scholar]
 Stepanov, A. V., Zaitsev, V. V., & Nakariakov, V. M. 2012b, Physics Uspekhi, 55, A04 [NASA ADS] [CrossRef] [Google Scholar]
 Thurgood, J. O., Morton, R. J., & McLaughlin, J. A. 2014, ApJ, 790, L2 [NASA ADS] [CrossRef] [Google Scholar]
 Tomczyk, S., McIntosh, S. W., Keil, S. L., et al. 2007, Science, 317, 1192 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Van Doorsselaere, T., Nakariakov, V. M., & Verwichte, E. 2007, A&A, 473, 959 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vasheghani Farahani, S., Van Doorsselaere, T., Verwichte, E., & Nakariakov, V. M. 2009, A&A, 498, L29 [Google Scholar]
 Verth, G., & Erdélyi, R. 2008, A&A, 486, 1015 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verwichte, E., Van Doorsselaere, T., White, R. S., & Antolin, P. 2013, A&A, 552, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zaitsev, V. V., & Stepanov, A. V. 1982, Sov. Astron. Lett., 8, 132 [NASA ADS] [Google Scholar]
 Zimovets, I. V., & Nakariakov, V. M. 2015, A&A, 577, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Table
120 coronal loop kink oscillations detected with SDO/AIA and their measured parameters.
All Tables
120 coronal loop kink oscillations detected with SDO/AIA and their measured parameters.
All Figures
Fig. 1 Active region from event number 40 from Table A.1. The three blue lines show some of the slits used to create timedistance maps for analysing the oscillations of the corresponding loops. 

In the text 
Fig. 2 Slit positions (x1, y1 from Table A.1) used to produce the sample of timedistance maps to analyse kink oscillations of coronal loops, plotted as blue asterisks. The overplotted red circles show the average slit position for each event. 

In the text 
Fig. 3 Three typical time distance maps, corresponding to loop 1 from event 32, loop 4 from event 40, and loop 2 from event 48. The fits correspond to the detrending polynomial added to the sinusoidal fit and multiplied by the exponential damping profile for panels a) and b). The red points in panel c) were taken by eye to map out the oscillation and used for the subsequent fitting. The vertical blue bars represent the measurement of the initial displacement (a_{0}) and the initial oscillation amplitude (a_{osc}) of the upper loop edge. 

In the text 
Fig. 4 Distribution of the measured initial displacement of 120 kink oscillations of coronal loops. The bin size is 2 Mm. 

In the text 
Fig. 5 Distribution of the measured initial oscillation amplitude of 120 kink oscillations of coronal loops, recorded from the first cycle of oscillations after the initial displacement. The bin size is 2 Mm. 

In the text 
Fig. 6 Initial oscillation amplitude of 120 kink oscillations of coronal loops, plotted against the initial displacement of the loop position, both measured in Mm. A linear fit of the data passing through the origin is shown by the solid black line with a gradient of 0.86 ± 0.01. 

In the text 
Fig. 7 Distribution of the measured periods of 120 individual kink oscillations of coronal loops. The bin size is 1 min. 

In the text 
Fig. 8 Distribution of the measured loop lengths for 118 individual coronal loops that undergo kink oscillations. The bin size is 20 Mm. 

In the text 
Fig. 9 Period plotted against loop length for 118 kink oscillations of coronal loops. The solid black line corresponds to an unweighted linear fit of the data. The best fitting linear function is P [min] = (0.025 ± 0.001) L[Mm], giving a kink speed of C_{k} = (1330 ± 50) km s^{1} from the gradient. The dashed lines correspond to kink speeds of 800 and 3300 km s^{1} for the upper and lower lines, respectively. 

In the text 
Fig. 10 Damping time plotted against period for 54 kink oscillations of coronal loops. The solid black line corresponds to a weighted linear fit of the data. The best fitting linear function is τ = (1.53 ± 0.03)P. The red circles correspond to damping envelopes best described by an exponential profile, and the blue squares correspond to those best described by a combination of nonexponential and exponential. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.