Free Access
Issue
A&A
Volume 654, October 2021
Article Number A112
Number of page(s) 14
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202141672
Published online 19 October 2021

© ESO 2021

1. Introduction

Type-II bursts are narrowband lanes of transient solar radio emission that appear to drift slowly toward lower frequencies over time in dynamic spectra (for detailed discussions about their properties, see the reviews by Nelson & Melrose 1985; Vršnak & Cliver 2008; Nindos et al. 2008; Pick & Vilmer 2008). The instantaneous bandwidth of type-II bursts may be as narrow as a few MHz (e.g., Nelson & Melrose 1985). They typically start below 150 MHz (e.g., Mann et al. 1996) although cases with starting frequencies equal to or higher than 500 MHz have been reported (e.g., Pohjolainen et al. 2008; Magdalenić et al. 2012; Cho et al. 2013). In the meter wavelength range type-II bursts last from less than 5 min to about 30 min and their drift rate lies between 0.1 and 0.4 MHz s−1 and increases with increasing starting frequency (e.g., Mann et al. 1995, 1996). These values of drift rates are consistent with magnetohydrodynamic shocks propagating upwards in the corona and driving radio emission near the plasma frequency and/or its harmonic via the plasma emission mechanism.

Only about 60% of coronal type-II bursts display fundamental-harmonic emission bands (Nelson & Melrose 1985), and the fundamental band is usually weaker than the harmonic band. It appears likely that the fundamental band is absorbed, especially when the radio emission occurs behind the limb and passes through dense regions of the corona. Occasionally, the fundamental and/or harmonic bands are divided into substructures, which are called multiple lanes when they show irregularly varying frequency ratios and split bands when they show an essentially constant frequency ratio, Δf/f ≈ 0.1 − 0.2. Multiple-lane events are interpreted in terms of emission from distinct source regions on the shock (e.g., Cairns & Knock 2006) that can have different plasma densities, propagation speeds, and geometry allowing the emissions to appear separate and drift at different rates (e.g., Robinson & Sheridan 1982; Vršnak et al. 2006; Zimovets & Sadykov 2015). A similar interpretation has been proposed for the split-band effect (McLean 1967; Schmidt & Cairns 2012). Alternatively, sometimes the band splitting is interpreted in terms of plasma emission from both upstream and downstream of the shock (Smerd et al. 1974) with the frequency difference allowing the Alfvén Mach number of the shock to be determined (e.g., Vršnak et al. 2001, 2002, 2004; Kouloumvakos et al. 2014). Although there is no accepted model for plasma emission from the downstream region (e.g., Cairns 2011) observational support for this interpretation has been provided by Zimovets et al. (2012), Zucca et al. (2014a), and Chrysaphi et al. (2018).

Although it is generally accepted that interplanetary (IP) type-II bursts are produced by coronal mass ejection (CME) driven shocks (e.g., Gopalswamy 2006), there is no consensus about the origin of coronal type-II bursts. Coronal shocks could be generated by different drivers: CME-related erupting structures, the pressure pulse of a flare, or flare-related small-scale ejecta. Both flare-related drivers and CMEs could act as pistons, with their main difference being the time over which energy is supplied by the piston (Vršnak & Cliver 2008). In the former case the piston is temporary and the shock freely propagates as a large-amplitude wave (blast shock), whereas in the latter case the wave is supplied by the energy provided by the CME.

Statistical studies (see Nindos et al. 2011, and references therein) show that most, if not all, coronal type-II bursts are observed during events with both flares and CMEs. Observations of metric type-II bursts without CMEs are rare and are usually attributed to source regions located within about 30° from the central meridian (e.g., Claßen & Aurass 2002), where the detection of CMEs is more difficult. On the other hand, the lack of association of a metric type-II burst with a flare is usually interpreted in terms of the occurrence of the flare behind the limb. In some events, there is a tight synchronization between the flare impulsive phase and the CME acceleration phase (e.g., Zhang et al. 2001, 2004; Temmer et al. 2008), which makes the identification of the shock driver difficult.

Support of the CME-driven scenario has been provided by the high correlation between metric type-II bursts and extreme ultraviolet (EUV) waves (Klassen et al. 2000; Biesecker et al. 2002), although each one may appear separately e.g., (Nitta et al. 2014). EUV waves (see the reviews by Patsourakos & Vourlidas 2012; Warmuth 2015; Long et al. 2017) are thought to be driven by the lateral expansion of CMEs. Furthermore, Gopalswamy et al. (2005, 2009) proposed that all metric type-II bursts are produced by CME-driven shock waves pointing out that the Alfvén speed profile may exhibit a local maximum in the corona, which could explain why metric type-II bursts usually cannot be followed through longer wavelengths even if they are produced by the same CME. However, Cane & Erickson (2005) have argued against a common CME driver for both coronal and IP shocks.

Type-II bursts limited to relatively high frequencies might be attributed to significant decay of the shock strength with height. Therefore they are usually associated with shocks driven by flare blast waves (e.g., Vršnak et al. 1995; Nindos et al. 2011; Magdalenić et al. 2012), although exceptions have also been reported (e.g., see Nindos et al. 2011; Kumari et al. 2017 for high-frequency type-II bursts associated with CMEs).

There are relatively few multiwavelength case studies of the origin of metric type-II bursts because (1) the metric emission comes from heights beneath the occulting disk of coronagraphs and (2) there are few interferometers that can provide routine observations of the Sun at metric wavelengths. Old Clark Lake and Culgoora imaging observations (see Gergely et al. 1983; Gary et al. 1984, and also the review by Nelson & Melrose 1985) show that in most cases the radio emission is located close to the CME front or its flanks.

More recent radio imaging observations have strengthened and refined the above conclusion because they have been accompanied by dynamic spectra with high temporal and spectral resolution, high-quality coronagraphic images, and imaging data in the EUV and soft X-rays (SXR). Probably the largest recent compilation of imaging type-II observations has been reported by Ramesh et al. (2012) who found that 38 of the 41 type-II sources that they observed at 109 MHz were associated with CMEs and were located near their leading edge. A similar conclusion for five events observed in the 164−435 MHz range was reached by Maia et al. (2000). The coronal shock radio sources may occur in front of an erupting flux rope (Bain et al. 2012; Zimovets et al. 2012; Zimovets & Sadykov 2015), above expanding soft X-ray loops (e.g., Klein et al. 1999; Dauphin et al. 2006), in association with an erupting jet during the progression of a CME (e.g., Zucca et al. 2014a; Maguire et al. 2021), in front of an EUV bubble in both radial and lateral directions (Kouloumvakos et al. 2014) or ahead of a CME that was deflected in the low corona (Pick et al. 2016). Other studies indicate that type-II sources are associated with the flanks of CMEs (e.g., Cho et al. 2007; Démoulin et al. 2012; Zucca et al. 2014b, 2018; Morosan et al. 2019) CME-streamer interactions in the low corona could also be important for the production of type-II bursts (e.g., Kong et al. 2012; Feng et al. 2012; Eselevich et al. 2015; Mancuso et al. 2019; Chrysaphi et al. 2020).

In a smaller number of events the type-II sources appear to be flare-related (e.g., White et al. 2007; Magdalenić et al. 2008, 2010, 2012; Nindos et al. 2011; Kumar et al. 2016). In these cases the type-II events were tightly synchronized with the related flares and either no CME was observed or the CME was not synchronized with the type-II burst. The above discussion shows that type-II bursts may occur under significantly different conditions.

In this article we study a complex event, SOL2010-02-012T11:25:00, whose dynamic spectrum contained at least five lanes of type-II emission and for which radio imaging data at ten frequencies were available. It is thus one of the most complex type-II events with adequate radio imaging data ever presented in the literature. We combine the radio spectral and imaging data with coronagraphic images and EUV and soft X-ray images from a variety of sources including Solar Terrestrial Relations Observatory (STEREO) Ahead (A) and Behind (B), in an attempt to associate the type-II-related shocks with disturbances in the corona. We present our data in Sect. 2 and our results in Sect. 3. In Sect. 4 we discuss the origin of the type-II bursts and in Sect. 5 we summarize our conclusions.

2. Observations

Observations with the Appareil de Routine pour le Traitement et l’Enregistrement Magnétique de l’Information Spectral-Jean Louis Steinberg (ARTEMIS-JLS) radiospectrograph (Kontogeorgos et al. 2006) and their analysis have been described in detail in previous articles (e.g., Nindos et al. 2011; Bouratzis et al. 2015; Alissandrakis et al. 2019). Here we only mention that the event was observed in high temporal resolution (10 ms) with the acousto-optic analyzer (Spectrograph Acousto-Optic, SAO) of the instrument in the frequency range of 265−470 MHz with a spectral sampling of 1.4 MHz. For an overview of the event we used data from the sweep-frequency Global Spectral Analyser (Analyseur de Spectre Global, ASG), which operated in the 100−700 MHz range with a time resolution of 100 ms and a spectral sampling of 1 MHz.

