[go: up one dir, main page]

A publishing partnership

Articles

KELVIN–HELMHOLTZ INSTABILITY OF A CORONAL STREAMER

, , and

Published 2013 August 26 © 2013. The American Astronomical Society. All rights reserved.
, , Citation L. Feng et al 2013 ApJ 774 141 DOI 10.1088/0004-637X/774/2/141

0004-637X/774/2/141

ABSTRACT

Shear-flow-driven instability can play an important role in energy transfer processes in coronal plasma. We present for the first time the observation of a kink-like oscillation of a streamer that is probably caused by the streaming kink-mode Kelvin–Helmholtz instability (KHI). The wave-like behavior of the streamer was observed by the Large Angle and Spectrometric Coronagraph Experiment C2 and C3 on board the SOlar and Heliospheric Observatory. The observed wave had a period of about 70–80 minutes, and its wavelength increased from 2 R to 3 R in about 1.5 hr. The phase speeds of its crests and troughs decreased from 406 ± 20 to 356 ± 31 km s−1 during the event. Within the same heliocentric range, the wave amplitude also appeared to increase with time. We attribute the phenomena to the MHD KHI, which occurs at a neutral sheet in a fluid wake. The free energy driving the instability is supplied by the sheared flow and sheared magnetic field across the streamer plane. The plasma properties of the local environment of the streamer were estimated from the phase speed and instability threshold criteria.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The Kelvin–Helmholtz instability (KHI), driven by a shear of the flow velocity, has been observed in various astrophysical environments. In the solar atmosphere, Berger et al. (2010) and Ryutova et al. (2010) discovered it in a prominence from Hinode Solar Optical Telescope observations. Foullon et al. (2011) and Möstl et al. (2013) found that it could develop at the flank of a fast coronal mass ejection (CME). Ofman & Thompson (2011) identified a vortex-shaped feature caused by the KHI on the boundary of a dimming area during the eruption of an active region. It has been shown by Chandrasekhar (1961) that the presence of a magnetic field tangential to the boundary surface and directed along the flow stabilizes the hydrodynamic KHI. For a symmetric shear layer, the instability can only grow if the velocity jump across the boundary between the two uniform regions exceeds twice the Alfvén velocity on either side.

The instability growth is modified if the field and density profiles across the boundary are asymmetric or if the shear layer has a finite width. Instability calculations and simulations have been performed to model the magnetospheric current sheet, astrophysical jets, and streamer sheets in the solar corona. The models adopted for these cases were two- and three-dimensional force-free or pressure-balanced current sheets with a jet or a wake layer centered at the current sheet (Lee et al. 1988; Wang et al. 1988; Dahlburg et al. 1998; Einaudi 1999; Dahlburg et al. 2001; Zaliznyak et al. 2003). Due to Galilean invariance, the wake, a layer of reduced speed embedded in a uniform background flow, is equivalent to the jet case for otherwise identical plasma parameters and produces the same instability growth rates but shifted phase speeds (Einaudi et al. 1999). More recently, 2.5 dimensional and three-dimensional simulations were performed to shed more light on the velocity and magnetic shear driven instabilities in the nonlinear regimes (Bettarini et al. 2006, 2009).

From all of these studies, three different modes could be distinguished during the initial phase of the instability. For low Alfvénic Mach numbers below about two, the instability is strongly influenced by the magnetic field and shows plasmoid formation as a result of the tearing mode instability (or resistive varicose mode; e.g., Einaudi et al. 1999). For Alfvénic Mach numbers above about two, the hydrodynamic behavior overwhelms the magnetohydrodynamic evolution. A symmetric oscillatory mode and an asymmetric oscillatory mode evolve in such a situation and are designated as kink (or sinuous) mode and sausage (or ideal varicose) mode.

