On Thermal-Pulse-Driven Plasma Flows in Coronal Funnels as Observed by the Hinode/EUV Imaging Spectrometer (EIS)

Using one-arcsecond-slit-scan observations from the Hinode/EUVImagingSpectrometer (EIS) on 5 February 2007, we find the plasma outflows in the open and expanding coronal funnels at the eastern boundary of AR 10940. The Doppler-velocity map of Fe xii 195.120 Å shows the diffuse closed-loop system to be mostly red-shifted. The open arches (funnels) at the eastern boundary of AR exhibit blue-shifts with a maximum speed of about 10 – 15 km s−1. This implies outflowing plasma through these magnetic structures. In support of these observations, we perform a 2D numerical simulation of the expanding coronal funnels by solving the set of ideal MHD equations in appropriate VAL-III C initial temperature conditions using the FLASH code. We implement a rarefied and hotter region at the footpoint of the model funnel, which results in the evolution of slow plasma perturbations propagating outward in the form of plasma flows. We conclude that the heating, which may result from magnetic reconnection, can trigger the observed plasma outflows in such coronal funnels. This can transport mass into the higher corona, giving rise to the formation of the nascent solar wind.


Introduction
The solar wind is the supersonic outflow of fully ionized gas from the solar corona streaming along the magnetic-field lines. It is well established that the polar coronal hole is the source of fast solar wind (Hassler et al., 1999;Wilhelm et al., 2000;Tu et al., 2005), while the slow solar wind originates from the boundary of active regions and along the streamers in the equatorial corona (Habbal et al., 1997;Sakao, Kano, and Narukage, 2007). Hassler et al. (1999) have found that the plasma can be supplied from the chromospheric heights in the network boundaries to the solar wind in polar coronal holes. More precise estimation of the formation of nascent wind in coronal funnels between 5 -20 Mm in coronal holes has been carried out by Tu et al. (2005). The bases of the polar coronal holes are mostly the sources of fast solar wind. It has also been suggested that plasma outflows observed at the edges of active regions are the source of the slow solar wind. Active-region arches, which can extend outward as the rays in the outer corona, may channel it (Slemzin et al., 2013). In addition to the large-scale origin, it has also been found that the small-scale outflows at the coronal hole boundaries (CHBs) can serve as the source of the slow solar wind (Subramanian, Madjarska, and Doyle, 2010). Recently, Yang et al. (2013) have presented a numerical model to describe the process of magnetic reconnection between moving magnetic features (MMFs) and the pre-existing ambient magnetic field that drives an anemone jet with inverted y-shape base and associated plasma blobs. They have found that an increase in the thermal pressure at the base of the jet is also driven by the reconnection, which induces a train of slow-mode shocks propagating upward resulting in plasma upflows. Their findings contribute to the formation of jets and small-scale flows in the quiet-Sun corona where MMFs undergo into the low atmospheric reconnection.
Two outstanding issues, however, remain unsettled: i) what are the drivers of these winds in the outer corona? and ii) what are the source regions and what drivers enable the mass supply to the lower solar atmosphere (chromosphere-TR, and inner corona)? There exist several studies advocating the role of Alfvén waves as a possible answer to the first question. The ion-cyclotron waves at kinetic scales are visualized as one of the possible candidates to provide momentum and heat the outer coronal winds both in theory and observations (e.g. Ofman and Davila, 1995;Tu and Marsch, 1997;Tu et al., 2005;Suzuki and Inutsuka, 2005;Dwivedi and Srivastava, 2006;Jian et al., 2009, and references cited therein). The second question is crucial at present. The consensus of solar-wind research, however, is unclear as to the origin of the mass supply to the supersonic wind. Tian et al. (2010) have found upflows in the open-field lines of coronal holes starting in the solar transition region and interpreted this as an evidence of the fast solar wind in the polar coronal holes. As far as the slow solar wind is concerned, the outflows at the boundaries of active regions can contribute at larger spatio-temporal scales to the mass supply as an expansion of the loops lying over these active regions (Harra et al., 2008). It has been shown recently that collimated jet eruptions can also contribute to the formation of the solar wind . The magnetic-field topology of structures such as jets (e.g. spray surges) may not contribute to the solar wind (Uddin et al., 2012). The contributions from the confined ejecta in the solarwind formation depend mainly on the local magnetic-field topology and plasma conditions. The question remains on whether the mass supply to the slow solar wind comes from the lower atmosphere expanding along the curved coronal fields. The question is what are the potential physical drivers? It has been found that the open-field lines at the boundary of active regions reconnect periodically with closed-field lines to guide the plasma motion in the form of solar wind (Harra et al., 2008). Similar examples are reported at the boundary of the coronal holes at small spatial scales as the source of slow solar wind (Subramanian, Madjarska, and Doyle, 2010). Therefore, an alternative option may be magnetic reconnection as a potential mechanism in the formation of slow solar wind. Apart from magnetic reconnection, the wave-heating scenario can shed light on the slow solar-wind source regions. Schmidt and Ofman (2011) have found that the energy stored in the slow magnetoacoustic waves propagating toward the higher atmosphere within expanding loops. This may be a potential candidate for the acceleration and formation of the slow solar wind. Wave activity at the bases of the fast (polar coronal holes) and slow (equatorial corona) solar wind can be important to power the energized plasma at greater heights up to the corona where it can be triggered supersonically in interplanetary space (e.g. Harrison, Hood, and Pike, 2002;Dwivedi andSrivastava 2006, 2008;De Pontieu et al., 2007;McIntosh et al., 2011;McIntosh, 2012, and references cited therein). Thus, there seems to be compelling evidence for the role of magnetic reconnection and wave phenomena in the solar-wind source region.
In the present article, we report the evidence of the outflowing magnetic arches acting as coronal funnels at the eastern boundary of an AR 10940 loop system observed on 5 February 2007. These coronal funnels seem to open up in the higher atmosphere to transport the outflowing plasma. Their footpoints are rooted in the boundary of the active region. They are the most likely heated regions that result in the activation of the outflowing plasma. We present a 2D MHD simulation of the open and expanding-funnel-type model atmosphere in which a rarefied and hot region is implemented near the footpoint that exhibits plasma perturbations similar to the observations. In Section 2 we summarize the observational results. We present the numerical model in Section 3. Discussion and conclusions are given in the last section.