The Nançay Radioheliograph (NRH; Kerdraon & Delouis 1997) is a synthesis instrument that provides two-dimensional images of the Sun with sub-second time resolution. For the event that we study here, the NRH provided data at ten frequencies (150.9, 173.2, 228.0, 270.6, 298.7, 327.0, 360.8, 408.0, 432.0 and 445.5 MHz) with a cadence of 250 ms. All ten frequencies were within the spectral range of the ASG, while the last eight were also within the range of the SAO.

The observed source position may oscillate due to ionospheric effects, which are stronger in the winter and for low solar elevation. Indeed, such oscillations affected the position of a noise storm continuum visible in the full-day summary NRH plots1. Our measurements showed that the rms amplitude of the ionospheric oscillations at 150.9 MHz was 74″ at the start of the daily NRH observations but dropped to 30″ before the start of our event, which corresponds to about 7% of the average source size (about 430″) at the same frequency. The shift was even smaller at high frequencies (6″, or about 3% of the source size at 408.0 MHz), and thus ionospheric effects could be safely ignored.

From the original NRH visibilities, we computed two-dimensional images with a resolution of 1.2′ by 1.9′ at 432 MHz. We also computed one-dimensional images using visibilities from the east-west (EW) and north-south (NS) arrays only; this improved the resolution by a factor of two, due to the fact that the NRH extension antennas have a very small contribution to the two-dimensional images (cf. Bouratzis et al. 2016).

In addition to the ARTEMIS-JLS and NRH data we used light curves from the Geostationary Operational Environmental Satellite (GOES) and the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), Hα images from Catania Observatory, soft X-ray images from the GOES Soft X-ray Imager (SXI) and the Hinode Soft X-ray Telescope (XRT), EUV images (195 and 304 Å) from STEREO Sun-Earth Connection Coronal and Heliospheric Investigation (SECCHI) Ahead and Behind (Howard et al. 2008), coronagraph images from STEREO Coronagraph1 (COR1) A and B (Howard et al. 2008), as well as white light images and magnetograms from the Michelson Doppler Imager, (MDI) before and after the event. No Solar and Heliospheric Observatory (SOHO) Extreme ultraviolet Imaging Telescope (EIT) data were available during the event, while the Transition Region and Coronal Explorer (TRACE) was pointed elsewhere.

We only had four Hα images and those after the start of the event, with one of them during the radio burst. The SXI provided 15 images during the event (9 during the burst), with various exposure times and filters; the pixel size was 5″. From the XRT we obtained a set of low resolution images with a pixel size of 8.2″ and a Ti poly filter covering the entire event with a cadence of 2 min and two sets of high resolution images (pixel size of 1″) with Al thick and Be thick filters with a 20 s cadence, starting in the early phase of the event. Among the SECCHI images, those from STEREO-A in the 195 Å band had the highest cadence, 150 s; the cadence from STEREO-B was 5 min. In the 304 Å band the cadence was 10 min, giving only one image during the metric burst. In the 171 and 284 Å SECCHI bands the event was not recorded due to the very low image cadence (2 h). We note that several EUV and soft X-ray images suffered from saturation effects.

COR1 gave images with a cadence of 5 min, with the first CME image beyond the occulting disk, located at 1.6 R, having been after the radio burst. We did not use any Coronagraph2 (COR2) images because the CME passed beyond the occulting disk well after the burst.

3. Results

3.1. Overview of the event

The event studied in this work was the third and strongest of a set of 3 homologous events, described by Chernov et al. (2014); it occurred around 11:25 UT on February 12, 2010, in active region 11046. It was a fairly strong (GOES class M9) and complex event, with an associated CME, the main phases of which are given in Table 1. Preliminary results were reported by Alissandrakis et al. (2011). The other two events on that day were weaker than the event studied here; they occurred around 07:21 UT (C7.9) and 09:40 UT (B9.6). Their dynamic spectra were rich in fine structure (see Chernov et al. 2014) but did not show any type-II emission. The 07:21 UT event was accompanied by a CME. We add that a noise storm continuum associated with the active region was observed by the NRH from February 10 to February 17, 2010.

Table 1.

Time line of the event of February 12, 2010.

Figure 1 shows the Dynamic Spectrum (DS) of the event, recorded with the ASG (sweep frequency) receiver of the ARTEMIS-JLS radio spectrograph, together with the GOES and RHESSI light curves. The DS shows that the event was rich in fast drift, type-III, U and other structures. Both the light curves and the spectrum indicate two phases, the first significantly stronger than the second. The second phase started at ∼11:28:15 UT with type-III bursts in the low-frequency part of the dynamic spectrum. A secondary maximum appeared around 11:30 in the GOES and RHESSI light curves. The bulk of the radio emission ended around 11:32 UT, but metric continuum emission appeared from 12:00 to 12:45 UT (see Fig. 45 of Bouratzis et al. 2015); this emission will not be considered in this work because it had no associated type-II bursts.

thumbnail Fig. 1.

Dynamic spectrum (negative) of the February 12, 2010 event (500 to 20 MHz), observed with the ASG receiver of the ARTEMIS-JLS radiospectrograph, together with light curves in soft X-rays from GOES and hard X-rays from RHESSI in arbitrary units and on a linear scale. Dashed horizontal lines mark the frequencies of the Nançay Radioheliograph. 1F and 1H mark the fundamental and harmonic emission of the first type-II burst, A, B1 and B2 mark the branches of the second type-II.

Two type-II bursts appeared, one in each phase of the event, in the metric part of the spectrum (> 100 MHz). The first type-II showed fundamental-harmonic (FH) structure, marked 1F and 1H in Fig. 1, as well as band splitting (see Sect. 3.4.1 for more details); no herringbone structures were detected. The second type-II differed from the first in that it had three branches with different drift rates, marked A, B1 and B2 in the figure, with B1 and B2 starting 20 s after A, with superposed pulsations (see also Sect. 3.4.2); there was no FH structure, no band splitting and no herringbones. At even lower frequencies, WIND/WAVES and STEREO/WAVES data showed no extension of the type-IIs, but only two strong type-III bursts, which were extensions of the type-III groups that appeared in the high frequency of the spectrum at the beginning of each phase.

Seen from the Earth, the event occurred near the central meridian (N25, E11; heliocentric distance 33°). Its evolution, through images from SXI, XRT, SECCHI (STEREO A and B), Hα, as well as white light images and magnetograms is presented in Fig. 2. From the position of STEREO A the flare was very close to the eastern limb, while from STEREO B it was seen near the western limb.

thumbnail Fig. 2.

Evolution of the event in Hα (Catania Observatory), soft X-rays (GOES/SXI and Hinode/XRT), 195 and 304 Å (STEREO A and B), together with a magnetogram and white light image from MDI. The white arc in the STEREO images marks the solar limb. The center of the field of view is at Carrington longitude 186.6, latitude 25.0 (solar x = −107″, y = 510″ for Earth images).

The left column of Fig. 2 shows the pre-flare situation, where a set of loops are visible in the SXR image (panel 29); the magnetogram is shown twice, saturated at ±200 G in panel 15. The first brightening appeared near 11:20:46 UT in the 195 Å STEREO-A image; it is well visible in the XRT image at 11:21:31 (panel 30 in Fig. 2), and was located near a compact negative polarity region adjacent to a positive polarity spot, marked as NP in panel 15 of the figure. We note that this happened about 3.5 min before the first type-III occurred.

The bright point evolved to a compact source; in the XRT high resolution image (panel 31 in the figure) it appears elongated, with approximate full width at half maximum of about 32″ × 15″, probably a compact loop. A second, smaller compact source appeared at the beginning of the second phase of the event in high resolution SXR images, located ∼8″ NE of the first (panel 32); by that time the peak of the previous bright point had shifted slightly to the north. In panels 33−35 a set of loops NW of the compact source is visible in the SXI images, with bright footpoints visible in Hα; these appeared during the second phase of the event, indicating a restructuring of the magnetic field. We note that, apart from these bright points, no extended flare ribbons were observed either in SXR or Hα.

Around 11:23 UT, ejecta in the form of a surge appeared in the STEREO images. The surge had both a hot component visible in emission at 195 Å and a cool component visible in emission at 304 Å and in absorption at 195 Å and Hα (Fig. 2); however, we could not identify any trace of the surge in the SXR disk images. From the Hα images (panels 25−27 in the figure) it is obvious that the surge did not originate in the compact source, but ∼30″ to its east: we note that surge activity in Hα was observed prior to and after the flare (panels 22 and 28, respectively), at different locations. In the SECCHI images the surge was visible well after the end of the metric burst, up to about 13:26 UT.

There was no flux rope proxy (e.g., hot channel, sigmoid, etc; see for example Table 1 in Patsourakos et al. 2020) in the employed limb and disk imaging EUV and SXR data. The surge itself, as observed in the Extreme Ultraviolet Imager (EUVI) 195 Å of STEREO-A contained seemingly intertwined bright and dark threads, in emission and absorption respectively, which might be an indication of the flux rope structure; in EUVI 195 Å of STEREO-B the surge appears almost exclusively in absorption. On the other hand, inspection of the 304 Å movies in both STEREO A and B suggests that the ejected plasma was rather laminar.

