PARALLEL EVOLUTION OF QUASI-SEPARATRIX LAYERS AND ACTIVE REGION UPFLOWS

, , , , , , , , and

Published 2015 August 11 © 2015. The American Astronomical Society. All rights reserved.
, , Citation C. H. Mandrini et al 2015 ApJ 809 73 DOI 10.1088/0004-637X/809/1/73

0004-637X/809/1/73

ABSTRACT

Persistent plasma upflows were observed with Hinode's EUV Imaging Spectrometer (EIS) at the edges of active region (AR) 10978 as it crossed the solar disk. We analyze the evolution of the photospheric magnetic and velocity fields of the AR, model its coronal magnetic field, and compute the location of magnetic null-points and quasi-sepratrix layers (QSLs) searching for the origin of EIS upflows. Magnetic reconnection at the computed null points cannot explain all of the observed EIS upflow regions. However, EIS upflows and QSLs are found to evolve in parallel, both temporarily and spatially. Sections of two sets of QSLs, called outer and inner, are found associated to EIS upflow streams having different characteristics. The reconnection process in the outer QSLs is forced by a large-scale photospheric flow pattern, which is present in the AR for several days. We propose a scenario in which upflows are observed, provided that a large enough asymmetry in plasma pressure exists between the pre-reconnection loops and lasts as long as a photospheric forcing is at work. A similar mechanism operates in the inner QSLs; in this case, it is forced by the emergence and evolution of the bipoles between the two main AR polarities. Our findings provide strong support for the results from previous individual case studies investigating the role of magnetic reconnection at QSLs as the origin of the upflowing plasma. Furthermore, we propose that persistent reconnection along QSLs does not only drive the EIS upflows, but is also responsible for the continuous metric radio noise-storm observed in AR 10978 along its disk transit by the Nançay Radio Heliograph.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The EUV Imaging Spectrometer (EIS: Culhane et al. 2007), on board the Hinode satellite (Kosugi et al. 2007), has provided observations of the presence of plasma flows in various solar environments, from coronal holes to active regions (ARs). One of the most remarkable EIS results was the finding of ubiquitous plasma upflows seen at the borders of ARs (Harra et al. 2008). These are located in regions of low electron density, low radiance, and over monopolar areas (Del Zanna 2008; Doschek et al. 2008; Harra et al. 2008). They have been observed to persist at nearly the same location from a day to at least one week (Doschek et al. 2008; Démoulin et al. 2013). The Fe xii 195.12 Å blueshifted line-of-sight velocities typically range from a few to 50 km s−1 and are faster in hotter coronal emission lines. Démoulin et al. (2013) carried out a detailed analysis of the evolution of upflows in an AR during its disk transit. They concluded that the global temporal variation of the velocities was consistent with a quasi-static flow subjected to a projection effect along the line of sight on upflows tilted from the vertical away from the AR core.

Several driving mechanisms have been proposed to explain their origin (see Baker et al. 2009 and references therein). In particular, noticing that in the analyzed examples the upflows appeared at locations where magnetic field lines with drastically different connectivities were anchored, Baker et al. (2009) proposed that in their case study (AR 10942) magnetic reconnection between closed field lines of the AR and either large-scale externally connected or "open" field lines was a viable mechanism for driving the upflows.

Del Zanna et al. (2011) analyzed EIS upflows in two ARs (AR 10961 and 10955) and found a null point high above in the coronal field in both ARs. These authors suggested that the continuous growth of the ARs maintained a steady reconnection process at the null point. In their view, interchange reconnection occurred between closed, high-density loops in the core of the AR and neighboring open, low-density flux tubes. In this way, magnetic reconnection created a strong pressure imbalance, which was the main driver of the plasma upflows (Bradshaw et al. 2011).

As in flares, upflows are not only related to reconnection at null points but also at quasi-separatrix layers (QSLs; see the examples studied by Baker et al. 2009; van Driel-Gesztelyi et al. 2012). QSLs are defined as thin volumes in which field lines display strong connectivity gradients (Démoulin et al. 1996a). When these gradients become infinitely large, a QSL becomes a separatrix. QSLs, like separatrices, are places where strong currents can form during the evolution of a magnetic field having a high Lundquist number (Aulanier et al. 2005; Büchner 2006; Effenberger et al. 2011; Janvier et al. 2014). Therefore, QSLs are natural locations where magnetic reconnection can take place (Priest & Démoulin 1995; Démoulin et al. 1996a). This was confirmed by MHD numerical simulations (Milano et al. 1999; Aulanier et al. 2006, 2010; Wilmot-Smith et al. 2010; Janvier et al. 2013), by kinetic numerical simulations (Wendel et al. 2013; Finn et al. 2014), and in laboratory plasmas (Lawrence & Gekelman 2009; Gekelman et al. 2012). Hesse & Schindler (1988) and Schindler et al. (1988) developed a general framework for three-dimensional (3D) reconnection based on the description of the magnetic field using Euler potentials and localized non-ideal regions. QSLs are the locations where these non-idealnesses can occur; therefore, the two approaches are complementary (Démoulin et al. 1996b; Richardson & Finn 2012).

Separatrices and null points are also embedded in QSLs in which reconnection complies with special properties, such as the slippage of field lines (Masson et al. 2009, 2012). Complex magnetic configurations involving both a null point and QSLs have been found to be associated with some observed upflows. In a quadrupolar configuration, formed by AR 10980 and a neighboring magnetic bipole, van Driel-Gesztelyi et al. (2012) showed that plasma upflows observed with EIS were co-spatial with QSL locations, including the separatrix of a null point for a fraction of the upflows. Global potential-field source-surface (PFSS) modeling indicated that part of the upflowing EIS plasma could access the solar wind along reconnected field lines, which extended up to the source surface, passing through the vicinity of the null point.

