Free Access
Issue
A&A
Volume 533, September 2011
Article Number L2
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201117605
Published online 11 August 2011

© ESO, 2011

1. Introduction

Solar flare emission is now observed over virtually the whole electromagnetic spectrum from radio frequencies as low as 10 − 100 MHz up to photon energies as high as a few hundred MeV. The bulk of this emission originates from the chromosphere, a small part of the solar atmosphere that is narrow in height and much smaller than the size of a flaring region. Being only about  < 3 Mm ( < 4′′) across, the study of the chromosphere is a formidable observational challenge in terms of both angular and temporal resolution especially for transient phenomena such as solar flares. During solar flares, large numbers of energetic electrons are accelerated in the solar corona and effectively release their energy into the dense chromosphere. These energetic particles are believed to be responsible (although not always directly) for the HXR, EUV, continuum “white-light” (WL), infrared, and radio emissions. On the basis of temporal correlation, the interconnection between these emissions was realized at an early stage (Najita & Orrall 1970; Švestka 1970). Owing to the lack of height-resolving observations, spectroscopic data have mostly been used to identify the structure of the emitting region (e.g. Canfield 1974). Therefore, the actual mechanisms of WL emission and their relation to HXR emitting electrons remains unclear (Neidig 1989; Neidig et al. 1993; Sylwester & Sylwester 2000; Matthews et al. 2003; Potts et al. 2010; Watanabe et al. 2010; Kretzschmar 2011). It is not even clear whether the emission is optically thick or thin or both. Thus the characteristic heights of continuum WL emission in various models are between the photosphere and the upper chromosphere. Observations supporting both (Fang & Ding 1995; Xu et al. 2006) a photospheric origin (e.g. Boyer et al. 1985; Ding et al. 1999; Chen & Ding 2006) presumably due to radiative back-warming, and the generation of WL emission in the higher temperature regions of the upper chromosphere (e.g. Machado & Rust 1974; Hudson 1972) have been made. The major observational challenge here is to have simultaneous HXR, UV, and WL images with sufficiently high angular resolution. The relative horizontal positions of HXR and WL flares have been investigated in detail and, although they have generally been found to agree (e.g. Krucker et al. 2011; Fletcher et al. 2007; Hudson et al. 2006; Metcalf et al. 2003), the uncertainty in the relative positions was large because of the limited absolute pointing accuracy of partial solar disk observations such as those from TRACE. Thus it was impossible to quantify the relative height of the sources above the photosphere in previous studies. With RHESSI (Lin et al. 2002) data, the height structure of HXR sources has become accessible for detailed observational studies (Aschwanden et al. 2002; Kontar et al. 2008; Kontar & Jeffrey 2010; Battaglia & Kontar 2011; Saint-Hilaire et al. 2010). These observations have demonstrated that the higher energy HXR footpoints originate from progressively lower heights in the chromosphere, in good agreement with collisional downward propagation of energetic electrons (see Brown et al. 2002, for details). The newly available SDO continuum images from HMI (Wachter et al. 2011) and AIA EUV (Lemen et al. 2011) observations allow us to scrutinize the spatial structure and more importantly the height structure of various types of emission associated with HXR footpoints. In this paper, we present the first spatially resolved simultaneous observations of HXR, EUV, and WL emission from a footpoint of a solar limb flare observed with RHESSI, SDO/AIA, and HMI. Using the recently developed hard X-ray visibility analysis technique (Hurfordet al. 2002; Schmahl et al. 2007), we determine the characteristic heights of the different types of emission using the method presented by (Battaglia & Kontar 2011; Kontar et al. 2010, 2008). The results are consistent with the EUV and continuum WL being produced above the HXR footpoints by lower energy electrons (~12 keV), hence suggest not a photospheric but an upper chromospheric origin of WL emission.

2. Height and density measurements

thumbnail Fig. 1

Top: RHESSI lightcurves of the main flare phase in different energy bands: 6–12 keV (green), 12–25 keV (red), 25–50 keV (blue), 50–100 keV (purple). Bottom: normalized AIA and HMI emission emanating from location of the southern HXR footpoint as function of time at several wavelengths (see legend).

