Investigation of the X-Ray Emission of the Large Arcade Flare of 2 March 1993

A large arcade flare, occurring on 2 March 1993, has been investigated using X-ray observations recorded by the Yohkoh and GOES satellites and the Compton Gamma Ray Observatory. We analyzed the quasi-periodicity of the hard-X-ray (HXR) pulses in the impulsive phase of the flare and found a close similarity between the quasi-periodic sequence of the pulses to that observed in another large arcade flare, that of 2 November 1991. This similarity helped to explain the strong HXR pulses which were recorded at the end of the impulsive phase as due to the inflow of dense plasma (coming from the chromospheric evaporation) into the acceleration volume inside the cusp. In HXR images a high flaring loop was seen with a triangular cusp structure at the top, where the electrons were efficiently accelerated. The sequence of HXR images allowed us to investigate complicated changes in the precipitation of the accelerated electrons toward the flare footpoints. We have shown that all these impulsive-phase observations can be easily explained in terms of the model of electron acceleration in oscillating magnetic traps located within the cusp structure. Some soft-X-ray (SXR) images were available for the late decay phase. They show a long arcade of SXR loops. Important information about the evolution of the flare during the slow decay phase is contained in the time variation of the temperature, T(t), and emission measure, EM(t). This information is the following: i) weak heating occurs during the slow decay phase and it slowly decreases; ii) the decrease in the heating determines a slow and smooth decrease in EM; iii) the coupling between the heating and the amount of the hot plasma makes the flare evolve along a sequence of quasi-steady states during the slow decay phase (QSS evolution).

The CGRO/BATSE HXR light curves for two flares: (a) 2 November 1991, (b) 2 March 1993. The red curves show the running mean of the original light curves; see the text for discussion.

Introduction
Quasi-periodic variations were observed in the hard X-ray (HXR) emission of many flares (Lipa, 1978; see also the review of Nakariakov and Melnikov (2009) and references therein). In our previous papers (Jakimiec and Tomczak, 2010 (Paper I), 2012 (Paper II)) we investigated quasi-periodic oscillations in flares with periods P = 10 -60 s, but in Paper I we have found three flares with periods P > 120 s. They turned out to be large arcade flares. Investigation of the quasi-periodic oscillations in such large flares is very important, since their large sizes allow us to investigate the structure of the oscillation volume more comprehensively. Unfortunately, appropriate observations of the X-ray oscillations in such large flares are very rare. In Paper III (Jakimiec and Tomczak, 2013) we investigated such a large arcade flare, one occurring on 2 November 1991. Its HXR light curve is shown in Figure 1a. We have found there direct observational evidence that the strong HXR pulse at 16:34 -16:35 UT is the result of dense plasma coming from the chromospheric evaporation and flowing into the acceleration volume located within the loop-top cusp structure.

Figure 2
The standard GOES soft X-ray light curves for the flare of 2 March 1993 (upper curve 1 -8 Å range, lower curve 0.5 -4 Å range, records from two GOES satellites are displayed). The hatched areas show the Yohkoh satellite nights.
In Paper III we have also found that i) the precipitation of accelerated electrons from the cusp structure is strongly asymmetric, i.e. there is a great difference in the precipitation into the northern and southern arms ("legs") of the flaring loop; ii) there are significant changes of the precipitation with time; iii) the properties of precipitation depend on the energy of the accelerated electrons.
It has been shown that these complicated properties of the precipitation can easily be explained in terms of our model of oscillating magnetic traps (see Paper III).
In the present paper we investigate a large arcade flare, one occurring on 2 March 1993. Its HXR light curve is shown in Figure 1b. We see a close similarity of the light curves in Figures 1a and 1b, which indicates that there is a close similarity in the impulsive-phase development of these two flares. Section 2 contains an analysis of the observations, Section 3 presents a discussion of the decay phase of the flare, and a summary of the paper is given in Section 4.