Figure 3 shows running difference images at 195 Å from STEREO A (top row) and from STEREO B (bottom row). The eruption produced a disturbance expanding in the corona, in the form of an EUV wave (e.g., Thompson et al. 1998). Essentially, the limit of the EUV wave separates the coronal region already affected by the eruption from the region that has not been affected yet. We note that, as the disturbance expands, it is better visible in the south of the event than in the north in the STEREO-A images. An interesting remark is that the top of the surge appears to be close or even beyond the limit of the EUV wave (panels d and e of the figure). The trace of the base of the EUV wave is clearly seen on the disk in panels d–g (STEREO-A) and i–k (STEREO-B) in Fig. 3. The EUV wave expanded outside the flaring active region and propagated over significant distances including surrounding quiet Sun areas. In addition, during its propagation, it caused deflections of ambient structures, best observed along its northern track in STEREO-A. We note in particular a disrupted set of loops N of the surge, better visible in the material presented in Sect. 3.3 below; these might be associated with the restructuring of loops seen in the SXR images during of the second phase of the event, as noted above.

thumbnail Fig. 3.

Running difference images at 195 Å from STEREO A and B. We note the formation of an expanding EUV wave. The white arc marks the limb.

The white light images from the STEREO-A coronagraph COR1 (Fig. 4, top row) show a streamer north of the active region that produced the flare; the CME was first seen above the occulting disk (1.6 R) at 11:35:21 UT. The difference images (middle row of the figure), computed by subtraction of the average of images at 11:25:21 and 11:30:21 UT, clearly show a deflection of this streamer to the north by about 150″ (marked SD on the figure), apparently a result of interaction with the expanding coronal disturbance.

thumbnail Fig. 4.

COR1 images and their base difference (top and middle row) and 304 Å images, saturated to show the surge better. We note the appearance of the CME above the occulting disk at 1.6 solar radii at 11:35:21 UT, the disruption of the streamer (marked SD) and the surge (E) seen both in the coronagraph and the 304 Å images. The dark circle marks the limb. The field of view is 2500 by 2500″.

The first image where the streamer deflection was recorded above the coronagraph occulting disk was after the end of the metric event. With a 5 min cadence and the occulting disk at 1.6 R it is hard to specify when this deflection started and at which height, it is thus quite possible that the deflection started earlier at lower heights. Indeed, the 195 Å STEREO-A image at 11:28:12 UT shows a disruption of loops at the base of the streamer, well visible in the material presented in Sect. 3.3 below at position angle of 40°. It appears that the streamer was an obstacle to the propagation of the disturbance to the north, and this is the reason that its trace was better visible south of the event in the 195 Å images, as noted in the previous paragraph. A final remark is that, in the last column of Fig. 4, the surge is visible both in the coronagraph and the 304 Å images.

3.2. Triangulation and velocity estimates

The cadence of STEREO images is not sufficient to give accurate velocities from triangulation. However, we could measure the apparent motion of the top of the surge in 195 Å STEREO-A images; we found that the radial component of the velocity was about 500 km s−1 and that the acceleration peaked around 11:25 UT (i.e. near the images shown in the first column of Fig. 3). This roughly coincides with the start of the bulk of the type-III emission, but it is 2 min before the start of the first type-II.

In addition, we were able to apply triangulation on several points along the surge to the 195 Å images at 11:30:46 UT; this gives a rough idea of the shape and, possibly, the trajectory of the surge. The maximum radial distance measured in this way was 1.27 R, that is, 265″ above the photosphere. By fitting a straight line to the positions of these points, we found that the projection of the surge on the photosphere was oriented 39° east of north, whereas its angle with respect to the vertical was 36°.

It is also possible to estimate velocities by measuring positions in the original images, in particular those of STEREO-A where the flare appeared very close to the limb. To this end we computed cuts of intensity along the radial direction, as a function of position and time, for a number of position angles. An example is shown in Fig. 5, where a radial cut of COR1-A intensity is shown, together with associated images.

thumbnail Fig. 5.

Coronagraph observations. Top: COR1A intensity as a function of position angle and distance from the center of the disk, r. Middle: same as top, with the pre-CME intensity subtracted. Bottom: cut of intensity as a function of time and r at a position angle of 61°. Horizontal lines are at the same r in all panels. Vertical lines in the top and middle panel mark the position angle of the cut. Bottom panel: vertical line marks the time of the other panels. The photospheric limb is at r = 972″.

Three features, moving at different radial velocities, are visible in the cut: the fastest is a blob in the CME, also visible in panel d of Fig. 4, moving at 680 ± 80 km s−1, not too different from the speed of the 195 Å surge mentioned in the previous paragraph. A slower feature, at 244 ± 16 km s−1, corresponds to the 304 Å surge seen in the images of the right column in Fig. 4; the lower speed of this feature, compared to the speed of the 195 Å surge, might be due to the fact that in the latter case the top of the surge was measured whereas here the average position of the blob, which apparently expanded radially, was recorded. The slowest feature, at 120 ± 8.5 km s−1, is associated with a second ejection seen in 304 Å images later during the event, visible at 60° position angle in the lower part of the middle panel of the figure. Velocity errors quoted here and later in the text correspond to the rms of values obtained at several position angles.

A linear extrapolation of radial positions placed the start of the CME around 11:26:30 UT. There was a data gap in SOHO/LASCO C2, whereas C3 recorded a halo CME centered around a position angle of 44° and moving at 509 km s−1; the estimated onset time was 11:18:22 UT. Taking into account the geometry of the event, the present estimates from COR1 should be more accurate than those from C3.

Applying the same technique to the CME associated with the earlier event of 07:21 UT (first CME), we detected a single feature moving at ∼740 km s−1 and no evidence of surge.

In a similar way we measured the lateral expansion of the southern front of the CME, using image cuts parallel to the limb. In the height range 1080−1530″ above the limb we found velocities of ∼320 km s−1, less than half the value of the radial expansion. The corresponding value for the first CME was 410 ± 20 km s−1.

The southern part of the expanding EUV wave in the 195 Å band shown in Fig. 3 and discussed in the previous section is well visible in image cuts parallel to the limb; the northern edge of the EUV wave is also detectable, though much weaker. At the radial distance of the flare, we measured lateral expansion velocities of about 400 km s−1 and 560 km s−1 for the southern and northern fronts of the disturbance, respectively.

The southern front of the expanding EUV wave can be followed above the limb, with its lateral velocity decreasing with r, down to ∼350 km s−1 at r = 1200″, approaching the lateral velocity of the CME. This is apparently a geometric effect, for a roughly spherical expanding disturbance.

Finally, the radial velocity of the EUV wave was measured from cuts in the radial direction, giving a value of about 460 km s−1. The velocities measured from white light (WL) and EUV images are summarized in Table 2.

Table 2.

Measurements of velocities from WL and EUV images.

3.3. Position of the type-II sources and other features with respect to the surge

In this section we will compare the information on the surge obtained by triangulation with the positions of the type-II bursts and other features seen from the Earth. Figure 6 shows the projection of the surge position on top of Hα and SXR images. We note that the extrapolation of the surge positions to lower heights passes very close to the Hα and SXR flares; we also note that the Hα surge is located at a different position and has a slightly different orientation, showing that it is not part of the same phenomenon.

thumbnail Fig. 6.

Projection of the surge on the sky plane as seen from the Earth (symbols), on top of an Hα image (left) and an XRT image (right). Different symbols refer to different parts of the surge, giving an idea of its lateral extent.

Figure 7 shows the positions of the NRH sources during the first and second type-II bursts. As the DS is rich in features, during the first phase of the event in particular, we took care to avoid emissions that were not related to the type-IIs. Before proceeding we should note that the observed positions may deviate from the true positions due to refraction and scattering effects; for a spherically symmetric corona, refraction moves the source closer to the disk center, whereas scattering (Steinberg et al. 1971; Bougeret & Steinberg 1977; Bastian 1994) moves it in the opposite direction and, at the same time, broadens the source. These effects are stronger for emission at the fundamental, where the refraction index deviates significantly from unity.

thumbnail Fig. 7.

Position of the type-II sources, together with the projected position of the surge (open circles) and the position of the Hα/SXR flare (large open circle). Left: first type-II, individual images; black is for 360.0 MHz, red for 327.0 MHz, green for 297.8 MHz and blue for 150.9 MHz. Right: average images for branches A (black), B1 (red) and B2 (blue) of the second type-II burst.

In recent works (e.g., Kontar et al. 2019; Zhang et al. 2021 and references therein) both refraction and anisotropic scattering are considered. Zhang et al. (2021) employed a spherically symmetric coronal model at 35 MHz and, for the heliocentric position of our flare, gave radial shifts of 0.08 to 0.24 R for the fundamental and −0.05 to 0.03 R for the harmonic, negative values denoting shifts toward the center of the disk (their Fig. 7). According to Kontar et al. (2019), the width and the decay time of the source scale approximately as f−1; if this scaling also holds for the shifts, they should be much smaller than the above values for our frequency range, a factor of four to 13 higher than 35 MHz, and probably negligible for the harmonic. Indeed, scattering computations (Kontar 2021, priv. comm.) showed that the source displacement for the harmonic at 300 MHz is 7″ and 1.1′ for the fundamental at 150 MHz, both considerably smaller than the NRH beam size (1.5′ and 4.4′, respectively); hence no corrections were applied.

