Shock wave energetics

The energetics of a global shock wave in the low solar corona


As the most energetic eruptions in the solar system, coronal mass ejections (CMEs) can produce shock waves at both their front and flanks as they erupt from the Sun into the heliosphere. However, the amount of energy produced in these eruptions, and the proportion of their energy required to produce the waves, is not well characterised. Here we use observations of a solar eruption from 2014 February 25 to estimate the energy budget of an erupting CME and the globally–propagating “EIT wave” produced by the rapid expansion of the CME flanks in the low solar corona. The “EIT wave” is shown using a combination of radio spectra and extreme ultraviolet images to be a shock front with a Mach number greater than one. Its initial energy is then calculated using the Sedov–Taylor blast–wave approximation, which provides an approximation for a shock front propagating through a region of variable density. This approach provides an initial energy estimate of 2.810 ergs to produce the “EIT wave”, which is approximately 10 % the kinetic energy of the associated CME (shown to be 2.510 ergs). These results indicate that the energy of the “EIT wave” may be significant and must be considered when estimating the total energy budget of solar eruptions.

Sun: coronal mass ejections (CMEs) — Sun: corona — Shock waves

1 Introduction

Typically observed as globally–propagating, diffuse features which can traverse the solar disk in under an hour, “EIT waves” were first observed in the solar corona using the Extreme ultraviolet Imaging Telescope (EIT; Delaboudinière et al., 1995) onboard the Solar and Heliospheric Observatory (SOHO; Domingo et al., 1995) spacecraft. These initial observations (Moses et al., 1997; Dere et al., 1997; Thompson et al., 1998) were quickly followed by many more, first from SOHO/EIT (e.g., Thompson & Myers, 2009), then subsequently from the Solar Terrestrial Relations Observatory (STEREO; e.g., Long et al., 2008; Veronig et al., 2008; Long et al., 2011a) and more recently the Solar Dynamics Observatory (SDO; e.g., Liu et al., 2010; Long et al., 2011b) spacecraft. Although originally considered to be the coronal counterpart of the chromospheric Moreton–Ramsey wave (Moreton, 1960; Moreton & Ramsey, 1960) as predicted by Uchida (1968), discrepancies in their kinematics and morphology raised difficulties with this assumption (e.g., Delannée, 2000).

Traditionally, “EIT waves” have been interpreted either as waves (fast–mode magnetoacoustic waves, shock waves or, alternatively, soliton waves e.g., Thompson et al., 2000; Wang, 2000; Wills-Davey et al., 2007; Grechnev et al., 2008) or as a brightening resulting from the restructuring of the coronal magnetic field (e.g., Delannée, 2000; Chen et al., 2002; Attrill et al., 2007; Delannée et al., 2008) during the eruption of a coronal mass ejection (CME). More recently, they have been interpreted using a hybrid approach, where the observed “EIT wave” is a freely–propagating magnetoacoustic wave initially driven by the rapid expansion of the CME as it erupts into the heliosphere (cf., Zhukov & Auchère, 2004; Downs et al., 2011). A more detailed discussion of the different theories proposed to explain “EIT waves” and the observations leading to these interpretations may be found in the recent reviews by Wills-Davey & Attrill (2009); Gallagher & Long (2011); Zhukov (2011) and Liu & Ofman (2014).

Figure 1: Running-difference images from the 193 Å passband showing the propagation of the global wave–pulse at 00:45:54 UT (panel a), 00:53:06 UT (panel b) and 01:01:54 UT (panel c). The EIS slit located at the southern polar coronal hole is indicated in each panel.

Although primarily observed as broad, diffuse features, there have been observations of sharply defined “EIT waves” which are usually identified as shock waves (e.g., Thompson et al., 2000; Grechnev et al., 2008, 2011; Veronig et al., 2010; Ma et al., 2011). These “S-waves” (e.g., Biesecker et al., 2002) have been noted to exhibit a drop in intensity, as well as a dramatic change in morphology with propagation (Warmuth et al., 2004a) and may be evidence for a shock wave decaying to an ordinary magnetohydrodynamic wave as it propagates away from its source (cf. Warmuth, 2007). This is consistent with recent observations by Carley et al. (2013), which showed evidence of acceleration of particles by a shock wave, driven in turn by an erupting CME. The shock wave was observed in the low corona as an “EIT wave” using EUV data from the Atmospheric Imaging Assembly (AIA; Lemen et al., 2012) onboard SDO. The paucity of shock wave observations compared to those of the more general “EIT wave” observations is consistent with the classification system proposed by Warmuth & Mann (2011), which groups “EIT waves” based on distinct kinematic classes. These include initially fast waves with pronounced deceleration, waves with moderate and nearly constant speeds and slow disturbances exhibiting erratic kinematics.

