Brought to you by:

Articles

OBSERVATIONS AND INTERPRETATION OF A LOW CORONAL SHOCK WAVE OBSERVED IN THE EUV BY THE SDO/AIA

, , , , , , , and

Published 2011 August 22 © 2011. The American Astronomical Society. All rights reserved.
, , Citation Suli Ma et al 2011 ApJ 738 160 DOI 10.1088/0004-637X/738/2/160

0004-637X/738/2/160

ABSTRACT

Taking advantage of both the high temporal and spatial resolutions of the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, we studied a limb coronal shock wave and its associated extreme ultraviolet (EUV) wave that occurred on 2010 June 13. Our main findings are: (1) the shock wave appeared clearly only in the channels centered at 193 Å and 211 Å as a dome-like enhancement propagating ahead of its associated semi-spherical coronal mass ejection (CME) bubble; (2) the density compression of the shock is 1.56 according to radio data and the temperature of the shock is around 2.8 MK; (3) the shock wave first appeared at 05:38 UT, 2 minutes after the associated flare has started and 1 minute after its associated CME bubble appeared; (4) the top of the dome-like shock wave set out from about 1.23 R and the thickness of the shocked layer is ∼2 × 104 km; (5) the speed of the shock wave is consistent with a slight decrease from about 600 km s−1 to 550 km s−1; and (6) the lateral expansion of the shock wave suggests a constant speed around 400 km s−1, which varies at different heights and directions. Our findings support the view that the coronal shock wave is driven by the CME bubble, and the on-limb EUV wave is consistent with a fast wave or at least includes the fast wave component.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

In 1997, the EUV Imaging Telescope (EIT; Delaboudinière et al. 1995) on board the Solar and Heliospheric Observatory (Domingo et al. 1995) observed globally propagating wave-like disturbances in the corona, known as EUV waves (also called "EIT" waves), for the first time (Moses et al. 1997; Thompson et al. 1998). Usually, EUV waves are observed as diffuse and irregular arcs of increasing coronal emission in the 195 Å channel. Compared to Moreton waves, EUV waves are a relatively frequent phenomenon and their speeds are relatively low. According to Biesecker et al. (2002), about 7% of the events in their large sample displayed sharp and bright wavefronts somewhat reminiscent of Moreton waves (Thompson et al. 2000). Such EUV waves are also called "bow waves" by Gopalswamy et al. (2000) or "S-waves" by Biesecker et al. (2002). In several events, the S-waves coincide spatially with Moreton waves observed at the same time (Khan et al. 2002; Warmuth et al. 2004), which would imply that at least S-waves are the long-sought coronal counterpart to Moreton waves. It is generally accepted that both Moreton waves and S-waves are shock waves.

However, the physical nature of the more common diffuse EUV waves is still being debated as to whether they are fast-mode MHD waves (e.g., Wills-Davey & Thompson 1999; Wang 2000; Wu et al. 2001; Ofman & Thompson 2002; Schmidt & Ofman 2010), slow-mode MHD waves (Wang et al. 2009), or non-waves related to a current shell or successive restructuring of field lines during the coronal mass ejection (CME; Delannée 2000; Delannée et al. 2008; Chen et al. 2002, 2005; Attrill et al. 2007). The observations from the Extreme Ultraviolet Imager (EUVI; Wuelser et al. 2004; Howard et al. 2008) instruments on board the twin Solar–Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008) spacecraft have advanced our understanding of the diffuse EUV wave (e.g., Long et al. 2008; Veronig et al. 2008; Gopalswamy et al. 2009; Attrill et al. 2009; Patsourakos et al. 2009; Ma et al. 2009; Kienreich et al. 2009; Zhukov et al. 2009; Veronig et al. 2010; Yang & Chen 2010). A hybrid wave/non-wave hypothesis was first suggested by Zhukov & Auchère (2004) and has been extended using computational simulation (Cohen et al. 2009; Downs et al. 2011). For detailed discussions of the different models and supporting observations, we refer to the recent reviews by Warmuth (2007, 2010); Vršnak & Cliver (2008); Wills-Davey & Attrill (2009), and Gallagher & Long (2010).

