The Effect Of Cooling On Driven Kink Oscillations Of Coronal Loops

Ever since their detection two decades ago, standing kink oscillations in coronal loops have been extensively studied both observationally and theoretically. Almost all driven coronal loop oscillations (e.g., by flares) are observed to damp through time often with Gaussian or exponential profiles. Intriguingly, however, it has been shown theoretically that the amplitudes of some oscillations could be modified from Gaussian or exponential profiles if cooling is present in the coronal loop systems. Indeed, in some cases the oscillation amplitude can even increase through time. In this article, we analyse a flare-driven coronal loop oscillation observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA) in order to investigate whether models of cooling can explain the amplitude profile of the oscillation and whether hints of cooling can be found in the intensity evolution of several SDO/AIA filters. During the oscillation of this loop system, the kink mode amplitude appears to differ from a typical Gaussian or exponential profile with some hints being present that the amplitude increases. The application of cooling coronal loop modelling allowed us to estimate the density ratio between the loop and the background plasma, with a ratio of between 2.05-2.35 being returned. Overall, our results indicate that consideration of the thermal evolution of coronal loop systems can allow us to better describe oscillations in these structures and return more accurate estimates of the physical properties of the loops (e.g., density, scale height, magnetic field strength).


Introduction
Standing kink oscillations were first observed in coronal loops by Aschwanden et al. [1999] and Nakariakov et al. [1999] using high-resolution imaging data collected by the Transition Region And Coronal Explorer (TRACE; Handy et al. 1999).Those authors found that a flare in the local Active Region (AR) caused the magnetic field guide of the coronal loop to shake from side-to-side in a manner analogous to oscillations of a guitar string.One of the most interesting aspects of kink oscillations in coronal loops was their rapid damping profiles, with many examples of flare-driven coronal loop oscillations damping to sub-resolution spatial scales within two or three periods.Such damping has been explained through a number of physical mechanisms, such as resonant absorption (e.g., Ruderman and Roberts 2002;Goossens et al. 2002), phase mixing, and foot-point damping.
Over the past decade, data from the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA; Lemen et al. 2012) has been used to perform statistical analyses of coronal loop oscillations (e.g., Zimovets and Nakariakov 2015;Goddard et al. 2016).Specifically, Goddard et al. [2016] analysed the damping profiles of 58 coronal loops finding that the amplitudes through time of the majority of oscillations could be modelled by exponential or quasi-exponential (close to but not exactly exponential) profiles.More recently, Pascoe et al. [2016] suggested that Gaussian profiles may model damping profiles of coronal loops as well, if not better, than exponential profiles in a large number of cases.It was, however, pointed out by Goddard et al. [2016] that some events were non-exponential (and also likely non-Gaussian) meaning, for some coronal loop oscillations, additional physics must be employed to explain the amplitude profiles.Interestingly, it has been shown analytically that it is possible that some driven coronal loops could oscillate and damp in a manner not adequately modelled by Gaussian nor exponential profiles through time.
As the majority of modelling of coronal loop oscillations is conducted using static background parameters (e.g., temperature, density), one opportunity to consider additional physics theoretically is to incorporate some form of time-dependence (see, for example, Dymova and Ruderman 2005;Al-Ghafri and Erdélyi 2013;Erdélyi et al. 2014).Morton and Erdélyi [2009] considered the damping of coronal loops due to cooling through time and found, that for typical oscillatory periods, cooling could play a key role in explaining observed damping profiles (this was shown further in Morton and Erdélyi 2010).The idea that amplification of coronal loop oscillations could occur due to cooling within coronal loops that contained flow was first suggested by Ruderman [2011a].This work was expanded upon in Ruderman [2011b] with the inclusion of a resonant layer, where it was found that cooling could cancel out the damping due to resonant absorption in some cases.The effects of loop expansion were considered for cooling coronal loops when no resonant layer was observed by Ruderman et al. [2017].Again it was found that amplification of the coronal loop oscillations could occur.
A more complete analysis comprising of loop expansion, resonant absorption, and cooling was recently conducted by Shukhobodskiy and Ruderman [2018] and Shukhobodskiy et al. [2018] for time independent and time dependent densities, respectively.Those authors found that in some cases, the combined effects of expansion and cooling could dominate over resonant absorption leading to a slowing down of the damping of the oscillations or even, in some cases, a brief amplification in kink mode coronal loop oscillations.Such amplification would only be short-lived and would be followed by the continued decay of the oscillation to sub-resolution levels.Additionally, assuming that both the external and internal densities were longitudinally stratified, the authors showed that the ratio of the frequency of a fundamental mode on the decrement of the kink oscillation is independent on the particular form of the density profile.A similar result was obtained previously by Dymova and Ruderman [2006] for kink oscillations in non-expanding magnetic flux tubes.It is possible, therefore, that accounting for cooling when modelling coronal loop oscillations may provide better fits for coronal loop oscillations, especially when increases in amplitude are evident.Additionally, important seismological information could be obtained from the system (such as the density ratio between the loop and the background) if one were to consider cooling which could help improve future numerical modelling.
The development of modelling of time-dependent coronal loop temperatures has a strong foundation in observations.Numerous authors have discussed the temperature evolution of coronal loop arcades, with cooling often being inferred through the observed progression of loops from hot to cold channels through time (see, for example, Winebarger and Warren 2005;López Fuentes et al. 2007;Aschwanden and Terradas 2008).It is well known that coronal loops can cool quickly, over the course of two or three oscillatory periods (e.g., Aschwanden and Terradas 2008), through processes such as the thermal instability, which can often lead to the occurrence of coronal rain (see Antolin et al. 2015 and references within).Such cooling means the application of theories which are not magnetohydrostatic, such as those recently developed by Shukhobodskiy et al. [2018], are important for further understanding and modelling coronal loop oscillations more generally.It should be noted that Morton and Erdélyi [2010] did conduct an application of theoretical work to observed kink oscillations, however, those authors did not include effects which could lead to amplification in their model.
In this article, we apply the theoretical models developed by Shukhobodskiy et al. [2018] to the amplitude profile of a kink oscillation in a coronal loop within AR 11598.We aim to showcase the seismological potential of cooling models through the inference of the density ratio between the loop and the background plasma.Our work is setup as follows: In Sect. 2 we describe the data studied here and the data analysis methods; In Sect. 3 we present our results including the application of the theoretical modelling and any evidence for cooling in the AR; In Sect. 4 we draw our conclusions and provide some suggestions about future work.