Observations and Their Analysis
In the present paper we investigate a large arcade flare, which occurred at the eastern limb on 2 March 1993. It was a long-duration event (LDE) of GOES class M5.1 (see Figure 2). The HXR light curves, recorded by the Yohkoh Hard X-ray Telescope (HXT; Kosugi et al., 1991) and the Compton Gamma Ray Observatory Burst and Transient Source Experiment (CGRO/BATSE; Fishman et al., 1992) are shown in Figures 3 and 4a, respectively. The nominal energy range of the BATSE observations is hν > 25 keV, but in Paper III we have found that the actual BATSE energy range was hν > 33 keV. The HXRs began to rise at The Yohkoh hard X-ray (HXR) light curves in four energy ranges. 21:01 UT (beginning of the impulsive phase). Unfortunately, there were no Yohkoh soft-X-ray (SXR) observations for the flare impulsive phase. Some SXR images were available only for the late phase of the flare decay (see Section 2.3).

Analysis of Quasi-periodicity of the HXR Pulses
We have investigated the quasi-periodicity of the HXR pulses using our standard method described in Papers II and III. We have calculated the normalized time series, where F (t) is the measured HXR flux andF (t) is the running average of F (t). The red curve in Figure 4a showsF (t) calculated with averaging time δt = 120 s. The normalized time series, S(t), is shown in Figure 4b. Next we have measured the time intervals, P i , between successive HXR peaks and calculated the period, P = P i , and its standard (r.m.s.) deviation, σ (P ). Our criterion for the quasi-periodicity is σ (P )/P 1. The values of P and σ (P ) are given in Table 1.  Figures 3 and 4, together with Table 1, show a clear quasi-periodicity of the HXR pulses during the impulsive phase rise (21:01 -21:12 UT). Table 1 also shows that the period P is  longer near the impulsive phase maximum than during the impulsive phase rise. We explain this effect as being due to a quick increase in the density inside the acceleration volume, which causes a significant decrease in the Alfvén speed, v A (see Section 3). The light curves in Figures 3 and 4a can be divided into two components: the HXR pulses and a "quasi-smooth" component (emission buried beneath the pulses in the figures). The full amplitude, Amp S, of the function S(t) [Equation (1)] for a HXR pulse is a measure of Figure 5 The Yohkoh HXT 23 -33 keV image for the HXR maximum of the 2 March 1993 flare. The contour levels are 12.5, 25, 37.5, 50, 62.5, 75, and 87.5 percent of I max . BPC is a triangular cusp structure at the top of flaring loop. F1 and F2 are the flare footpoints on the solar disc.
the ratio of the pulse intensity to the quasi-smooth component; the mean values of Amp S are given in Table 1. Our interpretation of the quasi-smooth component has been given in Paper II, where we explained it as a result of the superposition of the emission generated by many magnetic traps whose oscillations are shifted in phase.