The very high temporal and spatial resolution of SDO/AIA has greatly increased the number of observations of “EIT waves” (e.g., Nitta et al., 2013), and has contributed to a deeper understanding of the phenomenon (as evidenced by the detailed recent review by Liu & Ofman, 2014). However a number of questions remain, particularly with regard to the energy budget of the wave and how this is related to the energy of the associated flare and CME. Previous work, particularly by Emslie et al. (2004, 2012), has attempted to examine the partition of energy in a solar eruption, finding that the CME dominates the energy budget with an associated kinetic energy of 10 ergs for each of the events studied. However, neither of these discussions included the energy of an “EIT wave”.

Although not included in analysis of the fractionation of energy budgets, some work has been done on identifying the energy associated with the “EIT wave” itself. A lower bound on the energy of an “EIT wave” was proposed by Ballai et al. (2005) using the oscillation profile of a coronal loop which had been impacted by a propagating “EIT wave”; this was estimated to be 3.4 J (3.4 ergs). A more recent analysis of the conductive, kinetic and thermal energy of the “EIT wave” was performed by Patsourakos & Vourlidas (2012), whose generic estimate of the energy of the “EIT wave” using typical observational parameters is much higher, at 1.8 ergs, comparable to the energy of a small flare. The relationship between the energy of the propagating wavefront and the energy of the CME itself remains uncertain, although Liu et al. (2012) did use observations of a filament–cavity oscillation resulting from the impact of an “EIT wave” to estimate the energy of the wavefront at 10–10 ergs, much smaller than the kinetic energy of the associated CME (10 ergs). This estimate of an “EIT wave” having 1 % the energy of the associated CME is much lower than the 10 % proposed by Vourlidas et al. (2010), suggesting that more work is required on this topic.

In this paper, we present observations of an “EIT wave” propagating through the solar corona made by SDO/AIA. Section 2 presents the observations of the “EIT wave”, which are shown in Section 3 to be consistent with that of a shock wave initially driven by the CME. This allows the energy budget of the shock to be calculated in Section 4 before the implications of these observations are discussed in Section 5.

2 Observations and Data Analysis

The “EIT wave” studied here erupted from NOAA active region AR 11990 close to the East limb of the Sun on 2014 February 25 (as shown in Figure 1) and was associated with a GOES X4.9 class flare which began at 00:39 UT before peaking at 00:49 UT. The eruption was also associated with a CME and a Type \@slowromancapii@ radio burst, indicating the presence of a shock. Although the “EIT wave” pulse may be seen in all of the passbands observed by SDO/AIA, the 193 Å passband was primarily used in this analysis as it provided the clearest observations of the wave–pulse. Despite the intense nature of the erupting flare and the broad extent of the associated CME, the wave–pulse did not propagate equally in all directions across the Sun, being constrained by the locations of the neighbouring active regions. However, it was clearly observed to propagate south from the erupting active region through the quiet Sun towards the southern polar coronal hole.

The AIA data were prepared using the normal routines available from SolarSoft (Freeland & Handy, 1998). The intense nature of the erupting X4.9 flare meant that the data were strongly affected by the Automatic Exposure Control (AEC) algorithm which dynamically changes the exposure time of the images depending on the intensity of the most recent image. This can produce a constant-exposure time-series with a number of low-exposure-time images interleaved. As faint features in particular are difficult to discern in images with a low signal-to-noise, images with a lower than normal exposure time were not considered in this analysis.

At the time of the eruption, the Extreme ultraviolet Imaging Spectrometer (EIS; Culhane et al., 2007) onboard the Hinode (Kosugi et al., 2007) spacecraft was observing the southern polar coronal hole using a study with short exposure times (30 s), which reduced the signal-to-noise of the observations. Despite this, it was possible to use the Fe xii 195.12 Å emission line (which is the strongest line observed by Hinode/EIS) to examine the observed “EIT wave”. This line is also observed by the 193 Å passband on SDO/AIA, allowing a comparison between the two instruments.