However, upflowing plasma from other ARs may not have direct access to the solar wind. For example, AR 10978 was an isolated bipolar region that, according to a PFSS model, was completely covered by the separatrix surface of a helmet streamer (Culhane et al. 2014). Brooks & Warren (2011) and Brooks et al. (2012) analyzed EIS upflows and found that the abundance of Si was always enhanced over that of S by a factor of three to four (a classical value for FIP-bias enhancement in the corona). When the AR's western side was oriented in the Earth direction, the Si/S ratio, measured with the Solar Wind Ion Composition Spectrometer (Gloeckler et al. 1998) on board the Advanced Composition Explorer a few days later, was found comparable to the Si/S abundance ratio measured in the corona. This provided evidence to support a connection between the solar wind and the coronal plasma in the upflow region. Culhane et al. (2014) concluded that, even though AR 10978 was isolated and completely covered by closed streamer field lines, the coherent magnetic field, proton velocity, and density variation at L1, together with the matching FIP-bias evolution at the Sun and L1, were clear proof of the presence of AR plasma in the slow solar wind. Based on a global topology computation and analysis of noise-storm radio signatures, Mandrini et al. (2014b) proposed that the AR plasma could reach the solar wind via a two-step reconnection process. The first step was proposed to occur in the AR between closed AR loops and long externally connected loops, while the second step was shown to involve the large-scale global coronal field at a high altitude null point. Only this second step was analyzed in depth by Mandrini et al. (2014b), while the present work completes the previous study focussing on the role of QSLs during the first reconnection step.

The spatial relation between upflow and QSL locations and magnetic field line traces and connectivity, led Baker et al. (2009) and van Driel-Gesztelyi et al. (2012) to suggest that magnetic reconnection at QSLs was at the origin of EIS upflowing plasma. In this article, we put forward a proof of the concept. We analyze the temporal evolution of EIS upflows in relation to QSL evolution as AR 10978 crosses the solar disk during Carrington rotation (CR) 2064. Our results indicate that the evolution of the AR magnetic field leads to the evolution of QSL locations, which in turn leads to a spatial evolution of EIS upflows.

Our article is organized as follows. In Section 2.1, we describe the evolution of the photospheric magnetic and velocity fields of AR 10978. We briefly discuss the EIS upflow evolution in Section 2.2. In Section 3, we model the AR coronal field (Section 3.1) and search for the presence of magnetic null points (Section 3.2). Since reconnection at nulls cannot explain the observed upflows, we further analyze QSLs in Section 4. In particular, we demonstrate the spatial and temporal relation between upflow regions and QSLs (Section 4.2). Based on this analysis, combined with the photospheric magnetic field evolution, we unveil the characteristics of upflows originating from different AR locations (Section 4.3). Next, in Section 5, we find the presence of weak noise storms that remain located above the AR during its whole disk transit. We associate their origin to the magnetic reconnection process occurring at the AR QSLs. Finally, in Section 6, we summarize our results and draw our conclusions.

2. EVOLUTION OF AR 10978 DURING ITS DISK TRANSIT

2.1. Magnetic Field Evolution

AR 10978, observed with the Michelson Doppler Imager (MDI, Scherrer et al. 1995) on board the Solar and Heliospheric Observatory, rotated onto the disk on 2007 December 7. By this time, it was a mature AR. From its state of evolution and the separation of its opposite-sign polarity spots, it is likely that the AR was at least three to five days old when it appeared on the east solar limb (van Driel-Gesztelyi & Green 2015).