The most compact group of points in the left panel of Fig. 7 is that from 360 MHz; this is located about 80″ away from the flare and further away from the surge. Although the radio sources are big, from about 150″ by 90″ at 408.0 MHz to 400″ by 260″ at 150.9 MHz, their positions can be measured with much higher accuracy; hence the difference with respect to the surge position is real. Positions at 150.9 MHz show the largest scatter. In general the radio emission does not appear to be associated with the surge, although the radio sources are aligned roughly in the same direction. For the second type-II, the average position of the radio sources are at a different location, close to the flare and the lower part of the surge, but they are aligned along a direction quite different from that of the surge.

It is customary in solar metric radio astronomy to associate the emission frequency with the source height and thus obtain three-dimensional positions. There are some assumptions inherent to this, first of all that the emission is at the plasma frequency or its harmonic, which is true in the case of type-II emissions; moreover, an atmospheric model is required and, more serious, it is assumed that the corona is not disturbed, which is probably hardly the case in our situation. EUV waves are associated with small density and temperature enhancements (e.g., see the reviews of Patsourakos & Vourlidas 2009; Warmuth 2015; Long et al. 2017, and references therein). Our EUV data do not contain high-cadence images at multiple channels to perform differential emission measure analysis and unlock the thermal distribution of the plasma perturbed by the EUV wave. An idea about the degree of disturbance of the corona can be obtained from the properties of the EUV wave and the band splitting. As in almost all cases, our EUV wave is associated with weak (less than 20−30%) intensity enhancement and hence, possibly, even smaller density enhancements. On the other hand the compression ratio deduced from band splitting (see Sect. 3.4.4 below) is 1.5.

Having said the above, we obtained three-dimensional positions by assuming an isothermal corona at 2 × 106 K and a base electron density twice that of the Newkirk model. This gave radial distances from 1.21 R (408.0 MHz, harmonic) to 1.51 R (228.0 MHz, harmonic) and 1.35 R (150.9 MHz, fundamental), corresponding heights of 204″, 495″ and 340″ respectively.

Our next step was to compute the heliographic coordinates of the sources and from those their projected positions on the nearest 195 Å STEREO-A images. The results are shown in Fig. 8. Despite the scatter of the source positions during the first type-II (left panel of Fig. 7), here their projections are close to one another, near the edge of the 195 Å EUV wave, but far from the surge. For the second type-II the positions of all three branches line up, as if B1 and B2 which started 20 s after A, followed the same path. These are far from the surge, too; we note, however, that the extrapolation of branch A positions passes close to the disrupted loop system mentioned in Sect. 3.1, located north of the surge at a position angle of 56°.

thumbnail Fig. 8.

Projected positions of the type-II radio sources, on top of polar projections of running difference images from 195 Å STEREO-A. Top: type-II during the first phase of the event (11:27:00−11:28:05 UT); asterisks show the average positions from 360.8 MHz, “x” from 298.7 MHz, crosses from 298.7 MHz and open circles from 150.9 MHz. Bottom: type-II during the second phase of the event (11:29:21−11:30:57 UT); asterisks are for branch A, crosses for branch B1 and open circles for branch B2.

3.4. Detailed study of the type-II bursts

In addition to the overall evolution of the event described in the previous sections, some particular features are worth mentioning and analyzing; these will be discussed in the present section, together with an in-depth study of the burst properties.

3.4.1. An impressive fundamental-harmonic structure

The first interesting feature is a fundamental-harmonic emission (F/H) structure observed during the first phase of the event, already mentioned in Sect. 3.1, and shown in more detail in Fig. 9 where we have combined ASG and SAO data. We note that, in addition to the type-II emission, other emissions showed F/H structure as well.

thumbnail Fig. 9.

Fundamental-harmonic emission during the first phase of the event (composite ASG/SAO dynamic spectrum). The frequency range is from 106 to 450 MHz. Dotted horizontal lines mark the NRH frequencies. Channels around 200 and 310 MHz have been partly corrected for radio interference. The arrow points to a spectral feature discussed in the text.

The type-II itself showed band split, to be discussed further in Sect. 3.4.4. It had prominent fine structures at the fundamental, reminiscent of spike-like emissions reported by Armatas et al. (2019). The harmonic emission appeared more diffuse, the high-frequency split-band in particular, which showed some pulsating emission. Emission at the fundamental was not visible after 11:28:05 UT, while the harmonic emission extended up to 11:28:40 at least.

We employed two-dimensional cross-correlation to determine the time and frequency shift between the F/H emissions. For the full event we found a frequency ratio of 1.95 and a 0.6 s delay of the fundamental with respect to the harmonic. However, these values reflect mostly emissions outside the Type-II, which are the brightest. Isolating the type-II in the dynamic spectrum as much as possible, we obtained a similar frequency ratio and a smaller time delay: 0.23 s for the entire type-II and 0.5 s for the prominent feature at 11:28 UT (arrow in Fig. 9).

The measured F/H frequency ratio is consistent with previous works (Wild et al. 1954; Maxwell et al. 1963; Roberts 1959); it is expected to be slightly less than 2. Since radiation at the plasma frequency will not propagate, what is observed at the fundamental should originate in the short wavelength wing of the emission profile. We note, however, that in the presence of magnetic field the cutoff frequency of the extraordinary goes above the plasma frequency while that of the ordinary so that (cf. Eq. (2.122) in Zheleznyakov 1996):

(1)

where υ = (fpe/f)2 and (u = fce/f)2 are parameters expressing the electron density and the magnetic field through the plasma frequency and the electron gyrofrequency respectively. For example, with a magnetic field of 18 G (see Sect. 3.4.4 below) and an observing frequency of 150 MHz, we have and the cutoff frequency for ordinary radiation is about 1.15fp; this will facilitate the escape of ordinary mode radiation at the fundamental.

The time delay is attributed to the difference in the group velocity of the electromagnetic waves near the plasma frequency and its harmonic and to the different paths of the radiation due to refraction and scattering. Roberts (1959) was the first to measure time delays on the order of 1 s, from data with a time resolution of 0.5 s, whereas Maxwell et al. (1963) reported delays in the range of 2−150 s. Here, with an effective time resolution of 0.1 s and the use of cross-correlation, the delay can be more accurately measured.

The feature observed in the DS around 11:28 UT, marked by the arrow in Fig. 9, appears to be part of the type-II. However, one-dimensional NRH images (Fig. 10) show clearly that it was not, as it was located at a different position, west of the type-II source. It lasted for about 8 s and showed strong intensity fluctuations. Its polarization was of the same sense as that of the type-II, but weaker (on the order of 20%). As a whole it showed practically no frequency drift and had a stronger band split than the main event, visible only in the harmonic. Its role in the event is not clear, yet we considered that it was worth mentioning; its nature would not have been revealed without the imaging information.

thumbnail Fig. 10.

Expanded view of the spectral feature marked by the arrow in Fig. 9 (top panel, same for both columns). The lower rows show one-dimentional NRH images at three frequencies in the EW (left) and NS (right) directions.

3.4.2. Pulsations during the second type-II burst

Figure 11 shows the DS of the second type-II burst in more detail. As mentioned in Sect. 3.1, it consisted of three branches, a wide branch (A in the figures) and two narrow ones (B1 and B2 in the figures). Apart from that, the DS shows superimposed pulsations, an expanded view of which, with five times better time resolution, is shown in the top panel of Fig. 12. We note that the pulsations are seen in both the DS and the one-dimensional NRH images; most interesting, they are not limited inside the branches but extend over the continuum emission in between, as clearly shown in the bottom panel of the figure. Thus the bandwidth of the pulsations was about 200 MHz.

thumbnail Fig. 11.

Composite ASG/SAO dynamic spectrum of the second type-II burst.

We note that the period of pulsations decreased with time; a power spectrum analysis gave an initial period of 1.13 s, dropping progressively to 0.42 s. These values are within the range reported in the literature (e.g., Nindos & Aurass 2007; Bouratzis et al. 2015). It is also noteworthy that in the initial pulsations there is embedded temporal fine structure (right bottom panel of Fig. 12) and that the duration of their high intensity phase is longer than that of the low intensity phase.

thumbnail Fig. 12.

Expanded views of pulsations. Top: pulsations near the start of the second phase of the event. Bottom: DS of the first part of the pulsations recorded with the ASG (left) and the SAO (right), together with 1-dimensional EW images from the NRH placed at the corresponding DS frequencies (middle). The frequency range is 220−430 MHz; time tickmarks are every 2 s in the left and middle panels, every 1 s in the right panel.

A detailed examination of the one-dimensional NRH images revealed that, not only the pulsating source existed before the start of the second type-II, but there was also a second source of opposite circular polarization (Fig. 13), which started during the first phase of the event, around 11:28:46 UT. The two sources were visible in all NRH frequencies between 228.0 MHz and 408.0 MHz, but hardly visible in the DS due to their low intensity and phase differences in the intensity fluctuations which reduced the overall signal. The source with negative polarization disappeared near the start of the type-II, and it is tempting to speculated that the shock might be the result of the interaction of the two pulsating regions.

