A&A 413, 1125-1133 (2004)
DOI: 10.1051/0004-6361:20031497

Solar type II and type IV radio bursts observed during 1998-2000 with the ARTEMIS-IV radiospectrograph[*]

C. Caroubalos 1 - A. Hillaris2 - C. Bouratzis 2 - C. E. Alissandrakis 3 - P. Preka-Papadema 2 - J. Polygiannakis 2 - P. Tsitsipis 4 - A. Kontogeorgos 4 - X. Moussas 2 - J.-L. Bougeret 5 - G. Dumas 5 - C. Perche 5


1 - Department of Informatics, University of Athens, 15783 Athens, Greece
2 - Section of Astrophysics, Astronomy and Mechanics, Department of Physics, University of Athens, 15784 Panepistimiopolis Zografos, Athens, Greece
3 - Section of Astro-Geophysics, Department of Physics, University of Ioannina, 45110 Ioannina, Greece
4 - Department of Electronics, Technological Education Institute of Lamia, Lamia, Greece
5 - Observatoire de Paris, Departement de Recherche Spatiale, CNRS UA 264, 92195 Meudon Cedex, France

Received 2 April 2002 / Accepted 29 August 2003

Abstract
A catalogue of the type II and type IV solar radio bursts in the 110-687 MHz range, observed with the radio spectrograph ARTEMIS-IV operated by the University of Athens at Thermopylae, Greece from 1998-2000 is presented. These observations are compared with the LASCO archives of Coronal Mass Ejections and the Solar Geophysical Reports of solar flares (H$\alpha $ & SXR) and examined for possible associations. The main results are:

Key words: astronomical data bases: miscellaneous - catalogs - Sun: activity - Sun: corona - Sun: flares - Sun: radio radiation

1 Introduction

Type II bursts represent the passage of a MHD shock wave through the tenuous plasma of the solar corona (Wild & Smerd 1972); their radio emission is due either to energetic electrons accelerated at the shock front or plasma turbulence excited by the shock. Although they are associated with eruptive phenomena and are identified either with a flare blast wave (Vrsnak et al. 2001) or with a CME forward shock (Kahler et al. 1984; Maia et al. 2000; Aurass 1997) or with a shock driven by the flanks of a CME (Classen & Aurass 2002), their exact origin has not been unambiguously determined (Klassen et al. 1999); furthermore Pearson et al. (1989) question the type II association with the most energetic and impulsive flares, demonstrating an almost equal probability of type II association with smaller flares as well.

Recent observations place the origin in rapidly expanding or disrupting magnetic structures in the outskirts of flaring active regions, thus stressing the association of metric and decametric type II bursts with flares rather than CMEs (Klein et al. 1999); Classen & Aurass (2002) on the other hand, establish an almost equal probability for a type II to be driven by a flare blast, the front of a CME or its flank; they categorise accordingly type II as class 1, class 2 and class 3.

The dynamic spectra of type II events appear as slowly drifting bands, from high to low frequencies; they often exhibit a characteristic fundamental-harmonic structure (Maxwell & Thompson 1962), as well as significant fine structure (Mann et al. 1996; Slottje 1981) dubbed herringbone or backbone.

More often than not, the appearance of type II radio bursts is preceded by groups of bursts of the type III family. On the dynamic spectra the onset of the type III activity seems be in the same spectral range as the backward extrapolated type II lanes; Wild et al. (1963) named this type of event compound and concluded that both the type III and the type II bursts originated from a common source. Recently, Klassen et al. (1999) reported two more spectral features physically related to the type II exciter: the type II precursor consists of fast drift bursts or pulsations with a restricted bandwidth near the spectral range of the backward extrapolated type II lanes; the arc consists of a series of narrow-band bursts immediately preceding the type II onset; it has an inverted U-shaped envelope, hence its name.

The continua observed during periods of activity, on the other hand, represent the radiation of energetic electrons trapped within magnetic structures and plasmoids and they appear under the name type IV bursts (Kundu et al. 2001). The stationary type IV bursts (IV mB) emanate from magnetic structures usually located above active regions (Robinson 1985); they often present significant fine structure (Slottje 1981). The moving type IV bursts (Boischot 1957) are emitted from sources of meter wave continuum which are believed to move outwards at velocities of the order of 100-1000 km s-1; their spectrum is often featureless and sometimes last more than 10 min. Some of them appear following type II bursts and are possibly caused by energetic electrons produced in the wake of the type II shock. Others may originate from energetic electron populations trapped in expanding magnetic arches or plasmoids, i.e. blobs of dense plasma containing their own magnetic field (cf. Wild & Smerd 1972; Stewart 1985 and references therein). A number of moving type IV (IV mA) bursts are believed to originate within the densest substructures of CMEs (Aurass 1997; Aurass et al. 1999); the work by Klein & Mouradian (2002), Bastian et al. (2001) provides more evidence which corroborates this point of view, associating the type IV source with the erupting prominences within the CMEs.