Using linear theory, Lee et al. (1988) derived analytical expressions for the growth rates and phase speeds of the three modes of a jet, assuming a four-layer model with constant plasma parameters for each. They found that the kink mode dominates over the sausage mode when the value of β (ratio of thermal to magnetic pressure) in the ambient plasma exceeds unity. This value of β is, however, linked to the other equilibrium parameters of the current sheet by the requirement that the net pressure balance be maintained, e.g., the ratio of the density between the current sheet and its ambient plasma. In an accompanying paper, Wang et al. (1988) applied the wake model to the heliospheric current sheet and considered the influence of the difference between the width of the velocity shear layer and the width of the current sheet on the growth rate. Numerical simulations of the current carrying wake with various plasma parameters suggest that the kink mode often dominates for the large β parameters considered, except for long wavelengths of the perturbations (Dahlburg et al. 1998, 2001; Zaliznyak et al. 2003). Even though the instability growth rate decreases steadily with increasing sonic Mach number, i.e., with decreasing β, the simulations do not show a drastic threshold for β or for the sonic Mach number as exists for the Alfvénic Mach number.

The quasi two-dimensional current plane at the tip of a helmet streamer seems to become particularly close to the model geometry adopted in the above simulations. However, contrary to the ease with which the kink mode is excited in the simulations, it is not observed as frequently. Indeed, SOHO/LASCO observations show the almost continuous ejection of plasmoids from the tips of helmet streamers, (e.g., Sheeley et al. 1997). They were interpreted as the development of the tearing-type reconnection into the nonlinear varicose (sausage-like) fluid instability discussed above (Einaudi et al. 1999, 2001). On the other hand, to our knowledge, the equivalent kink-like instability along a streamer has not been observed.

In this paper, we report the rare observation of kink-like oscillations in the current sheet of a coronal streamer. On 2011 June 3, a narrow streamer was emerging from behind the occulter of SOHO/LASCO C2 and C3, which after about 7:30 exhibited a clear oscillating behavior.

Chen et al. (2010) also observed a wave-like motion of a streamer which they claim was excited by the eruption of a fast CME nearby. In their observations, the wave was excited impulsively and steadily decayed afterward. In the case of the streamer on June 3 reported here, no CME occurred close to the streamer in either time or space. Also, the wave amplitude in this event stayed constant for hours or even slightly increased with time. We therefore attribute this oscillation to an instability rather than to an impulsive excitation. To our knowledge, this is the first observation of a streamer wave that was probably caused by the streaming kink instability. In Section 2, the observations of the streamer wave by the Large Angle and Spectrometric Coronagraph Experiment (LASCO) coronagraph and the retrieval of the wave shape are described. In Section 3, we compare the potential of the streaming kink instability and other competing MHD modes of the KHI as sources of this wave. In the last section, the obtained results will be summarized and discussed.

2. OBSERVATIONS AND DATA REDUCTION

The oscillating feature was observed by the LASCO (Brueckner et al. 1995) C2 and C3 coronagraphs on board the Solar and Heliospheric Observatory (SOHO). LASCO C2 has a field of view (FOV) from 2 to 6 R and C3 from 3.7 to 30 R. The wave-like motion of the streamer first appeared in the C2 FOV at a distance of 3–4 R at 05:48 UT on 2011 June 3. The wave which developed then propagated almost radially to a distance of about 11 R at 10:30 UT. At larger distances, the phenomenon becomes impossible to observe due to the low signal-to-noise ratio. A movie showing the streamer wave can be found in the online journal. The end of the wave-like motion was accompanied by a disruptive eruption beginning around 13:25 UT. In Figure 1, a snapshot of the oscillation is presented. The streamer with the wave-like feature is located in the first quadrant of the image plane and is rooted at a latitude around 27fdg5.

Figure 1.

Figure 1. Oscillating streamer observed by C2 at 09:12 UT and by C3 at 09:06 UT. The white circle indicates the solar limb, and the plus sign represents the Sun's center. A movie showing the evolution of the oscillating steamer is available online. The units are in solar radius.(An animation of this figure is available in the online journal.)

Standard image High-resolution image