Investigation of the HXR Images
The sensitivity of the Yohkoh HXT was moderate and the counting rates were low in this flare (see Figure 3); therefore it was necessary to apply a rather long integration time in the reconstruction of HXR images. Figure 5 shows the 23 -33 keV image recorded during the maximum of the HXR emission. We see the high flaring loop above the eastern solar limb (its altitude was about 45 Mm). There was a triangular cusp structure, BPC, at the top and strong footpoint sources, F1 and F2, on the solar disc. The position of the footpoints indicates that the plane of the flaring loop was tilted from the plane of the image, i.e. that there is a significant geometrical foreshortening in the north-south dimension of the loop. We have seen in Papers I and III that the sources B and C were located at the places where the cusp structure was connected with the arcade channel. Figure 6 shows a sequence of 23 -33 keV images. This long sequence allowed us to investigate the asymmetry in the electron distribution and its changes in the precipitation in some detail. Between 21:06 and 21:08 UT (Figures 6a -c) there was a clear coupling between the top source P and the footpoint F1. During this time interval a gradual increase in the intensity of the P source was seen, which is certainly the result of the increase in density due to the chromospheric evaporation from the footpoint F1. About 21:10 UT ( Figure 6d) there had begun a strong precipitation toward the footpoint F2. This caused chromospheric evaporation, inflow of dense plasma into the acceleration volume, increase in the number of accelerated electrons, and generation of intense HXR pulses. This moment of time (≈ 21:10 UT) is analogous to 16:31 UT in the flare of 2 November 1991 (see Paper III), when the footpoint source F had appeared (see Figure 6 in that paper).   The maximum of the compression is different in different traps, i.e. their χ min values are different. The electrons inside the traps undergoing weaker compression (higher χ min values) achieved lower energies (≈ 15 keV) and the efficiency of their precipitation was low (see Figure 7). The electrons within the traps undergoing strong compression (low χ min values) reached higher energies (≈ 25 keV) and they efficiently precipitated toward the footpoints ( Figure 6). This means that the ensemble of the traps which was responsible for generation of most of the ≈ 15 keV electrons was different from the ensemble providing the ≈ 25 keV electrons. Figure 7 also shows that most of the ≈ 15 keV electrons emitted their energy and were thermalized within the BPC cusp structure. The strong source B in these images indicates that the traps which dominated in the generation of ≈ 15 keV electrons had a good connection with the arcade channel at B. Characteristic features of the HXR impulsive phase are the asymmetry in the precipitation of accelerated electrons from the cusp structure toward the footpoints and changes in the asymmetry with time and with the energy of the electrons (see Papers I -III). These features are clearly seen in Figure 6 of the present paper.  Figure 6f) the precipitation toward the footpoint F1 again dominated, but during the second strong HXR pulse (21:15 UT, Figure 6g) the precipitation toward the two footpoints was again of similar magnitude. In later images (after 21:16 UT, Figure 6h), the precipitation toward the footpoint F2 dominated.
These complicated variations in the asymmetry of precipitation of accelerated electrons can be adequately explained in terms of our model of oscillating magnetic traps.
i) The asymmetry in the precipitation for an HXR pulse may arise if the axis of the symmetry of the magnetic traps (approximately, it is the line joining source P with the middle of the segment BC) is not perpendicular to BC. A small deviation of the symmetry axis from perpendicularity to BC will introduce a small difference in the maximum compression (χ min ) at the opposite ends of the traps. These small differences in the χ min values lead to significant differences in the precipitation of electrons, since the efficiency of the precipitation steeply depends on χ min (see Section 3 in Paper III). Hence, a large asym- metry in precipitation may result from a moderate deviation of the axis of symmetry of the magnetic traps from perpendicularity to the line BC. The observation that for the footpoint intensities F1 > F2 means that (χ min ) 1 < (χ min ) 2 in the cusp, F1 ≈ F2 means (χ min ) 1 = (χ min ) 2 , and F1 < F2 means (χ min ) 1 > (χ min ) 2 . ii) Changes in the asymmetry of precipitation (F1/F2 ratio) mean that the direction of the axis of symmetry of the traps changes in time. Again, small changes in the direction can cause large changes in the asymmetry of precipitation, since small changes in χ min lead to large changes in the precipitation.
X-ray images investigated in Papers I and III have shown two important features of the large arcade flares: i) The triangular ("cusp") BPC structure was magnetically connected with the arcade channel, and therefore the accelerated electrons were able to penetrate into the channel and heat it. ii) The triangular cusp covered only a part of the length of arcade channel, i.e. the channel was significantly longer than the extension of the cusp measured along the channel (see an example in Figure 8). This indicates that the energy which penetrated into the arcade channel was efficiently transferred along the channel by thermal conduction. Unfortunately, the arcade channel could not be seen during the impulsive phase of the flare of 2 March 1993, since no SXR images are available for this phase.