Coronal Mass Ejections (CMEs) are solar transients characterised by the expulsion of solar plasma and frozen in magnetic field from the corona to the interplanetary space; their driving mechanism is probably magnetic (Amari et al. 2000; Vourlidas et al. 2000; Wu et al. 2001). The association of CMEs with other types of solar activity has been studied meticulously (cf. for example Hildner et al. 1976 on the CME-Sunspot Number association, Webb & Hundhausen 1987 and Munro et al. 1979 on the association with Eruptive Prominences and long duration events (LDEs) (also Gilbert et al. 2000; Kahler et al. 1989), flares (cf. also Svestka 2001), radio bursts (cf. also Maia et al. 2000)), flares and filament eruptions (Guiping Zhou et al. 2003), solar energetic electron events and coronal shocks (Klassen et al. 2002), but the results are not conclusive. Therefore both the CME origin and their association with active solar features are open research issues since our understanding of both is still uncertain. Catalogues of recorded CMEs, some including associated phenomena, have already been compiled and published (cf. for example Yashiro et al. 2001; Gil Moreno et al. 1998).

In this article we present a catalogue of all the type II and IV solar radio bursts observed during 1998-2000 by the ARTEMIS-IV Solar Radio spectrograph, including associated CMEs and flares from the LASCO archives and the Solar Geophysical Data (SGD) respectively. Throughout this report, type II bursts followed by type IV continua are referred to as type II/IV events. In Sect. 2 we describe our instrumentation and the observations used in the compilation of our catalogue; in Sect. 3 we present our results based on this catalogue. Discussion of results and conclusions appear in Sect. 4, while the catalogue is given in the appendix.

2 Observations

2.1 Instrumentation


  \begin{figure}