Due to the scarcity of direct observations of shock waves with high spatial resolution, the well-known signatures type II radio bursts (Wild & McCready 1950; Wild et al. 1954; Nelson & Melrose 1985) has long been used to indicate shock waves. More recently, coronagraphic observations in white light (Ontiveros & Vourlidas 2009) and UV spectra (Raymond et al. 2000; Mancuso & Bemporad 2009) have also been used to study coronal shock waves. Klassen et al. (2000) made a statistical study of type II bursts and found that 90% are associated with an EUV wave. However, the exact relationship between the EUV and shock wave is still unknown because of the low (or even no) spatial resolution of the radio data. For the same reason, the origin of the shock wave is also unknown. Two interpretations have been proposed: one suggesting a blast wave ignited by the pressure pulse of a flare and the other arguing for a piston-driven shock due to a CME (e.g., Reiner et al. 2001; Lin et al. 2006; Oh et al. 2007; Magdalenić et al. 2008; Vršnak & Cliver 2008).

Recently, Veronig et al. (2010) presented a study of a weak dome-like shock wave observed by EUVI which proved the feasibility of using EUV observations to study shock waves. The Atmospheric Imaging Assembly (AIA; Title & AIA team 2006; Lemen et al. 2011) on the Solar Dynamics Observatory (SDO) has seven EUV and two UV wavelengths, covering a wide range of temperatures, at high time cadence (12 s) and resolution (1farcs4, with 0farcs6 pixels). Taking advantage of these capabilities, Liu et al. (2010) studied a global coronal EUV wave and found multiple "ripples" for the first time; their results support the hybrid EUV wave hypothesis. Direct observations of shock waves with high spatial resolution in the low corona would greatly improve our understanding of the origin of shock waves and the nature of diffuse EUV waves (hereafter EUV wave).

In this paper, we present a detailed analysis of the 2010 June 13 event (cf. Kozarev et al. 2011), including time-dependent ionization and the nature of EUV wave near the solar surface. We briefly introduce the data used in this paper in the next section. The observations and results are presented in Section 3, and we present our discussion and conclusions in Section 4.

2. DATA AND OBSERVATION

On 2010 June 13, a limb solar eruption occurred in the active region (AR) NOAA 11079 (S24W91) from around 05:30 UT to 05:50 UT. This event involved a filament eruption, an M1.0 flare, and a CME with a strong and short-lived acceleration phase that has been studied in detail by Patsourakos et al. (2010). In particular, a distinct dome-like shock wave associated with this eruption was observed clearly by the AIA imagers in the coronal channels centered at 193 Å and 211 Å, allowing us to study the morphology and propagation characteristics of the shock wave in detail. In addition, a type II burst was also observed by some radio spectrographs, indicating the presence of a shock wave.