Data And Feature Selection
The data studied in this article were sampled by the SDO/AIA instrument on the 20th October 2012 between 17:50 UT and 18:50 UT.The 300 ′′ ×300 ′′ field-of-view (FOV) analysed here was focused on AR 11598 and was initially centred on co-ordinates of xc=−959 ′′ , yc=−85 ′′ .Four channels (namely the 304 Å, 171 Å, 193 Å, and 131 Å) were downloaded for analysis, using the ssw cutout service.proroutine, meaning we are able to make inferences about the thermal evolution of the loop system.The loop system is only weakly detectable in the 131 Å images and, as such, these data are not studied in detail in the remainder of this article.These data do, however, indicate that the coronal loop existed at non-flaring temperatures.Data in the UV sampled by the SDO/AIA instrument have a typical cadence of 12 s and a pixel scale of 0.6 ′′ .As the oscillations analysed here were driven by a flare in the same AR as the coronal loops, some frames were returned with a lower exposure time meaning the loop was difficult to detect due to the lower signal to noise ratio.For these frames, synthetic filling data was created with the intensity of each pixel being taken as the average intensity for that pixel in the frame before and the frame after.This should have no influence on our results as the loop analysed here oscillates with periodicities well above the cadence of that data.This coronal loop oscillation was previously discussed by Goddard et al. [2016] (Event 40, Loop 2) and its damping was classified as being best described by a combination of exponential and non-exponential fitting.This event was further analysed by Pascoe et al. [2016] and Pascoe et al. [2017] who suggested that the presence of multiple harmonics could explain the complicated amplitude profiles observed within this event.
The general evolution of the loop system analysed here is plotted in Fig. 1 for the 304 Å (top row), 171 Å (middle row), and 193 Å (bottom row) channels.At 17:50 UT, the coronal loop system can be detected in both the 171 Å and 193 Å channels, with the white line over-laid on the 171 Å images indicating the location of the slit studied here.The location of the flare in the AR, which occurred at approximately 18:08 UT, is pin-pointed by the typical disturbance patterns on the 193 Å image in the third column.In the southern part of the FOV (around 100 ′′ south of the oscillation analysed here), a large coronal loop arcade can be observed in all panels up until 18:19 UT.Over the subsequent 30 minutes (between 18:20 UT and 18:50 UT), however, the loop system completely fades from view in the coronal channels and large amounts of coronal rain are detected in the 304 Å data, draining material from the loops.The loop system completely fades from view by 18:49 UT.