\par\resizebox{18cm}{!}{\includegraphics[clip]{00208_01.EPS}}\par\resizebox{18cm}{!}{\includegraphics[clip]{00208_02.EPS}}
\end{figure} Figure 1: ASG dynamic spectrum of the 2000 February 8 type II/IV event (#23 in our catalog), preceded by a type III group with time resolution 1.0 s. The bursts of different types have been annotated on image. Upper panel: intensity spectrum, Lower panel: differential spectrum (time derivative of intensity).
Open with DEXTER


  \begin{figure}
\par\resizebox{18cm}{!}{\includegraphics[clip]{00208_25.EPS}}\par\resizebox{18cm}{!}{\includegraphics[clip]{00208_26.EPS}}
\end{figure} Figure 2: SAO dynamic spectrum of the type IV burst of the 2000 February 8 complex event shown in Fig. 1. Upper panel: intensity spectrum, Lower panel: differential spectrum, which enhances the fine structure (pulsating) of the type IV continuum.
Open with DEXTER


  \begin{figure}
\par\resizebox{12cm}{!}{\includegraphics[clip]{99630_C.eps}}
\end{figure} Figure 3: Composite dynamic spectrum of the 1999 June 30 type II/IV event (#18 in our catalog), preceded by a type III group, in the frequency range 687 MHz to 1 KHz with 18 s time resolution. Upper panel: WIND/WAVES spectrum, showing the type III in the interplanetary space; Middle panel: Nançay Decametric Array spectrum; Lower panel: ARTEMIS-IV ASG spectrum (cf. Fig. 4).
Open with DEXTER


  \begin{figure}
\par\resizebox{14cm}{!}{\includegraphics[clip]{99630_01.EPS}}\par\resizebox{14cm}{!}{\includegraphics[clip]{99630_02.EPS}}
\end{figure} Figure 4: Expanded ASG dynamic spectrum of the event shown in Fig. 3, with time resolution of 2.0 s. The bursts of different types have been annotated on the figure. Upper panel: intensity spectrum, Lower panel: differential spectrum (time derivative of intensity).
Open with DEXTER


  \begin{figure}
\par\resizebox{14cm}{!}{\includegraphics[clip]{99630_06.EPS}}
\end{figure} Figure 5: SAO differential spectrum of the 1999 June 30 type II/IV event (shown in Figs. 3 and 4), with time resolution 0.5 s. Intermediate drift bursts (fibers, see for example 11:29-11:31 UT) and pulsations have been enhanced by the differentiation which has suppressed the underlying continuum.
Open with DEXTER

Our radio data were obtained with the ARTEMIS-IV solar radio spectrograph (Caroubalos et al. 2001a), operating at Thermopylae, Greece (38$^\circ$49$^\prime$N, 22$^\circ$41$^\prime$E). The daily observations cover the time interval 05:30-15:00 UT. The instrument covers the frequency range of 110 to 687 MHz, using two receivers operating in parallel. The Global Spectral Analyser (ASG), a sweep frequency receiver and the Acousto-Optic Spectrograph (SAO), a multichannel acousto-optical receiver. The sweep frequency analyser (ASG) covers the full frequency band with a time resolution of 10 samples/s. The high sensitivity multi-channel acousto-optical analyser covers the 265-450 MHz range, with high frequency and time resolution (100 samples/s); its recordings are used, mostly, for the study of fine structures (cf. for example Caroubalos et al. 2001b).

2.2 Event List

In assembling the catalogue we included all the type II, IV and II/IV meter-wavelength events recorded by ARTEMIS-IV ASG during 1998-2000; during this three-year period the instrument was out of operation for technical reasons for about 12 months (see table in the appendix for the exact dates), thus the catalogue covers about twenty four months of observations. The dynamic spectra shown here are uncalibrated, but the minimum value has been subtracted from each channel. The SAO dynamic spectra in differential form were examined for type IV fine structure classification, as appropriate, when the recorded event, or part of it, was within the 265-450 MHz range. We note that in the fine structure detection and classification we use differential spectra, which suppress the underlying continuum, enhancing the structures with shorter time scales (examples of dynamic spectra appear in Figs. 1 to 5).

Since a number of events exceeded the spectral range of ASG, we supplemented the information in our list with recordings and spectra from the Izmiran Radio spectrograph[*], the Nançay Decametric Array[*] and AIP (Potsdam) which were available on the Internet; these were used in the classification of recorded type II bursts as fundamental or harmonic, or for the detection of type II bursts originating outside the ARTEMIS-IV spectral range, yet accompanying continua of the ARTEMIS-IV catalogue. In Fig. 3 we present an example of an extended dynamic spectrum in the 687 MHz-1 KHz range, using combined spectral data from ARTEMIS-IV, the Nançay Decametric Array and the WIND/WAVES[*] receivers (for a similar composite spectrum cf. Bougeret et al. 1998).

Most of the recorded type II bursts were in the first harmonic, as confirmed by the corresponding dynamic spectra of Izmiran and the 40-800 MHz sweep spectrometer of the (AIP) Astrophysical Institute of Potsdam (Mann et al. 1992). The type II shock velocities  $V_{\rm drift}$ were estimated from the observed frequency drift rate using the Newkirk (Newkirk 1961, 1967) density ($N_{\rm e}$, in cm-3) vs. height (R, in $R_{\hbox{$\odot$ }}$) coronal model:

\begin{displaymath}N_{\rm e} = 8.3 \times 10^4 \times 10^{4.32/R} = 8.3 \times 10^4 \exp \left(\frac{{9.947}}{R}\right)
\end{displaymath}

and the equation (Tlamicha & Karlicky 1976):

\begin{displaymath}\frac{1}{f}\frac{{{\rm d}f}}{{{\rm d}t}} = - V_{\rm drift} \frac{{[\ln (0.037f^2 )]^2 }}{{19.9}}
\end{displaymath}

where f is the frequency in MHz, $V_{\rm drift}$ is in $R_{\hbox{$\odot$ }}$/s (Solar radii per second) and we have assumed radial propagation of the MHD shock and emission at the first harmonic of the local Langmuir frequency. Under the same assumption (radial propagation, Newkirk Corona, radiation at the first harmonic) we estimated by extrapolation the type II launch time, using its velocity and the height corresponding to its appearance on our recordings.

Then the LASCO coronagraph (Brueckner et al. 1995) event list was examined. The CME velocities and start times were estimated from the LASCO movies using the linear regression by Yashiro & Michalek (Yashiro et al. 2001) and are available in the LASCO event list[*]. Following Classen & Aurass (2002), we have adopted an error of 100 km s-1 in the CME velocity estimate which seems to cover both projection and acceleration or deceleration effects.

Lastly, the NOAA Solar Geophysical Data catalogues[*] of GOES SXR and H$\alpha $ flares were examined for possible associations. The H$\alpha $ flares were, at times, reported by more than one observatory. In those cases we chose the observation with the largest duration.

As regards the compilation of the type II, IV and II/IV ARTEMIS-IV 1998-2000 catalogue we adopted spatial and temporal criteria for the selection of the associated CMEs and flares (SXR and H$\alpha $) as follows:

At this point we note that for a number of events the spatial criterion fails to lead to the acceptance or rejection of an association. This is due to either the lack of positional data for the radio burst or the flare or to the appearance of more than one accompanying flares as candidates for association. As regards the use of time window for temporal association, as already pointed out by Webb & Hundhausen (1987), it represents a compromise between rejection of significant associations and inclusion of spurious ones, as for every threshold-based selection process. To mediate this effect we use the above-mentioned combination of criteria.

3 Results

In the table of the appendix we include all solar type II and IV radio bursts recorded by ARTEMIS-IV in the 110-687 MHz range, during 1998-2000; we present the type II and IV burst parameters as well as the parameters of the associated CMEs, H$\alpha $ flares and SXR GOES flares.

Forty events were recorded by ARTEMIS-IV in the period 1998-2000, of which 8 type IV continua, 11 type II bursts and 21 type II bursts followed by a type IV continuum (Type II/IV). Among the type II/IV were the radio signatures of the July 14, 2000 major solar event (event 36; see Klein et al. 2001; Caroubalos et al. 2001b), the May 2, 1998 event (event 06; see Pohjolainen et al. 2001) and four radio rich CMEs (i.e. CMEs associated with a metric or decametric type II burst) from Gopalswamy et al. (2000) (their Table 1). The gross spectral characteristics and the associated CME and flare parameters are summarised in the appendix. We note that all but one of the forty events were associated with either an H$\alpha $ and/or a SXR (GOES) event. With regard to the CME-type II and CME-type IV association, our sample of event shows that:

In Fig. 6 we compare velocities of CMEs and associated type II drift velocities. To take into account all errors in the estimation of both, we adopted the criterion introduced by Classen & Aurass (2002), which states that, a velocity ratio in the 1/2-2 range may qualify the type II MHD shock as originating in the CME front (class 2 event); all but three of the CME associated type IIs in our data set seem to belong in this category.

In Fig. 7 we show histograms of velocities of CMEs associated with type II events as well as of CMEs not associated with type IIs. The former have an average velocity of ( $835 \pm 380$) km s-1while the latter of ( $500 \pm 150$) km s-1. Similar results are reported by Gopalswamy et al. (2000) who maintain that radio rich CMEs tend to be faster than average. The systematically lower velocity range, yet with some overlap, for the non type II associated CMEs implies that, although a relatively fast CME may be required for the formation of a MHD shock, still ambient medium parameters (density and magnetic field) are equally important since the CME needs exceed the Alfven velocity.


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[clip]{CMETII01.EPS}}
\end{figure} Figure 6: Velocities of CMEs in km s-1 ( $V_{\rm CME}$) versus drift velocities of associated type II bursts ( $V_{\rm drift}$) in km s-1. The lines mark: $V_{\rm CME}=2~V_{\rm drift}$, $V_{\rm CME}=V_{\rm drift}$, 2  $V_{\rm CME}=V_{\rm drift}$; a type II velocity in this range probably qualifies the MHD shock as piston driven by the CME front (see text).
Open with DEXTER


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[clip]{CMEV_01.eps}}
\end{figure} Figure 7: Velocity histograms of type II associated CMEs and CMEs without a type II association.
Open with DEXTER