Observational Results
The active region AR 10940 was observed by a one-arcsecond-slit scan of EUV Imaging Spectrometer (EIS: Culhane et al., 2006) onboard the Hinode spacecraft on 5 February 2007. The EIS is an imaging spectrometer of which 40-and 266-arcsecond slots are used for the image analyses using the light curves and emission per pixel. The one-and twoarcsecond slits are utilized for spectral and Doppler analyses using spectral-line profiles. EIS observes in two modes: i) scan; ii) sit-n-stare. The EIS observes high-resolution spectra in two wavelength intervals: 170 -211 Å and 246 -292 Å using, respectively, its Short-Wavelength (SW) and Long-Wavelength (LW) CCDs. The spectral resolution of the EIS is 0.0223 Å per pixel. The analyzed observations were taken on 5 February 2007 and the data-set contains spectra of various lines formed at chromospheric, transition region (TR), and coronal temperatures. The scanning observation started at 12:14:12 UT and ended at 13:31:21 UT on 5 February 2007. The scanning steps were without any offset in the region containing a coronal active region and its eastern boundary with open and expanding magnetic arches where we are interested in the present investigation. This provides us with an opportunity to understand the plasma activity along the open-field regions at the eastern boundary of the active region in between the diffused loop systems that reach up to a higher height in the corona (Figure 1). We refer to these structures as "coronal funnels". Such flow regions are heated at their base exhibiting outflows. To understand the approximate magnetic-field geometry associated with AR 10940 and its surrounding region, we perform a potential field source surface (PFSS) extrapolation. Fan lines over the Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) observations are shown in Figure 2. It is clear that the core loops of the active region are bipolar and connect the central opposite magnetic polarities (white loops shown by the black arrow). Another set of closed-field lines connect to the east-side weak negative polarity and the central positive  In order to obtain the velocity structures in such coronal funnels, we select the strongest EIS line [Fe XII 195.12 Å] in our study. We aim for the understanding of the impulsively generated plasma outflows in such funnels and associated physical processes. The slit started the scanning of the polar coronal hole with (X cen , Y cen ) ≈ (799.087 arcsecond, −19.185333 arcsecond). The observational windows acquired on the CCDs are 128 pixels high along the slit, and 111 pixels wide in the horizontal direction where spectra have a spectral resolution of 0.0222 Å per pixel.
We apply the standard EIS data-reduction procedures and calibration files/routines to the data obtained from the EUV-telescope, which is the raw (zeroth-level) data. The subroutines are found in the sswidl software tree under the IDL environment (www.darts.isas.jaxa.jp/pub/ solar/ssw/hinode/eis/). These standard subroutines reduce the dark-current subtraction, cosmic-ray removal, flat-field correction, hot pixels, warm pixels, and bad/missing pixels. The data are saved in the level-1 data file, and the associated errors are saved in the error file. We choose the clean and strong line Fe XII 195.12 Å to examine the spatial variations of the intensity and Doppler velocity in the observations scanned over AR 10940 on 5 February 2007. We co-align the Fe XII 195.12 Å map with respect to the long-wavelength CCD observations of He II 256.86 Å by considering it as a reference image and by estimating its offset. The orbital and slit-tilt are also corrected in the data using the standard method described in the EIS-software notes. We perform a double Gaussian fit for the removal of the weak blending of Fe XII 195.18 Å line as per the procedure described by Young et al. (2009), which is also outlined in the EIS Software Note 17. We constrain the weakly blended line 195.18 Å to have the same width as the line at 195.12 Å and to fix an offset of +0.06 Å relative to it. It is to be noted that we perform the fitting on 2 × 2 pixel 2 binned data to enhance the signal-to-noise ratio and also to be obtain a reasonable fit. Figure 1 displays the intensity (left) and Doppler-velocity (right) maps of the observed AR and open arches at its eastern boundary. The intensity map shows that the AR is made up of the core diffused loop systems lying at the lower heights in the plane perpendicular to the line-of-sight. In actuality, however, they may be tilted. These low-lying core loops exhibit high emission and downflows. At the eastern boundary, we have identified the open and expanding arches that exhibit coronal-funnel behavior to transport mass and energy into the higher corona. These expanding funnels may be the part of more quiescent and largescale loop systems opening higher up in the corona. It is seen in the enlarged Doppler map (see Figure 3) that four such regions are at least identified as blue-shifted and outflowing regions. As already noted, these could be the legs of higher and large-scale loops, which form the building blocks of the coronal funnels at relatively lower heights of the corona. They exhibit the slow and subsonic plasma outflows mostly spreading along their field lines with a maximum speed of about 10 -15 km s −1 . Figure 3 shows the selection of path 1 -4 (top panel). The first two paths are drawn over the large-scale open and magnetic-field arches (lines) extending towards QS -CH region in the North-East of this AR (not shown here). These two regions serve as the open and expanding funnel regions. Paths 3 and 4 are drawn to the expanding and blue-shifted regions that can also serve as coronal funnels. However, they are actually the lower parts of some core-loop systems. On this account, we have selected two coronal funnels (1 and 2), while the others (3 and 4) serve as funnels but are probably associated at larger heights with the curved-loop magnetic field. The middle and bottom panels (left to right), respectively, show the variation of projected L.O.S. Doppler velocity along paths 1 -4 with height in these funnels. For funnels 1 and 2, the footpoints exhibit stronger outflows with a maximum speed of 15 km s −1 . At heights beyond 10 Mm, the outflow speed weakens in these funnels. This signifies the start of the outflows due to heating near the footpoints of these funnels; the outflows weaken with the height as we move away from the heating source. The other two funnels (3 and 4) are likely the legs of the diffused core loops, exhibiting the outflows at certain lower heights with a maximum speed of 10 km s −1 . This indicates different locations of the heating. It is also noted that the outflows diminish at lower heights compared to the corresponding ones in funnels 1 and 2 in the form of open arches. This is because of their association with the curved loops at greater heights where plasma is trapped and flows downward in order to maintain a new equilibrium. It is also noted that the flow structure is more gentle in the open funnels 1 and 2. However, it is greater at the footpoint and decreases with height. In funnels 3 and 4, which are the lower parts of curved loops, the generation of the flow is rather impulsive. The outflow starts at a certain height above the loop's footpoint. It increases up to a certain distance and decreases thereafter. This shows that impulsive heating is at work near the loop footpoint, which causes enhanced upflows up to a certain distance. The downflowing plasma from the upper part of the loop may counteract with the upflows. The red-shifted apex of the core-loop system is clearly evident in Figures 1 and 3. Funnels 3 and 4 are the lower parts of this loop system.
In the next section, we outline the details of the 2D numerical simulation of such observed expanding coronal funnels and their plasma dynamics. We solve a set of ideal MHD equations in the appropriate VAL-III C initial temperature conditions and model atmosphere using the FLASH code.