thumbnail Fig. 13.

Pulsations before the start of the second type-II; Dynamic spectrum (top) and 1D NRH images at 270.6 MHz in Stokes I and V, EW and NS.

Pulsations are typical of type-IV bursts. In our case, there is certainly a type-IV continuum background; although this is not clearly visible in the DS due to the rich superposed fine structure, it is visible between the branches of the second burst as mentioned above. It appears that whatever caused the background pulsations induced pulsations in the emitting shock fronts as well; an alternative interpretation is that the continuum emission region is located in front of the shock front and is optically thin, so that its pulsations modulated the intensity of the type-II emission.

We used cross-correlation of the average intensity of selected DS spectra channels with all other channels in an effort to measure frequency drifts. We found a reverse relative drift of about 2 s−1; such a high drift cannot be interpreted in terms of motions, so that its origin needs further investigation.

3.4.3. Drift rates and velocity estimates

The drift rate in the early part of the first type-II was −7.7 × 10−3 s−1; around 11:27:20 UT the drift decreased and the first type-II appeared stationary until 11:28:05 UT (see Fig. 9), at which time a small jump to lower frequencies was observed, followed by another jump around 11:28:20 UT. Thus the value of the average drift over the ∼110 s duration of the burst dropped to −3.3 × 10−3 s−1 (Table 3).

Table 3.

Type-II drift rates and nominal velocities.

Assuming plasma emission, the frequency drift is related to the velocity of the exciter, υ and the scale length of the electron density, LNe along the path of the exciter, :

(2)

(3)

Classically, the drift is used to measure υ, assuming a constant value for LNe. However, in the present case of variable drift and the complex coronal environment, the observed decrease in the drift could be attributed to a decrease in the shock speed and/or to an increase in the density scale length. Moreover, the value of LNe is usually assumed to be equal to the scale height predicted by a particular spherically symmetric coronal model, so that the computed value of υ represents the radial component of the exciter velocity. With coronal models being practically isothermal at a temperature T, the scale height at a radial distance r, H(r), is given by:

(4)

where H is the scale height at r = R, μ is the average molecular weight and g the gravity at r = R. Thus the decrease in the drift rate might indicate that the shock encountered a high-temperature region with increased density scale length.

The above discussion shows that there are many uncertainties in computing the shock speed from the frequency drift. Nevertheless, we added the value of this nominal velocity, computed for a typical value of H = 105 km, in Table 3.

The bandwidth of branch A in the second burst increased with time, so that its low-frequency edge had a higher drift rate than its high-frequency edge (Table 3). Branches B1 and B2 were short (∼30 s) and narrowband, 7−14 MHz wide (relative bandwidth of 2.7−5.3%), with drifts between those of the two edges of branch A. No significant time variations in the drift rate were detected in this burst. We note that the drift rates of the components of the second type-II were up to a factor of 2 higher than the drift at the start of the first type-II and that all nominal velocities are considerably higher than the velocities derived from WL and EUV images (Table 2).

In the case of branch A of the second burst we had position measurements at a sufficient number of frequencies to obtain a three-dimensional trajectory by associating the emission frequency to the source height (Sect. 3.3 and Fig. 7). Using this information we computed a radial velocity of 2450 ± 260 km s−1. This value is very close to the average velocity of the low and high-frequency edges of the branch, derived from the frequency drift. From the same trajectory we computed a horizontal velocity of 2050 ± 230 km s−1, directed southward, and a total velocity of 3200 ± 230 km s−1.

3.4.4. Magnetic field from band splitting

Band splitting is widely accepted to be due to the density jump at the shock front, i.e. between the uncompressed plasma in front of the shock and the compressed plasma behind it. From the frequency ratio of the two bands the Alfvén Mach number, MA can be computed and from that, together with the density scale and the frequency drift, the intensity of the magnetic field can be estimated (Smerd et al. 1974; Vršnak et al. 2002, see also the review by Alissandrakis & Gary 2020).

Combining (3) with the definition of the Alfvén velocity we get for the magnetic field, B:

(5)

where MA is the Alfvén Mach number which, according to Vršnak et al. (2002), can be computed from:

(6)

with X being the square of the band-split frequency ratio.

During the initial phase of the first burst we measured a band-split frequency ratio of 1.23 ± 0.06, which corresponds to a density ratio of 1.5; applying the above method, we obtained MA = 1.18 ± 0.05 and a magnetic field strength of 18.4 ± 0.7 G for a 2 × 106 K isothermal coronal model with base density 2× the Newkirk model value. Here the errors refer to the rms of values measured in the DS at different frequencies and times. We note that the empirical relation of Dulk & McLean (1978) predicts B = 2.8 G at r = 1.32 R, which is the average radial distance to which our measurements refer, according to the model used. Our value is also higher than values reported in the literature (see Alissandrakis & Gary 2020). This difference could be attributed to the peculiarities of the event.

3.4.5. Polarization of the emission

We measured the circular polarization (Stokes parameter V) of the type-II burst from the NRH images. For the first type-II we found positive (right hand) polarization for both the harmonic and the fundamental emission. Taking into account the discussion in Sect. 3.4.1, the fundamental is polarized in the sense of the ordinary mode; the fact that the harmonic has the same sign of V implies that it is also polarized in the ordinary mode sense. The degree of polarization was measured to be about 35% for the fundamental and 65% for the harmonic. Normally the fundamental should be 100% polarized, but could be depolarized as a result of scattering and/or propagation effects (see Sect. 6.1 in Alissandrakis & Gary 2020). The second type-II burst, presumably emitting in the harmonic, was also polarized in the positive sense, the polarization degree ranging from 50% to 90%. We add that, with the exception of the pulsating source mentioned in Sect. 3.4.2, the entire event was polarized in the right hand sense.

In the literature there are diverse reports about type-II polarization, ranging from no polarization to strong polarization. Komesaroff (1958) and Roberts (1959) reported weak or no polarization, whereas Zlobec & Thejappa (1987) reported cases of both weak and strong polarization, noting that fine structures (herringbone) were more polarized than the burst continuum. Cairns & Robinson (1987) reported 15% polarization for herringbones, with the harmonic exhibiting weaker polarization than the fundamental. Hariharan et al. (2014) also reported low polarization (5−10%) in the harmonic for a single event, while Du et al. (2014) reported 30−60% for the fundamental in another single event. Thejappa et al. (2003) gave 50−90% polarization for a fragmented type-II burst at high frequencies, while strong polarization at the fundamental and weak for the harmonic was reported by Thejappa & MacDowall (2001) for IP type-IIs.

4. Discussion: Origin of type-II emission

The two type-II bursts were associated with an M9 flare and occurred in close time succession around the peaks of the two phases of the flare. Due to their timing with respect to the flare evolution, it is tempting to attribute their origin to pressure pulses generated by the flare. Our data did not allow us to perform a high-cadence emission measure analysis that could have directly revealed the possible existence of pressure pulses ignited by the flare and their relevance to the generation of the type-II bursts. Warmuth (2015) pointed out that flares igniting coronal waves are usually located at the periphery of their active regions and the waves are launched into directions away from the core of the active regions; this pattern facilitates the development and steepening of a pressure pulse since regions of low magnetic field strength and low Alfvén speeds are accessed. In our case the flare indeed occurred away from the core (defined as the location of the main polarity inversion line) of the active region (e.g., see panels 30−35 of Fig. 2) but the positions of the first type-II sources (see Fig. 8) are at odds with the suggestion of Warmuth (2015).

The positions of both type-II bursts are broadly consistent with the EUV wave whose generation and development were presented in Fig. 3. Other case-studies showing spatio-temporal association between EUV waves and metric type-II bursts include Kouloumvakos et al. (2014), Mancuso et al. (2019), Jebaraj et al. (2020), Koukras et al. (2020). There has been extensive literature about the drivers of EUV waves (see the reviews by Patsourakos & Vourlidas 2012; Warmuth 2015; Long et al. 2017, and references therein). As in most cases, the lateral propagation of our EUV wave makes it difficult to attribute its generation to the flare (e.g., Patsourakos & Vourlidas 2009). If we exclude the flare, we need an alternative low coronal driver in the form of mass motion. The only reliable candidate is the surge; from the timeline of the event it appears that the surge generated the disturbance that yielded the EUV wave. An example of an EUV wave linked to a surge has been presented by Zheng et al. (2013). Generally speaking, narrow eruptions (surges, jets etc.) could drive EUV waves (e.g., Warmuth 2015, and references therein).

This is somehow different from the scenario of Patsourakos & Vourlidas (2012) which identifies the lateral over-expansion of the CME cavity in the low corona as the EUV wave driver. However, our EUV data from both STEREO spacecraft do not show concrete evidence of a cavity. The same applies to the WL observations of the associated CME with COR1. It is possible that the line of sight was not favorable for revealing a complete CME cavity by, for example, not being aligned with its axis and/or not having enough mass along the line of sight to produce an appreciable signal in EUV or WL (see also Vourlidas et al. 2013, for a detailed discussion on how projection effects affect the appearance of CMEs). We therefore conclude that the scenario is not directly applicable to our event.