Tracking The Loop Displacement And Model Fitting
In Fig. 2, we plot a zoomed in FOV of the oscillation analysed here for five time-steps (top row).The white line over-laid on theses images indicates the location of the slit analysed in this article.The co-ordinates of this slit are the same as for Loop 40, Event 2 included in the table of Goddard et al. [2016].This event was studied subsequently by Pascoe et al. [2016] and Pascoe et al. [2017].Following the construction of the time-space diagram from 171 Å images (plotted in the bottom panel of Fig. 2), we applied a Canny edge-detection algorithm to obtain approximations for the coronal loop boundaries.The output for this edge detection routine is overlaid on the time-space diagram in Fig. 3 (blue lines).In order to model the oscillation of the loop axis, we assumed that the displacement is guided by midpoints between the boundaries of the upper loop as shown by the red dots in Figure 3.If we assume that the loop is cylindrical and radially symmetric, then the errors associated by this fitting can be assumed to be less than 1 pixel.Understanding whether such assumptions are justified would require much higher resolution data and so should be studied in the future.As the loop axis returned through time is qualitatively similar to that plotted in Fig. 1 of Pascoe et al. [2016] for the same event, in the following we assume that this oscillation profile is accurate and, therefore, neglect any further potential errors in our fitting.Ignoring any further errors introduced through the loop fitting is justified in this case, as our study aims to show, in principle, that including the effects of cooling could have important implications for coronal loop modelling.A larger statistical study conducted using a variety of fitting methods in the future will be required to fully understand the prevalence and importance of coronal loop cooling on kink-mode oscillations.
The oscillation extracted from this time-space diagram is plotted with red dots in Fig. 4. The approximate background trend of the loop was modelled by a polynomial of the 8th degree on all obtained data points (i.e., from frame 0 to frame 300) and is over-laid for the region of interest (between frames 130 and 230) as the green line in Fig. 4. The 8th order polynomial was selected as it best tracked the background of the amplitude profile throughout the entire timeseries (including the parts where no oscillation was detected).It should be noted that a number of background trends with different orders were tested with each returning similar seismological results.This is because the ratio of the internal and external densities, which we can calculate from this model, is dependent only on where, in this case, the increment to the amplitude is (which did not change for any background trend), not how large it is.This will be discussed in further detail later in the article.
The numerical fitting of the red points is done by considering the summation of four equations of Gaussian form: where Ai, µi and σi are variables to be fitted for each peak, t is time (in frames), and f (amp) is the final fitted function.The benefit of such fits over a typical sinusoidal fit is that no periodicity is assumed a priori.This fitting was completed using the Wolfram Mathematica 11.3 procedure NonlinearModelFit (which guarantees continuity of the summed functions) and is plotted as the blue line over-laid on Fig. 4. It is clear that this blue line accurately models the detected oscillation, however, it does slightly under-estimate the returned amplitude between frames 170 and 200.Therefore, we stress that the increase in amplitude detected during this oscillation (discussed subsequently) will be under-estimated by our model.