3 Characterisation of the wave–pulse

Figure 2: Top: Base-difference deprojected annulus image showing the eruption at 00:50:18 UT in the 193 Å passband. The vertical axis shows the height from Sun-centre in solar radii while the horizontal axis shows the angle clockwise from solar north. Bottom: Base-difference image showing the temporal variation at a height of 1.1 R (indicated by the dashed white line in the upper panel). The dashed white line shows the fit to the shock front at the leading edge of the bright feature, with the fitted initial velocity and acceleration given in the bottom right of the panel.

Considering the many interpretations for “EIT waves” and the different physics associated with each theory, it was important to characterise the true physical nature of the wave–pulse prior to any further analysis. The variations in Alfvén speed within the neighbouring active regions meant that the wave–pulse could not easily propagate through them and instead was forced to propagate primarily north and south along the solar limb as shown in Figure 1 and the associated online movie As a result, it was not possible to use the CorPITA technique (Long et al., 2014) to identify and characterise the wave–pulse. Instead, a deprojected annulus of 0.5 R width was used to study the evolution of the wave–pulse along the solar limb. This may be seen in the top panel of Figure 2 which shows the deprojected annulus image at 00:50:18 UT. Note that this is a base–difference image, with a pre–event image at 00:39:00 UT subtracted in order to highlight the wave–pulse. A similarly base–differenced stack plot is shown in the bottom panel of Figure 2, which gives the variation in the lateral extent of the wave–pulse with time at a height of 1.1 R (indicated by the dashed line in the top panel).

Figure 3: Radio spectra from the Learmonth Solar Radio Spectrograph (25–180 MHz) and the Siberian Solar Radio Telescope (SSRT; 180–400 MHz) showing the Type \@slowromancapii@ radio burst associated with this event. The Type \@slowromancapii@ burst is indicated by the dashed line.

The EUV observations of the wave–pulse as shown in Figure 1 show a clearly defined bubble with a sharp leading edge associated with the CME eruption. Although these observations are consistent with the interpretation of the wave–pulse as a shock wave, it was necessary to examine the kinematics of the wave–pulse to test this assumption. This was done by examining the variation in intensity over time, at a fixed height of 1.1 R, as shown in the stack plot in Figure 2. The leading edge of the bright feature was identified and fitted using a quadratic model (dashed line in the bottom panel of Figure 2), with this process repeated 10 times to minimise uncertainty. The dashed white line here indicates the fit to the leading edge of the shock front and has an estimated initial velocity of  km s and an estimated acceleration of  m s.

This approach was used to find the kinematics of the wave–pulse across a range of heights from 1 R to 1.3 R, giving an initial velocity range of  km s with a mean of  km s, and an acceleration range of  m s with a mean of  m s. These fitted initial velocities are larger than the 200–400 km s range estimated by Thompson & Myers (2009) using observations from SOHO/EIT and the mean velocity of 644 km s measured by Nitta et al. (2013) using higher–cadence observations made by SDO/AIA, indicating that this was a particularly fast “EIT wave” event.

Figure 4: The eruption from 2014 February 25 as observed by the SDO/AIA 193 Å passband. The eruption is shown close to the source using intensity and running-difference images in panels a & b respectively and over the South Pole using intensity and running-difference images in panels d & e respectively. Panels c & f show the intensity profiles from the white square regions shown in panels a & d. The Hinode/EIS field-of-view is shown in panels d & e, with the intensity profile from the EIS pixels within the white box shown in blue in panel f. The intensity images in panels a & d were processed for display here using the Multi-scale Gaussian Normalisation (MGN) technique of Morgan & Druckmüller (2014).