The GOES M3.5 limb flare occurred on 2011 February 24 with three well-pronounced main HXR (above  ~30 keV) peaks between 07:29 and 07:33 UT and an associated filament eruption to the south of the flaring site. Figure 1 (top) shows the RHESSI lightcurves in different energy bands of the thermal and non-thermal emission. The attenuator state was 1 during the main phase of the flare and the live-time was higher than 90% during the whole period of the observations, hence making the effects of pulse pile-up negligible. The flare appeared in the active region AR11163, was well observed in enhanced EUV emission visible in all AIA filters, and also displayed 6173 Å continuum emission observed by HMI. The continuum emission (with a maximum enhancement of  ~15% above the photospheric background) delineated two footpoints and followed the X-ray time profile. Unfortunately, AIA images became saturated at most wavelengths before the soft X-ray peak (2nd HXR peak) of the flare. Because of this saturation, we focus the analysis on the time of the first HXR peak. During the time around 07:30:10 UT, AIA was not saturated in the 94 Å, 193 Å, 304 Å, and 335 Å filters near the location of the southern footpoint.

thumbnail Fig. 2

Top: RHESSI soft (red) and hard X-ray emission (blue) superimposed on AIA 171 Å image. Bottom: AIA difference images at three wavelengths and HMI image (bottom right). Contours at 50, 70 and 90% of the maximum emission in RHESSI CLEAN images between 07:30:00 and 07:30:40 UT at 6–12 keV (red) and 25–50 keV (blue) are overlaid. The temporal evolution is shown in the movies available in the on-line edition. RHESSI light-curves (cf. top panel of Fig. 1) were added to the animation of the bottom panel to provide a link between the total X-ray flux evolution and the corresponding images.

To study the flare related enhancement in the EUV and continuum emission, the image from one minute before the analyzed time was subtracted. Figure 2 displays difference images in AIA 94, 193 and 304 Å wavelengths and the HMI continuum1. For the RHESSI analysis, a 40 s time interval at the first peak of the flare from 07:30:00 UT was chosen (see Fig. 1). The contours in Fig. 2 indicate the RHESSI emission between 07:30:00 and 07:30:40 UT at 50, 70, and 90% of the maximum in CLEAN images (Hurford et al. 2002). The images show a coronal (SXR) source above the limb and two footpoint (HXR) sources (north and south). During the first peak, the southern footpoint source was much more intense than the northern footpoint. Thus we focused the height analysis on the southern footpoint. We note that in the subsequent peaks the northern footpoint reached the same intensity as the southern footpoint. The RHESSI spectrum suggests emission well beyond 200 keV into the γ-ray range. During the analyzed time interval, the HXR spectrum is well fitted with a thermal component and a single thick-target power law with spectral index δ = 3.9. The flare was also observed by Fermi/GBM (Meegan et al. 2009) and the corresponding spectrum from FERMI data was found to be virtually identical to the RHESSI spectrum. The time evolution of the EUV emission at the position of the southern HXR footpoint can be seen in Fig. 1 (bottom). The total flux in a region of 5 by 10 arcsec extent around the position of the southern HXR footpoint is shown.