We checked data from the Solar-Terrestrial Relations Observatory spacecraft (STEREO-A and B) in order to constrain its emergence time using data from the Extreme-ultraviolet Imager (Wuelser et al. 2004). On the date of the AR appearance, STEREO A and B were separated from the Sun–Earth line by about 21° in each direction. In the 195 Å data of STEREO-B, the AR coronal loops could already be cearly seen over the east limb on December 5. STEREO-A data showed a potential first sign of flux emergence at the future location of AR 10978 on November 21–22 at the west limb, indicated by a significant increase of coronal emission from that location (http://stereo-ssc.nascom.nasa.gov/cgi-bin/images). Therefore, the first indications of the AR appearance could have been as old as 18 days when its magnetic field was first mapped with MDI. However, this early flux emergence (on November 21–22) might not have been a major episode and, perhaps, the flux decayed quickly.

AR 10978 had a peak mean magnetic flux of 3 × 1022 Mx (Mandrini et al. 2014b), which corresponds to the large AR category (see, e.g., van Driel-Gesztelyi & Green 2015). Such ARs may survive several months; indeed, AR 10978 returned in the following rotation as AR 10980 (see, e.g., van Driel-Gesztelyi et al. 2012) and its location was magnetically active for several more solar rotations.

From December 9 to 13, magnetic flux emergence continued within the AR (Figure 1). We were able to follow five significant flux emergence episodes, which are indicated in the figure with red ellipses and numbered as EFR1–EFR5 (emerging flux region; EFR). These occurred between the pre-existing polarities marked as PP and PF in the top left panel of Figure 1. Emerging bipoles are characterized by the divergence of opposite-sign polarity flux concentrations; therefore, these EFR sites are locations of fast photospheric motions. The evolution of the border of the supergranule to the north of PP is also noteworthy. This supergranule evolved as minor polarities emerged within it (not visible in Figure 1 because of the chosen high saturation level). By December 11, the supergranule border is no longer visible.

Figure 1.

Figure 1. Magnetic field evolution of AR 10978 for six days, from 2007 December 9 to 14. Five flux emergence episodes, indicated as EFR1–EFR5, are observed along this period. They occurred between the two pre-existing polarities marked as PF and PP, which emerged on the far side of the Sun. The saturation level is 800 G and the axes are in arcseconds. The vertical axis is the AR location in the solar north–south direction with the origin corresponding to the solar equator, while the horizontal axis indicates the east–west coordinates derotated to the time when the AR was close to the central meridian.

Standard image High-resolution image

We have computed the transverse flows of the photospheric magnetic field features employing a local correlation tracking technique (LCT, November & Simon 1988). The proper motions of the magnetic elements over the MDI sequence of magnetograms (spatial resolution 1farcs98) was computed using a Gaussian tracking window with an FWHM of 10''.

Different values for the correlation window were used, varying from 6'' to 14''. The general pattern of the horizontal velocity vectors did not vary significantly from case to case, in the sense that we obtained similar results in terms of their distribution by direct visual inspection of the computed flow maps. Therefore, the size of the correlation window (10'' ) was selected following a criterion based on reducing the level of noise and producing a coherent tracking of the magnetic features. Such window size is within the interval of sizes used in previous works (e.g., Nindos et al. 2003; Chae et al. 2004; Vemareddy et al. 2012).

The time cadence of the LCT, 96 minutes, is fixed by the observations. It filters the evolution of faster timescales and, in particular, the amplitude of the deduced velocity decreases (e.g., see Figure 7 of Chae et al. 2004). However, since the upflows that we analyze evolve on a timescale of days, we expect the 96 minute cadence to be sufficient to derive the photospheric flows relevant to understand these upflows. The time series for the analysis spans from 22:23 UT on 2007 December 9 to 06:23 UT on 2007 December 11. Prior to applying the LCT method, the sequence of images was aligned to eliminate the possible jitter and the solar rotation of the observed target within the field of view (FOV), which results in flow maps derotated to the AR central meridian passage (CMP) position. Finally, it should be noted that whereas the LCT results give a reliable characterization of the flow patterns, the absolute values of the transverse velocities should be taken with caution because this technique, in general, underestimates them (November & Simon 1988; Vargas Domínguez et al. 2008).

The maps of the transverse flows in Figure 2 display the velocities, using arrows, for pixels with magnetic field values over a threshold of 800 G. The computed flows show a global diverging pattern of the main photospheric magnetic polarities. Emergence episodes are also detected in the flow map, see, e.g., the one corresponding to EFR1 in Figure 2(a). Another noticeable flow pattern is that on the western portion of the northern negative polarity (Figures 2(b)–(d)). It indicates a motion of the magnetic features from the northeast to the southwest.

Figure 2.

Figure 2. Maps of the photospheric transverse velocities displayed over magnetograms. The flow maps, generated from a sequence of magnetograms previously derotated to central meridian passage position, are derived from an LCT analysis in the interval indicated in each map (top-left corner) employing an FWHM correlation window of 10''. Red and blue arrows represent velocities of up to 0.5 km s−1, the segment at the bottom left corner of panel (a) indicates the velocity scale. These velocities are averaged over the time intervals indicated in each panel. Notice that the velocity vectors show a global separation of the leading and following polarities as expected in a mature AR. The magnetic field background corresponds to the average image computed in the flow-map time interval. Positive (white) and negative (black) polarities are saturated at ±800 G. The field of view is the same as in Figure 1 with the origin set at the lower left corner.

Standard image High-resolution image

2.2. EUV Upflow Evolution

EIS is a raster-scanning instrument capable of constructing large FOV, up to ≈600'' in the dispersion direction and 512'' in the slit direction, using the 1'' and 2'' slits and the 40'' slot. It has a spectral resolution of 22.3 mÅ and 1'' pixels. EIS observes in two wavelength ranges 170–210 Å and 250–290 Å.

AR 10978 is selected for this study because, with the exception of a few other ARs (see, e.g., Del Zanna 2008; Del Zanna et al. 2011), it has the best EIS spatial and temporal coverage of an AR from limb to limb. This AR was tracked from 2007 December 6 to 19 using the full complement of EIS slits. In this article, we use a set of five large FOV slit rasters around the AR (CMP, see Figures 3 and 6), which occurred on December 11 at approximately 22:00 UT. The FOV cover both of the main AR polarities (see Table 1 in Démoulin et al. (2013) for more details). We selected EIS observations of the Fe xii emission line at 195.12 Å (T ≈ 1.4 MK) because it is the strongest line within EIS's wavelength ranges. We have also used EIS 40'' slot rasters of the AR that allow us to observe the bright coronal loops, which extend beyond the FOV of the slit rasters and are useful to constrain the free parameter of our magnetic field model (see Figure 4).

Figure 3.

Figure 3. Doppler velocity maps around AR 10978 CMP obtained from the Fe xii 195.12 Å EIS emission-line profiles. The color bar at the bottom indicates the velocity scale. The full pink circle surrounds low velocity regions, while the dashed pink lines separate upflow streams with different characteristics. Additional Doppler maps obtained with EIS have been published by Brooks & Warren (2011) and Démoulin et al. (2013). The vertical and horizontal axes indicate the position on the Sun in arcseconds.

Standard image High-resolution image
Figure 4.

Figure 4. (a) EIS Fe xii emission line intensity map of AR 10978 on December 11 at 10:25 UT overlaid with MDI contours with the same values as those in the right panel (drawn in white (black) for positive (negative) values). (b) Modeled MDI magnetic field contours (±100, ±500 G, positive (negative) shown in magenta (blue) color) on December 11 at 11:11 UT overlaid on the EIS slot image shown in grayscale. A set of computed field lines matching the global shape of the observed coronal loops has been added in a continuous line and orange color. The vertical and horizontal axes indicate the position on the Sun in arcseconds.

Standard image High-resolution image

EIS slit and slot raster data were processed using standard SolarSoft EIS routines to correct for dark current, cosmic rays, hot, warm, and dusty pixels and to remove instrumental effects of slit tilt and orbital variation in the line centroid position due to thermal drift. Doppler velocities for slit rasters were calculated by fitting a single-Gaussian function to the calibrated Fe xii spectra in order to obtain the line center for each spectral profile. Reference wavelengths were determined using the average wavelength value of a relatively quiescent Sun patch within each raster. Velocity maps follow the standard convention of blueshifts (redshifts) corresponding to negative (positive) Doppler velocity shifts along the line of sight (see Figures 3 and 6).

As an example, we show three Doppler velocity maps close to the AR's CMP in Figure 3. These are drawn using standard IDL routines. Following Démoulin et al. (2013), we have added guide marks in Figure 3 that help to visualize and track the main upflow structures. The pink circles surround weak velocity patterns (because of the projection effect), while the pink dashed lines separate flow streams.

When looking at the upflow pattern in Figure 3, a clear evolution is evident. This was interpreted by Démoulin et al. (2013) to be the signature of a projection effect on steady upflows that were inclined from the vertical at an angle that was larger to the east than to the west (see Figure 14 in Démoulin et al. 2013).

3. MAGNETIC FIELD CONFIGURATION

3.1. Magnetic Field Model of AR 10978

To compute the magnetic field topology of an AR, we first model its coronal field. We extrapolate the line-of-sight magnetic field of AR 10978 to the corona using the discrete fast Fourier transform method (as discussed by Alissandrakis 1981) under the linear force-free field (LFFF) hypothesis (${\boldsymbol{\nabla }}\times {\boldsymbol{B}}=\alpha {\boldsymbol{B}}$, with α constant). An example of extrapolation is shown in Figure 4(b) when the AR is at disk center on 2007 December 11. We use the MDI magnetogram at 11:11 UT as the boundary condition for the magnetic model. This magnetic map is the closest in time to the EIS slot image in Fe xii 195.12 Å (Figure 4(a)), which is large enough to identify the global shape of the coronal loops. The value of the free parameter of the model, α, is set to best match the observed loops following the procedure discussed by Green et al. (2002). The best-matching value is α = −3.1 × 10−3 Mm−1 (Figure 4(b)).

A region four times larger than that encompassed by AR 10978, and centered on it, is selected from each MDI full-disk magnetogram as the boundary condition for each extrapolation. This magnetic map is embedded within a region twice larger padded with a null vertical field component for two reasons. First, to decrease the modification of the magnetic field values since the method to model the coronal field imposes flux balance on the full photospheric boundary (i.e., the flux imbalance is uniformly spread on a larger area, so the removed uniform field is weaker as the area is larger). Second, to decrease aliasing effects resulting from the periodic boundary conditions used on the lateral boundaries of the coronal volume. The photospheric boundary condition is then written in a 1024 × 1024 horizontal grid to maintain the spatial resolution of the observations. This allows us to distinguish field lines that connect to the surrounding quiet-Sun regions from those that are potentially "open" lines because they leave the extrapolation box.

To compute the topology at different times during the AR transit, we use the same value of α with the corresponding MDI map as the boundary condition. We have checked that the large-scale loops observed in EIS slot images on the dates used in our analysis can be globally fitted with this α value. Indeed, the magnetic configuration of AR 10978 has a low shear and the coronal magnetic field configuration is mainly evolving due to the evolution of the photospheric boundary condition rather than to the change in the α value. We also carry out a transformation of coordinates from the local AR frame to the observed one (see Démoulin et al. 1997) to obtain a model for which the QSL locations can be compared to EIS upflows (Section 4).

3.2. Magnetic Null Points in the AR Neighborhood

The origin of some EIS upflows has been attributed, either directly or indirectly, to magnetic reconnection occurring in the vicinity of coronal magnetic null points (Del Zanna et al. 2011; van Driel-Gesztelyi et al. 2012; Mandrini et al. 2014b). Therefore, after summarizing the main properties of null points in the next paragraph, we investigate if the observed upflows in AR 10978 are due to magnetic reconnection at these null points.

The field connectivity around a null point is characterized by the presence of so-called spines and fans (see, e.g., Longcope 2005; Pontin 2011). A fan surface separates the coronal volume into two connectivity domains, while all field lines in the close vicinity of the fan converge to the associated spine line. The local field connectivity around a null point can be found using the linear term of the Taylor expansion of the magnetic field (see Démoulin et al. 1994; Mandrini et al. 2006, and references therein). From the diagonalization of the Jacobian matrix of the field, one finds three eigenvectors and the corresponding eigenvalues, which add up to zero to locally satisfy the field divergence-free condition. The eigenvalues are real for coronal conditions (Lau & Finn 1990). A positive null point has two positive fan eigenvalues and conversely for a negative null. The fan surface is defined by all of the field lines starting at an infinitesimal distance from the null in the plane defined by the two eigenvectors that correspond to the eigenvalues that have the same sign. In a similar way, the spines are defined by the field lines tangent to the third eigenvector.

We computed the location of magnetic null points for the MDI magnetograms closest in time to the EIS velocity maps described in Section 2.2. None of the null points found were located above the AR and most of the ones above quiet-Sun regions had no field lines connecting to the AR polarities. Figure 5 illustrates the location of all magnetic null points at heights above 20 Mm, which are related to AR 10978 for an MDI magnetogram on December 12. The location of the null points is indicated by the intersection of three segments that correspond to the direction of the three eigenvectors of the Jacobian matrix. These segments are color coded to indicate the magnitude of the corresponding eigenvalue. For a negative null, dark blue (light blue) corresponds to the highest (lowest) negative eigenvalue in the fan plane and red to the spine eigenvalue. All null points in Figure 5 are negative. Their spines, drawn in black in the figure, are connected to a quiet-Sun polarity at one end and to the main negative polarity of AR 10978 at the other. As an example, we have drawn in orange field lines belonging to the fan of one of the null points. Since these nulls are above the quiet Sun, they are associated with a magnetic field intensity, B, lower than in the AR (B lies in the range [20, 50] G). Such field values are present in small polarities with sizes ≈6 pixels, so they are not within the magnetogram noise. In fact, the height of these nulls, above 20 Mm, confirms that they are associated with real local magnetic polarities.

Figure 5.

Figure 5. Magnetic field model of AR 10978 and null points with a topological structure linked to its main polarities. The model corresponds to the MDI magnetic map at 03:11 UT on 2007 December 12. The locations of all magnetic null-points at heights greater than ≈20 Mm above the photosphere are shown as the intersection of three colored segments (see the text). We have drawn the null-point spines as black continuous lines and have added a set of field lines (orange continuous lines) in the vicinity of the fan plane below one of the nulls. The external western QSL trace, depicted as a thick black continuous line, overlays the Bn < 0 polarity (see Section 4). The axes in this panel are in Mm, with the origin set in the AR and the isocontours of the line-of-sight field correspond to ±100, ±500 G in continuous magenta (blue) style for the positive (negative) values.

Standard image High-resolution image

It is highly improbable that magnetic reconnection at these null points can drive all of the EIS upflows in AR 10978. On one hand, none of the null points has field lines linked to the upflow region on the eastern AR border. On the other hand, reconnection at these nulls could possibly result in upflows only in the neighborhood of the spines; then, this reconnection process can, at most, explain a fraction of the observed upflows (compare Figure 5 to Figure 6(d) using the QSL trace as a guide). Furthermore, reconnection at these nulls could only slightly affect the coronal loops since they are embedded in weak magnetic fields, which implies a low amount of energy release, and the structure of the pre- and post-reconnection loops is almost the same (except in the neighborhood of the nulls). Therefore, reconnection at these nulls, far away from the AR, is not expected to drive the upflows observed at the AR border.

Figure 6.

Figure 6. Evolution of the topology of AR 10978 for a set of magnetograms from 2007 December 10 to 13 (panels (a)–(e)). QSL traces related to the EIS upflows are marked with thicker black lines, while those not related are indicated with thinner ones. Most of these thinner black traces mark the drastic connectivity changes due to the presence of minor emerging bipoles in between the main AR positive and negative polarities. In panel (a), we have labeled the traces that we call outer and inner QSLs. The upflow (downflow) regions, for the Fe xii EIS emission line, are indicated as blue (red) shaded zones. The color bar below panel (d) indicates the velocity scale in km s−1. We have also added sets of computed field lines anchored in the vicinity of the outer QSL traces in panels (b)–(e) (see Section 4.2 and Figure 7 for their color convention). Field lines ending in a circle reach the borders of the box selected for the drawing. The convention for the magnetic field isocontours, in panels (a)–(e), and panel axes are the same as in Figure 5. Panel (f) depicts the line-of-sight magnetic field in grayscale with white (black) representing the positive (negative) field values. Two isocontours, with the same values as those in the other panels, have been added in gray continuous (dashed) line for the positive (negative) field. This panel has been added to better visualize the photospheric field distribution. The dates and times at the top of each panel correspond to those of EIS velocity maps. The axes are in Mm with the origin set in the AR.

Standard image High-resolution image

4. QSLS AND UPFLOWS

4.1. QSL characteristics

Based on our previous results (Baker et al. 2009; van Driel-Gesztelyi et al. 2012), we compute QSLs for a series of magnetic maps to confirm or refute the relationship between EIS upflow and QSLs during the AR evolution.

The method to compute QSLs was first described by Démoulin et al. (1996a). QSLs were defined using the norm, N, of the Jacobian matrix of the field-line mapping. This norm depends on the direction selected to compute the mapping; then, N has, in general, different values at both photospheric footpoints of a field line. To overcome this problem, Titov et al. (2002) proposed the introduction of a function called the squashing degree, Q, which is defined as N to the second power divided by the ratio of the vertical component of the photospheric field at the two opposite field-line footpoints. Q takes into account only the distortion of the field-line mapping, independently of the field strength, and is invariant along each field line.

To determine the QSL locations, we have to integrate a huge number of field lines. A key point is that a very precise integration method is used since derivatives of the mapping are needed to calculate N and Q. In order to decrease the computation time, we use an adaptive mesh, i.e., the mesh is iteratively refined only around the locations where the largest values of Q were found in the previous iteration. The fraction of points retained at each iteration controls the computation speed and how much the finally calculated Q-map will extend toward the lower values of Q. The iteration at a location is ended when the QSL is locally well resolved or, ultimately, when the limit of the integration precision is reached. Such computations can be performed at the photospheric level and also within the full coronal volume (Pariat & Démoulin 2012).

In complex magnetic configurations (e.g., quadrupolar or multipolar), the location of QSLs is strongly determined by the distribution of the magnetic field polarities at the photosphere (see Mandrini et al. 2014a, and references therein) and the maximum value of Q is typically very high (many orders of magnitude, up to infinity). Because the magnetic field configuration is more bipolar, the QSL locations become more influenced by the presence of magnetic shear and/or twist and the value of Q tends to be lower (see the discussion in Démoulin et al. 1997). AR 10978 is an isolated globally bipolar region within which several smaller bipoles emerged during its transit across the disk (see Figure 1 in Mandrini et al. 2014b and Figure 1 in this article); this increases the complexity of the QSL pattern and the Q values between the two main polarities.

4.2. QSL evolution

QSLs are the expected locations where magnetic reconnection can occur efficiently at coronal heights, releasing the stored magnetic energy (see references in Section 1). Then, the energy released is transported along field lines toward the chromosphere. Indeed, several studies have found that flare brightenings are located along chromospheric QSL traces (Démoulin 2007; Mandrini et al. 2010; Aulanier et al. 2011; Sun et al. 2013; Dalmasse et al. 2015; Vemareddy & Wiegelmann 2014; Savcheva et al. 2015).

In the case of EIS upflows, the relationship between QSLs and upflow regions is more indirect because the upflows are observed over a broad range of coronal heights. A direct comparison between upflow regions and QSL locations, considering the trace of reconnected field lines, is also complicated because the EIS velocity maps result from the integration of optically thin emission over a large depth along the line of sight; this integration is also affected by projection effects. However, taking into account the results from previous studies (Baker et al. 2009; van Driel-Gesztelyi et al. 2012), upflows are expected in the vicinity of QSLs where short- and high-density loops can reconnect with large-scale low-density closed or "open" loops. After reconnection, the higher pressure plasma from the initial short loops will be injected into newly formed large-scale loops creating a pressure gradient that drives the upflows (Bradshaw et al. 2011). This process is clearly illustrated in Figure 5 of van Driel-Gesztelyi et al. (2012).

Figure 6 shows the location of QSLs for a series of MDI magnetograms around AR 10978 CMP and illustrates all of the minor bipole emergences (compare to Figure 1). These magnetic maps are also the closest in time to the respective EIS velocity maps in Fe xii obtained during the AR disk transit. Extreme care was taken to coalign EIS data with MDI magnetograms. This was done not only using the image-header informations but also comparing the position of all structures (e.g., loop traces) with the location of the magnetic polarities. The QSL traces, black continuous thick and thin lines, have been overlaid on the photospheric magnetograms shown as isocontours of the field and the velocity maps indicated by blue (red) shaded regions corresponding to EIS upflows (downflows). The value of Q is above 400 for all of the traces shown; of course, at QSL locations where the null-point spines are anchored (see Figure 5) the value of Q is extremely high (Q ≥ 1012).

Figure 6 shows that the upflow regions are consistently located in the vicinity of QSLs drawn with thicker black lines on both AR main polarities. In panels (a) through (d), two main QSL traces, called outer and inner, extend along the negative western polarity. Sections of these traces (those drawn with thicker lines) are linked to the upflows. It is striking how well the inner QSL traces match the projected shape of the inner border of western upflow regions. Furthermore, the outer QSL shapes match the locations where the upflow velocities change in magnitude (see Démoulin et al. 2013, and Figure 3), indicating upflow regions with different characteristics. As the AR and its upflows evolve, see Figure 6(e), only the westernmost upflows associated to the outer QSL trace remain. On the eastern AR border, the outer QSL traces also match the projected shape of the upflows quite well. Next, as new minor bipoles emerge and evolve the shapes of the inner QSLs evolve as well and their complexity increases.

We have also added sets of field lines to all panels in Figure 6, except for panel (a) to avoid overloading it. These field lines have been computed starting integration in the vicinity of the outer QSLs, to the west (east) of the one on the main negative (positive) AR polarity. The projected shape of these field lines matches the spatial extension of the upflows at the borders of the AR, providing evidence of the close relationship between upflows and QSLs.

4.3. Detailed Connectivity Analysis on December 12

In this section, we present a scenario that, being consistent with the previously described observations and models, can explain why some sections of the QSL traces in AR 10978 are related to EIS upflows, while others are not and, also, why some upflow regions have different characteristics (as shown by Démoulin et al. 2013). As an example, we show in detail the field-line connectivity at a particular date and time. Because photospheric footpoint motions can provide forcing for the build-up of currents and induce reconnection along QSLs, we refer to the photospheric flows found in the AR (see Figure 2), which facilitate the sequence of reconnections described below.

Figure 7(a) shows two sets of field lines starting on the outer QSL trace on the negative main AR polarity. The northern field lines are anchored along the west side of the QSL trace and their opposite footpoints are located in quiet-Sun regions. These are large-scale loops probably filled with low-density plasma (low EUV emission). To the south of the same QSL trace, we have a set of shorter higher-density loops within the AR. This is the situation we envision before reconnection occurs.

Figure 7.

Figure 7. Magnetic field connectivity in the vicinity of the QSLs for the MDI magnetogram on 2007 December 12 at 3:26 UT. Panels (a) and (b) illustrate the pre-reconnection field lines in blue and the reconnected field lines in pink and red, respectively. These reconnected field lines would partially channel the EIS upflows located at the external eastern and western AR borders (Op and On; see also panel (d) in Figure 6). The corresponding EIS upflows are shown in the background of panel (b). The field lines that leave the drawn box correspond to large-scale loops connecting to positive quiet-Sun polarities (a circle is drawn at one of their ends). The dashed-line rectangle surrounds the region zoomed in panels (c) and (d). Panels (c) and (d) show the pre-reconnection and reconnected field lines using the same colors as in the top panels. These lines illustrate the reconnection process that would lead to the innermost EIS upflows. However, the red field lines in panel (d) do not follow the upflow spatial distribution; for this purpose, we have added a set of black lines in the vicinity of the western inner QSL (see Section 4.3 for an explanation). The convention for the magnetic field isocontours and axes is the same as in Figure 5. The spatial location of pre-reconnection and reconnected field-line footpoints on either side of the QSLs is typical of bipolar ARs (see examples in Démoulin et al. 1997).

Standard image High-resolution image

A large-scale photospheric flow pattern is present in the AR for several days, as shown for the period of December 10–12 (see Figures 2(b)–(d)). The main negative polarity is persistently moving to the southwest. As a result of this motion, the blue field lines anchored to the north of the outer QSL in Figure 7(a) may be forced to reconnect with the blue ones anchored to its south. After reconnection, we show the two sets of field lines in Figure 7(b). On one side, the process results in the long red field lines with footpoints at the western border to the south of the outer QSL. The injection of high pressure plasma from the shorter pre-reconnected loops could drive the observed upflows at the western south border of the AR. A similar process would drive, on the other side, the upflows on the eastern border of the main positive AR polarity along the reconnected pink field lines (Figure 7(b)). We remark that we can identify both pre- and post-reconnection field lines in the same magnetic field extrapolation because we have computed both QSLs and the driving photospheric flows. This is comparable to analyzing a snapshot in an MHD simulation knowing the velocity field direction. However, we have no way of using these observations to identify a pair of pre-reconnection field lines evolving into a pair of post-reconnection field lines.

The process discussed above would occur provided a large enough asymmetry in plasma pressure exists between the pre-reconnected loops. In this way, and for as long as a forcing is at work, we would observe the persistent EIS upflows in AR 10978. We speculate that reconnection would occur across the QSLs in the slipping mode analyzed by Aulanier et al. (2006) and further quantified by Janvier et al. (2013).

Comparison between upflow and QSL locations is not straightforward due to the changing viewing angle and the resulting projection effect, overlapping flows viewed in the optically thin corona, etc. However, if the reconnection process along QSLs lies at the origin of the EIS upflows, we expect that the projected spatial extension of the reconnected field lines at both AR borders compares to the spatial distribution of the observed Fe xii upflows. This is indeed the case for the pink reconnected field lines anchored in the main positive AR polarity (Op), but only partially for the red ones anchored in the main negative polarity (In; Figure 7(b)). However, field lines computed on the side and close to the QSL, thus in lower Q values such as the red field lines in Figure 6(d), have a projected extension that compares well with the extension of the upflow region at the western AR border, as do the pink ones in the same panel in relation to the upflows at the eastern AR border. Both of these sets would correspond to previously reconnected lines at the outer QSLs if we take into account that the QSL trace will shift to the east (west) on the main negative (positive) polarity as reconnection proceeds.

Though the upflow region to the west of the AR in Figure 7(b) looks in projection to be a single extended broad stream, it is in fact composed of at least two different streams (Figure 3). The easternmost border of this region lies in the vicinity of the inner QSL trace. In Figure 7(c), we show a set of blue field lines that have been computed starting integration from the eastern side of the inner QSL trace on the negative polarity. Magnetic reconnection between lines in this set and those that are anchored on the western border of the inner QSL on the positive polarity, also drawn in blue, would result in the red and pink lines shown in Figure 7(d). This reconnection process could be at the origin of the upflows located to the west of the inner QSL (at In). The projected extension of the reconnected red field lines does not match the apparent extension of the upflows; in fact, for this set of lines we would expect upflows to the east of this QSL trace. As done for the outer QSL, we integrated field lines anchored toward the west of the QSL trace. Such field lines, drawn in black in Figure 7(d), have a projected shape first directed to the west and then to the east. Plasma upflows along these lines would have the observed spatial distribution, i.e., we would observe stronger upflows close to the QSL trace fading toward the west. We also noticed that if we increased the absolute value of the α parameter in our LFFF model, keeping the same footpoints, the projected shape of the field lines matches the upflow extension better. Then, we attribute the discrepancy between upflows and field-line projected shapes to the limitation of our LFFF model that does not represent well the probably higher-sheared small-scale AR magnetic field.

Our results explain the origin of all the upflow regions, shown in Figure 6, for 2007 December 12 as being due to magnetic reconnection at QSLs. We have done similar detailed connectivity analyses for all of the upflows shown in that figure. They seem to originate either by reconnection within the outer QSLs, which are associated with the large-scale magnetic field configuration of the AR, or in the inner QSLs that develop as new bipoles emerge and evolve during the AR disk transit.

A two-step reconnection process was proposed by Mandrini et al. (2014b) to explain the way AR 10978 upflowing plasma could access open field lines and be observed by in situ instruments on board satellites at L1. In that article, we found evidence for the second step reconnection in the process, which opens the path for the closed-field confined plasma into the solar wind. Our present results offer proof of the first reconnection step, which is related to reconnection within the outer QSLs and provides a pathway for plasma originally confined along AR loops to flow into large-scale loops connecting to the quiet Sun. The latter is summarized in Figure 8. Furthermore, our findings of a coherent evolution between magnetic field, QSLs, and upflows prove the concept first proposed by Baker et al. (2009) to explain the origin of EIS persistent blueshifts. In this article, we go one step forward and suggest how, due to the photospheric field evolution, only sections of the QSLs in the AR are linked to persistent plasma upflows. Further evidence of a persistent reconnection process along QSLs in the AR can be revealed by radio observations of a persistent metric noise storm above AR 10978.

Figure 8.

Figure 8. Perspective view of the AR and the magnetic field connectivity sketching the main processes involved. Four field lines, corresponding to the four types of connectivity found in the coronal magnetic field model (Figure 7), are added on top of a magnetogram of the vertical field component. The photospheric flows shown in Figure 2, indicated by a yellow arrow, induce reconnection between over-pressured loops (dark-blue) located within the AR with under-pressured loops (light blue) connecting the AR to the quiet Sun by building up currents along QSLs in the corona. The pressure imbalance drives the upflows in the reconnected loops and the electrons accelerated during the quasi-steady sporadic reconnections generate a long-lived radio noise storm.

Standard image High-resolution image

5. MAGNETIC FIELD TOPOLOGY AND NOISE STORMS

5.1. Noise Storms: Characteristics

Following our study in Mandrini et al. (2014b), in this section, we search for evidence of energy release at QSLs in wavelengths typical of radio noise-storms. Noise-storms consist of a broadband continuum emission (bandwidth ∼200 MHz, around a central frequency of hundreds of MHz) lasting for hours to several days. Superimposed on this continuum, bursts of short duration with much smaller bandwidth (∼1 MHz) are observed. The continuum component exhibits a slow intensity variability. Mercier et al. (2014) analyzed radio observations in the range from 150 to 450 MHz with high spatial resolution and they concluded that noise storms show an internal structure with one or several compact cores embedded in a more extended and tenuous halo. The emission is due to the presence of a suprathermal electron population (energies ranging from one to a few tens of keV) injected and trapped in extended coronal structures, i.e., noise storms require a mechanism that quasi-continously accelerates electrons in the solar corona.

The onsets of noise-storms or their enhancements are often related to changes in the overlying corona (Kerdraon et al. 1983) and to energy release in the underlying AR (Raulin & Klein 1994; Crosby et al. 1996). These characteristics suggest that the plasma-magnetic field configuration is restructuring at the time and the place where the noise-storm is produced (see the next section).

5.2. Noise-storms: Evolution

We analyze the radio emission registered by the Nançay Radio Heliograph (Kerdraon et al. 1997) during the transit of AR 10978 across the solar disk. At higher observation frequencies (327, 408, 432 MHz), the emission indicates the presence of weak noise storms that remain almost unchanged during the whole observing period. Figure 9 illustrates the evolution of the noise storm above the AR. Radio emission contours at 80% of the maximum intensity are displayed over the nearest in time MDI magnetogram; these contours were built using 5 minute integrated data. They show no clear trend to strengthen (by increasing their sizes) or, conversely, to fade out (by diminishing their sizes).

Figure 9.

Figure 9. Radio noise storm sources over AR 10978 from 2007 December 8 to 12. The grayscale shows the magnetic field of the AR (MDI magnetograms). Superimposed, we show isocontours at 80% of the maximum intensity (5 minute integration) of the radio sources at 432 MHz (blue), 408 MHz (green), and 327 MHz (red). The selected radio observations show a typical distribution of the contours during the full observing period. The axes indicate the position on the Sun in arcseconds.

Standard image High-resolution image

At lower frequencies (150.9, 228 MHz), the radio emission contours embrace almost the full solar disk (not shown in Figure 9), indicating that the radiation corresponds to background coronal emission of thermal origin. Nançay Decameter Array (DAM, 10–80 MHz, Lecacheux 2000) did not observe Type III bursts, which are known to be caused by electrons that are accelerated outward along the open coronal field lines. DAM observes the upper corona from 0.7 to 3 R. The lack of Type III emission seems to indicate that accelerated particles were not injected into open coronal structures during the observing periods as the AR crossed the disk.

5.3. Noise Storms and Magnetic Field Connectivity

The global coronal magnetic field of CR 2064 was modeled in the PFSS approximation (Figure 10), assuming a current-free coronal field using a synoptic magnetogram as the photospheric boundary condition. To close the upper boundary, PFSS models assume that the field becomes purely radial at a given height, called the source surface. This is a free parameter usually set to the value 2.5 R. The PFSS model in this article was computed with the Finite Difference Iterative Potential-Field Solver code described by Tóth et al. (2011), using the corresponding MDI synoptic magnetogram for CR 2064 as the photospheric boundary condition.

Figure 10.

Figure 10. PFSS model of CR 2064 with AR 10978 located at the central solar meridian. The AR lies completely below the streamer belt and open field lines from two coronal holes are seen to the southeast and northwest of the AR. The field-line color convention is such that black indicates closed lines and pink (green) corresponds to open lines anchored in the negative (positive) polarity field. The magnetic field is saturated above (below) 30 G (−30 G). We have superimposed the NRH radio-emission isocontours at 432 MHz (blue) and 408 MHz (green) shown in Figure 9 on 2007 December 12. The radio emission is concentrated over the closed field lines anchored at the border of the AR and encloses it.

Standard image High-resolution image

As shown in Figure 10, the noise-storm radio emission is concentrated over the closed AR field lines with an extension to the south due to the projection effect. NRH isocontours fully enclose the AR. This spatial relation suggests that the same magnetic reconnection process that drives EIS upflows may accelerate the electrons that flow along the closed reconnected field lines originating from the radio noise-storm emission.

6. SUMMARY AND CONCLUSIONS

Since the discovery of ubiquituous plama upflows in EIS observations, several driving mechanisms were proposed. Among them were the impulsive heating at the footpoints of AR loops (Hara et al. 2008), "open" magnetic funnels explaining coronal plasma circulation (Marsch et al. 2008), chromospheric evaporation due to reconnection forced by flux emergence and/or braiding of field lines by photospheric motions (Del Zanna 2008), expansion of large-scale reconnecting loops (Harra et al. 2008), continual AR expansion (Murray et al. 2009), and, more recently, reconnection between over-pressure AR loops and neighboring under-pressure loops (Bradshaw et al. 2011). Baker et al. (2009) were the first to demonstrate the spatial relation between the location of upflows and QSLs at the border of a particular AR and to propose that magnetic reconnection at QSLs was at the origin of EIS upflowing plasma. van Driel-Gesztelyi et al. (2012) also found that EIS upflow regions at the border of an AR were cospatial with QSLs in another case study. However, an analysis of the spatial and temporal evolution of upflows and QSLs, which would provide strong support to the results found for individual examples, was still missing.

From an analysis of the evolution of the photospheric magnetic and velocity fields of AR 10978, as it transits the solar disk, combined with coronal magnetic field modeling and topology computations, we find that EIS upflow regions and QSLs evolve in parallel (Figure 6). Two sets of QSLs, called outer and inner (Figures 6 and 7), are found associated to EIS upflow regions with different characteristics (Démoulin et al. 2013). All of the EIS upflows in AR 10978 seem to originate either by reconnection within the outer QSLs, which are associated to the large-scale magnetic configuration of the AR, or within the inner QSLs that develop as new bipoles emerge and evolve within the AR during its disk transit (see an example in Figure 7). The reconnection process in sections of the outer QSLs is forced by a large-scale photospheric flow pattern which is present in the AR for several days. In our proposed scenario, which is summarized in Figure 8, upflows will be present provided a large enough asymmetry in plasma pressure exists between the pre-reconnected loops and for as long as a photospheric forcing is at work. A similar mechanism would be at work in sections of the inner QSLs, in this case forced by the emergence and evolution of the bipoles between the two main AR polarities. Furthermore, and within the limitations of our coronal field model, the projected extension of EIS upflows in AR 10978 match the projected shape of magnetic field lines computed in the vicinity of QSLs. Thus, our findings offer both observational and modeling support to the concept first put forward by Baker et al. (2009) and suggest how, due to the photospheric field evolution, only sections of the QSLs in the AR are linked to persistent plasma upflows.

Recent studies show that EIS upflowing plasma can gain access to open-field lines and be released into the slow solar wind via magnetic-interchange reconnection at magnetic null-points (see, e.g., van Driel-Gesztelyi et al. 2012). As shown in Figure 10, AR 10978 is completely covered by closed streamer field lines. Therefore, it seems unlikely that the upflowing plasma from AR 10978 can reach the solar wind; however, Culhane et al. (2014) found signatures of plasma with AR composition at 1 AU that apparently originated west of AR 10978. Based on a topological analysis of the global coronal magnetic field around AR 10978, Mandrini et al. (2014b) proposed a two-step reconnection process to explain the way a fraction of the AR upflowing plasma could access open field lines. Our present results describe the first reconnection step which occurs within the outer QSLs and brings the plasma originally confined along AR loops to flow into large-scale loops connecting to the quiet Sun. The second step reconnection in the process, modeled by Mandrini et al. (2014b), opens the path into the solar wind for the plasma confined into those large-scale loops. We also find further evidence of this first step in radio observations. Comparison of the large-scale global coronal field (Figure 10) to the location of a persistent metric noise-storm above AR 10978 observed by NRH, suggests that closed field lines reconnected within the outer QSLs may channel the accelerated electrons at the origin of the noise-storm.

A variety of magnetic configurations have been associated with EIS upflows located at AR borders. Some upflows have been related to magnetic reconnection at magnetic null-points and associated separatrices (e.g., Del Zanna et al. 2011), in other cases, no null points were present (e.g., Baker et al. 2009) or reconnection at null points could not explain the majority of EIS upflows (this article), while there are other examples in which reconnection explaining upflows could happen at both null points and QSLs (e.g.,van Driel-Gesztelyi et al. 2012). This is similar to what has been found for solar flares (e.g., Démoulin et al. 1994; Mandrini et al. 2010). Furthermore, solar flares can be either confined or associated with coronal mass ejections (CMEs). In the later case, they have a direct impact on the interplanetary space (e.g., Rouillard 2011). The same happens with upflows that can remain confined within coronal loops or become outflows and access the solar wind (e.g., van Driel-Gesztelyi et al. 2012; Culhane et al. 2014). On the other hand, solar flares are intrinsically impulsive events in which the magnetic energy, accumulated on the timescale of days, is released in a short time interval (typically, less than one hour) and the reconnection process is fast (Shibata & Magara 2011). Conversely, upflows involve magnetic reconnection on the timescale of days. Upflows are driven by the long-term slow evolution of the AR magnetic field (emergence, large-scale velocity patterns, diffusion), i.e., upflows are not driven by a global instability of the magnetic field like flares and CMEs, but rather by a gradual evolution of the magnetic field.

C.H.M., G.D.C., F.A.N., and A.M.V. acknowledge financial support from the Argentine grants PICT 2012-0973 (ANPCyT) and PIP 2012-0403 (CONICET). C.H.M. and G.D.C. are thankful for grant UBACYT 2013-0100321 (UBA). The research leading to these results has received funding from the European Commission's Seventh Framework Programme under grant agreement No. 284461 (eHEROES project). L.v.D.G. and D.B. acknowledge support from STFC Consolidated Grant ST/H00260/1. L.v.D.G.'s work was supported by the Hungarian Research grant OTKA K 109276. L.v.D.G. gratefully acknowledges support from the Leverhulme Trust under grant reference RPG-2014-051. This work was partially supported by a one-month invitation to P.D. to visit the Instituto de Astronomía y Física del Espacio and by a one-month invitation to C.H.M. to visit the Paris Observatory. C.H.M., G.D.C., and A.M.V. are members of the Carrera del Investigador Científico (CONICET). F.A.N. is a fellow of CONICET.

Please wait… references are loading.
10.1088/0004-637X/809/1/73