Model Equations
Our model system starts from a gravitationally stratified solar atmosphere, which can be described by the ideal two-dimensional (2D) MHD equations: Here , V , B, p = k B m T , T , γ = 5/3, g = (0, 0, −g) with its value g = 274 m s −2 , m, and k B are the mass density, flow velocity, magnetic field, gas pressure, temperature, adiabatic index, gravitational acceleration, mean particle mass, and Boltzmann's constant, respectively. It should be noted that we do not consider radiative cooling and thermal conduction in our present model for the sake of simplicity. We simulate only the dynamics of the plasma outflows to compare them with that of the observed coronal funnels.

Initial Conditions
We assume that the solar atmosphere is at rest [V e = 0] in equilibrium with a current-free magnetic field [∇ × B e = 0]. As a result, the magnetic field is force-free: The divergence-free constraint is satisfied automatically if the magnetic field is specified by the magnetic-flux function [A(x, y)] as B e = ∇ × (Aẑ).
Here the subscript "e" corresponds to equilibrium quantities. We set a curved magnetic field by choosing where S is the strength of a magnetic pole and (a, b) = (0, −10) Mm is its position. For such a choice of (a, b), the magnetic-field vectors are weakly curved to represent the expanding coronal funnels. As a result of Equation (5) the pressure gradient is balanced by the gravity force: −∇p e + e g = 0.
With the ideal-gas law and the y-component of Equation (8) where is the pressure scale height, and p 0 denotes the gas pressure at the reference level that we choose in the solar corona at y r = 10 Mm. We take an equilibrium temperature profile [T e (y)] (see top-left panel in Figure 4) for the solar atmosphere derived from the VAL-C atmospheric model of Vernazza, Avrett, and Loeser (1981) that is smoothly extended to the solar corona.
The transition region is located at y ≈ 2.7 Mm. There is an extended solar corona above the transition region, with the temperature minimum level located at y ≈ 0.9 Mm below the solar chromosphere. Having specified T e (y) (see Figure 4, left-top panel) with Equation (9), we obtain the corresponding gas-pressure and mass-density profiles.

