Acessibilidade / Reportar erro

Decimetric fine structures as a possible signature of chromospheric evaporation

Abstract

In this work, we concentrate in the analysis of radio slowly drifting fine structures associated with solar flares recorded in the frequency range of (1000 - 2500) MHz by the Brazilian Solar Spectroscope (BSS), in regular operation at National Institute of Space Research (INPE) - Brazil. The main morphological aspects of each fine structure are narrow-band of about 5 - 10 MHz and small duration of the order of 50 milliseconds. The majority of these fine structures are observed over a time interval of a couple of minutes before the maximum of the associated flare, in the impulsive phase. However, some observations during the gradual decay phase are also reported. They drift towards lower frequencies, with slow rates typically of about of 10 - 100 MHz s-1 Estimated velocities of the exciter from the negatively drifting structures suggest that the shock-like exciter is propagating in the higher chromosphere. Than, those fine structures are interpreted as a possible signature in decimetric emissions from the chromospheric evaporation phenomenon. Details of these observations and their interpretation in terms of the plasma emission produced by accelerated particles and the chromospheric evaporation front is presented.


Decimetric fine structures as a possible signature of chromospheric evaporation

F. C. R. Fernandes; J. R. Cecatto; M. C. Andrade; F. R. H. Madsen; L. C. P. Moraes; H. S. Sawant

Instituto Nacional de Pesquisas Espaciais, C.P. 515 12201-970, São José dos Campos, Brazil

ABSTRACT

In this work, we concentrate in the analysis of radio slowly drifting fine structures associated with solar flares recorded in the frequency range of (1000 - 2500) MHz by the Brazilian Solar Spectroscope (BSS), in regular operation at National Institute of Space Research (INPE) - Brazil. The main morphological aspects of each fine structure are narrow-band of about 5 - 10 MHz and small duration of the order of 50 milliseconds. The majority of these fine structures are observed over a time interval of a couple of minutes before the maximum of the associated flare, in the impulsive phase. However, some observations during the gradual decay phase are also reported. They drift towards lower frequencies, with slow rates typically of about of 10 - 100 MHz s-1 Estimated velocities of the exciter from the negatively drifting structures suggest that the shock-like exciter is propagating in the higher chromosphere. Than, those fine structures are interpreted as a possible signature in decimetric emissions from the chromospheric evaporation phenomenon. Details of these observations and their interpretation in terms of the plasma emission produced by accelerated particles and the chromospheric evaporation front is presented.

1 Introduction

Decimetric band is known to exhibit various types of fine structures [1-4]. However, observations in the decimetric band became stagnant till 80's. Around that time the Skylab satellite observations showed that energy to the flare is released where the electron density range is around 109 - 1010 cm-3, i.e., in the decimetric range of the solar radio spectrum [5 and references within]. Since then, various high frequency-time resolution spectroscopes were put into operation [6]. The Brazilian Solar Spectroscope (BSS) [7, 8] belongs to this new generation of solar radio spectroscopes. BSS is in regular operation since 1998 and it has recorded many different types of fine structures. Here, we analyze some particular groups of decimentric fine structures with slowly drifting high frequency edge. Karlický et al. [9] also reported the observations of slowly drifting structures in decimetric range.

Recently, the chromospheric evaporation has been invoked to explain the occurrence of the slowly drifting fine structures in radio wavelengths [10]. In fact, the first suggestion the chromospheric evaporation could be investigated by radio observations in decimetric wavelengths generated by electron beam going down in magnetic loop and crossing up-going evaporation front was given by [1].

The basic scenario of chromospheric evaporation is the following [12, 13]. During the very beginning of the flare, the turbulence in plasma of the active region begins and the flare energy is released up in the corona. This energy then accelerates particles and heats the plasma. The accelerated particles then move down toward the upper chromosphere along the magnetic lines, where they lose their energy by collisions with denser plasma. The heated plasma expands rapidly. However, it is confined by the magnetic loop structure, and then it flows upward with velocities of about a few hundreds of km s-1. The up-flow of the hot plasma is known as chromospheric evaporation [14]. The up-flowing hot plasma can be inferred by broadening and blue shifting of the lines as observed in soft X-ray [12, 13, 15-18].