In the LASCO images, the streaming feature appears as a white-light jet (Wang et al. 1998; Feng et al. 2012). Usually, a jet can be traced down to its footpoints at the solar surface. However, AIA/SDO did not observe any possible source for the jet at a lower altitude. Moreover, an elongated structure existed at the same position for many hours before the appearance of the wave. An alternative possibility is the leakage of plasma from the cusped magnetic field lines above a helmet streamer. The cusp is a hot region at the streamer top with low magnetic field strength where the magnetic boundary is concavely bent outward. In this region, the plasma β may easily have values above unity. Thus, the magnetic field confinement can be considered to be marginal and the streamer plasma may occasionally leak out of the cusp region to feed the slow solar wind. It is generally observed that the cusp altitude is below 2–2.5 R, which is very close to the inner boundary of the LASCO C2 FOV (Chen et al. 2010). In EUV observations, due to the relatively low density of the helmet streamer plasma and the restricted FOV, it is usually difficult to see the dome part of a streamer. For that reason, the Atmospheric Imaging Assembly (AIA) was probably not able to see the lower helmet counterpart associated with the observed streamer.

To clarify the spatial context of this event, we show an EUV image observed at λ = 28.4 nm by the leading spacecraft of the Solar Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008) in the upper panel of Figure 2. At the time of observation, STEREO A was located about 90° west of the Sun–Earth line and looks almost vertically down onto the longitude of the streamer. Two coronal holes, one above and one below the streamer latitude, are indicated.

Figure 2.

Figure 2. Upper panel: EUV image observed at 28.4 nm by STEREO A/EUVI. On 2011 June 3, STEREO A was located about 90° west of the Earth. Two coronal holes above and below the latitude of the wave-like feature are indicated. Lower panel: extrapolated coronal field lines from the HMI synoptic map for Carrington rotation 2110 using the PFSS model. The red lines indicate inward open field lines, and the green lines represent outward open field lines. The spherical source surface is set at 1.5 R.

Standard image High-resolution image

The lower panel of Figure 2 shows the field lines from a PFSS model emerging from a 20° wide longitude range around the limb and from latitudes of 27fdg5 ± 30°. The field extrapolation was produced from a synodic magnetogram observed by the Helioseismic and Magnetic Imager (HMI; Schou et al. 2012) on board the Solar Dynamic Observatory (SDO) for Carrington rotation 2110. In this magnetogram, the limb longitude is intersected by two neutral lines which, for a source surface height rSS < 1.6 R, result in two current sheets above the limb. The southern sheet very closely matches the latitude of the observed oscillating streamer. For larger rSS, the two current sheets merge and the PFSS field beyond the source surface height becomes unipolar. In this case, the observed streamer would have to be interpreted as a pseudostreamer (Wang et al. 2007). From the observations in Figure 1, we can state that the source surface height rSS for the oscillating streamer must be significantly less than the LASCO C2 occulter radius of about 2.2 R. The fact that a value of rSS < 1.6 R produces a current sheet within a few degrees of the observed streamer latitude gives us confidence that the observed streamer is associated with a current sheet and the streamer bulge is confined to low altitudes.

In Figure 3, we present some selected image frames that zoom in on the streamer wave. To make the streamer structure more prominent, a pre-event image immediately before the streamer's appearance was subtracted from the following sequence of C2 and C3 images, respectively. For a better visualization of the wave-like structure and its propagation, the streamer was rotated clockwise by 27fdg5 to horizontal. The shape of the streamer axis was traced by hand and marked by red plus signs. In each frame, if applicable, the crest and trough are indicated in yellow. The first crest is indicated by P1, the first trough by P2, and the second crest by P3. From the C2 observations, we can see that the streamer started to oscillate only above a distance of about 3 R from the Sun's center. After 08:12 UT, the streamer brightness intensified, which can be seen, e.g., from the large brightness gradient along the streamer near 4.6 R at 09:24 UT. This also modified the shape of the streamer, especially below a heliocentric distance of 5–6 R.

Figure 3.