Numerical Scheme and Computational Grid
Equations (1) -(4) are solved numerically using the FLASH code (Lee and Deane, 2009;Lee, 2013). This code uses a second-order, unsplit Godunov solver (Godunov, 1959) with various slope limiters and Riemann solvers, as well as Adaptive Mesh Refinement (AMR) Lee 2011, 2012). We set the simulation box of (−7, 7) Mm × (2.6, 51.6) Mm along the x-and y-directions ( Figure 5). The lower boundary, where we apply the heating pulse, is set at x 0 = 0 Mm, y 0 = 2.6 Mm. We set and hold fixed all of the plasma quantities at all boundaries of the simulation region to their equilibrium values, which are given by Equations (6) and (9). As the magnetic field is curved and plasma is stratified gravitationally, open boundaries would not be a perfect choice with regard to the fixed boundary conditions. We have verified this experimentally. As the FLASH code uses a third-order accurate Godunov-type method, a characteristic method is already built into this procedure. However, the Riemann problem at the boundaries corresponds to open boundaries, resulting in numerically induced reflections from these boundaries. On the other hand, the fixed boundaries lead to negligibly small numerical reflections. Therefore, we adopt this approach in the numerical simulations. In addition, we modify the equilibrium mass density and gas pressure at the bottom boundary as where and Here A and A p are the amplitudes of the perturbations, (x 0 , y 0 ) are their initial position, and w x , w y denote its widths along the x-and y-directions, respectively. The symbol τ denotes the growth time of these perturbations. We set and hold fixed A = 7, A p = 2, x 0 = 0 Mm, y 0 = 2.6 Mm, w x = 3 Mm, w y = 0.5 Mm, and τ = 10 seconds. These fixed boundary conditions perform much better than transparent boundaries, leading to only negligibly small numerical reflections of the wave signals from these boundaries.
In our modeling, we use an AMR grid with a minimum (maximum) level of refinement set to 3 (8) (see Figure 5). The refinement strategy is based on controlling numerical errors in mass density, which results in an excellent resolution of steep spatial profiles and greatly reduces numerical diffusion at these locations.
A standard procedure to check the magnitude of a numerically induced flow is to run the code for the equilibrium alone, without implementing any perturbation. We have verified that the numerically induced flow is of the order of 2 km s −1 in the solar corona, and the transition region is not affected by the spatial resolution which is about 20 km around the transition region. This resolution is much smaller than the width of the transition region (≈ 200 km) and the pressure scale height near the bottom, which is about 0.5 -1 Mm. We have also performed grid-convergence studies by increasing the spatial resolution by a factor of two at the transition region. As the numerical results have been found essentially similar results for the finer and coarser grid, we have limited our analysis to the latter.

