A&A 455, 1123-1133 (2006)
DOI: 10.1051/0004-6361:20065065
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
Abstract
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. |
Open with DEXTER |
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
displayed
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. |
Open with DEXTER |
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![]() ![]() ![]() |
Open with DEXTER |
![]() |
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![]() |
Open with DEXTER |
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
directed (
km s-1).
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. |
Open with DEXTER |
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![]() |
Open with DEXTER |
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
spectral signatures.
![]() |
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![]() ![]() |
Open with DEXTER |
![]() |
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![]() ![]() ![]() |
Open with DEXTER |
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. |
Open with DEXTER |
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). |
Open with DEXTER |
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.
Acknowledgements
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
,
where
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,
and
,
while for a neutral gas,
and
,
where
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.