Investigation of the Flare Decay Phase
There were no SXR observations for the impulsive phase of the investigated flare of 2 March 1993. Some SXR images, recorded with the thin Al.1 filter, were available only for the late decay phase of the flare (Figure 9). We see a long arcade of SXR loops at an altitude of about 95 Mm. The response function of the Yohkoh SXT observations with the Al.1 filter weakly depends on temperature in the wide range of temperatures T ≈ 2.5 -20 MK (Tsuneta et al., 1991). Therefore all plasma having T > 2.5 MK efficiently contributed to the recorded emission. The distribution of the intensity in the SXR images displays the distribution of the SXR emitting plasma.

Figure 9
The Yohkoh soft X-ray images for the late decay phase of the 2 March 1993 flare (reverse intensity scale in color). The dashed line shows the solar limb.

Figure 10
The Yohkoh soft X-ray image recorded with the Al.1 filter at 05:02:23 UT (reverse intensity scale in color). The contour levels show the 23 -33 keV HXR image recorded at 21:09 UT. See text for discussion. Figure 10 shows a comparison of an SXR arcade image with an impulsive-phase HXR image. The HXR image was recorded at 21:08 UT, when the arcade channel was at the level of the B and C sources (see Figure 5). Figure 10 shows that the BPC cusp structure covered only a small part of the arcade length (compare this figure with Figure 8).
In Figure 11 we show the time variation of the temperature, T , and the emission measure, EM, derived from the standard GOES observations (Figure 2). During the quick increase of the temperature (21:00 -21:20 UT; impulsive phase) the cusp structure and the arcade channel had been filled with the plasma coming from the chromospheric evaporation.
Between 21:20 and 22:00 UT the temperature decreased, but the emission measure continued to increase. This means that the energy release decreased, but it was sufficiently high to maintain the chromospheric evaporation. The decay phase of the investigated flare started at about 22:00 UT.
In steady-state loops the following condition is fulfilled: (E H is the heating rate per unit volume, E R is the radiative loss per unit volume, the integrals are over the whole volume of a loop), i.e. the total heating is balanced by the total radiative loss from the loop. The "scaling law" of Rosner, Tucker, and Vaiana (1978) is an analytical expression of the energy balance [Equation (2)]. It can be written in the following form (see Bąk-Stȩślicka and Jakimiec, 2005): where N is the mean electron number density in cm −3 , T is the temperature at the top of the loop in K, and L is the semilength of the loop in cm. It is very important that coronal loops have an efficient mechanism which allows them to satisfy the energy balance condition [Equation (2)], i.e. to reach a steady state. Here we briefly describe this mechanism. If E H dV > E R dV then the flux of energy which reaches the footpoints by thermal conduction is large and it generates chromospheric evaporation. This increases the density in the loop and increases E R dV , which allows the loop to achieve the energy balance [Equation (2)]. If, on the other hand, E H dV < E R dV , then the flux of energy which reaches the footpoints is low, since most of the energy is emitted above the footpoints. Therefore the pressure at the footpoints is too low to balance the weight of the plasma contained in the loop and some amount of plasma precipitates to the chromosphere (this is seen in numerical simulations described in Jakimiec et al. (1992)). The density in the loop and E R dV decrease and this allows the loop to achieve the energy balance [Equation (2)]. This mechanism of self-regulation of coronal loops is quick in comparison with the slow evolution of a flare during the slow decay phase. For example, numerical simulations have shown that a loop of semilength L = 20 Mm fits a change in heating rate in the time of only 5 min .
In previous papers we investigated the flare evolution in log T vs. log N or log T vs. log √ EM diagnostic diagrams (see Jakimiec and Bąk-Stȩślicka (2011) and references therein). We have found that during the slow decay phase flares evolve along the line of steady-state loops, i.e. the line described by Equation (3) with L = const. This indicates that during this evolution flares are close to the steady state [Equation (2)] with slowly decreasing heating, i.e. decreasing values of the integrals in Equation (2) (quasi-steady-state, or QSS, evolution of the loops). These results have been supported by numerical simulations of loops with slowly decreasing heating . These results can be summarized as follows: i) During the slow decay phase (after 22:30 UT in Figure 11) the loops seen in Figure 9 slowly evolved along a sequence of steady states (QSS evolution). ii) The loops were continuously heated to support this QSS evolution. The question arises: How was this continuous heating of the loops seen in Figure 9 brought about? Our proposed explanation is the following: i) Magnetic reconnection occurred at the tops of the loops and between the loops and arcade channel. This provided continuous heating of the loops. ii) New magnetic loops which were generated by the reconnection comprised only a small fraction of the volume of loops seen in Figure 9, and they quickly achieved the energy balance [Equation (2)] due to the self-regulation described above. Therefore the reconnecting loops did not much disturb the slow QSS evolution of the loops. This interpretation is supported by the close correlation between the temperature and emission measure during the slow decay phase (d log T /d log √ EM) ≈ 0.5, in agreement with Equation (2) and numerical simulations described in Jakimiec et al. (1992).
The two main processes of thermal energy loss from a flare kernel are thermal conduction and radiative losses. The rates of these losses, calculated per unit volume, are (see Bąk-Stȩślicka and Jakimiec, 2005) and where a is the radius of the X-ray kernel, L is the length of the flaring loop "legs" and (T ) is the radiative loss function. When the observed temperature changes are slow, i.e. the values of dT /dt are small, the loss of energy from the kernel is compensated for by the heating E H : We have applied Equations (3) -(6) to the beginning of the decay phase (22:30 UT in Table 2) and also to the temperature maximum (21:18 UT), since the change in the heating was slow also then (small dT /dt means small dE H /dt ). In Figure 5 we have measured the radius of the cusp BPC structure, a ≈ 7.9 Mm, and the length of the loop "legs", L ≈ 26 Mm. Using Equations (3) -(6) we have calculated the values of the physical parameters, which are given in Table 2. Table 2 shows that the energy release within the cusp (magnetic reconnection at the edges of the cusp, acceleration of the electrons and their thermalization) steeply decreased after the temperature maximum. Figure 11a shows that after 22:00 UT the decrease in the energy release was much slower, which allowed the flare to develop the long QSS decay phase.
We also see in Table 2 that E C > E R , i.e. the conductive loss of energy was larger than the radiative loss (see also Kołomański, 2007). Hence, we can use the following simple approximation in further estimates: E H ≈ E C ∼ T 3.5 . This gives T ∼ (E H ) 0.286 . The small value of the power index in the last formula implies that significant changes in E H induce only weak changes in T (see Table 2). Figure 9 indicates that the energy release during the decay phase occurred in the arcade of SXR loops. It is most probable that the weak energy release which maintained the temperature at T > 6 MK was mostly due to reconnection between the arcade loops and the arcade channel. The bright arcade loops in Figure 9 are the places where the reconnection was most efficient.
We have applied the above Equations (3) -(6) to one of the arcade loops. We have measured the radius in the southern loop-top kernel, a ≈ 11 Mm, and the length of the "leg" connecting the kernel with the loop footpoint, L ≈ 95 Mm. We have obtained E H ≈ 0.023 erg cm −3 s −1 and N ≈ 2.0 × 10 9 cm −3 . Hence, to maintain the long-duration decay phase, a weak energy release in the arcade loops was sufficient.
Important information is contained in the time variation of the emission measure, EM(t) (Figure 11b). This figure shows that during the impulsive phase a large amount of hot plasma had been accumulated within the flaring system (the flaring loop and arcade channel), and during the decay phase this amount of plasma slowly and smoothly decreased (the small peak in EM(t) about 00:25 UT is of instrumental origin; the peaks after 03:00 UT are due to other flares). The relationship between the physical parameters during the decay phase was the following: where (Here changes in the loop semilength, L, and changes in the volume, V , of the SXR emitting plasma have been assumed to be of minor importance.) Combination of the relationships [Equation (8)] gives This simple relationship stresses the fact that the slow decrease in EM during the decay phase is due to a slow decrease in the heating E H . In other words, the observed slow and smooth decrease of the emission measure, EM(t), indicates that the thermal energy release, E H (t), decreased slowly and smoothly during the long decay phase. It is most probable that this weak energy release was mostly due to the reconnection between the magnetic arcade loops and the arcade channel.