Results of the Numerical Simulation and Comparison with Observations
The results of the numerical simulations and their comparison with the observed plasma outflows are summarized as follows: Figure 6 (first column) displays the velocity vectors plotted over the total velocity maps for t = 50, 100,150,200,250,300,350,400,450 seconds. The diverging magneticfield lines of the model funnel is also over-plotted in these snapshots. It is clear from the t = 50 seconds snapshot that the implemented heating at the footpoint, just below the transition region, results in the alteration of the ambient plasma pressure, and the plasma starts flowing upward with a typical average speed of 40 km s −1 near the footpoint at this time. The outflow velocity is maximum near the heated region at the footpoint of the model funnel. This is higher than the observed line-of-sight outflow velocities (10 -16 km s −1 ) in various coronal funnels above their footpoints (see Figure 3). The simulated outflow velocities greatly depend on the initial conditions of the model funnel and magnitude of the heating pulse. The observed trend of the flows in funnels 1 and 2 (see Figure 3) and the same derived from the model match each other (see Figures 6 and 7). At higher altitudes, above the heating location in these funnels, the outflow weakens which qualitatively matches with the results obtained from our model (see Figures 6 and 7). Our model fits better the plasmaflow conditions in the open coronal funnels (e.g. funnels 1 and 2) where the gentle flows start and expand due to their footpoint heating. The flow becomes steady at each height after 300 seconds in the model funnel when the heated plasma reaches the height of 10 Mm. This also indicates that the heating pulse is at work for a certain duration, and after some time the generated flows reach a new equilibrium. The physical scenario is in agreement with the open coronal funnels 1 and 2 (see Figures 3 and 7). Comparison of these two scenarios supports the outflow of the plasma due to heating. The heating causes thermal flows guided in the magnetic-field lines of the open coronal funnels, while it subsides away from the source. The physical behavior of the velocity field and its spatial distribution as observed by Hinode/EIS along each funnel are consistent with the velocity field at a particular temporal span of the numerical simulation when the outflowing plasma rises to a maximum height of ≈ 10 Mm. For the funnels that are the lower parts of the curved-loop system (funnels 3 and 4; Figure 3), the energy release might occur where the blue-shift is enhanced at a certain height above the footpoint impulsively. Therefore, the locations in those observed funnels are identical to the energy-release site of the modeled funnel where outflows start as a result of heating. On the contrary, it seems that impulsive heating in funnels 3 and 4 causes the increment in the outflows up to a certain height and thereafter is balanced by the downflowing and trapped plasma from the loop apex.  Figure 6, second and third columns, respectively, display the temperature and density maps for t = 50, 100, 200, and 300 seconds. It is clear that during the heating, the plasma maintained at inner coronal/TR temperatures (sub-MK and 1.0 MK; the hot plasma envelops the cool one) and with somewhat higher density, starts flowing from the footpoint of the model funnel toward greater heights. The plasma is denser near the footpoint. It flows along the open funnels toward higher heights. Thermal perturbations created the slow and subsonic flows of the plasma; therefore, it reaches up to only the lower heights. This is a situation similar to that of the transition region within the funnel being pushed upward due to the evolution of the thermal perturbation underneath. Sub-mega-Kelvin and denser plasma from lower solar atmosphere move up and is enveloped by the hot coronal plasma. The denser plasma maintained at the TR and inner coronal temperatures is visible up to a height of 10 Mm in the model funnel. This is consistent with the observations. Only the lower parts of the coronal arches are associated with enhanced fluxes and thus with higher densities, while their higher parts are less intense and denser regions (see Figure 1, left panel).