Because of the limb location of the flare a height analysis is possible. However, SDO observations alone cannot infer the heights of the sources at different EUV wavelengths or for white-light emission. Here, utilizing RHESSI observations, we find the density structure of the chromosphere and importantly the reference height (projected radius of the point at the photospheric level beneath the HXR sources) as described by Battaglia & Kontar (2011) and Kontar et al. (2010, 2008). Following the method described by those authors, we use visibility forward fitting to find the position of the southern footpoint as a function of HXR energy in several energy bands (30–40 keV, 40–55 keV, 55–80 keV, 80–130 keV). The energy bands were chosen to be large enough to guarantee good count statistics for reliable fits. Both the southern and the northern footpoints are fitted with circular Gaussians. This assures that the emission of the northern footpoint is properly accounted for, even if it is not intense enough to permit fit parameters to be determined with acceptable errors. The flare morphology was such that the radial direction (direction of electron propagation along the loop) can be seen along the x-direction in the images (see Fig. 2). An exponential density profile (hydrostatic atmosphere at lower heights), n(h) = ncs + nphexp( − (h − rref)/h0), was fitted to the positions found by visibility forward fitting, where ncs = 4.6 × 1010 cm-3 is the (constant) loop density, nph = 1.16 × 1017 cm-3 the photospheric density (Vernazza et al. 1981), which corresponds to τ500   nm = 1, and h is the height above the photosphere2. The loop density ncs was assumed to be the same as the coronal source density determined from the emission measure given by the spectral fit and the size of the coronal source in RHESSI 6–12 keV images. Fitting the HXR energy-position relation with the above density model, we find a scale height of h0 = 211 ± 46 km and a reference height for the radial distance between the solar disk center and the point under the footpoints at the photosphere of rref = 929.4 ± 0.3 arcsec. The photospheric reference distance rref is used to find the absolute height of the EUV and WL sources above the photosphere. Using the HXR source size as the size of the magnetic flux tube, we selected a stripe along the x-direction of 5 arcsec width in the y-direction over the position of the southern RHESSI footpoint. At each horizontal distance x, the AIA and HMI counts were summed in the y-direction. The height of the emission is simply the x-position minus rref as found above. This is shown in Fig. 3. Proper co-alignment of the images is crucial for this study. The RHESSI disk center is known to an accuracy higher than 0.2 arcsec. Using SDO full disk images, the disk center of AIA and HMI can be determined. This position of solar-disk center is found to deviate by only –0.39/0.67 arcsec in the x/y direction from that of RHESSI, a deviation that, although small, was taken into account. The larger uncertainty is related to the roll-angle of the SDO images (rotation around the disk center). To investigate the effect, the HMI image was rotated by 0.1 degrees relative to the Sun center. This uncertainty translates mostly into a y-direction uncertainty of  <1.6 arcsec and  <0.5 for the x-direction for the limb event under study. For the height measurements of the EUV and HMI emissions relative to the RHESSI emission, we have taken this uncertainty into account as shown in Fig. 3. We find that the RHESSI sources at 35–100 keV originate from heights of between 1.7 and 0.8 Mm. The WL emission peaks at 2.5 Mm, which is  ~1 Mm above the 30 keV HXR emission, suggesting it originates from the upper chromosphere.

thumbnail Fig. 3

Top: height of HXR emission at the lowest (red vertical lines) and highest (blue vertical lines) imaged energies, EUV emission at different wavelengths (see legend), and WL relative enhancement where the black line indicates the height of maximum enhancement. Arrows indicate the corresponding uncertainties. Middle: energy deposition rate as a function of height for two different cases of electron cutoff energy. Bottom: fitted density model (thick solid line). The thin solid lines give the lower and higher limits of the density determined by the uncertainty in the fit of the scale height. The model from Vernazza et al. (1981) is given for comparison.

Using the density model found above and a power-law electron injection rate (electrons s-1 keV-1) found from the spatially integrated RHESSI X-ray spectrum, we can compute the energy deposition rate (ergs s-1 Mm-1) as a function of height (Brown 1973), (1)where F(E,h) is the electron rate spectrum (Brown 1971), K = 2πe4 where ΛeH and Λee are the Coulomb logarithms for electron-neutral and electron-electron collisions, and nH(h) and ne(h) are the hydrogen and electron number densities respectively. The resulting energy deposition rate as a function of height given by Eq. (1) is shown in Fig. 3 where we have assumed a step change in ionization from fully ionized to neutral at a height of 2 Mm. As expected, the deposition rate depends on the low energy cut-off of the injected electron spectrum F0(E0). From the thick target fit to the RHESSI spatially integrated HXR spectrum the cut-off energy is around Ec = 12 keV and the spectral index δ = 3.9. This leads to a total power of  ~2.7 × 1027erg   s-1. However, the cut-off is generally unknown as the low energy part of the X-ray spectrum is dominated by thermal emission, so lower cut-off values are possible and likely (e.g. Emslie 2003; Hannah et al. 2009). However, even with Ec = 12 keV, the energy deposited appears much larger than that emitted at other wavelengths. For the southern footpoint, an enhancement of  ~5% above the total solar irradiance of 3.4 × 1026   erg   s-1   arcsec-2 over a source area of  ~20 arcsec2 leads to a total power emitted of 3.4 × 1026   erg   s-1. The energy in the non-thermal electrons is therefore up to an order of magnitude larger than that needed to power the WL emission, which often constitutes the bulk of flare emitted energy (Neidig 1989). We note that the WL energetics obtained using the 6173 Å continuum is found to be lower than at TRACE wavelengths near the 1700 Å continuum, and assuming black body emission to deal with the TRACE response (e.g. Fletcher et al. 2007), but close to the values reported by Neidig (1989) based on spectroscopy above 2500 Å. Various optical bands indeed show different levels of enhancement during flares (Kretzschmar 2011). In addition, the maximum of 6173 Å continuum is reached later in the flare, after the analyzed time interval.