Figure 3. Some selection of the traced streamer in C2 and C3 images. The streamer was rotated clockwise by 27fdg5 to be horizontally located. The traced streamer at different instances is delineated by red plus signs. The first crest is indicated by P1, the first trough by P2, and the second crest by P3. For the context of the streamer, a strip from −1.5 R to 1.5 R in the north–south direction, from 2 R to 6.5 R in the east–west direction, was cropped in the C2 images. In the C3 images, the range in the east–west direction is from 4 R to 12 R.

Standard image High-resolution image

The complete data set of the traced streamer axis from the C2 and C3 images is displayed in Figure 4. The time cadences of the C2 and C3 images were both generally 12 minutes but occasionally had an observational gap of more than 12 minutes. The streamer traced in C2 is indicated by dashed lines and that in C3 by solid lines. Streamer traces observed with a time difference of less than six minutes are drawn in the same color. Although after 08:12 UT the full length of the streamer was traced out, in the rest of our analysis we neglect the segment behind the brightness intensification noted above.

Figure 4.

Figure 4. Traced streamer in all the frames observed from 07:12 UT to 10:30 UT by C2 and C3. The streamer within C2 FOV is indicated by the long dashed lines, and within C3 FOV by the solid lines. The streamers traced at different times are shown in red (R), green (G), and blue (B) color. Those with observational time by C2 and C3 within six minutes are displayed with the same color.

Standard image High-resolution image

The streamer positions traced in Figure 4 may be subject to some uncertainty. To investigate the amplitude variation with distance and time, in Figure 5, the streamer segments around the first and second crests P1 and P3 are compared after they propagated for slightly more than an hour from the C2 FOV into C3. The upper dashed horizontal line in Figure 5 marks the peak of P3 and the lower dashed line the peak of P1. These lines show that the upper boundary of crests P1 and P3 increased from 0.3 to 0.5 R when they propagated to larger distances from about 5.2 to 7.5 R, which indicates an increase in their amplitudes with distance. If we use the two short solid lines as the lower boundary of the streamer segment around its peak at 07:36 UT and 08:48 UT, then the amplitude, as calculated from half the average distance of the lower and upper boundaries, increases slightly from 0.14 to 0.17 R. It also increases slightly with time during our observations. Although these estimates of the amplitude increase are not very precise due to data noise, it is very unlikely that the streamer in our study experienced an amplitude decrease as was observed by Chen et al. (2010). Therefore, we rule out that the wave-like behavior of the streamer in our observations was caused by an interaction of the streamer with a CME as reported in their paper.

Figure 5.

Figure 5. Segments of the streamer around the crests P1 and P3 at different heliocentric distances. The dashed lines mark the upper boundary of the crests P1 and P3 and the the solid bars indicate their lower boundary. The X and Y axes are in units of R.

Standard image High-resolution image

Since the brightness in LASCO C2 and C3 is caused by Thomson scattering of electrons in the corona, the observed coronagraph brightness is proportional to the column density along the line of sight. Moreover, based on the fact that the streamer can distinctly be seen on the dark background sky and appears as a narrow ribbon, we conclude that we observe the streamer more or less edge-on. We therefore assume that the observed column densities can serve as a proxy of the volume density along a profile normal to the streamer axis. In the upper panel of Figure 6, the eight red lines mark the profiles along which the brightness measurements were performed. Note that the pre-event image was not subtracted from this image. Instead, a monthly minimal image was subtracted to remove the stray light and/or the F corona.

Figure 6.

Figure 6. Upper panel: cut of the C2 image at 07:24 UT in which the streamer started to show clear oscillations. The eight short red lines mark the positions where the brightness distribution across the streamer is measured. The X and Y axes are in units of R. Lower panel: magnetic field configuration of the streaming kink instability adopted from Figure 1 in Wang et al. (1988). The X and Y axes are in units of the half thickness of the current layer.

Standard image High-resolution image