Theoretical Modelling And Observed Amplitude Profiles
In this section, we fit the observed oscillatory profile plotted in Fig. 4 using the theoretical model proposed by Shukhobodskiy et al. [2018].The model proposed by those authors consisted of a straight magnetic flux tube with varying cross section along its length.The tube was modelled with three layers: the core which contained an arbitrary flow and oscillated as a solid body; a transitional layer (or annulus) with monotonically decaying density from the core to the external layer; and the surrounding background plasma.See Fig. 5 for a detailed schematic of this model.Ruderman et al. [2017] obtained the governing equation for such a model under the assumptions of a thin tube with thin boundaries, by considering jumps of displacement and pressure across the transitional layer.However, this equation was not closed as jump conditions were not defined in terms of the displacement.Later, Shukhobodskiy and Ruderman [2018] and Shukhobodskiy et al. [2018], closed that system for the time independent and time dependent densities, respectively.
In the model applied here, we consider a loop of half-circular shape, where the curvature affects only the density distribution.The temperature is modelled to be exponentially decaying in time (similar to the profiles studied in, for example, Aschwanden and Terradas [2008]; Morton and Erdélyi [2010]; Ruderman [2011b]; Ruderman et al. [2017]) and is approximated by: where T0 is the constant external temperature and t cool is the cooling time (assumed to be the total lifetime of the oscillation here).Additionally, we describe the variation of the loop crosssection with height, z, in a manner similar to that discussed by Ruderman et al. 2008Ruderman et al. , 2017;;Shukhobodskiy et al. 2018.This takes the form: (3) where R f is the radius of the magnetic flux tube at the foot-points, λ = R(0)/R f is the expansion factor and Lc is an arbitrary constant.Please note that Lc should be selected such that the expansion factor can be consistent with observations of coronal loops (1 < λ < 1.5).We also assume, similarly to Ruderman and Roberts 2002;Goossens et al. 2002, that the density in the transitional layer has a linear profile modelled by: where l is a constant determining the width of the transition layer and ρi and ρe are the internal and external densities, respectively.Using these equations, Shukhobodskiy et al. [2018] showed that the dimensionless amplitude, A where A(0) = 1, of the kink mode may be approximated by: where: and: Here, µ0 is the magnetic permeability of the free space, ω is the oscillation frequency, H0 is the external scale height, and ρ f and B f are the density of plasma and magnetic field strength at the foot-points inside the loop.Additionally, X is determined by the following boundary value problem: ) where X is the function of Z defined by the boundary values: In order to obtain results comparable with observed amplitude profiles, we set At = A(0)A Ob (0), where A Ob (0) is the measured initial amplitude of the observed oscillation, At is the factor by which the dimensionless amplitude is scaled, and set L/Lc = 6 (similar to the values used by Ruderman et al. 2008Ruderman et al. , 2017;;Shukhobodskiy et al. 2018).This allows us to obtain results for loop expansion factors in the range of 1-1.5, consistent with values discussed in the literature for coronal loops (see, for example, Klimchuk 2000;Watko and Klimchuk 2000).In Fig. 6, we plot the absolute displacement of the observed oscillation through time by subtracting the fitted displacement (blue line in Fig. 4) from the background trend (green line in Fig. 4).We define the initial amplitude of the oscillation, A Ob (0), by the first peak in Fig. 6.The observed amplitude profile is then constructed by calculating all values of the peaks in Fig. 6 and using the Interpolation procedure in Wolfram Mathematica 11.3.
The most interesting region of Fig. 6 (denoted by the dashed black box over-laid between frames 160-185) is plotted in the zoomed cut-out in the top right corner of the panel.It is immediately clear that peak four is slightly larger than peak three indicating an amplitude increase potentially occurs through time.Such an increase would be in agreement with the theoretical model of Shukhobodskiy et al. [2018] which suggests that cooling in the loop system can lead to wave amplification.We note that cooling does not strictly lead to an amplification of the oscillation, it only modifies the damping profile from Gaussian or exponential.Therefore, although the amplitude increase plotted here is sub-resolution, the fact that the oscillation is not damping in a manner consistent with a Gaussian or exponential profile is enough to suggest some sort of cooling may be occurring.It should also be noted that larger amplitude increases (greater than 1 pixel) were found for lower order polynomials, however, we only focus on the 8th order fitting here as that is sufficient as a proof of concept.The increase in the amplitude can be clearly seen in Fig. 7, where the blue line plots a fit between the measured peak amplitudes (green circles) through time.The red line over-laid on Fig. 7 indicates the line of best fit calculated by solving the system of equations ( 2)-(10) numerically.A maximum χ 2 value was found by looping over all variables, hence this is numerically expensive.It was assumed that the cooling starts at the first peak and ends in the last peak of Fig. 6 (i.e., that t cool is equal to the lifetime of the oscillation).
The software used to obtain these theoretical solution is Wolfram Mathematica 11.3.The numerical procedure for obtaining the the- ¢ t ( t, r, z)   i (t, r, z) £ e (t, r, z)   B (r, z) U (t, r, z) lR(z)