Based especially on X-ray lines observations, many estimations of the velocity of chromospheric evaporation have been done. The determined evaporation front velocity ranges from aproximately 100 km s-1, for the transition region [19] or of about 150 km s-1 [19, 20] for coronal lines up to 950 km s-1 for decimetric layer [10]. Also, Mariska et al. [21] estimated the velocity of 650 km s-1 and Silva et al. [13] inferred the interval of 350-450 km s-1 for the velocity of the evaporation front.

While the hot and dense plasma rises up it creates a discontinuity in temperature and density, in the loop. Beam of energetic electrons moving towards the foot-points of the loop interacts with this slowly upward moving "shock front" where a optical thickness (t) is reduced due to its temperature dependence (t µ n × T-3/2) allowing to escape of high frequency radio emission. As the beam travels further it encounters high densities and lower temperature which increases opacity and radio emission is absorbed. Thus, as the shock front moves up, the cut off in high frequency will be slowly decreased allowing to infer the velocity of shock front. Eventually, the shock front will come in equilibrium with loop plasma and enabling escape of radio emission. Based on this scenario, Aschwanden and Benz [11] have suggested that the decimetric fine structures with slow drifting variation can be observed as result of the chomospheric evaporation process. However, there is a lack of observations in radio wavelengths as an evidence of emission caused by the cromospheric evaporation.

Here, we present the decimetric fine structures associated with the solar flares recorded by BSS as a possible signature in decimetric waves of the chromospheric evaporation process. In section 2 we briefly describe the BSS instrument and its main characteristics. In section 3 the observations and the main properties of the fine structures recorded are presented. The shock velocity determinations and results are presented in section 4 and 5, respectively. The summary and conclusions are presented in section 6.

2 Instrumentation

The BSS was put into regular operation at INPE for systematic solar observations in 1998 [8, 22]. It operates in conjunction with a 9 meter diameter polar mounted parabolic antenna, over the frequency range of (1000 - 2500) MHz, with high time (10 - 50 ms) and frequency (3 MHz) resolutions. The data are recorded in up to 100 digital channels [23], with absolute timing accuracy of 3 milliseconds. Thus, the BSS system has the capability to detect fine structures with narrow bandwidth and short durations [24].

3 Decimetric observations

Since the BSS has been put into regular operation, in 1998, it recorded more than 350 groups of solar decimetric emissions, most of them showing fine structures in time and/or frequency [25-28]. Out of them, we have selected the emissions associated with 8 solar flares recorded by BSS showing decimetric fine structures with high frequency edge slowly drifting from high to low frequencies observed in the period of March, 1999 to July, 2002. The Figure 1 shows the dynamic spectra of the recorded fine structures selected for analysis.


In the Table 1 we present a list of these emissions. The associated soft X-ray flares (GOES satellite) were also included. For 6 flares the fine structures occur prior the maximum of the soft X-ray associated emission, in the impulsive phase. In 3 cases, they are associated with the gradual decay phase of the X-ray emission, and for 1 case they are observed during the peak time of the soft X-ray emission. Except for 1 flare (classified as C1.6), the soft X-ray flares reported by GOES satellite associated with the fine structures were very intense: 3 flares reached X class and others 4 the M class.

4 Data analysis

4.1 Frequency Drift Rate

The high frequency and time measurements for each individual fine structure have been done using the BSSView and BSSData software. Then, the drift rate of the high frequency edge of each group of fine structures was determined from a linear regression in the frequency-time plane as shown in the Figure 2. The drift rates determined are given in the Table 2.


4.2 Chromospheric Evaporation Shock Velocity

The frequency drift rate (df/dt) of the bursts group can be given by

where f is the plasma frequency, Ne is the electron density, h is the height above the photosphere, q is the angle between the direction of propagation of the front and the vertical direction and vce is the chromospheric evaporation front velocity.

For the plasma emission the plasma frequency can be given as function of electron density in the solar atmosphere as:

where s is the harmonic radiation number (s = 1 and s = 2 for fundamental and second harmonic emission, respectively).

We have used a density model for the solar chromosphere derived from the statistical relation between drift rate and frequency of decimetric type III bursts [29]:

where A = 0.09 and a = 1.35. The density model used [29], adapted from Aschwanden and Benz [11], assumes a power-law dependence for the low corona and an exponential shape for the upper corona as:

where h is the height in the solar atmosphere, h1 is the height of the transition in the density variation regime, where the power-law model is changed into exponential one, N1 is the density transition, and Nq is the electron density at the bottom of quiet solar corona obtained by barometric model, p is the power-law index and H is the density height scale. From the continuity conditions at the height h1 we have:

and

As shown by Aschwanden and Benz [11], the constants p and Nq can be expressed as:

and

where v is the beam velocity expressed in units of the speed of light, A and a are constants in the empirical type III drift rate [29]. Considering the second harmonic radiations (s = 2) and the following values for the constants H = 7 × 104 km, h1 = 3.5 × 105 km, Nq = 4.6 × 108 cm-3 and p = 5, we determined the density and height intervals covered by the fine structure groups.

5 Results

Considering the frequency drift rates for the high frequency edge of the decimetric emissions, the shock-front velocity of the source was estimated assuming a displacement in a vertical direction. The measurements of total duration of the global structure (Dt), frequency drift rate (df/dt), density interval (DNe), height interval (Dh) and shock-front velocity of the fine structure source (vce) are shown in Table 2.

6 Summary and discussions

We have presented some radio fine structures showing slowly drifting frequency associated with solar flares recorded Brazilian Solar Spectroscope (BSS), in the frequency range of (1000-2500) MHz. The common characteristics of those fine structures are the narrow-band and short duration. The high frequency edge of the fine structure groups slowly drifts towards lower frequency, with rates less than 100 MHz s-1.

Those fine structures were interpreted as a possible signature of chromospheric evaporation process in decimetric wavelengths. In this case, the evaporation shock front causes the frequency drift of layer with reduced optical thickness, from where the decimetric radiation can escape. This drift depends upon the evaporation front velocity. Thus, there is a drifting "window" in which the electromagnetic radiation of plasma emission can escape and can be observed. The scenario of the visibility "window" for decimetric emission (1000-2500 MHz) drifting towards lower frequencies is valid at the very beginning of the flare, during the pre-impulsive phase in which the accelerated electron beam is heating the dense plasma. However, we have observed slowly drifting fine structures during the pre-flare phase but also associated with the gradual phase. These results also reported previously by [30] is an evidence that, even in the late phases of some particular flares, the dynamic of the heated plasma into the magnetic configuration permits the escape of the electromagnetic radiation in the decimetric range.

We considered the plasma emission at second harmonic, which has higher possibility to escape from dense layer of the solar atmosphere (109 cm-3) and we have used a density model for the solar chromosphere based on frequency drift rates of type III bursts observations [29]. From the frequency drift rate of the high frequency edge of the fine structures we are able to estimate the front velocity of the evaporation shock in the range of 30-1700 km s-1. These velocities are comparable with those computed by Mariska et al. [21] and Karlický [10] and other authors as mentioned in the Introduction, and they agreed with the upper limit of the chromospheric evaporation velocities reported by Fisher et al. [31].

These results permit to interpret those fine structures as a possible signature of the chromospheric evaporation phenomenon in decimetric wavelengths. Complementary investigations of those emissions in other wavelengths are ongoing to confirm the results and will be reported elsewhere.

Acknowledgements

The BSS Project is supported by INPE, CNPq and FAPESP. F.C.R. Fernandes thanks MCT-CNPq for receiving scholarship (PCI-3849/2002-06). L.C.P. Moraes thanks CNPq for IC scholarship (PIBIC).