In Fig. 8 we present velocities of CMEs versus the total integrated flux of the associated SXR flare (1-8 A GOES channel); the velocity increases with increasing total SXR flare output, approximately following the empirical relationship $V_{\rm CME}=150~\ln(F_{\rm SXR})+1384$; this may imply a certain correlation between the extent of the magnetic restructuring linked to a CME launch and the extent of the reconnection releasing the free energy necessary for a flare. However, a cause and effect mechanism cannot be observationally established under the circumstances.

In Fig. 9 we summarize the total activity, including estimated CME take off time, SXR flare and type II, III and IV metric radio bursts, of the CME associated events of the ARTEMIS-IV catalogue. The origin of time coincides, for each event, with the time of the SXR flare maximum. We note that most of the CME launch times precede by 5-60 min (average 30 min) the SXR peak, and in most bursts the onset of the radio bursts. Type III groups precede most CME associated type II and type II/IV bursts.


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[clip]{CMEV_02.EPS}}
\end{figure} Figure 8: Velocities of CMEs in km s-1 ( $V_{\rm CME}$) versus integrated flux ( $F_{\rm SXR}$) of associated SXR flares in J/m2. The data are fitted with $V_{\rm CME}=150~\ln(F_{\rm SXR})+1384$.
Open with DEXTER


  \begin{figure}