3. Discussion and conclusions

Our observations suggest that the WL emission originates from either the upper chromosphere or lower transition region with a density of between 1011   cm-3 and 1013   cm-3. These upper and lower limits of the density are due to the uncertainties in the fit of the density scale height. Figure 3 (bottom) displays the fitted density model, as well as the upper and lower limits to the density. The arrows in Fig. 3 (top) illustrate the uncertainty in the RHESSI source heights resulting from the uncertainty in the density scale height. The peak of WL emission is found to be about 1 Mm higher than the position of the lowest observed RHESSI energy range (30–40 keV). Despite all the uncertainties in the height caused by the density fit or the relative pointing error, the WL source is still above the 30–40 keV RHESSI sources. Owing to the presence of the SXR coronal emission in the X-ray spectrum, the non-thermal HXR spectrum below  ≈ 20 keV is not known. However, extrapolating the height function found in Sect. 2 to lower energies, a source at 12 keV would be seen at about 2 Mm. This coincides with the height of the maximum energy deposition found from Eq. (1) for cut-off energies of 12 keV (Fig. 2). The observed position of the maximum of the WL emission would then correspond to energy deposition by the low energy part of the non-thermal electron distribution, in contrast to the assumption often made that the cut-off in electron energies is of the order of 20 keV. This also suggests that WL and HXR sources above 30 keV should be spatially separated by  <1 Mm with decreasing separation in events towards the disk center. However, we note that the height of the maximum energy deposition, especially in the case of low initial cut-off energies, is influenced by the loop density and the length of the loop. Assuming that the acceleration of the particles happens at the center of what is observed as the coronal source, Fig. 2 indicates a precipitation distance of about 13 arcsec from the site of acceleration to the footpoints. This value was used for the computation of the energy deposition rate. In this case, the maximum energy deposition of electrons with energy  ≲ 12 keV will be near the top of the loop. Should the loop density be lower or the acceleration region extended, the lower energy electrons would deposit their energy at greater depths. The other interesting aspect of this flare is that HXR emission above 20 keV is well fitted by a single power-law. Non-uniform ionization of plasma produces a break at the energy corresponding to the column depth of the transition region, which is observed in other flares (Kontar et al. 2002; Su et al. 2011). The absence of this break point to the transition region above the stopping depth of 20 keV electrons is again consistent with the height measurements. The observation of strong emission from heights  ≥ 2 Mm suggests that a substantial part of the WL continuum is formed in an optically thin or finite optically thick region with ionized or partially ionized plasma in the upper chromosphere favoring free-bound and free-free emission of the 6173   Å continuum. This is also supported by the height of the AIA emission at 94, 335, and 193 Å (representing temperatures  >106 K) originating from around 3 Mm and above the WL emission peak. Finally, it is likely that the WL is powered by lower energy electrons  ~12 keV, which do not penetrate deeply into the solar chromosphere.


1

Animated images of the flare are also available at http://www.astro.gla.ac.uk/users/mbattaglia/20110224_online_material/

2

In this work, we use a fit to the full expression for the HXR flux maximum and not the simplified version given by Eq. (7) of Brown et al. (2002). The latter is likely to underestimate the densities in a sharply changing chromosphere. See also Battaglia & Kontar (2011) and Kontar et al. (2010, 2008).