In Figure 7, we display the brightness distribution across the streamer for these eight profiles, extracted from the observation at 07:24 UT when the oscillation of the streamer started to become pronounced. Two dotted lines represent the maximum ρs and assumed background ρ of the brightness. Obviously, the ratio ρs declines approximately with heliocentric distance from ρs ≃ 2 to ρs ≃ 1.5, except at r = 4 R. The brightness distribution at r = 2.5 R was contaminated by a bright structure in the southern part of the streamer. Therefore, for this heliocentric distance we used the background estimate from r = 2.8 R instead. Interestingly, the streamer wave is discernible only at a distance beyond r = 3.5 R, when the density ratio has already settled to its lower bound of ρs ≃ 1.5. As another parameter, we can estimate from Figure 7 the width of the current sheet. For altitudes below r = 4 R the streamer is well concentrated with a half width of δ ≃ 0.2 R. Beyond 4 R, the streamer appears to be more diffusive, probably due to decreasing signal-to-noise ratio.

Figure 7.

Figure 7. Brightness distribution across the streamer at 07:24 UT along the short red lines in Figure 6 at different heliocentric distances. The distance and the corresponding ratio of the maximal to the background brightness are indicated in the upper part of each panel. Two different background models for the stray-light and F corona subtraction, that is, the monthly minimum and the yearly minimum, are applied. The ratio using the yearly minimum background is indicated in the parentheses in each panel. Two parallel dotted lines represent the maximal and the background brightness. The distance along the X axis is measured from the upper to lower positions on the short red line.

Standard image High-resolution image

There are two major sources of uncertainties involved in the estimate of the density ratio. One is the unknown line-of-sight depth of the streamer. We assumed the density ratio to be equal to the ratio of the observed column densities, which implies an infinite streamer depth. If the real depth of the streamer is much less than R, then we underestimate the density ratio. The other unknown is related to the background subtraction of the coronagraph image. To obtain an upper bound of the ratio, a yearly minimal image was utilized instead of the monthly minimum background. The corresponding density ratios are indicated in the parentheses in each panel. Although the ratios increase slightly, the trend of decreasing ratio with distance is clear.

Based on the retrieved wave-like structure in Figure 4, we have estimated parameters that characterize the streamer wave. The heliocentric distance of the crests and troughs, P1, P2, and P3, is plotted in Figure 8 as a function of time. Applying linear fits to their space-time positions yields a phase speed for P1, P2, and P3 of 406 km s−1, 374 km s−1, and 356 km s−1, respectively. The decrease of the phase speed from P1 to P3 can be directly concluded from the slight divergence of the fit lines in Figure 8. As seen in Figure 3, at 08:06 UT, the half wavelength $\overline{P1P2}$ was about 1 R; by 09:30 UT, the half wavelength $\overline{P1P2}$ had grown to 1.5 R. The period of the wave was measured from the time interval between the first and second crests, P1 and P3. We found a period of approximately 72 minutes at r = 5 R and 84 minutes at r = 9 R.

Figure 8.

Figure 8. Distance–time plot of P1, P2, and P3. Linear fits indicated by the solid lines are applied to derive their corresponding phase speeds, which are 406, 374, and 356 km s−1. The error bars are the assumed δ = ±0.2 R uncertainties. The propagation of this error to the phase speeds are indicated in the right corner as the 1σ uncertainty of the fitted slope.

Standard image High-resolution image

The major uncertainty involved in the above analyses comes from the delineation of the streamer profile in the LASCO observations (Figure 3). Once the profiles are determined, tracking the crests and troughs of the profiles and fitting a linear function to their distance–time plots only contributes a minor uncertainty. We assume that the uncertainty of the streamer identification corresponds to a half width δ of the streamer of about ±0.2 R. This is also the uncertainty of the positions of the crests and troughs in Figure 8. The linear fits bearing such uncertainties produce a 1σ error for the three velocity estimates of 16, 9, and 24 km s−1.

3. THE STREAMING KINK INSTABILITY

