A&A 455, 1123-1133 (2006)
L. Teriaca1 - A. Falchi2 - R. Falciani3 - G. Cauzzi2 - L. Maltagliati1,2
1 - Max-Planck-Institut für Sonnensystemforschung, Max-Plank Str. 2, 37191 Katlenburg-Lindau, Germany
2 - INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy
3 - Dipartimento di Astronomia e Scienza dello Spazio, Università di Firenze, Largo Fermi 2, 50125 Firenze, Italy
Received 22 February 2006 / Accepted 23 May 2006
Aims. We study the dynamics and the evolution of a C2.3 two-ribbon flare, developed on 2002 August 11, during the impulsive phase as well as during the long gradual phase. To this end we obtained multiwavelength observations using the CDS spectrometer aboard SOHO, facilities at the National Solar Observatory/Sacramento Peak, and the TRACE and RHESSI spacecrafts.
Methods. CDS spectroheliograms in the Fe XIX, Fe XVI, O V and He I lines allow us to determine the velocity field at different heights/temperatures during the flare and to compare them with the chromospheric velocity fields deduced from H image differences. TRACE images in the 17.1 nm band greatly help in determining the morphology and the evolution of the flaring structures.
Results. During the impulsive phase a strong blue-shifted Fe XIX component (-200 km s-1) is observed at the footpoints of the flaring loop system, together with a red-shifted emission of O V and He I lines (20 km s-1). In one footpoint simultaneous H data are also available and we find, at the same time and location, downflows with an inferred velocity between 4 and 10 km s-1. We also verify that the "instantaneous'' momenta of the oppositely directed flows detected in Fe XIX and H are equal within one order of magnitude. These signatures are in general agreement with the scenario of explosive chromospheric evaporation. Combining RHESSI and CDS data after the coronal upflows have ceased, we prove that, independently from the filling factor, an essential contribution to the density of the post-flare loop system is supplied from evaporated chromospheric material. Finally, we consider the cooling of this loop system, that becomes successively visible in progressively colder signatures during the gradual phase. We show that the observed cooling behaviour can be obtained assuming a coronal filling factor of 0.2 to 0.5.
Key words: Sun: flares - Sun: chromosphere - Sun: corona - line: profiles
Most solar flare models are based on the idea that flares arise from the sudden release of free magnetic energy stored in the corona in non-potential configurations (e.g., Klimchuk et al. 1988). The bulk of the energy is released through reconnection of the magnetic field lines in the corona (Shibata 1996; Kopp & Pneuman 1976; Petschek 1964) and then transported down to the chromosphere along the magnetic field lines by accelerated particles or by a thermal conduction front. Hydrodynamic simulations indicate that both non-thermal electrons (Mariska et al. 1989; Fisher et al. 1985a,b) and a thermal conduction front (Gan et al. 1991) should lead to the evaporation of chromospheric material. Models by Fisher et al. (1985a) result in upflows (gentle evaporation) or downflows (explosive evaporation), on both transition-region (TR) and chromosphere, depending on the energy flux of the impinging electrons, with downflows predicted for large energy fluxes. In the case of explosive evaporation the overpressure of evaporated plasma drives upward motions in the corona and downward motions in the chromosphere and it is expected that the momenta of the oppositely moving plasma should be balanced. The chromospheric evaporation injects mass into the corona to fill coronal loops that appear as the soft X-ray (SXR) flare and that eventually cool down radiatively and/or conductively forming cold H post-flare loops (see, e.g., Moore et al. 1980). Although this general picture seems well established, the details of some of the key processes assumed to take place remain to be verified. The small temporal and spatial scales involved, together with the need for comprehensive temporal and height coverage, conjure up to make "complete'' flare observations a very demanding task.
Blue-shifts of lines formed at flare temperatures during the impulsive phase (corresponding to upflow velocities in excess of several hundreds kilometre per second) have been observed since the early 1980s (e.g., Antonucci et al. 1982; Doschek et al. 1980) and confirmed with the Bragg Crystal Spectrometer aboard Yohkoh (BCS, Culhane et al. 1991). However, the majority of these observations revealed only a blue asymmetry, indicating the presence of a large, static component even in the early phases of the flare (e.g., Gan & Watanabe 1997; Bentley et al. 1994), rather than wholly blue-shifted profiles as predicted by hydrodynamic models (e.g., Gan et al. 1995; Emslie et al. 1992). To which extent the lack of spatial resolution might have influenced these observations is not easily quantifiable (neither BCS nor earlier spectrometers provided any effective spatial resolution on the solar surface). Since reconnection is expected to occur at very small spatial scales, it might well be that the associated chromospheric evaporation takes place in rather confined portions of the solar surface. This is indeed the case in eruptive flares, for which, at any given time, the evaporation in the reconnection scenario occurs in a very thin shell (1'') located on the outer part of the flaring loop system (Falchi et al. 1997; Forbes & Acton 1996), and the heating of individual field lines is supposed to last only 100 s or so (Forbes 2003). The smearing introduced by observing over the whole disk would then significantly alter the relevant signals. Hydrodynamic simulations modelling solar flare emission with multiple loops (Hori et al. 1998; Warren & Doschek 2005) have shown that it is possible to reproduce the line profiles dominated by the stationary component (as observed by BCS) with a simple model based on the successive independent heating of small scale threads. Only the simulated line profiles for individual threads do show strongly blue shifted line profiles.
Spatially and temporally resolved observations of flare lines have finally become available with the advent of the Coronal Diagnostic Spectrometer (CDS, Harrison et al. 1995) aboard SOHO, which provides stigmatic slit images in the 31 to 38 nm (NIS-1) and 51 to 63 nm (NIS-2) range with a spatial resolution of few seconds of arc and a cadence depending upon the adopted observing program. Various authors have since reported the presence of blueshifted emission in hot coronal lines during flares, signalling upflows of up to 200 km s-1 (Czaykowska et al. 1999; Brosius & Phillips 2004; Brosius 2003; Del Zanna et al. 2002; Teriaca et al. 2003), but the total number of events analysed up to now remains small. Moreover, the same authors report about intricate patterns of blue- and red-shifted emission in TR lines, depending both on the location within the flaring region and the phase of the flare, that still need to be consistently explained within the theory of chromospheric evaporation. Finally, spatially and temporally resolved observations of flare lines need to be combined with simultaneous chromospheric observations, in order to verify the spatial and temporal coincidence of chromospheric downflows with hot plasma upflows. In fact, while many authors have reported on the simultaneous appearance of redshifts in chromospheric lines and blue-shifts in flare lines during the early phase of flares (e.g., Wülser et al. 1994; Canfield et al. 1990b; Falchi et al. 1992), to date the only direct proof that such shifts arise from spatially coincident areas has been provided by Teriaca et al. (2003). However, due to the observational setup, their coronal and chromospheric data were not simultaneous, leading to the need of further observations.
In this context, we report here on comprehensive, multi-instrument observations of a small eruptive flare, covering its temporal evolution from the impulsive phase and related evaporation throughout the gradual phase and appearance of cold postflare loops. Spatially resolved CDS spectra allow us to study the dynamics of the upper atmosphere at various times during the flare, while simultaneous chromospheric observations provide a test for the chromospheric evaporation scenario. RHESSI data acquired during the gradual phase, combined with CDS measurements, provide a final element in support of the chromospheric evaporation. Other SOHO instruments and TRACE, provide various context information and help us in defining and constraining the flare evolution.
|Figure 1: TRACE image (logarithmic scale) in the 17.1 nm band showing AR 10061 at flare peak. The dark grey (red in the electronic version) box indicates the UBF FOV while the area rastered by CDS/NIS is represented by the light grey (blue) box. The white vertical line marks the position of the DST spectrograph slit. Here and in the following images, the white dashed line shows the apparent magnetic inversion line obtained from MDI data in the region around the filament.|
Space-born data - Spectroheliograms of a area (mostly overlapping the UBF FOV, see Fig. 1) were obtained with the Normal Incidence Spectrometer (NIS) of the CDS experiment (Harrison et al. 1995) starting at 14:26:27 UTC. Each raster was obtained in 210 s by stepping the 4'' wide slit eastward in 19 steps of 6''. On board binning over two rows yielded a pixel size of 3.4'' along the slit. Spectra were acquired in the He I 58.43 nm ( K), O V 62.97 nm ( K), Fe XVI 36.08 nm ( K) and Fe XIX 59.22 nm ( K) lines. Every ten rasters an additional step of westward was performed to compensate for the solar rotation. During the entire period of CDS observations AR 10061 was constantly monitored by the Transition Region and Coronal Explorer (TRACE, Handy et al. 1999) that provided images in the 17.1 nm band (Fe IX-X, K) with a spatial scale of 1'' pixel-1 (on board binning) and a cadence of s. The Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI, Lin et al. 2002) exited its night a few minutes after the peak phase of the flare, and provided information on hot ( K) thermal flare plasma as well as on non-thermal components. From RHESSI data in the energy band 3 to 10 keV, we derived images with spatial resolution up to 2'' and spectra of the spatially integrated flux with 1 keV spectral resolution. Unfortunately, due to particle precipitation events, only three different times could be analysed during the gradual phase of the flare. Finally, MDI (Scherrer et al. 1995) full-disk magnetograms taken at 14:27:34 UTC and 16:03:34 UTC were used to determine the position of the photospheric apparent magnetic inversion line.
CDS spectra - Spectra acquired after recovery of the SOHO spacecraft are characterised by broad and asymmetric line profiles that were fitted with an opportune template provided within the CDS software. The template consists of a Gaussian component (with all parameters free to vary) and a wing (non-Gaussian) component whose parameters are established for both NIS-1 and NIS-2 spectra.
The analysis of line positions reveals that no trends (either parallel or perpendicular to the slit direction) are present within a single raster. However, the analysis of the raster-averaged NIS-2 line positions versus time reveal a small linear trend that is equal in He I and O V. The shift amounts to 20 km s-1 over two hours. This trend has been accounted for by providing a corrected (shifted) wavelength vector for each NIS-2 line in each raster. No such a trend is visible in the Fe XVI NIS-1 line.
Flows in the He I 58.43 nm, O V 62.97 nm, and Fe XVI 36.08 nm lines have been measured fitting a single component to the line profiles. These lines have an obvious pre-flare component, so that a reference wavelength has been obtained averaging the central positions derived from the fitting for the whole dataset. Velocities have then been computed simply using the Doppler-shift of the flare profiles with respect to the reference wavelength.
The measurement of flows in the hot (K = 6.9) Fe XIX 59.22 nm line requires a more careful analysis. The line appears only during flares, so great caution is needed in choosing the reference profile (see discussions in Brosius 2003; Teriaca et al. 2003). The line profiles obtained averaging all CDS pixels with a radiance greater than 0.16 W m-2 sr-1 are identical in width and position in three successive scans during the late phase of the flare (from 14:53:10 to 15:00:10 UTC, times at raster centres), so their average was used as our reference profile. A multi-component fitting was hence applied to the profiles in various pixels and times, constraining the rest component by imposing width and position of the reference profile, while the background (characterised by the Fe XII 59.26 nm line, see Del Zanna & Mason 2005) was constrained through the average non-flare spectrum. Uncertainties on best-fit parameters are evaluated accounting for data noise (photon statistics, pulse height distribution and readout noise) and the average 3- uncertainty amounts to 15 to 18 km s-1 over the active region for single component fitting and to 25 to 30 km s-1 for double component fitting.
Since RHESSI was in the Earth shadow during the initial phase, no hard
X-ray observations were available
to define the impulsive phase of the flare. Moreover, the low temporal
resolution attained with the CDS rastering mode
prevents the use of EUV radiance peaks as a proxy for hard X-ray (HXR)
bursts, as was done, e.g., by
Brosius (2003). Hence we
used the time-derivative of the GOES SXR flux
as a proxy for the HXR emission (Neupert effect; Neupert 1968).
Both the 0.1 to 0.8 nm (1.6 to 12.4 keV) flux and its derivative are
around the time of
the flare development in Fig. 2.
The main episode of energy release indicated by the Neupert peak
is around 14:41:20 UTC, but significant
signal is present up to the GOES maximum, i.e., up to around 14:48 UTC.
We assume the gradual phase starts after this time.
|Figure 2: Soft X-ray full disk flux at 1 UA in the 1 to 8 Å band around the time of the flare (GOES 8 data). The time derivative (dashed line) represents a proxy for the hard X-ray curve (Neupert effect) and shows its maximum around 14:41:20 UTC.|
Figure 3 outlines the main phases of the flare evolution as seen in the TRACE 17.1 nm channel and H centre. A filament located on the magnetic neutral line is clearly visible in both spectral ranges before the flare. The activation of a portion of this filament is observed around 14:13 UTC in the H images and lasts up to 14:37 UTC, when it becomes clearly visible also in the TRACE images. The eruption was detected by CDS in the O V line, where upward speeds around 30 km s-1 were detected at 14:36:33 UTC.
Few minutes later the footpoints of a large loop system (spanning over km on the solar surface) start to brighten, reaching their maximum around 14:42 UTC. The flare ribbons are clearly visible and nearly identical in both TRACE and H images. Although the east-most ribbon is outside of the FOV of the UBF H images, we verified its existence in the images acquired by the HASTA telescope (see Fig. 1 in Maltagliati et al. 2006). During the impulsive phase of the flare, the ribbons are associated with large coronal upflows and chromospheric and TR downflows (see next section). The brightenings last several minutes and fade off around 14:50 UTC in the 17.1 nm band while their decay is much longer in H.
|Figure 3: Evolution of the flare: TRACE images (logarithmic scale) in the 17.1 nm band are shown in the first two rows and H images in the last two rows. The white box in TRACE images indicate the H FOV. The images at 14:34 UTC show the region just before the flare. The areas (labelled A and B) over which data were integrated to obtain the light and velocity curves shown in Fig. 4 are indicated in black on the TRACE image taken at 14:42:00 UTC. The arrow in the last H image indicates the cool post-flare loops.|
|Figure 4: Light and velocity curves at the footpoints (areas A and B in Fig. 3) of the loop system. For clarity reasons uncertainties are only shown for the Fe XIX speeds. O V, Fe XVI, and TRACE radiances are normalised to pre-flare values while Fe XIX radiances are normalised to the average value reached after the flare (around 15:32 UTC). The average uncertainties in the normalised radiance are of about 0.15. On footpoint A, the average uncertainties in speed are 12 km s-1 (O V) and 20 km s-1 (Fe XVI). On footpoint B, the average uncertainties in speed are 9 and 12 km s-1 for O V and Fe XVI. The shaded area represents the time during which downflows between 4 and 10 km s-1 were observed in H. (We remind that only for footpoint B we have chromospheric observations). On both panels, the scale of downward (positive) velocities has been expanded for sake of clarity.|
Immediately after the GOES maximum, the first RHESSI spectrum (14:50 UTC) between 4 and 10 keV shows a non-thermal component indicating that electron acceleration is still present. During the gradual phase of the flare, two cooling loop systems connecting the ribbons appear in the TRACE images, one around 15:40 UTC, and the other around 16:03 UTC. Figure 3 shows them at the time of maximum radiance. Finally, these loops appear also as dark features in the H images at 16:21 UTC, revealing that the plasma has cooled down to temperatures around 104 K.
In the following, we will concentrate on the footpoints of the second loop system, contained in the CDS FOV, where most of the brightenings and velocity episodes are visible in our dataset. Figure 4 summarises the radiance and velocity evolution in these footpoints, labelled A and B in Fig. 3 (TRACE image of 14:42:00 UTC). A and B are optimised to obtain a high S/N ratio in the smallest region still showing relevant coronal velocities. The area A consists of three CDS pixels (6 ), while the area B, due to the weak emission of the Fe XIX line in the westward ribbon, consists of eleven CDS pixels (224 arcsec2).
At this time the Fe XIX spectrum obtained integrating over A shows a blue-shifted emission with a velocity between -170 km s-1 (single component fitting) and -250 km s-1 (double component fitting, see Fig. 5, top-right). However, given the noise level in the spectrum, it is difficult to establish whether a two component fitting is justified. We thus can assume as a lower limit the value obtained with a single component fit. The same region also shows upward motions around -70 km s-1 in the Fe XVI line. Due to the broad NIS-1 profiles, we have not attempted to separate the strong Fe XVI background component, resulting in a likely underestimate of the upward velocity.
These upward motions of hot plasma seem to
originate in the region between the flow patterns of
opposite sign observed
in TR, but from an area where O V motions are still downward
Furthermore, if we consider a semi-circular loop with radial legs,
the projection effect would amount to about
in the x-direction
and thus suggest a link
between the upward motions in Fe XIX and
the even stronger downflows seen in the TR (the target area is on the west hemisphere,
the loop direction is east-west, and we assume the Fe XIX
emission from a region 5000 km above that of O V and He I).
|Figure 5: Eastward ribbon (A) at 14:40:45 UTC. The top left panel shows a TRACE image in the 17.1 nm band. The bottom left panel shows a map of the O V LOS velocity within the area encompassed by the black box on the TRACE image. On both panels isocontours of the O V radiances at a level of 5 W m2 sr-1 are shown with white thick lines while the thick black isocontours on the velocity map indicate Fe XIX radiances at a level of 90 mW m2 sr-1. The Fe XIX profile shown in the top-right panel is obtained by integrating over three pixels along the CDS slit, centred at the location indicated by a white x on the velocity map (area A defined in Fig. 3). The single pixel O V profiles on the bottom right panel refer to the locations indicated by a black + (thick histogram) and a white + (thin histogram) on the velocity map.|
Around 14:44 UTC, CDS scanned again this flaring kernel (Fig. 6). The resulting Fe XIX profile (bottom panel of the figure) shows a strong blue-shifted component dominating the whole profile. The corresponding velocity is of about -200 km s-1. Upflows of -70 km s-1 are still measured in the coronal Fe XVI. The same velocity pattern seen in O V and He I during the previous passage is still discernible, with the downflows now weaker ( km s-1) and unchanged upflows.
Finally, strong coronal upflows are still visible
in the kernel A when CDS scans the area a third time during the
impulsive phase, around 14:47:30 UTC (-170 km s-1 with a double component
fitting and -100 km s-1 with a single component fitting).
O V and He I do not show any downflow anymore, while the nearby
upflow is still present.
F OOTPOINT B: This footpoint is within the DST FOV. Starting from around 14:41 UTC, the difference of H images obtained at both pm and 150 pm shows a red-wing emission excess in correspondence of the flaring kernels, resulting in a downflow velocity of 4 to 10 km s-1 (cf. Sect. 2.3).
Around 14:42 UTC this area was rastered over by CDS. Velocities in O V and He I indicate only downwards motions (20 km s-1) outlining the tiny bright features visible at the higher spatial resolution afforded by H data. Chromospheric downflows are measured at the same time in the same area, and are indicated with thin black contours in the upper panels of Fig. 7. Upflows of more than 200 km s-1 were measured in Fe XIX by averaging over the whole footpoint area. The corresponding profile of the Fe XIX line is shown in the bottom panel of the same figure. The fitting procedure reveals a blue-shifted component of the same strength of the stationary component, with a speed of -240 km s-1. A single component fit provides a velocity of about -125 km s-1.
During the next raster at 14:45:23 UTC, CDS still measures strong upward speeds in the Fe XIX line (Fig. 4), cospatial and cotemporal with downflows in the low chromosphere (H). The radiance of the blue-shifted component is comparable to that of the stationary component also at this time. At later times, neither the Fe XIX nor the H data show any relevant motions. Hence in this footpoint, the oppositely directed motions measured in low-chromosphere and in hot flare lines, are verified to be co-spatial and co-temporal and to persist for at least 200 s.
|Figure 6: Eastward ribbon (A) at 14:44 UTC. The top left panel shows a TRACE image in the 17.1 nm band while the H image closest in time is shown in the top right panel. On both images the thin black contours indicate chromospheric downflows of 4 km s-1. The Fe XIX line profile shown on the bottom panel is obtained by integrating over three pixels along the CDS slit (area A defined in Fig. 3).|
Comparing the two footpoints, we notice that the Fe XIX blue-shifted component is clearly dominating over the stationary one only in the eastward ribbon (A), possibly due to the smaller area used to obtain reliable line profiles. This fact may then support the assumption that multiple thin loops are present in the flaring ribbons at any time, and that they are heated successively, as in the simulations of Warren & Doschek (2005). This is further confirmed by TRACE and H images that show very clearly the existence of multiple footpoints within the flaring ribbons, brightening at different times. Note also that during the gradual phase (see Sect. 5) post-flare loops of 1 to 2 diameter clearly emerge as independent features in the TRACE images.
Finally, for both footpoints, the large upflows measured in the Fe XIX line are spatially related to patches of downflows in the TR, themselves closely associated to the flaring kernels. However, a significant upflow pattern lasting at least 600 s is also present nearby footpoint A, but we lack of an obvious explanation for it.
Hydrodynamical simulations of explosive chromospheric evaporation predict the equality of momenta between the hot plasma moving upward and the dense cold plasma moving downward during the impulsive phase of a flare (Fisher et al. 1985b). If the oppositely directed flows that we measure during the impulsive phase indeed signal chromospheric evaporation, we should be able to verify such equality. Earlier attempts have been performed, using BCS/SMM coronal data, by Zarro et al. (1988) and Canfield et al. (1990a). These authors find agreement between momenta of the oppositely directed plasma (integrated over the whole impulsive phase) within one order of magnitude, although the lack of spatial resolution in coronal data leaves open the issue of co-spatiality of the flows. We have shown that chromospheric downflows and coronal upflows are co-spatial and co-temporal at least in the case of footpoint B, so we can compare the "instantaneous'' momenta at 14:42 UTC in this footpoint.
The calculation of the downward chromospheric (
and the upward coronal
momenta is described in Appendix A,
where the uncertainties of the measured or
assumed physical parameters are taken into account. We also note that
the calculations of both momenta involve totally
independent sets of physical quantities. For the downward
momentum we obtain
g cm s
g c m s-1.
The range of values depends from the pre-flare hydrogen
density, the chromospheric velocity and the condensation thickness
that we can estimate only within a factor of two.
For the upward momentum, we obtain
g cm s
g cm s-1,
where the lower and upper limits refer to
the velocity and density values obtained assuming a single or a
double-component fit of the Fe XIX line, an iron abundance varying
between photospheric and coronal values and the minimum and maximum
value of the coronal filling factor. Possible departures from ionisation
equilibrium can also contribute to the uncertainty of the result.
In principle we could also attempt to check the momentum balance at 14:45:23 UTC,
when we still measure co-spatial coronal and chromospheric flows on
the area analysed above.
However, at this time the chromospheric density
has already increased to flare values, difficult to estimate without an
ad hoc flare model, which we cannot obtain for lack of suitable
|Figure 7: Westward ribbon (B) at 14:42 UTC. Panels are as in Fig. 6. On both images the thin black contours indicate chromospheric downflows of 4 km s-1 and 7 km s-1. The Fe XIX line profile shown on the bottom panel is obtained by integrating over the area within white solid lines in the H image (area B defined in Fig. 3). A single component fit (plus background) would result in an upflow speed of () km s-1.|
|Figure 8: Gradual phase at a time when the CDS slit was crossing the western footpoint. The top left panel shows a TRACE image in the 17.1 nm band while the H image closest in time is shown in the top right panel. On the TRACE image isocontours of the Fe XIX radiance at levels of (0.11, 0.22) W m-2 sr-1 are indicated by thick white solid lines. On the H image, instead, the white thick contours refer to a RHESSI image obtained between 14:50 and 14:51 UTC in the 6 to 12 keV band. The contours are traced at levels of 40% and 75% of the maximum. The Fe XIX line profile shown on the bottom panel is obtained by integrating over the areas within white solid thin lines in the H image (area B). In this case the radiance of the rest component (dashed line) was set to zero. The dotted line indicates the background (lines plus constant background) defined from pre-flare spectra. The thick solid line indicates the fitted profile.|
The equality within one order of magnitude of the "instantaneous'' chromospheric and coronal momenta supports the explosive chromospheric evaporation model, even though the flare analyzed is quite small and thus unlikely to provide the large energy flux requested by the simulations. More in general, we note that even with simultaneous, spatially resolved observations we cannot obtain a verification of the momentum balance better than one order of magnitude due to inherent uncertainties in the physical parameters.
RHESSI data were available from 14:50 UTC. Images were reconstructed in the 3 to 12 keV band, at 14:50, 14:52 and 15:13 UTC, with an integration time of 60 s using the Pixon algorithm of the standard RHESSI software (Maltagliati et al. 2006). In Fig. 8 we show the TRACE and H images at 14:48 UTC, with overplotted, respectively, the contours of Fe XIX radiance and of the RHESSI image in the 6 to 12 keV band at 14:50 UTC (thick solid white line). The RHESSI emission encompasses the western ribbon and extends eastward outlining a loop structure that is also partially visible from the Fe XIX radiance isocontours (note that the loop top lies outside the CDS FOV), but no emission is registered on the eastern footpoint (A). Thus, also RHESSI data outline a different behaviour of the two flaring ribbons.
The corresponding integrated RHESSI spectrum mainly describes the emission of an area closely overlapping kernel B. This spectrum (Fig. 9, left panel) reveals the broadened emission line feature around 6.6 keV, due to lines of highly ionised iron (Fe XXIV-XXVI) observable in RHESSI data at temperatures above K (Phillips 2004). Indeed, the best fit to the spectrum indicates a thermal component of K, as well as a non-thermal one with a spectral index . values of the same order have already been observed in RHESSI data during the decay phases of small flares (Krucker et al. 2002; Hannah et al. 2004). The presence of non-thermal electrons at these low energies, however, is not obviously linked to any evaporation process, as no upflows are measured in the coronal signatures (bottom panel of Fig. 8). Around 23 min. later (Fig. 9, right panel) the spectrum exhibits only a thermal component with a temperature K.
A coronal plasma heated above 107 K in kernel B might explain the very different behaviour displayed by the two footpoints, both in radiance and velocity. In A, after the end of the impulsive phase, both coronal signatures undergo a rapid decay, while in B the radiance of Fe XIX continues its increase up to 15:05 UTC, i.e., well into the gradual phase. The Fe XVI shows a similar trend, but with a more significant delay, reaching its maximum radiance in B almost one hour after the GOES peak. This seems consistent with the rapid cooling of "hot'' plasma (>107 K), evaporated during the impulsive phase, to Fe XIX temperatures, and successively, to values characteristic of Fe XVI emission.
In this same footpoint, downward directed motions are measured in both
Fe XIX and Fe XVI around the time of their peak emission, although
the amplitudes are close to the error level (for Fe XIX,
(20 to ) km s-1 for a double-component fitting
and (6 12) km s-1 for a single-component fit; for Fe XVI,
() km s-1). Downflows of
km s-1, are further observed in O V to outline
the B footpoint around the time of maximum O V emission during the
gradual phase (15:42 UTC). No obvious downflows are instead measured
in A during the gradual phase.
|Figure 9: RHESSI spectra in the 4 to 10 keV energy range at 14:50 UTC ( left panel) and at 15:13 UTC ( right panel). Notice the presence of a broadened emission feature at about 6.6 keV corresponding to a group of emission lines due to Fe XXIV-XXVI.|
Finally, combining gradual phase
RHESSI data with CDS measurements up to 14:50 UTC, we evaluate
the contribution of the evaporated chromospheric plasma in filling the flare
loops. Using the Fe XIX data, we can estimate the total number of
electrons transported with evaporation during the impulsive phase,
assuming a symmetric evaporation on the two footpoints:
where is the electron density of the upflowing plasma, f the coronal filling factor, the area of the observed upflows (both footpoints), the velocity, and the duration of the evaporation. The electron density and the upward velocity are the same as in Sect. 4.2 (see also Eq. (A.3)) and depend on whether we use one ( km s-1) or two ( km s-1) components to fit the Fe XIX line profile at 14:42:03 UTC and upon the assumed iron abundance. We assume such densities and velocities as typical throughout the evaporation. cm2 accounts for the fact that the plasma is flowing from both footpoints. Finally, s is taken equal to the duration of chromospheric downflows, for which we have a better temporal resolution.
With these values, for the two cases of single and double fitting and for both iron abundances of 7.5 and 8.1, we obtain the minimum and maximum estimate of = [3.8, 10.] electrons.
The above values can be compared with the total number of electrons in the
entire loop system after the end of the evaporation phase.
This can be written as:
|Figure 10: Cooling curve obtained from the time of maximum appearance of progressively colder features. RHESSI temperature at 14:50 UTC refers to the thermal component (see left panel of Fig. 9). Each set of curves (thin and thick) refer to three different values of the volumetric filling factor f (0.1, dot-dashed; 0.3, dashed; 0.5, dotted). The solid part of the curves indicates the times when conduction dominates, while the non-solid part indicates the time when the plasma cools radiatively. The thick curves are obtained assuming constant density only during the conductive phase, while the thin curves are obtained assuming constant density through the entire cooling (see text).|
We compare here the temporal evolution defined by the data with the results of the analytical formulae provided by Cargill et al. (1995) for the cooling time of a flare plasma. The first RHESSI spectrum at 14:50 UTC still shows a non-thermal component (see Sect. 5), suggesting further energy input at the beginning of the gradual phase. Hence we begin the analysis of the cooling at 15:13 UTC, when RHESSI shows only a thermal component. The first data point corresponds to a K, and an average electron density of cm-3 (for f=1). We used the formulae provided for the static case since after 15:13 UTC no evaporative motions are visible anymore. The calculations were performed following two different approaches.
In the first case a constant density was assumed throughout the entire cooling period and the complete radiative loss function of Rosner et al. (1978) was used in the calculations. For this purpose Eqs. (7S) and (8S) of Cargill et al. (1995) have been generalised to any value of the exponent of the power law. Using such equations, we computed the time until when conduction dominates over radiation. The temperature evolution of the cooling plasma was obtained through Eq. (3a) of Cargill et al. (1995) before this time and Eq. (4) after, with the appropriate values for each piece of the radiative loss function of Rosner et al. (1978). The thin lines in Fig. 10 show the results of these calculations for three values of the filling factor f (0.1, dot-dashed; 0.3, dashed; 0.5, dotted). The solid part of the curves indicates the times when conduction dominates, while the non-solid part indicates the time when the plasma cools radiatively.
In the second approach we assumed that the density is constant only during the conductive phase while in the radiative phase (Serio et al. 1991). We also assumed the simplified radiative loss function that fits, within a factor two, the CHIANTI radiative loss function for coronal abundances for . During the time dominated by conduction, the calculations are similar to the previous case. After this time, being , we can use Eq. (6) of Cargill et al. (1995) to calculate the temperature evolution during the time when radiative cooling dominates. The results are represented by the thick curves in Fig. 10, for the same values of f. Again, the solid part represent the conductive cooling while the non-solid part is the radiative one.
Both approaches show that the observed curve can be explained by the cooling of the hot flare plasma injected into the loop system only if a volumetric filling factor of 0.2 to 0.5 is assumed, a value well in agreement with the estimated filling factors of cooling loops given by Aschwanden et al. (2003). We notice that the estimated cooling time in the temperature range is much shorter than the observed one. This is expected because our initial assumptions are no longer valid in this range, in particular due to recombination of hydrogen.
During the impulsive phase of the flare, large coronal upflows of -200 km s-1 or more are measured from the Doppler shift of the Fe XIX line profiles on the ribbons at the footpoints of the flaring loop system. Interestingly, the two footpoints show quite a different behaviour. Footpoint A shows a fast (within our time resolution) rise of the radiances of all our signatures, followed by a slower decay of the Fe XVI and Fe XIX emission. On this footpoint, during the impulsive phase, we can select a small area (3 CDS pixels) where a blue-shifted component clearly dominates the Fe XIX profile. Upward motions of about -70 km s-1 (likely underestimated) are also measured from the Fe XVI CDS data. These coronal upward motions occur in the location where downflows are observed in the O V and He I lines for at least 200 s but, unfortunately, no H Dopplergrams are available for this footpoint.
On footpoint B instead, the radiances of Fe XIX and Fe XVI lines increase very slowly and reach their maxima respectively 20 min and 60 min after the start of the flare. RHESSI data at 14:50 UTC show the emission in the 3 to 12 keV band to be mainly associated to this footpoint, with the spectrum indicating a temperature of the thermal component of 10.4 MK, i.e. larger than the Fe XIX formation temperature (8 MK). This would explain the slow rise of Fe XIX and Fe XVI radiances as due to progressive cooling of this "hot'' plasma. During the impulsive phase, the blue-shifted and the stationary component of the Fe XIX line profile on this area are of comparable strength, and no clear motion is measured from the Fe XVI line. The O V and He I lines show downward ( km s-1) motions outlining the flaring kernels. For this footpoint, within the FOV of the UBF imager, we observed simultaneous H downflows that we estimated to be about () km s-1. To our knowledge, this is the first time that co-spatial, co-temporal, oppositely directed flows during the impulsive phase of a flare are measured in low chromosphere and corona, supporting the model of explosive chromospheric evaporation.
Such direct measurement gave us the possibility to compare the "instantaneous'' momenta of the coronal and chromospheric moving plasma, whose equality is required in the hypothesis of explosive evaporation. We found that such momenta are equal within one order of magnitude. Previous works (Canfield et al. 1990a; Zarro et al. 1988) verified the equality of the momenta within the same limits using data without spatial resolution and integrating over the whole impulsive phase. Although we can utilise spatially resolved data, it remains very difficult to obtain more accurate estimates of the momentum balance, due to the inherent uncertainties in physical parameters such as coronal electron density and abundances, chromospheric hydrogen density and condensation height. To improve the knowledge of these parameters a larger number of spectral lines both in the EUV (including density diagnostics) and in visible/near-infrared (to build a semi-empirical model of the chromosphere) are necessary. However, this is likely to jeopardise the temporal and spatial resolution needed to study the impulsive phase of a flare. New space instruments with larger effective areas, such as the forthcoming EIS on Solar B, should improve the situation on the EUV side.
A further element in support of the chromospheric evaporation was given by RHESSI data acquired during the gradual phase, from which we obtained an estimate of the mass in the coronal loop system. Combining the velocity and radiance data from CDS, we could also estimate the mass evaporated during the impulsive phase of the flare, and found that it provided a large contribution to the density of the X-ray emitting coronal material.
Finally, during the extended gradual phase, the magnetic loops filled by evaporated plasma become visible in all of our spectral lines and band-passes, with progressively cooler features appearing at later times. The observed temporal evolution can be reproduced using the analytical formulae of Cargill et al. (1995). Initially the cooling is dominated by conduction, while radiative losses take over later on. A cooling time comparable to the observed one can be obtained only assuming a volumetric filling factor of 0.2 to 0.5. This in turn is consistent with the the presence of multiple thin loops heated successively, suggested by the observed dominant blue-shifted component in the Fe XIX line profile only when we can use a small area to obtain reliable data.
To end this paper, we would like to make some remarks about the complexity of the TR behaviour. While our observations show downflows in the flaring kernels also from TR lines, in qualitative agreement with some hydrodynamical simulations of explosive evaporation (Fisher et al. 1985a,b), we note that recent hydrodynamic models of chromospheric evaporation (Allred et al. 2005) always predict upflows for the TR plasma. Several observations obtained with CDS during the impulsive phase of flares, with programs privilegeing either the spatial or the temporal resolution, offer a wide range of results. Upflows have been reported by Teriaca et al. (2003), using an observational sequence similar to the one presented in this paper, and Brosius (2003), that adopted a very high cadence program, but with a narrow field of view and limited spatial resolution. Upflow is also measured in our data during the impulsive phase, in the eastern footpoint nearby the area A but not appear to be obviously connected to the coronal upflows.. Downflows are instead reported by Brosius & Phillips (2004), using the same observing program as Brosius (2003), and Kamio et al. (2005) that analyse four flares in a larger FOV but with intermediate temporal resolution. Although at much lower spatial and temporal resolution, the same contradictory indications were obtained with SMM and OSO-8 instruments (see the discussion in Teriaca et al. 2003). It is likely that the behaviour of the TR is strongly related to the structure of the atmosphere depending upon the local magnetic topology. Space observations with substantially higher spatial resolution than now () are most likely necessary to better understanding the dynamics of the TR and provide constraints to theoretical models. Further, we would like to point out that hydrodynamical simulations are generally computed for energetic flares (GOES class X-M) with a much harder spectrum for the impinging electrons ( from 4 to 6) than normally observed with RHESSI ( from 7 to 10) in small flares such as the one studied here.
The authors thank Prof. E. Priest for suggesting calculating the amount of plasma evaporated into the loop and Dr. K. Wilhelm for carefully reading the manuscript. We are also grateful to Dr. V. Andretta and Dr. U. Schühle for useful comments and suggestions. Many thanks are due to the NSO and CDS staffs for their help and assistance in acquiring the data here analysed. We would also like to thank the RHESSI staff for their help in analysing the data. SOHO is a mission of international cooperation between ESA and NASA. NSO is operated by AURA, Inc., under cooperative agreement with the NSF.
The mass density can be written as
is the mean particle weight,
g is the
proton mass, and
the total particle number density.
For a fully-ionised plasma of hydrogen with 10% helium,
while for a neutral gas,
is the hydrogen number density.
Downflow momentum in the chromosphere (neutral gas) can then be written as:
At coronal level, the momentum of the upflowing plasma is given by:
The radiance of the Fe XIX 59.2 nm line is proportional to for and it is, hence, suitable for emission measure calculation. The iron abundance is 8.1 in the corona (Feldman et al. 1992). However, because the plasma is ablated from the chromosphere, the photospheric abundance of 7.6 (Grevesse & Sauval 1998) may also apply. Thus, we will assume an abundance varying between 7.5 and 8.1.
Coronal upflows are measured in eleven effective CDS pixels, so cm2. L is taken to vary between cm ( ) and cm (obtained estimating the geometrical depth of the flaring loop along the line of sight, assuming the loop is perpendicular to the solar surface). The electron density, , of the evaporating plasma (from the emission measure) and the coronal velocity, , depend on whether the Fe XIX spectra are fitted with one or two components. The computed values of are between and cm-3 for f=1. Finally, we note that in the extreme case of magnetic flux tubes not expanding with height, the ratio yields a lower limit of f=0.2, so we assume .
For both values of f, is calculated with the values of EM and obtained in both cases of single and double component fitting and for both estimates of L and . The maximum and minimum of the sixteen obtained values are taken to give g cm s g cm s-1.