Acknowledgments

The authors thank the anonymous referee for helpful comments and suggestions. This work is supported by the Leverhulme Trust (M.B., E.P.K.) and a STFC rolling grant (E.P.K.). Financial support by the European Commission through the FP7 HESPE network is gratefully acknowledged.

References

  1. Aschwanden, M. J., Brown, J. C., & Kontar, E. P. 2002, Sol. Phys., 210, 383 [NASA ADS] [CrossRef] [Google Scholar]
  2. Battaglia, M., & Kontar, E. P. 2011, ApJ, 735, 42 [NASA ADS] [CrossRef] [Google Scholar]
  3. Boyer, R., Sotirovsky, P., Machado, M. E., & Rust, D. M. 1985, Sol. Phys., 98, 255 [NASA ADS] [CrossRef] [Google Scholar]
  4. Brown, J. C. 1971, Sol. Phys., 18, 489 [NASA ADS] [CrossRef] [Google Scholar]
  5. Brown, J. C. 1973, Sol. Phys., 31, 143 [NASA ADS] [CrossRef] [Google Scholar]
  6. Brown, J. C., Aschwanden, M. J., & Kontar, E. P. 2002, Sol. Phys., 210, 373 [NASA ADS] [CrossRef] [Google Scholar]
  7. Canfield, R. C. 1974, Sol. Phys., 34, 339 [NASA ADS] [CrossRef] [Google Scholar]
  8. Chen, Q. R., & Ding, M. D. 2006, ApJ, 641, 1217 [NASA ADS] [CrossRef] [Google Scholar]
  9. Ding, M. D., Fang, C., & Yun, H. S. 1999, ApJ, 512, 454 [NASA ADS] [CrossRef] [Google Scholar]
  10. Emslie, A. G. 2003, ApJ, 595, L119 [NASA ADS] [CrossRef] [Google Scholar]
  11. Fang, C., & Ding, M. D. 1995, A&AS, 110, 99 [NASA ADS] [Google Scholar]
  12. Fletcher, L., Hannah, I. G., Hudson, H. S., & Metcalf, T. R. 2007, ApJ, 656, 1187 [NASA ADS] [CrossRef] [Google Scholar]
  13. Hannah, I. G., Kontar, E. P., & Sirenko, O. K. 2009, ApJ, 707, L45 [NASA ADS] [CrossRef] [Google Scholar]
  14. Hudson, H. S. 1972, Sol. Phys., 24, 414 [NASA ADS] [CrossRef] [Google Scholar]
  15. Hudson, H. S., Wolfson, C. J., & Metcalf, T. R. 2006, Sol. Phys., 234, 79 [NASA ADS] [CrossRef] [Google Scholar]
  16. Hurford, G. J., Schmahl, E. J., Schwartz, R. A., et al. 2002, Sol. Phys., 210, 61 [NASA ADS] [CrossRef] [Google Scholar]
  17. Kontar, E. P., & Jeffrey, N. L. S. 2010, A&A, 513, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Kontar, E. P., Brown, J. C., & McArthur, G. K. 2002, Sol. Phys., 210, 419 [NASA ADS] [CrossRef] [Google Scholar]
  19. Kontar, E. P., Hannah, I. G., & MacKinnon, A. L. 2008, A&A, 489, L57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Kontar, E. P., Hannah, I. G., Jeffrey, N. L. S., & Battaglia, M. 2010, ApJ, 717, 250 [NASA ADS] [CrossRef] [Google Scholar]
  21. Kretzschmar, M. 2011, A&A, 530, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Krucker, S., Hudson, H. S., Jeffrey, N. L. S., et al. 2011, ApJ, in press [Google Scholar]
  23. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2011, Sol. Phys., in press [Google Scholar]
  24. Lin, R. P., Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Phys., 210, 3 [Google Scholar]
  25. Machado, M. E., & Rust, D. M. 1974, Sol. Phys., 38, 499 [NASA ADS] [CrossRef] [Google Scholar]
  26. Matthews, S. A., van Driel-Gesztelyi, L., Hudson, H. S., & Nitta, N. V. 2003, A&A, 409, 1107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Meegan, C., Lichti, G., Bhat, P. N., et al. 2009, ApJ, 702, 791 [NASA ADS] [CrossRef] [Google Scholar]
  28. Metcalf, T. R., Alexander, D., Hudson, H. S., & Longcope, D. W. 2003, ApJ, 595, 483 [NASA ADS] [CrossRef] [Google Scholar]
  29. Najita, K., & Orrall, F. Q. 1970, Sol. Phys., 15, 176 [NASA ADS] [CrossRef] [Google Scholar]
  30. Neidig, D. F. 1989, Sol. Phys., 121, 261 [NASA ADS] [CrossRef] [Google Scholar]
  31. Neidig, D. F., Wiborg, P. H., & Gilliam, L. B. 1993, Sol. Phys., 144, 169 [NASA ADS] [CrossRef] [Google Scholar]
  32. Potts, H., Hudson, H., Fletcher, L., & Diver, D. 2010, ApJ, 722, 1514 [NASA ADS] [CrossRef] [Google Scholar]
  33. Saint-Hilaire, P., Krucker, S., & Lin, R. P. 2010, ApJ, 721, 1933 [NASA ADS] [CrossRef] [Google Scholar]
  34. Schmahl, E. J., Pernak, R. L., Hurford, G. J., Lee, J., & Bong, S. 2007, Sol. Phys., 240, 241 [NASA ADS] [CrossRef] [Google Scholar]
  35. Su, Y., Holman, G. D., & Dennis, B. R. 2011, ApJ, 731, 106 [NASA ADS] [CrossRef] [Google Scholar]
  36. Sylwester, B., & Sylwester, J. 2000, Sol. Phys., 194, 305 [NASA ADS] [CrossRef] [Google Scholar]
  37. Švestka, Z. 1970, Sol. Phys., 13, 471 [Google Scholar]
  38. Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 [NASA ADS] [CrossRef] [Google Scholar]
  39. Wachter, R., Schou, J., Rabello-Soares, M. C., et al. 2011, Sol. Phys., 19 [Google Scholar]
  40. Watanabe, K., Krucker, S., Hudson, H., et al. 2010, ApJ, 715, 651 [NASA ADS] [CrossRef] [Google Scholar]
  41. Xu, Y., Cao, W., Liu, C., et al. 2006, ApJ, 641, 1210 [NASA ADS] [CrossRef] [Google Scholar]