Analytical and MHD simulations have shown that streaming sausage, kink, and tearing instabilities can occur in the presence of sheared velocity and magnetic field systems. At the tip of a helmet streamer, the magnetic field has opposite directions on two sides with a neutral sheet inside. The field strength inside the streamer is probably well below the value outside and the pressure deficit is compensated for by an enhanced plasma density. The plasma velocity assumes a minimum in the neutral current sheet that can be considered the source of the slow solar wind and it increases toward both sides where the dilute plasma is accelerated into the fast solar wind. Plasma blobs, which were probably formed by the streaming sausage and tearing instabilities and convected away by the surrounding solar wind flow, have been observed at the tip of the helmet streamer (e.g., Sheeley et al. 1997). So far, however, kink-mode waves have not been detected. In view of the amplitude increase of the wave-like structure that we observe, we tend to believe that it is instability-driven. The increase of the perturbation is probably fed by the kinetic energy of the surrounding solar wind plasma.

A very robust threshold for the streaming kink and sausage instabilities found in many simulations concerns the velocity difference of the shear flow. In the simulation of Chen et al. (2009), they assumed an Alfvén speed of 20 km s−1 inside and 75 km s−1 outside of the current sheet. With these parameters, Chen et al. (2009) could excite both modes with a velocity difference in the shear flow of 100 km s−1 for the kink and about 150 km s−1 for the sausage mode.

If we take the measurements by the Ulysses spacecraft (McComas et al. 2008) as a typical estimate of the solar wind speed and its gradients, then these threshold values for the velocity shear can easily be achieved. Note, however, that the Ulysses measurements were made beyond 1 AU where the solar wind has been accelerated to its asymptotic speed, which has probably not yet been reached at 10 R. Also, the Alfvén speed assumed by Chen et al. (2009) is relatively low. A typical value beyond 2 R is 200 km s−1 (e.g., Warmuth and Mann 2005). However, in cusp and streamer regions, the field strength is probably reduced with respect to this average.

Lee et al. (1988) found from an analytic solution of a four-layer model that the kink mode dominates over the sausage mode when β exceeds unity in the ambient plasma, or equivalently, when ρs is less than about two. Their model consists of two layers of width δ with adjacent half spaces. The magnetic field is tangential to the boundaries and asymmetric to the central boundary. The flow is parallel to the magnetic field and symmetric with respect to the central boundary. If we keep the plasma temperature constant throughout, then there are six parameters aside from the temperature: magnetic field, flow speed, and plasma density, both inside the inner layers and in the half space, which can be varied. The unperturbed pressure balance of the current sheet reduces the number of independent parameters by one.

Considering this number of model parameters, we have re-examined Lee et al.'s dispersion relation (Equations (34) and (35) of Lee et al. 1988) in the parameter range that we find most appropriate for our observations. We found a half width of the streamer of about δ ≃ 0.2 R and a typical wave number of $k=2\pi /2.5 \, R_\odot ^{-1}$ for the wave structure, which yields kδ = 0.5. Assuming the magnetic field in the streamer was half of the outside value, the temperature was about 1 MK, and taking into account the pressure balance (Equation (10) in Lee et al. 1988)

the observed ρs = 1.5 results in VA ≈ 100 km s−1, and the plasma β then adjusts to values of 1.5 outside and 9 inside the streamer.

In Lee et al.'s (1988) work, large constant growth rates are obtained for kδ > 2 for both the kink and sausage modes. They are probably due to the sharp step-like gradients between the piecewise constant model profiles assumed. Numerical simulations with more realistic profiles (e.g., Zaliznyak et al. 2003) yield maximum growth rates well below kδ ≃ 1. From these dispersion relations, we find that the phase speed does not vary drastically with wave number and in the limit kδ → asymptotically becomes

where ρ and V are the density and flow speed outside the streamer, and ρs and V + ΔVs are the same parameters inside the streamer. Note that different from Lee et al. (1988), we here assume a wake model for the streamer so that ΔVs is negative. In order to reach the observed phase velocities of Vϕ ≃ 400, ..., 350 km s−1, the choice of values for V and ΔVs is considerably constrained. The set of values for which we find a reasonable agreement with the observations is V = 500 km s−1 and ΔVs = −150 km s−1 = −1.5VA.