In addition to the EUV observations, a Type \@slowromancapii@ radio burst associated with this event was observed by both the Learmonth Solar Radio Spectrograph (Lobzin et al., 2010) and the Siberian Solar Radio Telescope (Smolkov et al., 1986) in the period 00:46–01:06 UT (see Figure 3). Although quite complex and fragmented, the burst shows clear characteristics of a Type \@slowromancapii@, with a slow frequency drift from 400 to 100 MHz. The height of the shock was estimated by assuming that we observe plasma emission and can therefore convert directly to electron density using . The radial stratification of density with height was then estimated using density measurements derived from all of the EUV filters of AIA as described by Zucca et al. (2014). This approach allows the height that the radio emission was produced at to be inferred from its frequency. In this case, the Type \@slowromancapii@ burst starting frequency of 350 MHz placed it at an initial height of 1.1 R(assuming the emission is the 2nd harmonic of the plasma frequency). Furthermore, a velocity of 550 km s could be derived from its frequency drift assuming radial propagation. This identification of a shock at low heights in the corona is consistent with the assertion of the “EIT wave” being a shock. However, while it is possible that the “EIT wave” and Type \@slowromancapii@ burst are from the same MHD disturbance in the corona, they may be at two spatially distinct locations and propagating in different directions. We are therefore cautious not to compare the Type \@slowromancapii@ and “EIT wave” velocities as they are not generally expected to be the same (Klassen et al., 2000), instead using the presence of a Type \@slowromancapii@ burst as a positive identification of a shock in the low corona.

The high measured velocities of the “EIT wave” and the occurrence of a Type \@slowromancapii@ burst are all consistent with the presence of a shock front, indicating that it may be possible to estimate the magnetosonic Mach number M of the disturbance. This can be achieved using the Rankine–Hugoniot relations for a perpendicular magnetosonic shock (), described as,


where is the density compression ratio , the plasma- and an adiabatic index has been assumed (Priest, 1982; Vršnak et al., 2002). The density compression ratio can be estimated from EUV images by examining the variation in image intensity as the shock front passes through the field of view as,


where and are the pre–event image intensity and density respectively while and are the image intensity and density measured during the passage of the wave–pulse. This relationship was used to calculate the variation in image intensity close to the source of the eruption and also at the southern polar coronal hole as shown in Figure 4. This approach has previously been employed by Muhr et al. (2011) and Zhukov (2011) and allows an estimate to be made of Mach number assuming no strong variation in temperature and a constant density along the depth of emission (assumed to be the pressure scale height). Since the total emission will be due to the emission measure along the line of sight, only a fraction of which is occupied by the plasma experiencing the shock, Equation 2 is therefore an underestimate of the true density enhancement.

A clear jump in intensity is apparent in panel c of Figure 4, which gives the intensity profile for the white box shown in panels a and b. This enhancement in image intensity is mirrored in the corresponding density compression ratio of 1.07 as given in the legend. The resulting magnetosonic Mach number of 1.10 confirms that the observed “EIT wave” was shocked at this stage in its propagation.

This process was repeated close to the southern polar hole as shown in panels d, e & f of Figure 4, with the intensity variation measured in a location at the bottom of the Hinode/EIS slit as shown in panel d. The temporal variation in the measured intensity was much lower here, likely reflecting the distance from the source. The density compression ratio of 1.04 and the resulting magnetosonic Mach number of 1.07 for the 193 Å passband are lower than close to the source but still exceed unity, indicating that while the wave–pulse remains shocked, the shock has decayed with propagation (cf. Kantrowitz & Petschek, 1966; Mann, 1995; Warmuth, 2007, and references therein).

The position of the EIS slit allowed the intensity ratio and hence density compression ratio and Mach number to be compared between the different instruments. The intensity of the EIS 195.12 Å emission line was averaged along the EIS pixels corresponding to the white box shown in Figure 4, producing the blue intensity profile shown in panel f of Figure 4 (note that the EIS intensity has been scaled to match the AIA intensity in the plot). Both intensity profiles exhibit similar variations, indicating that the intensity increase is a true feature. Although the EIS intensity exhibits more random variation with time compared to AIA, this is most likely due to the lower signal-to-noise in the narrow-band (monochromatic) EIS signal as opposed to the higher signal-to-noise in the broadband AIA signal. The density compression ratio and resulting Mach number were then calculated for the EIS intensity profile using Equations 1 and 2 and were found to be 1.05 and 1.08 respectively. These are comparable to those measured with AIA and provide independent verification of both the approach and the results.

Figure 5: The Alfvén speed at a height of 1.1 R as derived using the technique of Zucca et al. (2014). The fitted kinematics of the shock front are indicated by the dashed line.

