Different Patterns of Chromospheric Evaporation in a Flaring Region Observed with Hinode/EIS
We investigate the chromospheric evaporation in the flare of 2007 January 16 using line profiles observed by the EUV Imaging Spectrometer (EIS) onboard Hinode. Three points at flare ribbons of different magnetic polarities are analyzed in detail. We find that the three points show different patterns of upflows and downflows in the impulsive phase of the flare. The spectral lines at the first point are mostly blue shifted, with the hotter lines showing a dominant blue-shifted component over the stationary one. At the second point, however, only weak upflows are detected; in stead, notable downflows appear at high temperatures (up to 2.5–5.0 MK). The third point is similar to the second one only that it shows evidence of multi-component downflows. While the evaporated plasma falling back down as warm rain is a possible cause of the redshifts at points 2 and 3, the different patterns of chromospheric evaporation at the three points imply existence of different heating mechanisms in the flaring active region.
Chromospheric evaporation refers to the drastic mass motions in flaring loops caused by rapid energy deposit in chromospheric layers (or probably higher) by non-thermal electrons or by thermal conduction. When the flare energy is transported toward lower layers, it can produce a local overpressure that drives both upward and downward mass motions. These motions can be detected through Doppler shift measurements in chromospheric and coronal lines. Antonucci et al. (1982, 1985), Antonucci & Dennis (1983), Canfield et al. (1987), Zarro & Lemen (1988), Wülser et al. (1994), Ding et al. (1996), and Doschek & Warren (2005) obtained blueshifts of 200–400 km s in the Ca XIX line using the Bent and Bragg Crystal Spectrometer (BCS) onboard the Solar Maximum Mission (SMM; Acton et al. 1980) data and Yohkoh/BCS (Culhane et al., 1991) data, respectively. Similar measurements using data from the Coronal Diagnostic Spectrometer (CDS; Harrison et al. 1995) onboard the Solar and Heliospheric Observatory (SOHO) revealed upflow velocities of 60–300 km s in the Fe XIX line (Teriaca et al., 2003, 2006; Brosius & Phillips, 2004; Del Zanna et al., 2006; Milligan et al., 2006a, b). On the other hand, in chromospheric and transition region lines, redshifts of about tens of km s were observed, implying downward motions (Wülser et al., 1994; Czaykowska et al., 1999; Teriaca et al., 2003, 2006; Brosius, 2003; Kamio et al., 2005; Del Zanna et al., 2006). Recently, Milligan & Dennis (2009) and Chen & Ding (2010) measured Doppler velocities in multi-lines at flare regions to study the evaporation process using Hinode/EIS (Culhane et al., 2007) data . They detected both upflows and downflows in different emission lines.
Chromospheric evaporation can be classified into two types, the explosive evaporation and the gentle one (Fisher et al., 1985a, b, c; Milligan et al., 2006a, b; Brosius, 2009). When the lines formed in the upper chromosphere and transition region are red shifted and the hotter lines are blue shifted, the case is considered as an explosive evaporation. When all the lines appear to be blue shifted, it is considered as a gentle one. In the hydrodynamic simulations of a flaring atmosphere subject to thick-target electron beam heating, Fisher et al. (1985a) found that an energy flux of about 10 ergs cm s can be served as an effective threshold between the gentle and explosive evaporation. Note that in observations, explosive evaporation may appear in both large flares and microflares (Milligan et al., 2006a; Brosius, 2009; Veronig et al., 2010; Brosius & Holman, 2010; Chen & Ding, 2010).
Although many observational results support the evaporation model, there are still some inconsistencies between the observations and models. Some hot lines, such as Ca XIX, Fe XIX, and Fe XXIII, show a dominant stationary component with, however, a relatively weak blue-shifted component indicative of hundreds of km s upflows (Antonucci et al., 1982; Ding et al., 1996; Milligan et al., 2006a; Milligan & Dennis, 2009); while the flare dynamic models predict that these lines should be mostly blue shifted. In the case of explosive evaporation, there is a conversion from redshifts to blueshifts. The models (Fisher et al., 1985c) predict that the downflows occur only at temperatures of 1 MK, which was supported by many obsevations (Kamio et al., 2005; Milligan et al., 2006a; Del Zanna et al., 2006). However, downflows were also detected at much higher temperatures, e.g., 2.0 MK (Milligan, 2008). These pose challenges to the theoretical models.
With the high spatial and spectral resolution data of Hinode/EIS, we investigate in detail the process of chromospheric evaporation in a flaring region using multi-lines with different temperatures. Our work focuses on the two problems as shown above. Besides, we also find that there appear opposite line shifts in regions of different magnetic polarities. In the following, we describe the observations and data reduction in §2. The results are shown in §3. We make a discussion in §4.
2 Observations and data reduction
The observations presented here are for a GOES C4.2 class flare that started at 02:22 UT on 2007 January 16. The flare is located in the core of NOAA AR 10938. Figure 1 shows the Hinode/EIS Fe XII 195 Å intensity map and the magnetogram measured by the Michelson Doppler Imager (MDI; Scherrer et al. 1995) onboard SOHO for the active region. This region presents a bipolar magnetic structure, which is favorable for production of flares. The 1 slit of EIS was used to raster over the AR of 240240 with an exposure time of 5 s requiring a total duration of about 26 min. EIS scanned this region 3 times from 01:54:11 UT, corresponding to the preflare, impulsive, and post-impulsive phases of the flare. Our main interest focuses on the second one between 02:20:30 and 02:46:49 UT that covers the impulsive phase. The area with a FOV of 8050, marked by the black box in Figure 1, contains the flare and is studied in detail here. Figure 2 shows the Ca II H images observed by the Solar Optical Telescope (SOT; Tsuneta et al. 2008) onboard Hinode, the SOT filter (FG) magnetogram, the Fe XII 195 Å intensity map, and the Doppler velocity of the Fe XII line for this region. The flare has two ribbons as seen from the Ca II H images. The right ribbon started to brighten at 02:30 UT; the left one appeared at 02:32 UT and disappeared later than the right one. In particular, we select three points showing significant upflows or downflows in most spectral lines for study. From the Ca II H images and the magnetogram, these points are located at the flare ribbons and close to strong magnetic patches. Point 1 lies in the positive polarity region and at the western part of the right ribbon, while points 2 and 3 lie in the negative polarity region and at the left ribbon. The time for this scan (the second one) is just before the GOES soft X-ray flux reached its maximum at 02:42 UT. Figure 3 shows the GOES 1–8 Å light curve of the flare and the SOT Ca II H intensity evolution at the three points in the impulsive phase. We also mark the EIS scanning time ranges and the time when EIS scanned over the three selected points. Note that the GOES light curve shows two peaks, which may correspond to brightenings of the two flare ribbons.
The EIS spectrum at each location in the raster contains 17 spectral lines. We select 11 lines among them spanning a temperature range of 0.05–16 MK. Details of the lines are shown in Table 1. The majority of these lines are well resolved with no blends or only trivial blends that can be safely ignored in the active region (Young et al., 2007b). We find that most of the observed line profiles have symmetric Gaussian shapes that can be fitted using a single Gaussian function. However, in some areas, the line profiles are asymmetric and can be well fitted by double Gaussian components.
We reduce the data using the standard EIS software data reduction package. This includes the correction of detector bias and dark current, as well as hot pixels and cosmic ray hits, resulting in absolute intensities in ergs cm s sr Å. We also make a correction for a slight tilt of the slit on the CCDs. An additional effect that is corrected for is a variation of line positions over the Hinode orbit due to temperature variations in the spectrometer. Such an orbital variation is obtained by averaging the centroid positions over a length of slit for which the underlying solar region is mostly a quiet region. We choose the bottom 50 rows in the EIS Fe XII raster (see the red box in Figure 1) as a quiet region. The result is then subtracted from the line center positions measured in all wavelength windows.
EIS does not have an absolute wavelength calibration. We adopt the observed line centers averaged over the quiet region (the red box in Figure 1) as the rest wavelengths. Because of the EIS effective area that only peaks at about 195 and 275 Å in the short wavelength (SW) and long wavelength (LW) bands, respectively, and the weak emission in the quiet region, we cannot obtain the accurate line centers for some lines, especially the high temperature lines. Therefore, we use the method by Brown et al. (2007) and the CHIANTI package (Dere et al., 1997, 2009) to determine the reference wavelengths for those lines. To check the reliability of the methods, we use four strong lines, e.g., Fe XII, Fe XIII, Fe XIV, and Fe XV, for test. We find that the wavelengths of the line centers determined using the above methods are nearly the same with a deviation being within 0.003 Å, which induces a velocity uncertainty of no larger than 5 km s.
3.1 Point 1: Upflows-dominated in the positive polarity region
EIS scanned point 1 at 02:33:58 UT in the impulsive phase of the flare (see Figure 3). This point shows upflows in all emission lines except for the He II line. Because of line blending we do not use the Ca XVII and Fe XXIV lines. We also ignore the Fe XXIII line since it is very weak in this region. We plot the line profiles and their fitting curves in Figure 4. Note that in the wavelength window of the Fe X 184.54 line, there exists the Fe XI 184.41 line; however, the latter is quite distinguishable from the former and therefore does not affect the fitting result (see also Figures 5 and 6). Some of the line profiles are fitted by double Gaussian components. We also calculate the intensity ratio of the blue component to the stationary one. The ratio and the Doppler velocity for the blue component are listed in Table 2. From the results we can see that most of the lines show obvious blueshifts. The upflow velocity increases with temperature from several tens of km s to the highest one of 116 km s; the intensity ratio of the two components, derived from the double Gaussian fitting, has also an increasing tendency. For the Fe XIII line, the intensity of the blue component is just 1.72 times that of the stationary component. However, the ratio increases to 9.95 for the Fe XVI line. This indicates that the blue components are dominant over the stationary components for most of the lines, especially the lines with higher formation temperatures.
We also measure the average Doppler velocity over the 9 EIS pixels around point 1 in the preflare and post-impulsive phases and plot the value against the temperature in Figure 7. The same is done for point 2. Note that the velocities corresponding to the blue components from the double Gaussian fitting are marked with the plus symbols. We find that most of the lines show significant blueshifts or blue-shifted components in the impulsive phase; while the shifts are trivial in the preflare and post-impulsive phases.
3.2 Points 2 and 3: Downflows-dominated in the negative polarity region
In contrast to the case of point 1, point 2 shows obvious downflows. EIS scanned point 2 at 02:38:14 UT, slightly before the GOES soft X-ray peak time at 02:42 UT. We fit all the line profiles at point 2 using a single Gaussian function (shown in Figure 5). Different from point 1, point 2 shows strong emission in high temperature lines (e.g., the unblended Fe XXIII line). Young et al. (2007b) reported that the Ca XVII line completely dominates the other lines in large flares. Therefore, we can ignore the blendings of the Ca XVII line with the O V line and the Fe XI line. The blending of the Fe XXIV line with the Fe XI line can also be ignored for the same reason. The average Doppler velocity over the 9 pixels around point 2 is also plotted in Figure 7. We find that most of the lines are red shifted, including the Fe XVI line (2.5 MK) and even the Ca XVII line (5.0 MK); the velocity decreases with temperature in the impulsive phase. Only the Fe XXIII and Fe XXIV lines show blueshifts with speeds of tens of km s. This means that the temperature division between upflows and downflows is about 5 MK. Note that in previous studies, the highest temperature division is 2 MK as reported by Milligan (2008). Our result poses a challenge to the flare dynamic models.
For point 3, scanned by EIS at 02:39:33 UT, we also detect obvious downflows in most of the lines. It shows a similar behavior to point 2 through the entire flare process. One of the differences between them is the temperature division between upflows and downflows in the impulsive phase of the flare. It is about 2 MK at point 3. Another difference is that not all the line profiles at point 3 can be fitted with a single Gaussian function. Some line profiles are asymmetric that can be well fitted with two components (see Figure 6), especially for the low temperature lines. We notice that line blending can cause asymmetric line profiles, such as the possible blending of the Fe XII 195.12 and Fe XII 195.18 lines. However, the Fe X line is not blended with others. The blending in the Fe VIII line can be ignored at the flare ribbons (Young et al., 2007b). Therefore, the apparent asymmetries, as well as the large widths, of the Fe X and Fe VIII lines imply that there exist multi-component downflows in the line formation layers (Dere et al., 1984; Li & Ding, 2009).
4 Discussions and conclusions
We have presented the line profiles and Doppler velocities at flare ribbons of different magnetic polarities to study the chromospheric evaporation process. The key findings of our study are as follows: (1) The flare ribbons with different magnetic polarities can show different patterns of upflows or downflows in the impulsive phase of the flare; (2) the line profiles at one point of positive magnetic polarity are mostly blue shifted, with the blue-shifted component dominating over the stationary one for especially hotter lines; (3) downflows are detected at the points of negative magnetic polarity in lines of relatively high temperatures (up to 2.5–5.0 MK); (4) there exist multi-velocity components in some lines, either cooler or hotter, depending on the points in the flaring regions.
The upflows at point 1 in the positive magnetic polarity region are basically consistent with the scenario of a gentle evaporation process. The most significant finding here is the dominance of blue-shifted components in the line profiles, which is well consistent with the prediction of theoretical dynamic models. In previous observations, high temperature lines (e.g., Fe XIX, Ca XIX, and Fe XXIV) were found to be dominated by a stationary component (Antonucci et al., 1982; Ding et al., 1996; Milligan et al., 2006a; Milligan & Dennis, 2009). It is difficult to explain such observational results in terms of the basic dynamic models. Doschek & Warren (2005) argued that the stationary component is from the top of the flare loop where the evaporated mass is accumulated, or that the mass is moving perpendicular to the line of sight. They also suggested a possibility that the instrumentation is not sensitive enough to detect the earliest blueshifted emission; by the time the emission level has risen sufficiently, the flare loops have already been filled. Falewicz et al. (2009) took into account the geometrical dependence of the line-of-sight velocities of the plasma motions along the loops inclined toward the solar surface as well as a distribution of the flare sites over the solar disc. They concluded that the stationary component can be observed for all flares during their early phases of evolution; on the other side, the blue-shifted component may be undetectable even for plasma moving along the flaring loop with a very high velocity. In spite of these explanations, searching for a dominant blue-shifted component in high temperature lines during flares is still an interesting task to reduce the discrepancy between observations and models. In this work, although we cannot accurately determine the Doppler velocities of the Fe XXIII (12.5 MK) and Fe XXIV (15 MK) lines at point 1 because of the EIS sensitivity limitation and line blending, we confirm that some line profiles are dominated by blue-shifted components. Our result shows further that the intensity ratio of the blue-shifted component to the stationary one increases with temperature.
The counterpart of upflows is downflows in the flaring loop subject to a rapid energy deposition and momentum balance. Downflows were mostly detected in chromospheric lines before the SOHO era (e.g., Ichimoto & Kurokawa 1984; Ding et al. 1995). This is known as the phenomenon of chromospheric condensation (Fisher, 1989). With space instruments, downflows are also seen in some coronal lines. Most previous observations showed that the downflows occur in temperatures of 1 MK (Kamio et al., 2005; Milligan et al., 2006a; Del Zanna et al., 2006). Quite recently, Milligan & Dennis (2009) detected downflows at a temperature of 1.5 MK. Moreover, Milligan (2008) reported a redshift of 14 km s in the Fe XV line (2 MK) at flare footpoints. The fact that the downflows exist in such a high temperature provides a constraint to the dynamic model. Our measurements presented here show that redshifts can appear in lines of even higher temperatures (Fe XVI, Ca XVII) at the points in the negative magnetic polarity region. Most of the line profiles at point 2 show redshifts while only very few are blue shifted. This case can be regarded as an explosive evaporation (Milligan et al., 2006a; Brosius, 2009). If we take the division temperature between upflows and downflows as the site of energy deposition by the electron beam, appearance of the unusually high division temperature at point 2 implies a rather high energy deposition site. A possible reason is that the flare loop has been filled with enough mass before the impulsive phase. A high coronal density can prevent the electron beam from penetrating deeper. Here we do not detect blueshifts at temperatures higher than 15 MK due to the limitation of the EIS dynamic range, but we find that the Fe XXIII and Fe XXIV line intensities at this point are still strong and their profiles are well Gaussian-shaped in this event. This seems to support the hypothesis of a high coronal density. In addition, Brosius (2003) detected redshifted emission that persisted at least 20 min after the cessation of blueshifts and suggested that flare plasma, heated and accelerated upward during the impulsive phase, subsequently cooled and fell back down in what may be thought of as “warm rain”. For the present event, scanning the active region by EIS lasted about 26 min; therefore, we do not have an enough high temporal resolution here. However, from Figure 7, we can find that the downflows still exist in the post-impulsive phase at point 2, contrary to the blueshifts at point 1, which are diminished to nearly zero in the post-impulsive phase. From the SOT Ca II H movie, we find that the flare ribbons are not brightened simultaneously. Therefore, we cannot exclude the warm rain as a possible cause of the downflows.
The three points discussed here are located in different magnetic polarity regions. They show different patterns of mass flows. Point 1 lies in the positive polarity region while points 2 and 3 in the negative polarity region. We check further the magnetic connectivities between the positive and negative magnetic polarities. The available data cannot provide a definite conclusion regarding whether point 1 is magnetically connected to point 2 or 3. However, even if they belong to different magnetic loops, the different evaporation patterns shown at them suggest that in one flaring region, the heating mechanisms and atmospheric conditions may vary from point to point.
- Acton, L. W., et al. 1980, Sol. Phys., 65, 53
- Antonucci, E., et al. 1982, Sol. Phys., 78, 107
- Antonucci, E., & Dennis, B. R. 1983, Sol. Phys., 86, 67
- Antonucci, E., Dennis, B. R., Gabriel, A. H., & Simnett, G. M. 1985, Sol. Phys., 96, 129
- Brosius, J. W. 2003, ApJ, 586, 1417
- Brosius, J. W., & Phillips, K. J. H. 2004, ApJ, 613, 580
- Brosius, J. W. 2009, ApJ, 701, 1209
- Brosius, J. W., & Holman, G. D. 2010, ApJ, 720, 1472
- Brown, C. M., et al. 2007, PASJ, 59, 865
- Canfield, R. C., Metcalf, T. R., Strong, K. T., & Zarro, D. M. 1987, Nature, 326, 165
- Chen, F., & Ding, M. D. 2010, ApJ, 724, 640
- Culhane, J. L., et al. 1991, Sol. Phys., 136, 89
- Culhane, J. L., et al. 2007, Sol. Phys., 243, 19
- Czaykowska, A., de Pontieu, B., Alexander, D., & Rank, G. 1999, ApJ, 521, L75
- Del Zanna, G., Berlicki, A., Schmieder, B., & Mason, H. E. 2006, Sol. Phys., 234, 95
- Del Zanna, G., Mitra-Kraev, U., Bradshow, S. J., Mason, H. E., & Asai, A. 2010, A&A, in press
- Dere, K. P., Bartoe, J.-D. F., & Brueckner, G. E. 1984, ApJ, 281, 870
- Dere, K. P., Landi, E., Mason, H. E., Monsignori Fossi, B. C., & Young, P. R. 1997, A&AS, 125, 149
- Dere, K. P., Landi, E., Young, P. R., Del Zanna, G., Landini, M., & Mason, H. E. 2009, A&A, 498, 915
- Ding, M. D., Fang, C., & Huang, Y. R. 1995, Sol. Phys., 158, 81
- Ding, M. D., Watanabe, T., Shibata, K., Sakurai, T., Kosugi, T., & Fang, C. 1996, ApJ, 458, 391
- Doschek, G. A., & Warren, H. P. 2005, ApJ, 629, 1150
- Falewicz, R., Rudawy, P., & Siarkowski, M. 2009, A&A, 508, 971
- Fisher, G. H., Canfield, R. C., & McClymont, A. N. 1985a, ApJ, 289, 414
- Fisher, G. H., Canfield, R. C., & McClymont, A. N. 1985b, ApJ, 289, 425
- Fisher, G. H., Canfield, R. C., & McClymont, A. N. 1985c, ApJ, 289, 434
- Fisher, G. H. 1989, ApJ, 346, 1019
- Guo, Y., Ding, M. D., Jin, M., & Wiegelmann, T. 2009, ApJ, 696, 1526
- Harrison, R. A., et al. 1995, Sol. Phys., 162, 233
- Ichimoto, K., & Kurokawa, H. 1984, Sol. Phys., 93, 105
- Kamio, S., Kurokawa, H., Brooks, D. H., Kitai, R., & UeNo, S. 2005, ApJ, 625, 1027
- Li, Y., & Ding, M. D. 2009, Research in Astronomy and Astrophysics, 9, 829
- Milligan, R. O., Gallagher, P. T., Mathioudakis, M., Bloomfield, D. S., Keenan, F. P., & Schwartz, R. A. 2006a, ApJ, 638, L117
- Milligan, R. O., Gallagher, P. T., Mathioudakis, M., & Keenan, F. P. 2006b, ApJ, 642, L169
- Milligan, R. O. 2008, ApJ, 680, L157
- Milligan, R. O., & Dennis, B. R. 2009, ApJ, 699, 968
- Neupert, W. M. 1968, ApJ, 153, L59
- Scherrer, P. H., et al. 1995, Sol. Phys., 162, 129
- Teriaca, L., Falchi, A., Cauzzi, G., Falciani, R., Smaldone, L. A., & Andretta, V. 2003, ApJ, 588, 596
- Teriaca, L., Falchi, A., Falciani, R., Cauzzi, G., & Maltagliati, L. 2006, A&A, 455, 1123
- Tsuneta, S., et al. 2008, Sol. Phys., 249, 167
- Veronig, A. M., et al. 2010, ApJ, 719, 655
- Wülser, J.-P., et al. 1994, ApJ, 424, 459
- Young, P. R., Del Zanna, G., Mason, H. E., Doschek, G. A., Culhane, L., & Hara, H. 2007a, PASJ, 59, 727
- Young, P. R., et al. 2007b, PASJ, 59, 857
- Zarro, D. M., & Lemen, J. R. 1988, ApJ, 329, 456