Discussion
In Paper III we have found clear observational evidence that the strong HXR pulses at the HXR maximum (see Figure 1a) were the result of the inflow of dense plasma (coming from the chromospheric evaporation) into the acceleration volume inside the cusp structure. For the investigated flare of 2 March 1993 we were not able to monitor the inflow of plasma into the cusp, since we had no SXR imaging observations for the impulsive phase of the flare. We can, however, obtain simple estimates for the increase in density inside the cusp using only the HXR light curves.
In Table 1 we see that the time interval, P i , between the pulses increased during the HXR maximum (P 2 /P 1 = 1.42, where P 1 is the period before 21:10 UT, and P 2 is the period during the HXR maximum). We assume that this increase in the period is the result of an increase in density and we consider two extreme cases: i) We assume that B 2 /8π p inside the traps, where B is the magnetic field strength and p is the pressure. Then the magnetic field does not change significantly during the increase in pressure in the traps and we have where v A is the Alfvén speed. ii) We assume that B 2 /8π ≈ p. Then the increase in pressure causes broadening of the traps and the magnetic field in the traps decreases. For this estimate we assumed that where B 1 is the magnetic field strength before the HXR maximum and B 2 is the field strength during the maximum. Then Combination of Equations (12) and (13) gives and therefore i.e. ρ 2 /ρ 1 = 1.42.
On the other hand, we can estimate the ratio of the densities directly from the ratio, J 2 /J 1 , of the HXR fluxes (J 1 is the flux between 21:05 -21:08 UT and J 2 is the flux at the HXR maximum). For the 23 -33 keV flux the ratio was J 2 /J 1 ≈ 2.0 (see Figure 3). During these time intervals the emission came mainly from the flare footpoints (see Figure 6); hence it is proportional to the number of precipitating electrons per second, which, in turn, is proportional to the electron number density, N , within the cusp. Hence, The last estimate is independent of the previous two. Therefore, we consider these estimates to be a confirmation that the density has quickly increased by a factor between 1.4 and 2.0 (the mean value is 1.8 ± 0.2) and the longer period, P i , at the HXR maximum is the result of this increase in density.