These numbers may be further examined by comparing the velocity of the wave–pulse to the local Alfvén speed, which can be taken as a proxy measure of the fast-mode speed. As described by Zucca et al. (2014), the local Alfvén speed may be calculated using a combination of emission measure and polarisation brightness. The local Alfvén speed of the plasma that the wave–pulse propagated through (at a height of 1.1 R) is shown in Figure 5 to decrease from an initial peak of 900 km s near the source to under 300 km s in the quiet corona at a position angle of 200 clockwise from solar north. It then increases to 600 km s at the southern coronal hole. The speed of the wave–pulse (indicated by the dashed line in Figure 5) exceeds the local Alfvén speed at all locations, indicating that the wave–pulse was a shock. However, the resulting Mach number where the intensity ratio was measured near the source and close to the polar coronal hole was quite low, implying that it was very weakly shocked in these locations and that the assumption of negligible temperature variation is valid in this case.

4 Energetics of the Eruption

With the observed wave–pulse identified as a propagating shock wave, it was possible to estimate its energy using the Sedov–Taylor relation (Sedov, 1946; Taylor, 1950a). This technique was originally developed during the Second World War to estimate the energy of an atomic bomb, and assumes a spherical blast wave emanating from a point source. It should be noted at this point that this interpretation may not be strictly correct in this case as the shock wave may be initially driven by the lateral expansion of the associated CME rather than originating from a point source such as a flare (cf. Patsourakos et al., 2010). However, it serves as a good first–order approximation, allowing an initial estimate of the energy associated with the shock front. The Sedov–Taylor relation may be derived using dimensional analysis and, assuming that the shock front is passing through a region of constant density, can be written in terms of the radius at time as,


where is the energy of the blast wave and is the density. This may then be used to obtain an estimate of the energy E by plotting log  against log  and fitting the result with a line of slope (cf. Taylor, 1950b). This is shown in Figure 6, which shows log  plotted against log  at a heliocentric distance of 1.1 R. It is clear that the line of slope (the dashed–dotted line) is not a good fit to the data, indicating that the assumption of constant density is not valid in this case.

An alternative approach is to consider a variable-density medium where (with being the distance from the source). This approach was outlined numerically by Sedov (1959) and was previously employed in a solar context by Grechnev et al. (2008). In this case, the radius of the shock front at time is defined as,


where the shock decelerates if and accelerates if . This allows to be estimated by fitting a line of the form,


to the data shown in Figure 6. This provides a much better fit to the data (indicated by the dashed line in Figure 6) for a value , which is consistent with the deceleration given by the fitted kinematics in Figure 2. The assumption of a variable-density medium approximates the conditions through which the shock propagates as it passes from an active region through the quiet corona towards a coronal hole.

Figure 6: Plot of log  versus log  for the distance of the shock front from the source at a height of 1.1 R. The constant density fit is shown by the dot–dash line with the dashed line indicating the fit assuming a variable-density medium ( with ).

By assuming a typical coronal density of  cm, it is possible to make an estimate of the initial energy required to produce the shock front. The dashed line in Figure 6 shows Equation 5 fitted to the data with . This returns an energy estimate of  ergs for this eruption, which is much larger than the energy previously estimated using typically observed parameters by Patsourakos & Vourlidas (2012) ( ergs) and the minimum energy threshold estimated by Ballai et al. (2005) ( ergs).

Although the energy value estimated here is much larger than previous estimates, it has been suggested by Vourlidas et al. (2010) that the energy of an “EIT wave” is 10 % that of the associated CME. This suggestion may be compared with our results by estimating the energy of the associated CME using the equations


where and E are the kinetic and potential energy of the CME respectively, is the gravitational constant, is the solar mass, is the height in solar radii at which the CME front first appeared in the LASCO field of view, is the mass of the CME and is the instantaneous velocity of the CME (see e.g., Carley et al., 2012, for more details). These can then be combined to estimate the total mechanical energy of the CME.

The velocity of the CME was derived from a point-and-click trace of the CME front in running-difference images produced from the LASCO C2 and C3 coronagraphs during the period 01:20–03:20 UT. There was no acceleration in the C2 and C3 fields of view (as is expected of a fast CME), and a linear fit to the heliocentric distance vs. time data gave a velocity of 1550 km s. The mass of the CME was derived using the Thomson scattering approach first outlined by Billings (1966) and the methods detailed in a variety of studies such as Vourlidas et al. (2002, 2010) and Carley et al. (2012). In this case, base difference images from C2 and C3 were used and it was assumed that the CME propagated in the plane of the sky. This resulted in a mass of 2.210 g, which is on the larger end of the scale of CME masses (Vourlidas et al., 2010).