Evidence for the propagation of the disturbance was also provided by the streamer disruption while its evolution was detected as a white-light CME in the coronagraph data. The loop system located at a position angle of 56° north of the surge (see Sect. 3.3) was most probably disrupted by the wave but was not the driver of the wave. This conclusion is also supported by the following arguments: (1) At 11:28 UT (that is, around the start of the second phase of the event) the STEREO A images indicate that the EUV wave has propagated well beyond these loops to the north. (2) The disrupted loops expand along a direction that is markedly different from the direction defined by the positions of the second shock (see Fig. 8, bottom). (3) The amplitude of the loop oscillations appears too small to give rise to a major disturbance such as the one associated with the second shock.

Our results directly link the same EUV wave with the origin of both type-II bursts. This conclusion is different from the traditional interpretation of the origin of successive type-II bursts in the absence of successive CMEs which invokes a CME driver for one burst and a flare-related origin for the other (e.g., Shanmugaraju et al. 2005; Subramanian & Ebenezer 2006). On the other hand, the origin of multilane type-II bursts have occasionally been interepreted as coming from distinct locations of a single CME-driven shock (e.g., Feng et al. 2015; Lv et al. 2017).

In the light of the association of the EUV wave with the type-II bursts a few remarks are in order: (1) The association of two separate type-II bursts with the same EUV wave indicates that each type-II is related with different parts of the wave; supporting evidence for this argument is primarily provided by the positional offsets between the sources of the two type-II bursts and by the different rate of expansion across the wave. The three-dimensional structure of the Alfvén speed above the active region as well as the magnetic field configuration that would favor the acceleration of electrons at distinct locations of the shock front might have also contributed to the appearance of the type-II sources at the observed positions. Unfortunately our data were not suitable for a detailed study of these effects. (2) The speed of both type-II bursts computed from either the drift rates measured in the dynamic spectra or from the projected positions of the radio sources are systematically higher than the driver (that is, the surge which generated the EUV wave) speed (see Sect. 3.4.3). Combining the speed at the start of the first type-II burst (see Table 3) and the value of the Alfvén Mach number found in Sect. 3.4.4 we find an Alfvén speed of about 1300 km s−1. Therefore the driver propagates at sub-Alfvénic speeds. This situation is indicative of a piston-driven shock (e.g., see Vršnak & Cliver 2008), during the first stages of the event. However, the movies suggest that after some point in time the EUV wave attained large lateral distances from the surge and, therefore, the wave (shock) was then possibly freely propagating.

5. Summary and conclusions

In this work we used all available data, from kilometric radio waves to hard X-rays, in order to study the event SOL2010-02-012T11:25:00 of February 12, 2010 in conjunction with dynamic spectra from ARTEMIS-IV/JLS and images from the NRH. Our emphasis was on the origin of the two metric type-II bursts.

The low cadence of the non-radio data made difficult the association of the fast evolving metric emission with the phases of the event in other parts of the electromagnetic spectrum. Still, the NRH observations near the central meridian reduced the influence of refraction and scattering on the source positions, whereas observations from STEREO A and B provided three-dimensional information about the event. These gave to our study some important advantages compared to other similar works.

Based on computations for homogeneous models, we did not apply any refraction and scattering corrections to the position of the radio sources because the resulting corrections were too small to affect our conclusions. However, this issue warrants further investigation as has been done for the decametric wavelength range (e.g., Zhang et al. 2021); it is also important to investigate what happens in the (normal) inhomogeneous case. In addition to refraction and scattering, we ignored ionospheric effects, which, according to our measurements, were considerably smaller than the size of the radio sources.

The complex event originated in a very compact flare and was accompanied by a surge, a CME, a streamer disruption, and an EUV wave. Our analysis, discussed in detail in Sect. 4, showed that the two associated type-II bursts were not related either with the surge (see Fig. 14) or with the streamer disruption but were related with the EUV wave which was probably driven by the surge. We could not identify any nonthermal radio signatures of the surge or the streamer disruption; it is possible, however, that such signatures might be buried deep into the rich structure of the dynamic spectrum. Emission from the surge was detected at a later phase of the event, and this will be presented in a subsequent publication.

thumbnail Fig. 14.

Three-dimensional rendering of the surge and NRH source positions for branch A of the second type-II burst (blue and red crosses respectively), together with their projection on the basic planes, xy, xz and yz (blue and red lines). The black dot marks the position of the flare.

It is also important to discuss the possible effect of the adopted coronal model to the derived NRH source positions, through the model dependence of the emission height. To this end we computed heights for two more isothermal models, one (model 2) with a coronal temperature of 2 × 106 K and a base density equal to that of the Newkirk model and another (mode 3) with T = 1.4 × 106 K and a base electron density twice that of the Newkirk model. Although both models gave smaller source heights, they did not bring the radio sources any closer to the surge. As for the speed derived from the frequency drift for the second type-II, model 2 gave a value close to that of the adopted model whereas model 3 gave a 40% lower value, still more than a factor of two above the WL/EUV measurements.

In addition to the origin of the type-II emission, our analysis of the radio observations gave a number of interesting results that are worthy of further investigation:

Using cross-correlation we measured the frequency ratio and the time delay between emission at the plasma frequency and its harmonic during the first type-II burst; we found values of 1.95 and 0.23−0.6 s, respectively. This is a powerful technique that can be used systematically to derive precise values of these parameters and thus provide more accurate information on the formation of the fundamental as well as on wave propagation effects.

The first burst started with a relatively high frequency drift; subsequently the drift diminished to increase again later on. Chrysaphi et al. (2020) reported a case of transition of a type-II from stationary to drifting, but we are not aware of any report of drifting-stationary-drifting transition. Such a transition might be due to the shock encountering regions of different density scale lengths.

For the same burst we inferred a magnetic field of 18 G from band splitting, a rather high value compared to those reported in the literature. We note that the derived magnetic field is model-dependent, dropping to 14 G for model 2 (see above) and to 11 G for model 3.

We observed pulsations superposed on the second type-II burst. The pulsations were broad-band with a bandwidth of about 200 MHz and they started before the type-II. They had a measurable (reverse) drift of about 2 s−1. Moreover, thanks to the NRH images, we identified a second pulsating source of opposite circular polarization, which disappeared near the start of the burst. This is a potentially important finding, since we are not aware of any other report of multiple pulsation sources, but its full investigation is outside the scope of this work.

Finally, we would like to stress that this work confirms the importance of imaging spectroscopy for understanding the physics of solar radio emission.


Acknowledgments

The authors are grateful to the other members of the ARTEMIS group, in particular to C. Caroubalos, P. Preka-Papadema and X. Moussas (University of Athens), as well as to A. Kontogeorgos and P. Tsitsipis (University of Thessaly). We also wish to thank L. Klein and the entire group of the Nançay Radioheliograph for the original NRH visibility data, G. Chernov for interesting discussions and E. Kontar for the scattering computations. Data from STEREO (SECCHI and coronagraphs), SXI, Hinode, the Catania observatory, RHESSI, and GOES were obtained from the respective data bases; we are grateful to all those who contributed to the operation of these instruments and made the data available to the community. Finally we wish to thank the Onassis Foundation for financial support (Grant 15153) for the continued operation of the ARTEMIS-JLS radio spectrograph.