Summary
The sequence of HXR images for the impulsive phase of 2 March 1993 flare has allowed us to investigate the asymmetry in the precipitation of the accelerated electrons from the acceleration volume within the cusp toward the flare footpoints. According to our model of the acceleration of the electrons in oscillating magnetic traps, the precipitation was most efficient during the maximum compression of the traps (see Papers I and III). The asymmetry of the precipitation shows that the maximum of compression was different at the opposite ends of the traps, i.e. the axis of the symmetry of the traps was slightly inclined toward one of the ends of the traps (see Paper III). The changes in the asymmetry from one HXR pulse to another indicate that the inclination of the axis of symmetry of the traps changed. SXR images, which were available for the late decay phase, show a long arcade of SXR loops (Figure 9). In Figure 10 a contrast between the "slim" shape of the flaring loop and the long arcade is seen. This contrast is enhanced by the fact that there were no SXR images for the flare impulsive phase and therefore the arcade channel was not seen during this phase. In other flares such SXR images show a connection between the cusp and arcade channel and the heating of the channel by the cusp (see Figure 8). Similar cases like in Figure 8 will be shown in our next paper. Important information as regards the evolution of a flare during slow decay phase is contained in the time variation of the temperature, T (t), and emission measure, EM(t). This information is the following: i) weak heating occurs during the slow decay phase and it slowly decreases; ii) the decrease in the heating determines a slow and smooth decrease in EM; iii) the coupling between the heating and the amount of the hot plasma makes the flare evolve along a sequence of quasi-steady states during the slow decay phase (QSS evolution).