To study the shock wave and associated phenomena, we mainly used the following data.

  • 1.  
    The AIA on SDO provides multiple simultaneous high-resolution full-disk images of the corona and transition region up to 0.5 R above the solar limb with 1.5 arcsec spatial resolution and 12 s temporal resolution (Lemen et al. 2011). Seven narrow EUV bandpasses centered on specific lines: Fe xviii (94 Å), Fe viii,xxi (131 Å), Fe ix (171 Å), Fe xii,xxiv (193 Å), Fe xiv (211 Å), He ii (304 Å), and Fe xvi (335 Å) have been employed. The temperature diagnostics of the EUV emissions cover the range from 6 × 104 K to 2 × 107 K (Lemen et al. 2011). Here, we mainly use the channels centered at 171 Å, 193 Å, 211 Å, and 335 Å (level 1.5 images). To show the shock wave more clearly, we employed an improved radial filter technique first developed by S. Cranmer and A. Engell following Ma et al. (2010). The radial filter technique involves dividing the Sun into different concentric "rings" and calculating the minimum (Imin) and maximum (Imax) intensities in each radial ring. The scaled intensity at each pixel in the ring is then obtained by using the equation Iscaled = (IImin)/(ImaxImin).
  • 2.  
    The data from EUVI on board the STEREO Ahead (STA) spacecraft. STA is 73fdg6 ahead of AIA on its orbit during the eruption, therefore this event is a limb event for AIA and an on-disk event for STA. The two different points of view allow us to construct the three-dimensional structure of the coronal wave and pursue the nature of coronal waves. EUVI observed the chromosphere and low corona in four spectral channels (304 Å, 171 Å, 195 Å, and 284 Å) out to 1.7 R with a pixel-limited spatial resolution of 1farcs6 pixel−1 (Wuelser et al. 2004). Here we only use the 195 Å images from EUVI A with a time cadence of 2.5 minutes.
  • 3.  
    The data from the solar radio spectrograph of San Vito observatory which sweeps the frequency range 25–180 MHz every 3 s. It monitors solar radio emissions originating mainly in the solar corona. It has a low band (25–75 MHz) antenna (non-tracking semi-bicone) and a high band (75–180 MHz) antenna (tracking log-periodic). The radio data used in this paper was downloaded from the USAF Radio Solar Telescope Network (RSTN) solar radio fixed frequency and spectral data at the NOAA National Geophysical Data Center (NGDC) Web site (http://www.ngdc.noaa.gov).

3. RESULTS

3.1. Shock Wave in EUV Wavelengths

3.1.1. Physical Characteristics

Figure 1 shows the AIA observations (processed using the radial filter technique) of this event in the 171 Å, 193 Å, and 211 Å channels. The top row images in Figure 1 show the morphology of the corona before the shock wave appeared and the bottom row images show the morphology during the shock wave propagation. During the eruption, a sharp spherical CME bubble (indicated by the black arrows in Figure 1) can been seen in all three channels. For the first time, a spherical shock wave (indicated by white arrows in Figure 1) was observed by AIA at EUV wavelengths (see the online animations associated with Figure 1). Note that this event was also studied by Kozarev et al. (2011). The shock wave was ahead of the CME bubble from its first appearance in AIA images at about 1.23 R (from the Sun center) at ∼05:38 UT. The thickness of the layer of shocked gas is ∼2 × 104 km based on the measurement from AIA images along the radial direction (meridional 116 degree) in the 193 Å channel.

Figure 1.

Figure 1. AIA images showing the morphology of the shock wave (reversed color table). The AIA images were contrast enhanced using a radial filter technique described in Section 2.(Animations [A, B, C] of this figure are available in the online journal.)

Standard image High-resolution image

Interestingly, the shock wave appeared as a distinct bright feature only in the 193 Å and 211 Å channels. No clear signature of the shock wave can be identified from the 171 Å image in Figure 1, but the online animation Figure 1A shows that the shock wave appeared as a dark feature in the 171 Å channel. The 193 Å and 211 Å channels are dominated by the Fe xii lines (log T ∼ 6.2) and Fe xiv lines (log T ∼ 6.3) for AR observations, respectively. The 171 Å channel observes the Fe ix line (log T ∼ 5.85) for AR plasma (see O'Dwyer et al. 2010 for detailed information). We show below that the temperature of the dominant plasma in the shock wave is around 2.8 MK.

The white box in Figure 1 marks the area where we measured the intensity change before and during the passing of the shock wave. The intensity tracking in the white box is placed in Figure 2, with the intensities in 211 Å and 335 Å scaled by multiples 5 and 250, respectively. The plot shows that after the shock wave passed, the intensity in the 193 Å, 211 Å, and 335 Å channels increased rapidly and then went through a rapid decrease followed by a slow one. However, the intensity in the 171 Å channel decreased when the shock wave passed and increased slowly until the CME bubble arrived. The intensity jumps (maximum intensity divided by the pre-shock intensity) are 1.2, 1.5, and 1.7 in the 193 Å, 211 Å, and 335 Å bands, respectively.

Figure 2.

Figure 2. Intensity flux tracking in the white box in Figure 1. The green, pink, blue, and olive curves refer to the intensity flux in 171 Å, 193 Å, 211 Å (×5), and 335 Å (×250). The dashed black line and dash-dotted line refer to the time when the shock and CME bubble arrived in the white box.

Standard image High-resolution image

3.1.2. Upward Propagation

The three dotted lines in the bottom panels in Figure 1 indicate the radial directions at 115°, 116°, and 117° clockwise from the solar north pole, respectively, along which we analyzed the upward propagation of the shock wave and associated CME.

Figure 3 contains time–distance slit images showing the shock wave formation and propagation in different wavelengths along the three radial directions marked by the dashed black lines in Figure 1. The images in the left, middle, and right columns are taken from the 171 Å, 193 Å, and 211 Å channels, respectively. For each column, the top, middle, and bottom panels are obtained by placing the slit along the radial directions at 115°, 116°, and 117°. The dashed curves in each panel of Figure 3 are the GOES flux in 1–8 Å passband showing the associated flare. A filament associated with this eruption rose sharply beginning at 05:32 UT (panels (d) and (g) in Figure 3), indicating that the filament initiated this solar eruption.

Figure 3.

Figure 3. AIA slit images (reversed color table) in 171 Å (left column), 193 Å (middle column), and 211 Å (right column) showing the time–distances of the CME bubble and shock wave along 115° (top row), 116° (middle row), and 117° (bottom row), respectively. The contoured dash curves are the GOES 1–8 Å flux.

Standard image High-resolution image

The CME bubble is most apparent in the 171 Å channel, so we use it to study the propagation of the CME bubble. In each radial direction, the CME bubble shows similar propagation: going through a rapid acceleration phase and then moving with a nearly constant or slightly decreasing speed. The slit image at 117° in the 171 Å channel (panel (c) of Figure 3) is then used to measure the distance, speed, and acceleration of the CME bubble. The result is given in the top panel of Figure 4 with the red, blue, and orange lines with plus symbols indicating the distance, third degree polynomial fit speed, and acceleration of the CME, respectively. The upward speed of the CME initially increased from 0 to 500 km s−1 and then decreased to 250 km s−1 before leaving the AIA field of view. The corresponding acceleration of the CME decreased from around 6 to −4 km s−2. (For a detailed analysis of the CME, refer to Patsourakos et al. 2010.)

Figure 4.

Figure 4. CME dynamics and GOES flux. The red, blue, and orange curves with plus (diamond) symbols indicate the time–distance, time–speed, and time–acceleration profiles of CME bubble (shock), respectively.

Standard image High-resolution image

The shock wave can be identified in the 193 Å and 211 Å channels, although its appearance differs slightly in each. In the 193 Å channel, the shock wave front is relatively sharp and distinctly separated from the CME bubble (see panels (d)–(f) of Figure 3). However, in the 211 Å channel, the intensity in the region between the shock wave front and the CME bubble appears evenly increased (see panels (g)–(h) of Figure 3). As discussed in Section 3.3, this is because it takes a finite time to ionize the plasma from Fe xii, which dominates the 193 Å channel, to Fe xiv, which dominates the 211 Å channel. The shock wave first appeared in AIA images at 1.23 R from the surface of the Sun at around 05:38 UT (dot-dashed line in panel (e) of Figure 3), two minutes after the CME showed a distinct rise (dotted line). For measuring convenience, we chose panel (e) of Figure 3 to analyze the shock wave propagation along the radial direction. Rough estimates showed that the linear fit speed of the shock wave is around 600 km s−1, larger than the CME bubble's linear fit speed of 410 km s−1 (not including the rapid rising phase).

We identified the brightest point of the shock wave in each image and used a semi-automatic tracking method to obtain the distance of the shock wave along the radial direction. The result is shown in the top panel of Figure 4, with the red, blue, and orange lines with diamond symbols giving the distance, third degree polynomial fit speed, and acceleration of the shock wave front. The third polynomial fit speed of the shock wave showed a slight decrease from 600 km−1 to 550 km s−1 with a monotonically decreasing acceleration less than −1 km s−2.

The bottom panel of Figure 4 displays GOES soft X-ray flux. The orange curve, the purple curve, and the green curve are the time derivatives of the 1–8 Å, 0.5–4 Å, and 1–8 Å flux. The three dotted vertical lines in each of the panels from left to right in Figure 4 mark the times 05:36 UT, 05:37 UT, and 05:38 UT, respectively. The figure shows that the flare began at around 05:35:30 UT and peaked at 05:39 UT. The distinct disturbance of coronal loops which later became part of the CME bubble first appeared at 05:36 UT (left line), while the sharp CME bubble came into being at 05:37 UT (middle line). The shock wave appeared at 05:38 UT (right), which is later than the first appearance of both flare and CME.

3.1.3. Lateral Propagation and EUV Wave

Figure 5 gives base difference images showing the evolution of the shock wave and the associated EUV wave. The EUV wave here refers to the disturbance propagating laterally along the solar surface which normally appears as a projected bright circular feature propagating on the solar disk in EUV difference images. Panels (a)–(h) of Figure 5 are AIA base difference images in the 193 Å channel obtained by subtracting the 193 Å image at 05:32:08 UT from the present images. Panels (i)–(l) of Figure 5 are EUVI A images in the 195 Å channel obtained by subtracting the 195 Å image at 05:33:00 UT from the current images. The AIA data (panels (a)–(h), see also the online animation (A) associated with Figure 5) show a close relationship between the shock wave and the EUV wave: (1) they appeared at the same time, ∼05:38:08 UT and (2) the front of the laterally expanding shock wave is tightly connected to the position of the EUV wave front. In EUVI observations, we could not identify the coronal shock wave, probably because the shock wave is weak and the column depth is small. The EUV wave can be identified in EUVI as a circular disturbance propagating from AR 11079 (see also the online animation (B) associated with Figure 5).

Figure 5.

Figure 5. AIA 193 Å (top two rows) and EUVI A 195 Å (bottom row) observations show the propagation of the shock wave and EUV wave.(Animations [A and B] of this figure are available in the online journal.)

Standard image High-resolution image

The top panel of Figure 6 is a projected base difference image (see the online animation associated with Figure 6 for more information) where the horizontal coordinate is the polar angle (clockwise from north pole of the Sun) and the vertical coordinate is radial distance in solar radii. The five horizontal lines (which are parallel to the solar surface) from bottom to top indicate the layers located at 0.98, 1.01, 1.04, 1.07, and 1.11 R, respectively. By placing the slit along these layers, we obtained five different slit images from the base difference images in 193 Å. These slit images are placed in panels (b)–(f) of Figure 6. The horizontal coordinate of the slit images is the meridional angle in degrees while the vertical coordinate is time in minutes. The black area indicates the dimming area changing with time and the white features indicate the EUV brightening. Taking the propagating white feature as the coronal wave front, we estimate the EUV wave speed. The dotted lines and the numbers around them in Figure 6 indicate the linear fit distances of the EUV wave fronts and the estimated speeds. The propagation of the EUV wave toward the north (left) and south (right) is not symmetric; the speed in the southward direction appears slower than the northward direction.

Figure 6.

Figure 6. AIA 193 Å slit images showing the coronal wave propagating at different heights parallel to the solar surface. The base time of the event is 05:32:08 UT.(An animation of this figure is available in the online journal.)

Standard image High-resolution image

A bright feature was also observed propagating in the reverse direction, implying the reflection of the EUV wave (indicated by the "RF" in panels (d) and (e) of Figure 6). However, a detailed analysis of this observation is left for future work.

3.2. Shock Wave in Radio Observation

According to the NGDCs event listing,5 a type II burst associated with this event was observed by several radio spectrographs. In Table 1 we list the extracted dates, observatory stations, event start time, event end time, spectral class (SC), lower frequency (LF), upper frequency (UF), and estimated shock speeds in km s−1. The start time of the type II burst fits well with the AIA observation of the shock wave and there is no other solar eruption at the same time. Therefore, the shock wave observed at radio and EUV wavelengths should be the same one.

Table 1. Record of the Type II Burst that Occurred on 2010 June 13

Station Start Time End Time SC LF UF Shock Speed
        (MHz) (MHz) (km s−1)
BLEN 0537.0 0541.1 II 175X 332  
CULG 0537.0 0543.0 II 57X 150  
HIRA 0537.0 0550.0 II 50 310  
CULG 0538.0 0551.0 II 57X 200 700
SVTO 0538.0 0552.0 II 35 180 621
LEAR 0539.0 0548.0 II 55 180 665

Download table as:  ASCIITypeset image

Figure 7 shows the metric type II radio emissions from San Vito radio (SVTO) spectrograph. The top panel is the dynamic spectrogram. The type II burst occurring from 05:38 UT to 05:53 UT can be easily identified. The black curve indicates the harmonic frequency (2fp) in the type II burst along which we measured the local plasma frequency (fp) drift. The result showed that 2fp drifts from 165 MHz to 40 MHz. Considering the relationship between plasma frequency and density,

Equation (1)

we determine the density of the shock wave front (middle panel of Figure 7). The density of the shock wave dropped from 8.6 × 107 cm−3 at around 05:38 UT to 4.5 × 106 cm−3 at about 05:52 UT. The inset in the middle panel is the coronal plasma density model from Sittler & Guhathakurta (1999) as used by Lin et al. (2006),

Equation (2)

where y is the height above the solar surface in solar radii and n0 is the electron number density at the solar surface, which is chosen as 109 cm−3 here. The density model was used to estimate the position of shock wave front (the black line in the bottom of Figure 7), showing that the shock wave formed at around 1.25 R. The derived speeds are displayed by the dotted and dashed lines in the bottom panel of Figure 7. The dotted line shows the linear fit speed (527 km s−1) and the dashed line shows the second degree polynomial fit speed which decreased from around 600 km s−1 to 400 km s−1. The speed of the shock wave estimated using the model of Sittler & Guhathakurta (1999) is smaller than the speeds listed in Table 1. As the speed derived from radio emissions is strongly dependent on the coronal density model used, the error of the speed is quite large. Considering this situation, the speed of the shock wave measured from radio emission is consistent with that obtained from AIA observation.

Figure 7.

Figure 7. Dynamic spectrogram from San Vito observatory in the range of 25–180 MHz and the derived density and speed of the shock wave.

Standard image High-resolution image

The two black stars mark the frequencies used to measure the density jump at the shock using the equation

Equation (3)

As described in Vršnak et al. (2002), the band-split frequencies in type II emission map the electron densities behind and ahead of the shock front. In front of the shock, the plasma is characterized by the electron density n1 and emits radio waves at the frequency fL (LF branch), while the plasma behind the shock is characterized by the electron density n2 and emits radio waves at the frequency fU. At 05:40 UT, the lower and upper frequencies of the harmonic bands in type II burst are around 132 ± 5 MHz and 165 ± 15 MHz (marked by the black stars in the top panel of Figure 7), respectively. Using Equation (3), we obtained the density jump X = 1.56 ± 0.1 at 05:40 UT. Gopalswamy et al. (2011) obtained a similar result using radio data from HiRAS. For an oblique shock (taking the adiabatic index γ = 5/3), the Alfvén Mach number MA and the density jump X are related by

Equation (4)

(Vršnak et al. 2002). For β → 0, in the case of the perpendicular shock (θ = 90°), $M_{A}=\sqrt{X(X+5+5\beta)/2(4-X)}=1.45$ and in the case of longitudinal shock (θ = 0°), $M_{A}=\sqrt{X}=1.25$. Considering MA = v/vA and MA = 1.35  ±  0.1 (the actual error maybe larger) this produces an estimated Alfvén speed 450 ± 30 km s−1. The estimated sound speed $C_s=\sqrt{(\gamma \kappa T)/m}$ is 156 ± 30 km s−1 for a pre-shock temperature of (1.8 ± 0.4) × 106 K. According to the differential emission measure solutions for regions 2 and 3 in Kozarev et al. (2011), the pre-shock temperature (peak value) is around log T = 6.25 ± 0.1 (1.8 ± 0.4 MK). The area focused on in this paper differs slightly from the region 2 and 3 studied by Kozarev et al. (2011), so there may be a small discrepancy in the temperature, although it should not be too different. In the following context, we use 1.8 ± 0.4 MK as the pre-shock temperature T1. The resulting fast magnetosonic speed $c_f=\sqrt{\vphantom{A^A}\smash{\hbox{${v_A^2+c_s^2}$}}}$ is 476 ± 38 km s−1. In the solar atmosphere with $\tilde{\mu }=0.6$ and γ = 5/3, the magnetic field strength can be derived (from the Equation (2.48b) of Priest 1982) as B = 3.57 × 10−4n1/20VA. For a pre-shock density of 6 × 107 cm−3 from the type II frequency at around 05:40 UT, the value of vA implies B ∼ 1.3 Gauss. Gopalswamy et al. (2011) obtained an estimate for B that is consistent with this using a completely different technique based on the standoff distance between the flux rope and the shock.

3.3. Verifying the Shock Wave Interpretation

Coronal shocks have been identified from white light coronagraph images by morphology and density enhancement (e.g., Vourlidas et al. 2003; Ontiveros & Vourlidas 2009) and from UV spectra (Raymond et al. 2000; Mancuso & Bemporad 2009) by measuring oxygen kinetic temperatures. AIA shows brightening ahead of the bright CME loops, but that could in principle be part of an ejected streamer rather than a shock. The typical attributes of a shock are density compression, gas heating, and subsequent gas ionization. The intensities obtained from AIA imagers combine the electron density and ionization state of the plasma. To test the relation between compression and heating by an MHD shock, we employ Vs = 600 km s−1 from the observed speed of the shocked front and ρ21 = 1.56 from the band splitting in the radio. If the shock was purely gas dynamic, with a sound speed Cs = 156 km s−1, and a Mach number M = 600/156 = 3.8, the compression would be ρ21 = 3.3, instead of the observed 1.56. Therefore the magnetic field B is important. According to the jump condition for a perpendicular shock (Equation (2.19) of Draine & McKee 1993), the compression ratio is

Equation (5)

where D = (γ − 1) + (2M−2 + γM−2A). Assuming M = 3.8 and ρ21 = 1.56, we find MA = 1.55.

For a perpendicular shock, the shocked plasma parameters (v2, ρ2, p2, B2) are related to those of the unshocked plasma (v1, ρ1, p1, B1) by the equations for conservation of momentum and energy:

Equation (6)

(Priest 1982). Considering ρ21 = X, B2/B1 = X, v2/v1 = 1/X, p = nkT, c2s = γp11, v2A = B12/μρ1, M = v1/cs, and MA = v1/vA, the relationship between the pre-shock temperature and post-shock temperature can be obtained as

Equation (7)

Taking X = 1.56, M = 3.8, MA = 1.55, and γ = 5/3, T2/T1 is around 1.57 according to Equation (7). Assuming the pre-shock temperature of T1 = (1.8 ± 0.4) × 106 K as we explained in Section 3.2, the post-shock temperature T2 is around (2.8 ± 0.6) × 106 K.

Using the ionization rates qi at T = 2.8 × 106 K, we can estimate the ionization timescales and then compare them with the observed timescales of the band ratios. The ionization timescale to reach Fe xiv (211 Å) and Fe xvi (335 Å) can be obtained from the following equations:

Equation (8)

Equation (9)

where $q_{{\rm Fe\,\mathsc{xii}}}$, $q_{{\rm Fe\,\mathsc{xiii}}}$, $q_{{\rm Fe\,\mathsc{xiv}}}$, and $q_{{\rm Fe\,\mathsc{xv}}}$ are ionization rate coefficients in units of cm3 s−1. Given the temperature of plasma the ionization rate coefficients can be easily obtained from CHIANTI by using the code ioniz_rate.pro. In Table 2, we list the ionization rate coefficients at 2.4 MK, 2.6 MK, 2.8 MK, 3.0 MK, 3.2 MK, respectively, in units of cm3 s−1. Taking ne = 6 × 107 × 1.56 = 9.4 × 107 cm−3, the ionization timescales can be obtained by using the Equations (8) and (9). The corresponding results are placed in the right two columns in Table 2.

Table 2. Ionization Rate Coefficients (in units of cm3 s−1) at Different Given Temperatures and the Corresponding Estimated Ionization Timescales

T $q_{{\rm Fe\,\mathsc{xii}}}$ $q_{{\rm Fe\,\mathsc{xiii}}}$ $q_{{\rm Fe\,\mathsc{xiv}}}$ $q_{{\rm Fe\,\mathsc{xv}}}$ t211 t335
(MK)         (s) (s)
2.4 4.0982035e-10 2.4158748e-10 1.4844586e-10 1.0106319e-10 110 386
2.6 4.6829123e-10 2.8404318e-10 1.7771562e-10 1.2294273e-10  94 323
2.8 5.2553271e-10 3.2686128e-10 2.0770653e-10 1.4567346e-10  83 277
3.0 5.8123741e-10 3.6961965e-10 2.3806008e-10 1.6894433e-10  74 242
3.2 6.3520065e-10 4.1199241e-10 2.6848624e-10 1.9249588e-10  66 215

Download table as:  ASCIITypeset image

Figure 8 displays time intensity ratio profiles, which are derived from the intensity changes in different wavelengths (Figure 2) for the region marked by the white box in Figure 1. The rise phase of the intensity ratio profile indicates the ionization time. The ionization timescales of 211 Å and 335 Å obtained from observation are around 100 ± 12 s (marked by the top two vertical lines in Figure 8) and 275 ± 12 s (marked by the bottom two vertical lines in Figure 8), respectively. In general, the ionization times obtained from observation are consistent with that derived from the theory of shock waves, 85 ± 13 s for t211 and 288 ± 53 s for t335. In other words, the brightening ahead of the CME bubble is consistent with the shock jump conditions, while it would be purely coincidental if the structure were not a shock. The theoretical estimates of the ionization times are limited by the uncertainty in the pre-shock temperature, and this might be improved by further analysis of the AIA data. We also note that the intensity ratios of the raw band intensities can be misleading, because they convolve the ionization times with the timescale for the increase of emission from the shocked gas compared to the intensity from the foreground and background plasma outside the shock. We have estimated the foreground and background emission based on the variation of the 171 Å band, and in this case the timescales for the intensity ratios are similar to those obtained from the raw band intensities. Moreover, a more detailed analysis would include the pre-CME density profile and the adiabatic cooling as the shocked gas expands. It should be possible to obtain accurate values of the electron temperature at different positions around the shock to study electron–ion temperature equilibration and energy losses to particle acceleration.

Figure 8.

Figure 8. Time intensity ratio profiles from the white box in Figure 1. The thick black curve is for the intensity ratio of 211/193 and the thin black curve is for the intensity ratio of 335/193 (×30).

Standard image High-resolution image

4. DISCUSSION AND CONCLUSIONS

On 2010 June 13, a CME-associated shock wave was directly observed by AIA. Combining these observations, EUV data from STEREO A and radio data from San Vito observatory, we present a detailed study of the shock wave and its low coronal reactions.

4.1. Main Findings

Our main findings are as follows.

  • 1.  
    The shock wave front first appeared at 05:38 UT as a dome-like bright feature, which is clearly distinguishable from its associated CME bubble. The top of the dome-like shock wave originated from about 1.23 R with the thickness of the shocked layer around 2 × 104 km. A type II radio burst occurred at the same time, further confirming the presence of the shock wave. The electron density jump X (the ratio of the plasma density behind the shock and that in front of the shock) is around 1.56, indicating that the shock wave is a weak shock wave. The pre-shock density and the density jump from the radio data, along with a shock speed of 600 km s−1, imply a post-shock temperature of 2.8 MK. The ionization timescales of 211 Å and 335 Å obtained from the AIA observation are consistent with the ionization timescale derived from shock wave jump relation in general, which proves that the dome-like bright feature is a shock wave.
  • 2.  
    The shock wave appeared distinctly only in the channels centered at 193 Å and 211 Å with the estimated relative intensity of the shock wave front increasing by 1.2, 1.5, and 1.7 in 193 Å, 211 Å, and 335 Å, respectively. All of the EUVI channels showed an increasing intensity with the exception of the 171 Å channel which showed a clear intensity decrease.
  • 3.  
    The upward speed of the shock wave shows a slight decrease from about 600 km s−1 to 550 km s−1 with a deceleration less than 1 km s−2. The lateral speed of the shock wave shows no acceleration, but varies according to both height and direction. The linear fit lateral speeds toward the north range from 246 km s−1 to 397 km s−1 while the speeds toward the south range from 342 km s−1 to 486 km s−1. The lateral speed of the shock wave seems to increase with the height below 1.1 R. In general, the lateral speed of the shock wave is less than its upward speed.

4.2. Origin of the Shock Wave

The high time cadence of the AIA observations makes the pursuit of the origin of the low coronal shock waves possible. AIA data show that a filament rising was observed at 05:32 UT, at the very beginning of this solar eruption. About four minutes after the filament rising, a flare and CME bubble appeared nearly simultaneously. At 05:38 UT, two minutes after the flare and CME bubble appeared, the shocked front came into being. Considering the time sequence, it is likely that the filament is a trigger for the whole eruption. Because the flare and CME bubble appeared simultaneously and earlier than the shock wave, neither of them can be eliminated as the driver of the shock wave according to the time sequence alone. However, considering the strong similarity of the shock wave and the CME in both morphology (dome-like) and kinematics (slightly decreasing), we would suggest that the CME bubble played the role of a piston driving the shock wave.

4.3. Nature of the EUV Wave

The direct observation of the low coronal shock wave by AIA strongly supports the existence of a dome-like fast-mode wave observed in EUV wavelengths. The fastest front of the EUV wave coincides with the lateral expansion of the shock wave. Therefore, our result supports the hypothesis that the EUV wave is a fast wave or at least includes a fast wave component. Considering the projection effect of the on-disk observation as discussed in Ma et al. (2009), the on-disk observation may include a non-wave component. However, due to the weakness of the on-disk wave front observed by EUVI, we did not get strong evidence for non-wave component as proposed by Liu et al. (2010). More recently, Chen & Wu (2011) showed an on-disk event and found both fast and slow propagation fronts. However, there are still inconsistencies between the observed phenomena and the different theories proposed to explain the EUV wave. The high cadence observations from SDO should allow these issues to be resolved.

We sincerely thank the anonymous referee for very helpful and constructive comments that improved this paper. We acknowledge the AIA team for the easy access to calibrated data. We are grateful to San Vito Observatory, NOAA's NGDC, and RSNT for providing the radio data and corresponding analysis software. The STEREO/SECCHI data are produced by an international consortium: NRL, LMSAL, NASA, GSFC (USA); RAL (UK); MPS (Germany); CSL (Belgium); and IOTA, IAS (France). S.M. thanks Kelly Korreck, Mark Webb, Ed Deluca, and Meredith Wills for helpful discussions. S.M. is also grateful to Steven Cranmer, Alexander Engell, and Henry Winter for help with techniques. At SAO, this work was supported by subcontract SP02H1701R from Lockheed Martin, by NASA contract NNM07AB07C, and by NASA grant NNX08BA97G. J.C.R. acknowledges support from NASA grant NNX11AB61G. J.L.'s work was supported by Program 973 grant 2011CB811403 and NSFC grant 10873030.

Footnotes

Please wait… references are loading.
10.1088/0004-637X/738/2/160