With the CME first appearing in the LASCO field-of-view at a height of 5 R, a constant velocity and constant mass were assumed, ignoring any mass accretion by the CME. This analysis produces an estimated potential energy of 5.010 ergs and kinetic energy of 2.510 ergs, giving a total mechanical energy estimate for the CME of 2.610 ergs. This indicates that the observed “EIT wave” had an energy 10 % that of the associated CME (which had quite a large energy in terms of the overall distribution of CME energies) consistent with the results of Vourlidas et al. (2010). The total mechanical energy of the eruption (excluding the energy released by the flare) can therefore be estimated at 2.810 ergs.

5 Discussion and Conclusions

In this paper, we use high–cadence SDO/AIA observations of an eruption on 2014 February 25 to estimate the energy of the associated “EIT wave”. The “EIT wave” was initially produced by the erupting CME before propagating freely, although neighbouring active regions meant that it propagated primarily north and south along the limb rather than across the disk visible from Earth. This allowed the lateral motion of the wave–pulse to be used to estimate its energy using the blast wave approximation of Sedov (1946) and Taylor (1950a, b).

The wave–pulse was identified as a shock wave using a combination of radio spectra, EUV intensity ratios, and the measured kinematics of the wave–pulse. A Type \@slowromancapii@ radio burst was observed which indicated the presence of a shock, with the high frequency of the burst suggesting a low coronal origin. The measured kinematics of the wave–pulse indicated that the observed “EIT wave” was particularly fast and exhibited deceleration as it propagated, suggesting that it may have propagated faster than the local fast–mode speed and was therefore a shock front. This was confirmed by examining the intensity ratio of 193 Å EUV images using the technique previously discussed by Muhr et al. (2011) and Zhukov (2011). This approach allows an estimate to be made of Mach number assuming no strong variation in temperature and a constant density along the depth of emission.

The primarily lateral motion of the wave–pulse along the limb away from the source meant that it was not possible to use the CorPITA technique of Long et al. (2014) to examine variations in the kinematics of the pulse. However, it was possible to use this lateral motion to estimate the initial energy budget of the wave–pulse using the blast wave approximation developed by Sedov (1946) and Taylor (1950a, b).

It was found that the equation of a blast wave in a region of constant density does not provide a good fit to the observed propagation of the wave–pulse. However, the equation of a blast wave in a variable-density medium provides an excellent fit to the data. This model is consistent with the observed propagation of the wave–pulse as it travels from the high density active region where it was produced through the lower density quiet corona towards the very low density coronal hole at the south pole. In addition, the degree of the variation in density was consistent with the observations of pulse deceleration, strongly indicating the validity of this approach. The resulting energy of the wave–pulse was found to be 2.810 ergs, approximately 10 % of the mechanical energy of the associated CME, which was estimated to be 2.610 ergs.

This consistency between the observations and the Sedov–Taylor relation is interesting and requires further investigation. The Sedov–Taylor relation is strictly valid for a spherical blast–wave originating from a point source, which is not the case here, and as a result can be considered as providing a first–order approximation. However, a growing number of both observations and theories suggest that “EIT waves” are consistent with the concept of a freely–propagating shocked simple wave formed by the rapid lateral expansion of a CME in the low corona (e.g., Vršnak & Cliver, 2008; Patsourakos et al., 2010). In this interpretation, the timescale over which the driver acts to produce the shock wave would have to be very short to allow it to be interpreted as a blast wave. The ability of the Sedov–Taylor relation to approximate the observations here suggests that this interpretation is valid and the timescale over which the CME acts to produce the “EIT wave” is very short.

This interpretation is supported by the excellent fit provided by the variable density medium formulation of the Sedov-Taylor relation. The accuracy of this model compared to the constant density formulation indicates that the wave was strongly affected by density variations in the medium that it propagated through. This would not be observed if the shock front were continuously driven by the erupting CME, suggesting that the wave front was freely–propagating.