References

  1. Alissandrakis, C., & Gary, D. E. 2020, Front. Astron. Space Sci., 7, 77 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alissandrakis, C., Nindos, A., Patsourakos, S., Hillaris, A., & Artemis Group 2011, EPSM-13 Conference, Rhodes, Greece [Google Scholar]
  3. Alissandrakis, C. E., Bouratzis, C., & Hillaris, A. 2019, A&A, 627, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Armatas, S., Bouratzis, C., Hillaris, A., et al. 2019, A&A, 624, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bain, H. M., Krucker, S., Glesener, L., & Lin, R. P. 2012, ApJ, 750, 44 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bastian, T. S. 1994, ApJ, 426, 774 [Google Scholar]
  7. Biesecker, D. A., Myers, D. C., Thompson, B. J., Hammer, D. M., & Vourlidas, A. 2002, ApJ, 569, 1009 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bougeret, J. L., & Steinberg, J. L. 1977, A&A, 61, 777 [NASA ADS] [Google Scholar]
  9. Bouratzis, C., Hillaris, A., Alissandrakis, C. E., et al. 2015, Sol. Phys., 290, 219 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bouratzis, C., Hillaris, A., Alissandrakis, C. E., et al. 2016, A&A, 586, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Cairns, I. H. 2011, in Coherent Radio Emissions Associated with Star System Shocks, eds. M. P. Miralles, & J. Sánchez Almeida, 4, 267 [NASA ADS] [Google Scholar]
  12. Cairns, I. H., & Knock, S. A. 2006, Planetary Radio Emissions VI, 419 [Google Scholar]
  13. Cairns, I. H., & Robinson, R. D. 1987, Sol. Phys., 111, 365 [Google Scholar]
  14. Cane, H. V., & Erickson, W. C. 2005, ApJ, 623, 1180 [NASA ADS] [CrossRef] [Google Scholar]
  15. Chernov, G. P., Fomichev, V. V., Gorgutsa, R. V., et al. 2014, Geomagnet. Aeron., 54, 406 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cho, K.-S., Lee, J., Moon, Y.-J., et al. 2007, A&A, 461, 1121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cho, K.-S., Gopalswamy, N., Kwon, R.-Y., Kim, R.-S., & Yashiro, S. 2013, ApJ, 765, 148 [NASA ADS] [CrossRef] [Google Scholar]
  18. Chrysaphi, N., Kontar, E. P., Holman, G. D., & Temmer, M. 2018, ApJ, 868, 79 [NASA ADS] [CrossRef] [Google Scholar]
  19. Chrysaphi, N., Reid, H. A. S., & Kontar, E. P. 2020, ApJ, 893, 115 [NASA ADS] [CrossRef] [Google Scholar]
  20. Claßen, H. T., & Aurass, H. 2002, A&A, 384, 1098 [Google Scholar]
  21. Dauphin, C., Vilmer, N., & Krucker, S. 2006, A&A, 455, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Démoulin, P., Vourlidas, A., Pick, M., & Bouteille, A. 2012, ApJ, 750, 147 [Google Scholar]
  23. Du, G., Chen, Y., Lv, M., et al. 2014, ApJ, 793, L39 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dulk, G. A., & McLean, D. J. 1978, Sol. Phys., 57, 279 [CrossRef] [Google Scholar]
  25. Eselevich, V. G., Eselevich, M. V., Sadykov, V. M., & Zimovets, I. V. 2015, Adv. Space Res., 56, 2793 [NASA ADS] [CrossRef] [Google Scholar]
  26. Feng, S. W., Chen, Y., Kong, X. L., et al. 2012, ApJ, 753, 21 [NASA ADS] [CrossRef] [Google Scholar]
  27. Feng, S. W., Du, G. H., Chen, Y., et al. 2015, Sol. Phys., 290, 1195 [NASA ADS] [CrossRef] [Google Scholar]
  28. Gary, D. E., Dulk, G. A., House, L., et al. 1984, A&A, 134, 222 [NASA ADS] [Google Scholar]
  29. Gergely, T. E., Kundu, M. R., & Hildner, E. 1983, ApJ, 268, 403 [CrossRef] [Google Scholar]
  30. Gopalswamy, N. 2006, Washington DC American Geophysical Union Geophysical Monograph Series, 165, 207 [NASA ADS] [Google Scholar]
  31. Gopalswamy, N., Aguilar-Rodriguez, E., Yashiro, S., et al. 2005, J. Geophys. Res.: Space Phys., 110, A12S07 [NASA ADS] [CrossRef] [Google Scholar]
  32. Gopalswamy, N., Thompson, W. T., Davila, J. M., et al. 2009, Sol. Phys., 259, 227 [NASA ADS] [CrossRef] [Google Scholar]
  33. Hariharan, K., Ramesh, R., Kishore, P., Kathiravan, C., & Gopalswamy, N. 2014, ApJ, 795, 14 [NASA ADS] [CrossRef] [Google Scholar]
  34. Howard, R. A., Moses, J. D., Vourlidas, A., et al. 2008, Space Sci. Rev., 136, 67 [NASA ADS] [CrossRef] [Google Scholar]
  35. Jebaraj, I. C., Magdalenić, J., Podladchikova, T., et al. 2020, A&A, 639, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Kerdraon, A., & Delouis, J. M. 1997, in Coronal Physics from Radio and Space Observations, ed. G. Trottet (Berlin: Springer-Verlag), Lect. Notes Phys., 483, 192 [Google Scholar]
  37. Klassen, A., Aurass, H., Mann, G., & Thompson, B. J. 2000, A&AS, 141, 357 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Klein, K.-L., Khan, J. I., Vilmer, N., Delouis, J.-M., & Aurass, H. 1999, A&A, 346, L53 [Google Scholar]
  39. Komesaroff, M. 1958, Aust. J. Phys., 11, 201 [NASA ADS] [CrossRef] [Google Scholar]
  40. Kong, X. L., Chen, Y., Li, G., et al. 2012, ApJ, 750, 158 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kontar, E. P., Chen, X., Chrysaphi, N., et al. 2019, ApJ, 884, 122 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kontogeorgos, A., Tsitsipis, P., Caroubalos, C., et al. 2006, Exp. Astron., 21, 41 [NASA ADS] [CrossRef] [Google Scholar]
  43. Koukras, A., Marqué, C., Downs, C., & Dolla, L. 2020, A&A, 644, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Kouloumvakos, A., Patsourakos, S., Hillaris, A., et al. 2014, Sol. Phys., 289, 2123 [NASA ADS] [CrossRef] [Google Scholar]
  45. Kumar, P., Innes, D. E., & Cho, K.-S. 2016, ApJ, 828, 28 [Google Scholar]
  46. Kumari, A., Ramesh, R., Kathiravan, C., & Gopalswamy, N. 2017, ApJ, 843, 10 [Google Scholar]
  47. Long, D. M., Bloomfield, D. S., Chen, P. F., et al. 2017, Sol. Phys., 292, 7 [CrossRef] [Google Scholar]
  48. Lv, M. S., Chen, Y., Li, C. Y., et al. 2017, Sol. Phys., 292, 194 [NASA ADS] [CrossRef] [Google Scholar]
  49. Magdalenić, J., Vršnak, B., Pohjolainen, S., et al. 2008, Sol. Phys., 253, 305 [Google Scholar]
  50. Magdalenić, J., Marqué, C., Zhukov, A. N., Vršnak, B., & Žic, T. 2010, ApJ, 718, 266 [Google Scholar]
  51. Magdalenić, J., Marqué, C., Zhukov, A. N., Vršnak, B., & Veronig, A. 2012, ApJ, 746, 152 [Google Scholar]
  52. Maguire, C. A., Carley, E. P., Zucca, P., Vilmer, N., & Gallagher, P. T. 2021, ApJ, 909, 2 [NASA ADS] [CrossRef] [Google Scholar]
  53. Maia, D., Pick, M., Vourlidas, A., & Howard, R. 2000, ApJ, 528, L49 [CrossRef] [Google Scholar]
  54. Mancuso, S., Frassati, F., Bemporad, A., & Barghini, D. 2019, A&A, 624, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Mann, G., Classen, T., & Aurass, H. 1995, A&A, 295, 775 [Google Scholar]
  56. Mann, G., Klassen, A., Classen, H. T., et al. 1996, A&AS, 119, 489 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Maxwell, A., Hughes, M. P., & Thompson, A. R. 1963, J. Geophys. Res., 68, 1347 [NASA ADS] [CrossRef] [Google Scholar]
  58. McLean, D. J. 1967, PASA, 1, 47 [Google Scholar]
  59. Morosan, D. E., Carley, E. P., Hayes, L. A., et al. 2019, Nat. Astron., 3, 452 [Google Scholar]
  60. Nelson, G. J., & Melrose, D. B. 1985, in Type II Bursts, eds. D. J. McLean, & N. R. Labrum, 333 [Google Scholar]
  61. Nindos, A., & Aurass, H. 2007, in Pulsating Solar Radio Emission, eds. K. L. Klein, & A. L. MacKinnon, 725, 251 [NASA ADS] [Google Scholar]
  62. Nindos, A., Aurass, H., Klein, K.-L., & Trottet, G. 2008, Sol. Phys., 253, 3 [NASA ADS] [CrossRef] [Google Scholar]
  63. Nindos, A., Alissandrakis, C. E., Hillaris, A., & Preka-Papadema, P. 2011, A&A, 531, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Nitta, N. V., Liu, W., Gopalswamy, N., & Yashiro, S. 2014, Sol. Phys., 289, 4589 [NASA ADS] [CrossRef] [Google Scholar]
  65. Patsourakos, S., & Vourlidas, A. 2009, ApJ, 700, L182 [Google Scholar]
  66. Patsourakos, S., & Vourlidas, A. 2012, Sol. Phys., 281, 187 [NASA ADS] [Google Scholar]
  67. Patsourakos, S., Vourlidas, A., Török, T., et al. 2020, Space Sci. Rev., 216, 131 [Google Scholar]
  68. Pick, M., & Vilmer, N. 2008, A&ARv, 16, 1 [NASA ADS] [CrossRef] [Google Scholar]
  69. Pick, M., Stenborg, G., Démoulin, P., Zucca, P., & Lecacheux, A. 2016, ApJ, 823, 5 [NASA ADS] [CrossRef] [Google Scholar]
  70. Pohjolainen, S., Pomoell, J., & Vainio, R. 2008, A&A, 490, 357 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Ramesh, R., Lakshmi, M. A., Kathiravan, C., Gopalswamy, N., & Umapathy, S. 2012, ApJ, 752, 107 [NASA ADS] [CrossRef] [Google Scholar]
  72. Roberts, J. A. 1959, Aust. J. Phys., 12, 327 [Google Scholar]
  73. Robinson, R. D., & Sheridan, K. V. 1982, PASA, 4, 392 [NASA ADS] [Google Scholar]
  74. Schmidt, J. M., & Cairns, I. H. 2012, J. Geophys. Res.: Space Phys., 117, A04106 [NASA ADS] [Google Scholar]
  75. Shanmugaraju, A., Moon, Y. J., Cho, K. S., et al. 2005, Sol. Phys., 232, 87 [NASA ADS] [CrossRef] [Google Scholar]
  76. Smerd, S. F., Sheridan, K. V., & Stewart, R. T. 1974, in Coronal Disturbances, ed. G. A. Newkirk, IAU Symp., 57, 389 [NASA ADS] [CrossRef] [Google Scholar]
  77. Steinberg, J. L., Aubier-Giraud, M., Leblanc, Y., & Boischot, A. 1971, A&A, 10, 362 [NASA ADS] [Google Scholar]
  78. Subramanian, K. R., & Ebenezer, E. 2006, A&A, 451, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Temmer, M., Veronig, A. M., Vršnak, B., et al. 2008, ApJ, 673, L95 [Google Scholar]
  80. Thejappa, G., & MacDowall, R. J. 2001, J. Geophys. Res., 106, 25313 [NASA ADS] [CrossRef] [Google Scholar]
  81. Thejappa, G., Zlobec, P., & MacDowall, R. J. 2003, ApJ, 592, 1234 [NASA ADS] [CrossRef] [Google Scholar]
  82. Thompson, B. J., Plunkett, S. P., Gurman, J. B., et al. 1998, Geophys. Res. Lett., 25, 2465 [Google Scholar]
  83. Vourlidas, A., Lynch, B. J., Howard, R. A., & Li, Y. 2013, Sol. Phys., 284, 179 [NASA ADS] [Google Scholar]
  84. Vršnak, B., & Cliver, E. W. 2008, Sol. Phys., 253, 215 [Google Scholar]
  85. Vršnak, B., Ruzdjak, V., Zlobec, P., & Aurass, H. 1995, Sol. Phys., 158, 331 [Google Scholar]
  86. Vršnak, B., Aurass, H., Magdalenić, J., & Gopalswamy, N. 2001, A&A, 377, 321 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Vršnak, B., Magdalenić, J., Aurass, H., & Mann, G. 2002, A&A, 396, 673 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Vršnak, B., Magdalenić, J., & Zlobec, P. 2004, A&A, 413, 753 [Google Scholar]
  89. Vršnak, B., Warmuth, A., Temmer, M., et al. 2006, A&A, 448, 739 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Warmuth, A. 2015, Liv. Rev. Sol. Phys., 12, 3 [Google Scholar]
  91. White, S., Mercier, C., Bastian, T., & Bradley, R. 2007, NRAO Newslett., 110, 9 [Google Scholar]
  92. Wild, J. P., Murray, J. D., & Rowe, W. C. 1954, Aust. J. Phys., 7, 439 [NASA ADS] [CrossRef] [Google Scholar]
  93. Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452 [Google Scholar]
  94. Zhang, J., Dere, K. P., Howard, R. A., & Vourlidas, A. 2004, ApJ, 604, 420 [NASA ADS] [CrossRef] [Google Scholar]
  95. Zhang, P., Wang, C., & Kontar, E. P. 2021, ApJ, 909, 195 [NASA ADS] [CrossRef] [Google Scholar]
  96. Zheleznyakov, V. V. 1996, in Radiation in Astrophysical Plasmas, Astrophys. Space Sci. Lib., 204 [Google Scholar]
  97. Zheng, R., Jiang, Y., Yang, J., et al. 2013, ApJ, 764, 70 [NASA ADS] [CrossRef] [Google Scholar]
  98. Zimovets, I. V., & Sadykov, V. M. 2015, Adv. Space Res., 56, 2811 [NASA ADS] [CrossRef] [Google Scholar]
  99. Zimovets, I., Vilmer, N., Chian, A. C. L., Sharykin, I., & Struminsky, A. 2012, A&A, 547, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Zlobec, P., & Thejappa, G. 1987, Hvar Obs. Bull., 11, 111 [Google Scholar]
  101. Zucca, P., Pick, M., Démoulin, P., et al. 2014a, ApJ, 795, 68 [NASA ADS] [CrossRef] [Google Scholar]
  102. Zucca, P., Carley, E. P., Bloomfield, D. S., & Gallagher, P. T. 2014b, A&A, 564, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Zucca, P., Morosan, D. E., Rouillard, A. P., et al. 2018, A&A, 615, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1.