L
Figure 5: Equilibrium configuration of theoretical model.The black lines with labels r, φ and z correspond to the cylindrical polar coordinates system.B(r, z) is the background magnetic field, U(t, r, z) is the background flow parallel to the magnetic field, R(z) is the radius of the tube, l is a constant determining the thickness of the transitional layer, ρi(t, r, z) is the density in the internal region, ρe(t, r, z) is the density in external region, ρt(t, r, z) is the monotonically decaying density in the transitional layer from ρi(t, r, z) to ρe(t, r, z), L is the length of tube for the standing waves and the characteristic value of the wavelength for propagating waves.
oretical results may be summarised as follows: To obtain the solution to the boundary value problem, Eqs. ( 9) and (10), we used the NDEigensystem procedure; Eq. ( 6) was then integrated numerically using the NIntegrate procedure subject to the GlobalAdaptive method (which uses various numerical integration methods and chooses the most accurate and fastest version for the particular problem).Finally we obtained the amplitude of the oscillation by using the NDSolve procedure.
The application of this model revealed several interesting effects with regards seismology of the loop system.Firstly, it was found that ζ (the density ratio at the foot-points of the loop at t = 0) and κ (the scale height) determine the position of the turning points of the amplitude profile (i.e., where the gradient of the amplitude profile is zero).However, the effect of κ on the position of the turning point reduces as the value of κ itself increases.As a result, we can neglect this effect for sufficiently long loops with κ > 1.6.Therefore, by minimising the difference between the turning points of the theoretical and observed amplitude profiles we can determine an approximate value of ζ for this system without knowing any other background parameters.For the example studied here, fitting a density ratio of between 2.05 and 2.35 provides sufficiently good approximations for the position of the local amplitude increase.As this work aims to provide a proof-of-concept of the application of cooling theory to coronal loop oscillations and as the density ratio is only dependent on where the amplitude deviates from a Gaussian or exponential decay, we neglect any errors in the measured amplitude (green dots) here; however, the effects of any errors on our seismological should be analysed using a larger statistical sample in the future.The differences in nature between the red and blue curves are due to the spline fitting used here.Secondly, the values of κ, α, and λ determine the speed of the amplitude decay.Unfortunately, for the example studied here, various sets of these parameters return similar shapes and amplitude profiles meaning we are unable to provide accurate seismological estimates at this stage.A larger statistical sample of cooling events and further numerical work will be required to attempt to estimate such values in the future.

Evidence For Cooling In The AR
Evidence of cooling in the coronal loop system can be inferred through analysis of the intensities within SDO/AIA imaging channels.For the loop studied here, no evidence of coronal rain formation was found in the 304 Å filter indicating that catastrophic cooling likely did not take place within this loop.The loop intensity within the 171 Å and 193 Å filters did appear to decrease slightly (potentially below any level of significance) during the oscillation, however, whether this was due to cooling, some line-of-sight effects (e.g., the supposition of multiple loops within one pixel), or purely noise is unknown.The flare within the local AR lead to large swathes of saturation and fringing patterns on the 193 Å data meaning we were unable to study the temporal evolution of the loop in detail in this wavelength during the period of interest for cooling for this event.Future work should aim to analyse a larger statistical sample of coronal loop oscillations with non-Gaussian and non-exponential damping profiles in order to detect whether direct signatures of cooling can be found within the loops themselves.
Significant cooling can be inferred elsewhere in this AR during the course of this oscillation, however.Large amounts of coronal rain can be observed in the 304 Å channel in the loop system approximately 100 ′′ south of the loop analysed here potentially indi- The application of the theories tested here on loops which are catastrophically cooling would be an interesting project for the future.