Received on 6 February, 2004; revised version received on 26 April, 2004

  • [1] F. Droege, Astron. Astrophys. 57, 285 (1977).
  • [2] M. Guedel and A. O. Benz, Astron. Astrophys. Supp. Series 75, 243 (1988).
  • [3] M. A. F. Allaart, J. Nieuwkoop, C. Slottje and L. H. Sondaar, Solar Phys. 130, 183 (1990).
  • [4] H. Isliker and A. O. Benz, Astron. Astrophys. Supp. Series 104, 145 (1994).
  • [5] R. Moore et al., Report of NASA Skylab Workshop on Solar Flares, Colorado (1980).
  • [6] A. Kruger and W. Voight, Solar Phys. 161, 393 (1995).
  • [7] F. C. R. Fernandes, PhD Thesis, INPE-6396-TDI/612, INPE, Săo José dos Campos, Brazil (1997).
  • [8] H. S. Sawant, K. R. Subramanian, C. Faria, F. C. R. Fernandes, J. H. A. Sobral, J. R. Cecatto, R. R. Rosa, H. O. Vats, J. A. C. F. Neri, E. M. B. Alonso, F. P. V. Mesquita, V. A. Portezani and A. R. F. Martinon, Solar Phys. 200, 167 (2001).
  • [9] M. Karlický, F. Farnik and H. Meszarosova, Astron. Astrophys. 395, 677 (2002).
  • [10] M. Karlický, Astron. Astrophys. 338, 1084 (1998).
  • [11] M. J. Aschwanden and A. O. Benz, Astrophys. J. 438, 997 (1995).
  • [12] E. Antonucci, A. H. Gabriel and B. R. Dennis, Astrophys. J. 287, 917 (1984).
  • [13] A. V. R. Silva, H. Wang, D. E. Gary, N. Nitta and H. Zirin, Astrophys. J. 481, 978 (1997).
  • [14] P. A. Sturrock, Bulletin of the American Astronomical Society, 5, 28 (1973).
  • [15] E. Antonucci, B. R. Dennis, A. H. Gabriel and G. M. Simnett, Solar Phys. 96, 129 (1985).
  • [16] K. Tanaka, Astrophysics and Space Science 118, 101 (1986).
  • [17] T. Watanabe, Solar Phys. 126, 351 (1990).
  • [18] E. Antonucci, M. A. Dodero and R. Martin, Astrophys. J. Supp. Series 73, 147 (1990).
  • [19] L. Teriaca, A. Falchi, G. Cauzzi, R. Falciani, L. A. Smaldone and V. Andretta, Astrophys. J. 588, 596 (2003).
  • [20] G. Cauzzi, A. Falchi, L. Teriaca, R. Falciani, L. A. Smaldone and V. Andretta, S. Vito Conference Ser., 561 (2002).
  • [21] J. T. Mariska, A. G. Emslie and P. Li, Astrophys. J. 341, 1067 (1989).
  • [22] F. C. R. Fernandes, J. R. Cecatto, J. A. C. F. Neri, C. Faria, A. R. F. Martinon, R. R. Rosa,F. P. V. Mesquita, V. A. Portezani, M. C. Andrade, E. M. B. Alonso, H. O. Vats and H. S. Sawant, Boletim da SAB 20, 33 (2000).
  • [23] H. S. Sawant, K. R. Subramanian, C. Faria, S. Stephany, F. C. R. Fernandes, J. R. Cecatto, R. R. Rosa, V. A. Portezani, F. P. V. Mesquita and E. M. B. Alonso, Astron. Soc. Pacific Conf. Ser. 206, 347 (2000).
  • [24] F. C. R. Fernandes, J. R.; Cecatto, J. A. C. F. Neri, M. J. B. Silva,R. R. Rosa, C. Faria, M. C. Andrade, E. M. A. Alonso, A. R. F. Martinon, K. M. Dobrowolski, T. M. Barbosa, D. C. Freitas, V. F. Nascimento, A. L. Reis and H. S. Sawant, Boletim da SAB 21, 49 (2001).
  • [25] F. C. R. Fernandes, Scientific Report INPE-9654-RPQ/740. INPE, Săo José dos Campos (2003a).
  • [26] F. C. R. Fernandes, Scientific Report INPE-9653-RPQ/739. INPE, Săo José dos Campos (2003b).
  • [27] F. C. R. Fernandes, Scientific Report INPE-9652-RPQ/738. INPE, Săo José dos Campos (2003c).
  • [28] F. C. R. Fernandes, Scientific Report INPE-9881-RPQ/742. INPE, Săo José dos Campos (2003d).
  • [29] J. L. Melendez, H. S. Sawant, F. C. R. Fernandes and A. O. Benz, A. O., Solar Phys. 187, 77 (1999).
  • [30] A. Czaykowska, B. De Pontieu, D. Alexander and G. Rank, Astrophys. J. 521, L75 (1999).
  • [31] G. Fisher, R. C. Canfield and A. N. McClymont, Astrophys. J. 281, L79 (1984).

Publication Dates

  • Publication in this collection
    01 Mar 2005
  • Date of issue
    Dec 2004

History

  • Accepted
    26 Apr 2004
  • Received
    06 Feb 2004
Sociedade Brasileira de Física Caixa Postal 66328, 05315-970 São Paulo SP - Brazil, Tel.: +55 11 3091-6922, Fax: (55 11) 3816-2063 - São Paulo - SP - Brazil
E-mail: sbfisica@sbfisica.org.br