Online material

Movie of Fig. 2 (top) Access here

Movie of Fig. 2 (bottom) Access here

All Figures

thumbnail Fig. 1

Top: RHESSI lightcurves of the main flare phase in different energy bands: 6–12 keV (green), 12–25 keV (red), 25–50 keV (blue), 50–100 keV (purple). Bottom: normalized AIA and HMI emission emanating from location of the southern HXR footpoint as function of time at several wavelengths (see legend).

In the text
thumbnail Fig. 2

Top: RHESSI soft (red) and hard X-ray emission (blue) superimposed on AIA 171 Å image. Bottom: AIA difference images at three wavelengths and HMI image (bottom right). Contours at 50, 70 and 90% of the maximum emission in RHESSI CLEAN images between 07:30:00 and 07:30:40 UT at 6–12 keV (red) and 25–50 keV (blue) are overlaid. The temporal evolution is shown in the movies available in the on-line edition. RHESSI light-curves (cf. top panel of Fig. 1) were added to the animation of the bottom panel to provide a link between the total X-ray flux evolution and the corresponding images.

In the text
thumbnail Fig. 3

Top: height of HXR emission at the lowest (red vertical lines) and highest (blue vertical lines) imaged energies, EUV emission at different wavelengths (see legend), and WL relative enhancement where the black line indicates the height of maximum enhancement. Arrows indicate the corresponding uncertainties. Middle: energy deposition rate as a function of height for two different cases of electron cutoff energy. Bottom: fitted density model (thick solid line). The thin solid lines give the lower and higher limits of the density determined by the uncertainty in the fit of the scale height. The model from Vernazza et al. (1981) is given for comparison.

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.