Discussion And Conclusions
In this article, we have studied a flare-driven kink oscillation in a coronal loop observed in AR 11598 by the SDO/AIA instrument (Fig. 1).This oscillation corresponded to Loop 40, Event 2 from Goddard et al. [2016].The specific slit analysed here is indicated by the white line over-laid on the top row of Fig. 2, with the returned time-space diagram being plotted in the bottom panel of Fig. 2. In Fig. 3 the time-space diagram is plotted with the output from a Canny edge-detection algorithm over-laid in blue.The midpoints between the edge-detection outputs (red dots) were assumed to track the displacement of the centre of the loop as it oscillated.The returned oscillation was qualitatively similar to that returned by Pascoe et al. [2016] for the same event indicating that our fitting method was sound.This oscillation was then fitted with a summation of four Gaussian functions (corresponding to the four peaks in the oscillation) of the form Eq. 1 and an 8-th order background trend removed (the blue and green lines in Fig. 4, respectively).Again, it should be emphasised that numerous background trends were tested in order to assure that we were not introducing important effects in our subsequent analysis.
Once the oscillation and the background had been fitted, we removed the background trend from the oscillation and plotted the absolute values of the kink mode (see Fig. 6) through time.The cut-out over-laid on Fig. 6 clearly highlights a deviation in oscillatory amplitude from a typical Gaussian or exponential damping between frames 160 and 185 (18:22 UT -18:27 UT).Indeed, there are even hints that an amplitude increase could be present, however, this is small.It should be noted that some background trends returned profiles with amplitude increases of over 1 pixel.Such an amplitude increase would obviously not be expected from a typical Gaussian or exponential decay (as has been considered previously by, for example, Goddard et al. 2016;Pascoe et al. 2016).However, models proposed by Shukhobodskiy and Ruderman [2018] and Shukhobodskiy et al. [2018] which consider cooling within the system can account for amplitude increases.Therefore, in addition to the multiple harmonics scenario suggested by Pascoe et al. [2017], future work should consider cooling when trying to understand complex amplitude profiles in coronal loop oscillations.
Although there were some hints that the loop faded slightly in the coronal 171 Å and 193 Å channels over the course of the hourlong dataset analysed here, it is unclear whether this is an effect of cooling, due to line-of-sight effects, or purely noise.Little other evidence for cooling could be obtained through analysis of SDO/AIA time-series.Additionally, no evidence of coronal rain was observed in the 304 Å channel.More obvious cooling could be observed within the local AR, however, in a loop system located approximately 100 ′′ south of the loop analysed here.This loop contains large amounts of coronal rain in the 304 Å channel potentially hinting at the occurrence of the thermal instability, or catastrophic cooling.The effects of such catastrophic cooling on coronal loop oscillations would be an interesting topic for future study.
The key seismological result obtained here is that the location in time at which the amplitude begins to increase (an effect of the cooling on the oscillation) is dictated solely by the ratio between the internal and external densities for sufficiently large coronal loops.Therefore, fitting the model such that the difference between the observed turning point and the theoretical fitting point is minimised (as is shown in Fig. 7) allows us to calculate an estimate of the density ratio, ζ, an important parameter for modelling.For this loop, inversions of the theoretical model provide a good fit to the observed amplitude profiles when ζ is in the range 2.05-2.35(i.e., the loop foot-point is initially more than twice as dense as its surroundings).
In order to follow on from this work, we aim to complete two further studies.Firstly, we will conduct a statistical analysis of oscillations within potentially cooling coronal loops in the solar corona.This will provide important constraints on density ratios of loops for future modelling.Secondly, through further numerical work, it should prove possible to conduct further seismology in order to return values such as the scale height of the loop and the annulus thickness.These values will, again, provide further constraints for future modelling.Overall, the theoretical work of Shukhobodskiy et al. [2018] has proved adept at modelling the oscillations of the loop analysed here and should be considered by authors in the future when analysing coronal loop oscillations.

Figure 1 :
Figure 1: The evolution of the coronal loop system during the hour around the oscillations analysed here.Plotted are the SDO/AIA 304 Å (top row), 171 Å (middle row), and 193 Å (bottom row) channels.The white slit over-laid on the middle row indicates the location of the slit analysed in this article.

Figure 2 :
Figure 2: (Top row) The coronal loop system analysed here plotted using 171 Å data at five times during the oscillation.The white lines indicate the slit studied in this article, which corresponds to Loop 40, Event 2 from Goddard et al. [2016].(Bottom row) Time-space diagram constructed from the slit over-laid on the 171 Å data.The dashed vertical lines indicates the time steps plotted in the top row (from left to right).The bright chord starting and finishing at 12 ′′ is the oscillation studied in this article.

Figure 3 :
Figure 3: Zoomed time-space diagram from Fig. 2 with the Canny edge-detection output (blue lines) and the returned oscillation profile (red line) over-laid.The returned oscillation is qualitatively similar to the oscillation returned by Pascoe et al. [2016] for the same event.

Figure 4 :
Figure4: The position-time dependence of the coronal loop.Red dots correspond to data points collected from the 171 Å time-space diagram plotted in Fig.3.The green line corresponds to the fitted 8th-order background trend and the blue line is the Gaussian fit of the observed data points computed of the form Eq. 1.The two vertical dashed lines indicate the boundaries of the area over which numerical fitting took place.The Gaussian fit under-estimates the oscillation amplitude between frames 170 and 200 meaning the amplitude at this time will be under-estimated.

Figure 6 :
Figure6: The absolute value of the displacement of the Gaussian fitted position profile (blue line from Fig.4) from the background 8th order polynomial trend (green line from Fig.4) through time.The cut-out in the top right corner provides a magnified view of the region on the graph bounded by the black dashed box.This cut-out clearly shows the deviation of the damping profile from a typical Gaussian or exponential fit during this time due to a slight increase in the measured amplitude (compared to the background trend).Similar results are found for other background polynomial trends.

Figure 7 :
Figure 7: The dependence of the amplitude of the coronal loop oscillation on time.Green points correspond to maxima of the peaks plotted in Fig. 6, the blue line is the fitted observed amplitude-time profile, and the red line represents the theoretical amplitude-time profile with ζ = 2.2, λ = 1.05, κ = 2.4 and α = 1.4.