\par\resizebox{8.8cm}{!}{\includegraphics[clip]{TOTAL_01.eps}}
\end{figure} Figure 9: Schematic representation of the total activity of the CME associated events. It depicts durations of SXR flares, Type III, II, IV activity and extrapolated CME launch time; all times are measured from the SXR flare peak.
Open with DEXTER

4 Summary and conclusions

The Artemis-IV catalogue of type II, IV and II/IV events recorded during 1998-2000 includes information on associated CMEs and flares; the examination of these data has provided us with some interesting statistical results.

All but one of the catalogued events were flare-associated; 68% were also CME-associated. Most (92%) of the CME take off times precede the associated SXR peak flux (GOES Channel 1-8 A) by 30 min on the average, while an important fraction (72%) also precedes the flare onset. Previous observations (Andrews 2001) maintain that the estimated CME launch precedes or follows, within 40 min, the SXR flare onset which does not contradict our statistic. Furthermore, our results are in accordance with a previous case study on four CME associated flares by Zhang et al. (2001), who stated that the initial acceleration phase of a CME coincides with the rise phase of the accompanying SXR flare, thus supporting a CME triggering of the flare scenario.

Of the type II bursts in our list, 67% are associated with CMEs. This corroborates the observations of Sheeley et al. (1984), Classen & Aurass (2002) who reported that 70% of type II bursts are associated with CMEs; they also pointed out that majority of the type II associated CMEs have velocities greater than 400 km s-1 which is the case with our data set; this is also in accordance with Gopalswamy et al. (2001) whose data point out that the faster CMEs drive type II busts, resulting in radio rich events.

Lastly, our results imply that type IV and composite type II/IV bursts manifest an even better association with CMEs (77% and 79% respectively) than type II bursts alone, thus providing a more reliable proxy for the CME onset.

Acknowledgements
This work was financially supported by the Research Committee of the University of Athens. The LASCO CME catalogue is generated and maintained by the Center for Solar Physics and Space Weather, The Catholic University of America in cooperation with the Naval Research Laboratory and NASA. SOHO is a project of international cooperation between ESA and NASA. The Solar Geophysical Data Catalogue is compiled and maintained by the US Department of Commerce. The authors wish to thank the anonymous referee for a number of extremely useful suggestions and comments on the first draft of this article.

References

  
Online Material

Appendix A: ARTEMIS-IV, LASCO, GOES and H$\alpha $ Observations

Table A.1

The table provides a concise summary of metric bursts (type II, IV and II/IV events), with their associated CMEs and SXR-H$\alpha $ flares, in the 1998-2000 period.

In Col. 1 the table gives the type of activity. Under the heading "associated bursts" we included all bursts of the type III family preceding type II, IV and II/IV events, as well as "arcs" and "precursors". In parenthesis, next to H$\alpha $, is the reporting station. Column 2 gives the date and start time of the events; for CMEs and type II bursts recorded by ARTEMIS-IV we have included the extrapolated launch time as specified in the "event list section", together with the time of the first detection; for the CMEs the time of first detection is the time of the first appearance in the field of view of the LASCO/C2 Coronagraph. In Col. 3 we give the peak time of SXR and H$\alpha $ flares and in Col. 4 the end times of bursts and flares.

Column 5 gives the flare class (B, C, M, X for the SXR; SF, SN, ..., 3B for H$\alpha $, from SGD) and Col. 6 the SXR integrated flux. In Col. 7 we give the location of the flare on the disk, in heliographic coordinates, as well as the central position angle and the width (entry underneath the position angle) of the CME. Column 8 shows the frequency range of radio bursts in MHz, while Col. 9 gives the logarithmic frequency drift of the type II bursts. Column 10 gives the drift velocities of type II bursts, calculated from the frequency drift rate and a Newkirk coronal density model assuming radial propagation (cf. Sect. 2.2) and the CME velocities from the LASCO lists, estimated by linear regression (Yashiro et al. 2001). In Col. 11 we give the CME acceleration from the same LASCO lists.

In the remarks column we provide information on the fine structure of type IV continua, from the SAO dynamic spectra; we also mark the type IIs which were outside the ARTEMIS-IV spectral range and the relevant information has been taken from the AIP, Izmiran or Nançay Decametric Array archives. In addition, we give information about the solar active region associated with each burst. For the radio bursts with imaging data from the Nançay radioheliograph, ensuring an "excellent" association with the accompanying H$\alpha $ flare, we record the NOAA number of the region of origin. Bursts without a NOAA region entry in this column, for lack of positional information, are only deemed in "good'' association with the accompanying activity.



Copyright ESO 2004