In Figure 9, for different values of the density ratio ρs we show the threshold values for ΔVs, growth rate, and phase speed for V = 500 km s−1 and ΔVs = −150 km s−1. From the threshold diagram, we conclude that waves with kδ ≃ 0.3–0.5 can readily be excited by a velocity shear ΔVs not much above VA. We have to keep in mind that the Alfvén velocity is reduced inside the streamer layer by a factor of two. The value of |ΔVs| = 1.5VA for the growth rate and phase speed that we adopt below is just above the threshold. If the threshold is exceeded, then the largest growth rates for the kink mode occur at kδ > 0.5 and are achieved for the density ratio ρs between 1.5 and 3. For larger density ratios, the sausage mode grows faster. For the observed density ratio of ∼1.5, the sausage mode is stable and has a zero growth rate. Therefore, only the kink mode was observed. For kδ much larger than 0.5, we believe that the analytical dispersion based on a model with sharp, step-like gradients overestimates the growth rate.

Figure 9.

Figure 9. Normalized velocity shear threshold, growth rate, and phase speed vs. wave number for different density ratios ρs from the dispersion equation for the kink(black) and sausage(red) modes of Lee et al. (1988). The curves are labeled according to the density ratio ρs. For the ratios 1.5 and 1.1, the sausage mode is stable in the wavelength range shown. The growth rate and phase speed were calculated for ΔVs = −1.5 VA. The other parameters are given in the text.

Standard image High-resolution image

The phase speed of the kink mode for the observed density ratio ∼1.5 and the observed wave number kδ ≃ 0.5 is Vϕ = 3.8 VA = 380 km s−1. Again, this comes close to what we observed. Note that each unstable kink and sausage mode is a coalescence of an up- and downshifted propagating MHD boundary wave, so that the phase speed splits for low wave numbers where the mode becomes stable. When unstable, the two waves merge into a pair with complex conjugate frequencies and a unique phase speed.

The growth rates rise quite steeply once the threshold is exceeded and a value of $\gamma =0.2 kV_{A\inf }$ does not seem unrealistic. However, we have to keep in mind that the instability occurs in an inhomogeneous environment with fixed boundary conditions on the Sunward side of the streamer. At this side, the wave amplitude can only have noise level. Under these conditions, the instability becomes convective instead of showing temporal growth.

The growth scale that we can expect is Vϕ/γ = Vϕ/(0.2 kVA). With the parameters used in our model, Vϕ = 380 km s−1 and VA = 100 km s−1, we estimate the e-folding growth scale to be 3.8/0.2 · λ/2π = 3λ. From the increase of the observed wave amplitude that is apparent from the streamer traces in Figure 4, this number does not seem unrealistic.

4. CONCLUSIONS AND DISCUSSION

We present for the first time the observation of an oscillating streamer probably driven by the streaming kink instability, a mode of the MHD KHI. The streamer presented a period of about 70–80 minutes, and its wavelength increased from 2 to 3 R in about 1.5 hr. The phase speed of the oscillation varied from 406 ± 20 to 356 ± 31 km s−1 for different crests and trough. These derived properties of the wave structure place constraints on the plasma environment in which the instability took place. We find a low Alfvén speed of about 100 km s−1 and a high flow speed of about 500 km s−1 outside the streamer. It also turns out that the flow speed within the streamer was 150 km s−1 lower.

In previous studies, oscillations of a streamer have been observed following an interaction with a fast CME (Chen et al. 2010, 2011). The phenomenon was termed a "streamer wave" by the authors. A distinct characteristic of their observation was the decrease of the oscillation amplitude with time (see their Figure 4) due to the convection of the energy with the outward propagation of the wave. In the event reported here, no CME was detected in the LASCO FOV. Also, the coronagraph observations by STEREO A and STEREO B did not show any CME activity close to the streamer. Another distinction from the previously reported oscillating streamer, within the same distance range, was the increase in the amplitude of the wave with time. The streamer also displayed a pronounced increase in amplitude when it convected to larger heliocentric distances. This is evidence that the wave motion was excited by an instability.