Time line of the event of February 12, 2010.

Table 2.

Measurements of velocities from WL and EUV images.

Table 3.

Type-II drift rates and nominal velocities.

All Figures

thumbnail Fig. 1.

Dynamic spectrum (negative) of the February 12, 2010 event (500 to 20 MHz), observed with the ASG receiver of the ARTEMIS-JLS radiospectrograph, together with light curves in soft X-rays from GOES and hard X-rays from RHESSI in arbitrary units and on a linear scale. Dashed horizontal lines mark the frequencies of the Nançay Radioheliograph. 1F and 1H mark the fundamental and harmonic emission of the first type-II burst, A, B1 and B2 mark the branches of the second type-II.

In the text
thumbnail Fig. 2.

Evolution of the event in Hα (Catania Observatory), soft X-rays (GOES/SXI and Hinode/XRT), 195 and 304 Å (STEREO A and B), together with a magnetogram and white light image from MDI. The white arc in the STEREO images marks the solar limb. The center of the field of view is at Carrington longitude 186.6, latitude 25.0 (solar x = −107″, y = 510″ for Earth images).

In the text
thumbnail Fig. 3.

Running difference images at 195 Å from STEREO A and B. We note the formation of an expanding EUV wave. The white arc marks the limb.

In the text
thumbnail Fig. 4.

COR1 images and their base difference (top and middle row) and 304 Å images, saturated to show the surge better. We note the appearance of the CME above the occulting disk at 1.6 solar radii at 11:35:21 UT, the disruption of the streamer (marked SD) and the surge (E) seen both in the coronagraph and the 304 Å images. The dark circle marks the limb. The field of view is 2500 by 2500″.

In the text
thumbnail Fig. 5.

Coronagraph observations. Top: COR1A intensity as a function of position angle and distance from the center of the disk, r. Middle: same as top, with the pre-CME intensity subtracted. Bottom: cut of intensity as a function of time and r at a position angle of 61°. Horizontal lines are at the same r in all panels. Vertical lines in the top and middle panel mark the position angle of the cut. Bottom panel: vertical line marks the time of the other panels. The photospheric limb is at r = 972″.

In the text
thumbnail Fig. 6.

Projection of the surge on the sky plane as seen from the Earth (symbols), on top of an Hα image (left) and an XRT image (right). Different symbols refer to different parts of the surge, giving an idea of its lateral extent.

In the text
thumbnail Fig. 7.

Position of the type-II sources, together with the projected position of the surge (open circles) and the position of the Hα/SXR flare (large open circle). Left: first type-II, individual images; black is for 360.0 MHz, red for 327.0 MHz, green for 297.8 MHz and blue for 150.9 MHz. Right: average images for branches A (black), B1 (red) and B2 (blue) of the second type-II burst.

In the text
thumbnail Fig. 8.

Projected positions of the type-II radio sources, on top of polar projections of running difference images from 195 Å STEREO-A. Top: type-II during the first phase of the event (11:27:00−11:28:05 UT); asterisks show the average positions from 360.8 MHz, “x” from 298.7 MHz, crosses from 298.7 MHz and open circles from 150.9 MHz. Bottom: type-II during the second phase of the event (11:29:21−11:30:57 UT); asterisks are for branch A, crosses for branch B1 and open circles for branch B2.

In the text
thumbnail Fig. 9.

Fundamental-harmonic emission during the first phase of the event (composite ASG/SAO dynamic spectrum). The frequency range is from 106 to 450 MHz. Dotted horizontal lines mark the NRH frequencies. Channels around 200 and 310 MHz have been partly corrected for radio interference. The arrow points to a spectral feature discussed in the text.

In the text
thumbnail Fig. 10.

Expanded view of the spectral feature marked by the arrow in Fig. 9 (top panel, same for both columns). The lower rows show one-dimentional NRH images at three frequencies in the EW (left) and NS (right) directions.

In the text
thumbnail Fig. 11.

Composite ASG/SAO dynamic spectrum of the second type-II burst.

In the text
thumbnail Fig. 12.

Expanded views of pulsations. Top: pulsations near the start of the second phase of the event. Bottom: DS of the first part of the pulsations recorded with the ASG (left) and the SAO (right), together with 1-dimensional EW images from the NRH placed at the corresponding DS frequencies (middle). The frequency range is 220−430 MHz; time tickmarks are every 2 s in the left and middle panels, every 1 s in the right panel.

In the text
thumbnail Fig. 13.

Pulsations before the start of the second type-II; Dynamic spectrum (top) and 1D NRH images at 270.6 MHz in Stokes I and V, EW and NS.

In the text
thumbnail Fig. 14.

Three-dimensional rendering of the surge and NRH source positions for branch A of the second type-II burst (blue and red crosses respectively), together with their projection on the basic planes, xy, xz and yz (blue and red lines). The black dot marks the position of the flare.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.