These results suggest that the Sedov–Taylor blast wave approximation may be used to estimate the energy of an “EIT wave” shock front propagating through the solar atmosphere. In addition, the energy of the “EIT wave” shock is comparable to the energy of a flare and is not a negligible fraction of the total energy budget (cf. Emslie et al., 2004, 2012). These results indicate that the energy of the “EIT wave” is significant and must be considered when estimating the total energy budget of a solar eruption.

The authors wish to thank the anonymous referee whose comments helped to improve the paper. DML received funding from the European Commission’s Seventh Framework Programme under the grant agreement No. 284461 (eHEROES project). DB acknowledges support by STFC Consolidated Grant ST/H00260/1. PZ is supported by a TCD Innovation Bursary. EC is supported by ELEVATE: Irish Research Council International Career Development Fellowship – co-funded by Marie Cure Actions. SDO/AIA data is courtesy of NASA/SDO and the AIA science team. Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with NAOJ as a domestic partner, NASA and STFC (UK) as international partners. Scientific operation of the Hinode mission is conducted by the Hinode science team organized at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner countries. Support for the post-launch operation is provided by JAXA and NAOJ (Japan), STFC (U.K.), NASA, ESA, and NSC (Norway).


  1. Attrill, G. D. R., Harra, L. K., van Driel-Gesztelyi, L., & Démoulin, P. 2007, ApJ, 656, L101
  2. Ballai, I., Erdélyi, R., & Pintér, B. 2005, ApJ, 633, L145
  3. Biesecker, D. A., Myers, D. C., Thompson, B. J., Hammer, D. M., & Vourlidas, A. 2002, ApJ, 569, 1009
  4. Billings, D. E. 1966, New York: Academic Press, —c1966,
  5. Carley, E. P., McAteer, R. T. J., & Gallagher, P. T. 2012, ApJ, 752, 36
  6. Carley, E. P., Long, D. M., Byrne, J. P., et al. 2013, Nature Physics, 9, 811
  7. Chen, P. F., Wu, S. T., Shibata, K., & Fang, C. 2002, ApJ, 572, L99
  8. Culhane, J. L., Harra, L. K., James, A. M., et al. 2007, Sol. Phys., 243, 19
  9. Delaboudinière, J.-P., Artzner, G. E., Brunaud, J., et al. 1995, Sol. Phys., 162, 291
  10. Delannée, C. 2000, ApJ, 545, 512
  11. Delannée, C., Török, T., Aulanier, G., & Hochedez, J.-F. 2008, Sol. Phys., 247, 123
  12. Dere, K. P., Brueckner, G. E., Howard, R. A., et al. 1997, Sol. Phys., 175, 601
  13. Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1
  14. Downs, C., Roussev, I. I., van der Holst, B., et al. 2011, ApJ, 728, 2
  15. Emslie, A. G., Kucharek, H., Dennis, B. R., et al. 2004, Journal of Geophysical Research (Space Physics), 109, 10104
  16. Emslie, A. G., Dennis, B. R., Shih, A. Y., et al. 2012, ApJ, 759, 71
  17. Freeland, S. L., & Handy, B. N. 1998, Sol. Phys., 182, 497
  18. Gallagher, P. T., & Long, D. M. 2011, Space Sci. Rev., 158, 365
  19. Grechnev, V. V., Uralov, A. M., Slemzin, V. A., et al. 2008, Sol. Phys., 253, 263
  20. Grechnev, V. V., Afanasyev, A. N., Uralov, A. M., et al. 2011, Sol. Phys., 273, 461
  21. Kantrowitz, A. & Petschek, H. E.: in Plasma in Theory and Application, ed. W. B. Kunkel (McGraw Hill, Oxford 1966), 148
  22. Klassen, A., Aurass, H., Mann, G., & Thompson, B. J. 2000, A&AS, 141, 357
  23. Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3
  24. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17
  25. Liu, W., Nitta, N. V., Schrijver, C. J., Title, A. M., & Tarbell, T. D. 2010, ApJ, 723, L53
  26. Liu, W., Ofman, L., Nitta, N. V., et al. 2012, ApJ, 753, 52
  27. Liu, W., & Ofman, L. 2014, Sol. Phys., 67
  28. Lobzin, V. V., Cairns, I. H., Robinson, P. A., Steward, G., & Patterson, G. 2010, ApJ, 710, L58
  29. Long, D. M., Gallagher, P. T., McAteer, R. T. J., & Bloomfield, D. S. 2008, ApJ, 680, L81
  30. Long, D. M., Gallagher, P. T., McAteer, R. T. J., & Bloomfield, D. S. 2011, A&A, 531, A42
  31. Long, D. M., DeLuca, E. E., & Gallagher, P. T. 2011, ApJ, 741, L21
  32. Long, D. M., Bloomfield, D. S., Gallagher, P. T., & Pérez-Suárez, D. 2014, Sol. Phys., 66
  33. Ma, S., Raymond, J. C., Golub, L., et al. 2011, ApJ, 738, 160
  34. Mann, G. 1995, Journal of Plasma Physics, 53, 109
  35. Moreton, G. E. 1960, AJ, 65, 494
  36. Moreton, G. E., & Ramsey, H. E. 1960, PASP, 72, 357
  37. Morgan, H., & Druckmüller, M. 2014, Sol. Phys., 289, 2945
  38. Moses, D., Clette, F., Delaboudinière, J.-P., et al. 1997, Sol. Phys., 175, 571
  39. Muhr, N., Veronig, A. M., Kienreich, I. W., Temmer, M., & Vršnak, B. 2011, ApJ, 739, 89
  40. Nitta, N. V., Schrijver, C. J., Title, A. M., & Liu, W. 2013, ApJ, 776, 58
  41. Patsourakos, S., Vourlidas, A., & Stenborg, G. 2010, ApJ, 724, L188
  42. Patsourakos, S., & Vourlidas, A. 2012, Sol. Phys., 281, 187
  43. Priest, E. R. 1982, Dordrecht, Holland ; Boston : D. Reidel Pub. Co. ; Hingham,, 74P
  44. Sedov, L. I. 1946, Journal of Applied Mathematics and Mechanics, 10, 241
  45. Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics, New York: Academic Press, 1959,
  46. Smolkov, G. I., Pistolkors, A. A., Treskov, T. A., Krissinel, B. B., & Putilov, V. A. 1986, Ap&SS, 119, 1
  47. Taylor, G. 1950, Royal Society of London Proceedings Series A, 201, 159
  48. Taylor, G. 1950, Royal Society of London Proceedings Series A, 201, 175
  49. Thompson, B. J., Plunkett, S. P., Gurman, J. B., et al. 1998, Geophys. Res. Lett., 25, 2465
  50. Thompson, B. J., Reynolds, B., Aurass, H., et al. 2000, Sol. Phys., 193, 161
  51. Thompson, B. J., & Myers, D. C. 2009, ApJS, 183, 225
  52. Uchida, Y. 1968, Sol. Phys., 4, 30
  53. Veronig, A. M., Temmer, M., & Vršnak, B. 2008, ApJ, 681, L113
  54. Veronig, A. M., Muhr, N., Kienreich, I. W., Temmer, M., & Vršnak, B. 2010, ApJ, 716, L57
  55. Vourlidas, A., Buzasi, D., Howard, R. A., & Esfandiari, E. 2002, Solar Variability: From Core to Outer Frontiers, 506, 91
  56. Vourlidas, A., Howard, R. A., Esfandiari, E., et al. 2010, ApJ, 722, 1522
  57. Vršnak, B., Magdalenić, J., Aurass, H., & Mann, G. 2002, A&A, 396, 673
  58. Vršnak, B., & Cliver, E. W. 2008, Sol. Phys., 253, 215
  59. Wang, Y.-M. 2000, ApJ, 543, L89
  60. Warmuth, A., Vršnak, B., Magdalenić, J., Hanslmeier, A., & Otruba, W. 2004, A&A, 418, 1101
  61. Warmuth, A. 2007, Lecture Notes in Physics, Berlin Springer Verlag, 725, 107
  62. Warmuth, A., & Mann, G. 2011, A&A, 532, A151
  63. Wills-Davey, M. J., DeForest, C. E., & Stenflo, J. O. 2007, ApJ, 664, 556
  64. Wills-Davey, M. J., & Attrill, G. D. R. 2009, Space Sci. Rev., 149, 325
  65. Zhukov, A. N., & Auchère, F. 2004, A&A, 427, 705
  66. Zhukov, A. N. 2011, Journal of Atmospheric and Solar-Terrestrial Physics, 73, 1096
  67. Zucca, P., Carley, E. P., Bloomfield, D. S., & Gallagher, P. T. 2014, A&A, 564, A47
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description