We assume that the streamer is a one-dimensional layer with a plasma flow and magnetic field mutually parallel and normal to the layer gradient according to the model invoked by Lee et al. (1988). The plasma velocity is symmetrically distributed, the magnetic field varies asymmetrically, and the plasma temperature is constant throughout. We do not distinguish between the widths of the velocity shear and the magnetic field shear, and the scale of the density variation. With these model assumptions, we could reproduce many features of our observation.

In order to reproduce the observed phase speed, we had to invoke a high flow speed and a low Alfvén speed outside of the streamer. Typically assumed values for the slow solar wind in the distance range of 8–10 R are VA ≃ 200 km s−1 (Warmuth and Mann 2005) and V ≃ 200 km s−1 (e.g., Sheeley et al. 1997). The latter authors determined V from the velocity of density fluctuations drifting off streamer tips. They also show that the speed of individual density blobs scatters enormously around the above quoted average and may reach more than 600 km s−1 in some cases. The value of V = 500 km s−1 that we obtain from our observations is therefore not unrealistic. Moreover, at the time of our observations, the streamer was sandwiched between two nearby coronal holes to the north and south of the streamer (see Figure 2). It does not seem impossible that the fast solar wind from these coronal holes had some impact on the flow speed V in the vicinity of the streamer.

In our analysis, we also had to assume a low Alfvén speed VA in order to reproduce the observed phase speed of 350–400 km s−1. Our value of VA = 100 km s−1 is low, but still slightly higher than the value taken by Chen et al. (2009) in their simulations of a coronal streamer. Higher values for VA require even larger flow speeds V and also a larger velocity shear magnitude |ΔVs| inside the streamer. From our estimates, ΔVs = −1.5VA just exceeds the threshold condition. For a larger ΔVs and also for density ratios ρs > 3, the kink instability would be superseded by the sausage mode. We believe that these special conditions are why the kink mode is only seldom observed with coronal streamers compared to the sausage mode. This alternative mode manifests itself in the formation of density fluctuations and plasmoids that are more commonly observed in coronal streamers.

The streamer is also special in that its cusp is probably lower than for conventional streamers. We have tried to construct the magnetic field near the streamer by a PFSS model. For mathematical reasons, these models invoke a spherical outer boundary at some distance rSS where the field becomes entirely radial. Typically, rSS is chosen to be between 2 and 2.5 R, but for the real coronal magnetic field the streamer cusps do not need to have rSS at the same height. Modeling the field in the vicinity of the streamer, we found that its current sheet is mapped into the heliosphere at exactly the observed latitude if rSS is given a value of less than 1.6 R.

The alternative possibility is a pseudostreamer configuration (Wang et al. 2007) for the observed oscillating streamer. In the case of the pseudostreamer, the initial magnetic field has the same direction in different layers. We have gone through the formula in Lee et al. (1988) in order to check how the dispersion equation and the shape of a streamer would change when the configuration of the magnetic field varies from the case of the helmet streamer to the case of the pseudostreamer. Although the details of the derivation are beyond the scope of this paper and are not presented here, surprisingly, we find that the dispersion equation and the related instability threshold, growth rate, phase speed, and shape of the streamer do not change. The only difference in the resulting magnetic field is its direction. For helmet streamers, its direction is opposite on two sides of the streamer axis, while for pseudostreamers, the direction of the magnetic field is the same.

We thank the anonymous referee for insightful comments and a careful reading of our paper. We also benefited from valuable discussions with Yao Chen, Jiansen He, and Hui Tian. SOHO and STEREO are projects of international cooperation between ESA and NASA. The SECCHI data are produced by an international consortium of NRL, LMSAL, and NASA GSFC (USA), RAL, and U. Birmingham (UK), MPS (Germany), CSL (Belgium), IOTA, and IAS (France). SDO is a mission of NASAs Living With a Star Program. L.F. and W.Q.G. are supported by MSTC Program 2011CB811402, NSFC under grants 11003047, 11233008, 11273065 and by grants BK2012889. L.F. also acknowledges the Key Laboratory of Dark Matter and Space Astronomy, CAS, for financial support. The work of B.I. is supported by DLR contract 50 OC 1301.

Please wait… references are loading.
10.1088/0004-637X/774/2/141