Discussion and Conclusions
We have presented observations of outflowing coronal arches (coronal funnels 1 and 2) lying at the eastern boundary of the AR 10940 observed on 5 February 2007, 12:15 -13:31 UT. The scanning observations show that these arches open up in the nearby quiet-Sun corona and exhibit plasma outflows maintained at coronal/TR temperature of around 1.0 MK. They serve as the expanding coronal funnels from which the plasma moves to the higher corona and serves as a source of the slow solar wind (Harra et al., 2008). The plasma outflows may be generated in such open-field regions because of the low atmospheric reconnections between the open-and closed-field lines (Subramanian, Madjarska, and Doyle, 2010). Episodic heating mechanisms are now well observed and interpreted as the drivers of outflowing plasma in the curved coronal loops (Klimchuk, 2006;Del Zanna, 2008;Brooks and Warren, 2009). The steady heating may generate the hot-plasma upflows in the corona (Tripathi et al., 2012).
There has so far been little effort made to model the plasma outflows in the corona and that too in an entirely different context of the large-scale evolution of coronal magnetic field. Murray et al. (2010) have modeled in 3D the origin and driver of coronal outflows, and found that outflows are the result of the expansion of an active region during its development. Harra et al. (2012) have modeled the AR coronal outflows as a consequence of compression during the creation and annihilation of the magnetic-field lines.
However, our objective of the present investigation is to implement and test our 2D numerical-simulation model of the localized expanding coronal funnels with Hinode/EIS observations of the outflowing open field arches (i.e. funnels) at the boundary of AR 10940. We solve the set of ideal MHD equations in the appropriate VAL-III C initial temperature conditions and model atmosphere using the FLASH code. The key ingredient of our model is the implementation of realistic ambient solar atmosphere, e.g. realistic temperature, presence of the TR, implementation of expanding coronal fields, stratified atmosphere in the initial equilibrium, which significantly affect such plasma dynamics. We have implemented a rarefied and hotter region at the footpoint of the model funnel that triggers the evolution of the slow and subsonic plasma perturbations propagating outward in the form of plasma flows similar to the observed dynamics. We implemented the localized heating below the transition region at the footpoint of the funnel. The plasma is considered rarefied in the horizontal direction, mimicking the structured, open, and expanding magnetic funnels. The heated plasma evolves and exhibits the plasma perturbations as has been observed as outlined in the Hinode/EIS observations. The outflows start at the base of the funnel that further weakens with the height, which is suggestive of the plasma dynamics due to heating near the footpoint. A similar physical scenario is observed in the selected outflowing magnetic arches mimicking the coronal funnels 1 and 2 at the eastern boundary of AR 10940.
We conclude that the implemented episodic heating can excite plasma outflows in expanding coronal funnels residing at boundary of solar active regions as well as in the quiet-Sun. These slow, subsonic plasma outflows may not be launched at higher altitudes in the corona. We have examined the presence of hot and denser plasma up to 10 Mm height as triggered by thermal perturbations in our model. Observations also show that the plasma only rises significantly in the lower parts of these funnels up to inner coronal heights with significant intensity and velocity distributions. However, even if such flows reach up to the inner coronal heights of 10 Mm in these funnels, they may also contribute to the mass supply to the slow solar wind. Therefore, our model and observations invoke the dynamics of the plasma in the localized coronal funnels, which may be important candidates to transport mass and energy into the inner corona. However, more observational studies should be performed with new spectroscopic data (e.g. Interface Region Imaging Spectrograph (IRIS)) to compare with our proposed 2D model, specifically under the physical conditions of different types of localized flux tubes in the solar atmosphere, which can serve as plasma outflowing regions due to episodic heating.