Issue 
A&A
Volume 518, JulyAugust 2010
Herschel: the first science highlights



Article Number  A37  
Number of page(s)  9  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201014394  
Published online  31 August 2010 
Linear Alfvén waves in the solar atmosphere
K. Murawski^{1}  Z. E. Musielak^{2,3}
1  Group of Astrophysics, Institute of Physics, UMCS,
ul. Radziszewskiego 10, 20031 Lublin, Poland
2  Department of Physics, University of Texas at Arlington, Arlington, TX 76019, USA
3  KiepenheuerInstitut für Sonnenphysik, Schöneckstr. 6, Freiburg 79104, Germany
Received 10 March 2010 / Accepted 11 May 2010
Abstract
Aims. We aim to analytically and numerically explore smallamplitude Alfvén waves in the solar atmosphere.
Methods. We transform the wave equations to obtain the
cutofffrequency and wave travel time for strictly linear Alfvén waves.
The wave equations are solved numerically to find out spatial and
temporal signatures of the waves.
Results. The analytical predictions are verified by numerically
solving the wave equations for linear Alfvén waves. The waves are
impulsively generated and their characteristics and behavior in the
solar atmosphere are investigated by the numerical simulations. The
derived cutofffrequency is used to determine regions in the solar
atmosphere where strong reflection occurs for Alfvén waves of different
frequencies.
Conclusions. The numerical results reveal that impulsively
generated smallamplitude waves exhibit characteristic spatial and
temporal signatures which agree with the predictions of the analytical
theory.
Key words: magnetohydrodynamics (MHD)  Sun: atmosphere
1 Introduction
The Sun is a compressible medium and so is able to support acoustic waves. However, the presence of a magnetic field introduces additional elasticity to the natural compressibility of the plasma. A magnetic field imparts magnetic pressure and tension forces, which sustain additional wave motion. Among others Alfvén waves are of great interest as they may participate in the heating of the solar corona (e.g., Hollweg 1991) and the acceleration of the solar wind (e.g., Parker 1991). As a result, the evolution of Alfvén waves in the solar atmosphere is a fundamental problem which requires further exploration. This problem is formidable though, because the solar atmosphere is a strongly inhomogeneous and structured medium (e.g., Zaqarashvili & Murawski 2007).
We consider an atmospheric model for which the basic conditions for the propagation of linear Alfvén waves can be determined by an analytic theory. To verify the prediction of this theory, we numerically solve the resulting wave equations and compare the results of numerical simulations with those obtained analytically.
Some Alfvén waves originate from the dense photosphere where they can be reflected back into the photosphere, but some of them get through the cutoffregion and reach the solar corona (Musielak & Moore 1995). Through the coupling provided by the magnetic field the solar chromosphere and corona host a large amount of these waves. Indeed, observational evidence of Alfvéniclike fluctuations in the corona and chromosphere (DePontieu et al. 2007; Tomczyk et al. 2007; Bonet et al. 2008; Fujimura & Tsuneta 2009; Jess et al. 2009; WedemeyerBöhm & van der Voort 2009) indicate frequencies on the order of 39 min. However, some of these waves were interpreted by Van Doorsselaere et al. (2008a,b) as fast magnetoacoustic kink waves.
While there are several studies of the Alfvén wave propagation and reflection in stratified atmospheres that relies on the Fourier decomposition of an incident signal (e.g. Velli 1993), very few of them considered the problem for impulsively generated waves. For instance, Hollweg & Isenberg (2007) showed that for impulsively generated Alfvén wave the dynamics could be understood in terms of the first two reflected signals. The cutofffrequency criterion presented here offers an alternative approach to the usual Fourier analysis.
It has been shown that Alfvén waves may become nonlinear already in the chromosphere, exerting a ponderomotive force on magnetoacoustic waves (Suzuki & Inutsuka 2005; Verdini & Velli 2007; Verdini et al. 2009; Matsumoto & Shibata 2010). Numerical simulations of Alfvén waves launched from the chromosphere reveal that highfrequency modes are attenuated before they reach the corona by coupling with magnetoacoustic waves (Matsumoto & Shibata 2010). For our purpose here we only discuss smallamplitude Alfvén waves. Largeamplitude waves are left for future studies.
This paper is organized as follows. A physical model of the solar atmosphere is presented in Sect. 2. The following part of the paper discusses linear wave equations for Alfvén waves. Smallamplitude Alfvén waves are studied by numerical means in Sect. 4. The analytical and numerical results are discussed in Sect. 5. We conclude with a summary of the main results.
2 A physical model
2.1 MHD equations
Our model system is taken to be composed of a
gravitationallystratified solar atmosphere that is described by the
ideal magnetohydrodynamic (MHD) equations
Here is the mass density, is the flow velocity, is the magnetic field, p is the gas pressure, is the adiabatic index, is a vector of gravitational acceleration with its value g=274 m s^{2}, T is the temperature, m is the mean particle mass and is Boltzmann's constant.
2.2 A model of the solar atmosphere
We assume that the solar atmosphere is in static equilibrium (
)
with a forcefree magnetic field, i.e.
Henceforth the subscript corresponds to equilibrium quantities. As a result, the pressure gradient is balanced by the gravity force
Using the ideal gas law of Eq. (6) and the ycomponent of the hydrostatic pressure balance indicated by Eq. (8), we express the equilibrium gas pressure and mass density as
Here
is the pressure scaleheight, and denotes the gas pressure at the reference level that is set and hold fixed at Mm.
We adopt a smoothed stepfunction profile for the plasma temperature
where and denote the chromospheric temperature at its lower part. The symbol corresponds to the temperature of the solar corona that is separated from the chromosphere by the transition region that is located at Mm. A width of the transition region is km. We set and hold fixed .
According to Eq. (9) the temperature profile of Eq. (11) uniquely determines equilibrium mass density and gas pressure profiles (Fig. 1). Both and experience a sudden falloff below the transition region.
Figure 1: Equilibrium mass density (top panel) and gas pressure (bottom panel) in units of 10^{12} kg m^{3} and 1 Pa, respectively. The distance y is expressed in 1 Mm. 

Open with DEXTER 
We consider a uniform vertical magnetic field such as
Here is a unit vector along the ydirection. This choice of results in Eq. (7) being satisfied and the Alfvén speed, , is given as
whose profile is displayed in Fig. 2 (top panel). Note that in the low chromosphere km s^{1}. It raises up abruptly through the transition region, reaching a value of km s^{1} in the solar corona. The growth of results from the falloff of with height. For comparison purposes the sound speed,
(15) 
is shown in Fig. 2 (bottom panel). The sound speed reaches a value of about 7 km s^{1} below the altitude y=1 Mm and it grows to about 105 km s^{1} in the solar corona, for y> 3 Mm.
As a result of the Alfvén and sound speed profiles the plasma
(16) 
is illustrated in Fig. 3. Note that for y=1 Mm and it falls off abruptly at high altitudes reaching a value of at y=10 Mm.
Figure 2: Equilibrium profile of Alfvén speed (top panel), (middle panel) and sound speed (bottom panel). These speeds, and the distance y are expressed in units of 1 Mm s^{1}, 1 s^{1} and 1 Mm, respectively. 

Open with DEXTER 
3 Linear Alfvén waves
3.1 Wave equations for Alfvén waves
We now discuss smallamplitude perturbations of the equilibrium that is
described in Sect. 2.2. Assuming that plasma evolution occurs
only along the ydirection, the zcomponents of momentum Eq. (2) and induction Eq. (3) can be linearized, leading to the following equations:
Here and denote the perturbed zcomponents of the plasma velocity and magnetic field, respectively. Hence we arrive at the following wave equations:
Note that V_{z} and B_{z} are described by different wave equations, which means that there is a phase shift between the wave variables and, as a result, they behave differently. The derived wave equations have nonconstant coefficients and become the classical wave equations with constant coefficients when a homogeneous plasma with is considered.
Figure 3: Equilibrium profile of the plasma vs. height y (in 1 Mm). 

Open with DEXTER 
3.2 Cutofffrequency with approximate wave travel time
We rewrite the wave equations given by Eqs. (18) and (19)
in their standard form, i.e., without the first derivative term; the wave equations
written in their standard forms are often referred to as KleinGordon equations (e.g., Rae & Roberts 1982; Musielak et al. 1987, 2006).
Because there is no term with the first derivative in the wave equation for V_{z}, this
equation is already written in its standard form. However, in order to cast the wave
equation for B_{z} in its standard form, we must use the transformation
The resulting standard wave equation can be written as
where is the socalled critical frequency (e.g., Musielak et al. 1992, 2006) given by
Note that the critical frequency for V_{z} is zero.
We make Fourier transforms in time
and
,
and obtain
and
The transformed wave equations are in the form of Euler's equation , where is the second derivative of u with respect to x. This equation has the oscillatory and nonoscillatory solutions when and , respectively, and it has the turningpoint that separates these solutions when (e.g., Murphy 1960).
We define the turningpoint frequencies
and
for
and
,
respectively, by taking the nonconstant coefficients in Eqs. (23) and (24) to be equal to the
term 1 / 4 y^{2} from Euler's equation (see Musielak et al. 2006).
The results are
and
Note that the obtained turningpoint frequencies depend on the approximate wave travel time (Routh et al. 2010), which is defined as
Here denotes an atmospheric height at which waves are initially (at t=0 s) launched. To distinguish between and the actual wave travel time, the latter will be introduced in the next subsection.
Using the oscillation theorem (e.g., Kahn 1990), it is easy to show that the conditions and give oscillatory (or propagating wave) solutions for the corresponding wave variables. Obviously, the solutions are nonoscillatory (or nonpropagating waves) when and .
After obtaining the turningpoint frequencies for both wave variables,
we now follow Musielak et al. (2006) and select the cutofffrequency
from the condition
3.3 Cutofffrequency with actual wave travel time
To obtain a cutofffrequency that involves the actual wave travel time
,
we follow Musielak et al. (1992) and introduce
which allows us to write Eqs. (18) and (19) in the following forms:
Here
Using the transformations
and
we cast the wave equations given by Eqs. (18) and (19) in their standard forms
and
where the critical frequencies are defined as
and
The critical frequencies can also be expressed in terms of the Alfvén speed, , and its derivatives. Using Eq. (32), we obtain
and
It is interesting to note that the above expressions for the critical frequencies are the same as those obtained by Routh et al. (2010), who studied the propagation of torsional waves along solar magnetic flux tubes. This clearly implies that the physical properties of torsional tube waves and the Alfvén waves investigated here are the same.
After obtaining the critical frequencies, we make Fourier transforms in
time
and
,
and obtain
and
We use these equations and follow the approach described in the previous subsection to define the turningpoint frequencies and for the transformed wave variables and , respectively. The results are
It must be pointed out that both turningpoint frequencies are expressed in terms of , which is the actual wave travel time in an inhomogeneous medium
Here is the integration constant to be determined from a given model of the background medium.
Based on the oscillation theorem (e.g., Kahn 1990), the wave propagation
conditions are
and
.
These conditions require that both wave variables
are described by the propagating wave solutions and can be simplified
by introducing the cutofffrequency
given by
It is important to point out that the resulting cutofffrequency depends on the actual wave travel time and that can only be determined once a model of the background medium is specified.
3.4 The cutofffrequency
In the previous two subsections, we derived two different cutofffrequencies and . We obtained the cutofffrequency from the transformed wave equations that describe the wave propagation using the variable y, and we demonstrated that this cutoff explicitly depends on the approximate wave travel time. On the other hand, the cutofffrequency was obtained from another set of transformed wave equations that involved the variable , which represents the actual wave travel time. There were also significant differences between the critical frequencies for these two cases.
With two different cutofffrequencies, we must now determine which one
is the true cutofffrequency that correctly describes the propagation
conditions for the considered linear Alfvén waves. Because
correctly accounts for the actual wave
travel time, we choose it to be the true cutofffrequency. Hence, we
have
or
Once the true cutofffrequency is determined, we may convert into y and give the explicit form of . The result is
where the turningpoint frequencies are given by
and
with the actual wave travel time defined by Eq. (45). In addition, the explicit forms of the critical frequencies are
Now, the cutofffrequency can be obtained when we know which turningpoint frequency is higher. This can only be determined when a model of background atmosphere is specified. This model is given in Sect. 2, and we will use it below to obtain the cutofffrequency for linear Alfvén waves propagating along uniform magnetic field embedded in an nonisothermal atmosphere.
3.5 Physical model considered in this paper
We consider a stratified and nonisothermal medium
with a uniform background magnetic field (for details, see Sect. 2.2).
Let us introduce the density scale height
,
which is defined as
Note that for the background medium considered here, the density and pressure scale heights are different (see however Sect. 3.6).
Using the above definition, we obtain
and
Hence, the critical frequencies (see Eqs. (51) and (52)) can be written as
and
where
The above equations clearly show that in the considered model either and are both positive, or at least one of them is positive and the other one is negative. To demonstrate this, we introduce and consider two cases: and . In the first case, we have and when , and when , and and when . In the second case, we always have and . An interesting result is that no temperature distribution in this model can lead to both critical frequencies being negative.
For the temperature distribution given by Eq. (11), the density
scale height can be written as
where the pressure scale height is given by Eq. (10) and , with the sound speed that corresponds to the fixed temperature (see Eq. (11)). In addition, the derivative of the density scale height with respect to y is
Using Eqs. (59) and (60), we obtain the critical frequencies and then determine the cutofffrequency with the help of Eq. (48). Hence we get the Alfvén cutoffwave period,
This cutoffperiod is displayed in Fig. 4 (solid line) and its approximate analog is illustrated by a dashed line. In the low chromosphere s at Mm and at higher levels it falls off abruptly with y reaching a minimum s around the transition region. At higher altitudes it raises up, reaching a value of s at y=8.5 Mm. Note that lies close to for y<2 Mm, but it departs from in the solar corona.
Note here that a different approach in literature is based on a use of the Elsässer variables (Velli 1993; Hollweg & Isenberg 2007; Verdini & Velli 2007; Verdini et al. 2009). For instance, Velli (1993) used a Fourier approach, demonstrating that wave reflection is a global phenomenon. It is interesting that a rough estimate of the cutofffrequency, i.e. comparing the wave frequency with the maximum of the Alfvén speed gradient (e.g. Eq. (4) of Velli 1993), yields a period of about 7 s, which is not so different from the value of about 4 s obtained from Fig. 4 at the transition region.
Figure 4: Alfvén cutoffwave period of Eq. (61) (solid line) and its approximate analog of Eq. (28) (dashed line) in units of 1 s. Distance y is expressed in units of 1 Mm. 

Open with DEXTER 
The Alfvén cutoffperiod plotted in Fig. 4 can be used to determine the height in the solar atmosphere at which wave reflection becomes dominant. The condition for this height is , where is the wave period. For example, the reflection of Alfvén waves with a period of 80 s becomes important in the solar chromosphere at y = 1 Mm and in the solar corona at y = 6 Mm. The situation is different for short period (5 s) Alfvén waves. According to the results of Fig. 4, these waves are reflected in the solar transition region, but once they reach the solar corona, their reflection is negligible. It must be noted that neither the reflection of the 80 s waves nor the reflection of the 5 s waves is total; actually, the wave reflection in the atmosphere considered here is always partial (see Sect. 5 for details).
3.6 Special case of isothermal medium
A special case of a stratified and isothermal medium with the
uniform background magnetic field was considered by Musielak &
Moore (1995). There the density and pressure
scale heights are identical, which means that
(see Eq. (60)). Using Eqs. (56) and (57), we obtain
which are the same critical frequencies as those obtained by Musielak & Moore (1995).
Since
,
we get
the relationship
for the turningpoint frequencies, and obtain the cutofffrequency
Using this cutofffrequency, we write the condition for the propagation of linear Alfvén waves along the uniform magnetic field oriented along the yaxis in a stratified and isothermal medium.
4 Numerical simulations of the wave equations
As analytical solutions of Eqs. (16) and (17) are found so far for only a few specific choices of (Didenkulova et al. 2009), we adopt CLAWPACK which is a software package designed to compute numerical solutions to hyperbolic partial differential equations (LeVeque 2002). The simulation region 0.5 Mm <y< 20 Mm is covered by a uniform grid of 10^{4} numerical cells. We verified by convergence studies that this grid well represents the simulation region. We set open (with zerogradient) boundary conditions at the bottom and top boundaries of the simulation region.
The system of Eqs. (16) and (17) was initially (at t=0 s) perturbed by a Gaussian pulse as
We set and hold fixed the initial pulse width w_{y}=0.3 Mm and consider two cases: (a) waves launched below the transition region at Mm; and (b) waves launched above the transition region at Mm.
Figure 5: Wave profiles at (from top to bottom) t=10 s, t=60 s, t=160 s and t=500 s for the linear Eqs. (16) and (17) for Mm. Solid (dashed) line corresponds to V_{z} (B_{z}), which is expressed in units of 1 Mm s^{1} (10^{3} T). 

Open with DEXTER 
Figure 6: Timesignatures of V_{z} (solid line) and B_{z} (dashed line) collected at y=10 Mm, resulting from numerical solutions of Eqs. (16) and (17) for Mm. Time, V_{z} and B_{z} are expressed in units of 1 s, 1 Mm s^{1} and 10^{3} T, respectively. 

Open with DEXTER 
Figure 7: Normalized Fourier power spectrum. Waveperiod P is expressed in 1 s. 

Open with DEXTER 
Figure 8: Actual wave (solid line) and approximate (dashed line) travel times to the detection point Mm for Mm. Time and distance are expressed in units of 1 s and 1 Mm, respectively. 

Open with DEXTER 
Figure 5 illustrates the wave profiles for the case (a) at t=10 s (top panel), t=60 s (top middle panel), t=160 s (bottom middle panel), and t=500 s (bottom panel). Note that ( ) at the leading (trailing) part of the initial pulse of Eq. (64). According to Eq. (17) we get ( ) there. As a result B_{z} becomes negative (positive) at the wavefront (wavetail), which is clearly seen in Fig. 5 (top panel). The initial pulse splits into counterpropagating pulses. Downward propagating waves are evident at t=60 s (top middle panel). The initial pulse already spread into counterpropagating pulses. The upward propagating pulse partially reflects back into the low chromosphere and during this reflection the wave phase is not reversed. At t=500 s (bottom panel) both the reflected and originally propagating downward pulses almost merged into a single quantity. As Eqs. (16) and (17) are coupled any initial signal in V_{z} excites perturbations in B_{z}. Indeed, these perturbations are obvious in Fig. 5.
Note that Eq. (19) can be rewritten as
(66) 
From Fig. 2 we infer that the gradient of Alfvén speed attains its maximum at the transition region. As a result of that downward propagating perturbations in B_{z} (Fig. 5, dashed line) are much larger than the perturbations in V_{z}. Note that in order to compare magnitudes of V_{z} and B_{z} the latter was expressed in 10^{3} Teslas.
Timesignatures which correspond to the wave profiles of Fig. 5 are displayed in Fig. 6. These timesignatures are made by collecting wave signals at the detection point y=10 Mm. These waves passed through the transition region and reached the solar corona. Note that the wave signals grow fast, reaching their maximum at s. Later, the wave signals fall off.
Figure 7 shows a Fourier power spectrum, which reveals that most power is concentrated at the long wave period waves with a maximum located at s, which falls within the range of 39 min oscillations detected in the solar corona (Jess et al. 2009; DePontieu et al. 2007; Tomczyk et al. 2007; Bonet et al. 2008; WedemeyerBöhm & van der Voort 2009; Fujimura & Tsuneta 2009).
Note that the approximate arrival time to the level y=10 Mm is s, which is shorter than s (Fig. 8) to which the numerical data fits very well.
Figure 9: Wave profiles at t=2 s ( top panel), t=5 s ( middle panel) and t=10 s ( bottom panel) for the linear Eqs. (16) and (17) for Mm. Solid (dashed) line corresponds to V_{z} (B_{z}), which is expressed in units of 1 Mm s^{1} (10^{6} T). 

Open with DEXTER 
Figure 10: Timesignatures of V_{z} (solid line) and B_{z} (dashed line) collected at y=1.3 Mm, resulting from numerical solutions of Eqs. (16) and (17) for Mm. V_{z} and B_{z} are expressed in units of 1 Mms^{1} and 10^{6} T, respectively. 

Open with DEXTER 
Figure 11: Normalized Fourier power spectrum vs. wave period P (in 1 s) for the timesignature of Fig. 10. 

Open with DEXTER 
Figure 12: Actual wave (solid line) and approximate (dashed line) travel times (in units of 1 s) to the detection point y = 1.3 Mm for Mm. Distance y is expressed in units of 1 Mm. 

Open with DEXTER 
Figure 10 shows timesignatures which are made by collecting wave signals at the detection point y=1.3 Mm. These signatures consist of short period waves with the maximum Fourier power spectrum amplitude located at s (Fig. 11). These waves have passed into the bottom layer through the transition region. Note that the signal arrives at the detection point at s, which agrees both with the approximate ( ) and accurate ( ) arrival times of Fig. 12.
5 Discussion of the analytical and numerical results
We analytically and numerically studied the propagation of linear Alfvén waves in the solar atmosphere. Our model of the solar atmosphere is nonisothermal and its temperature distribution approximates that observed in the solar chromosphere, transition region, and corona. The results of our analytical and numerical studies extend those previously obtained by Musielak & Moore (1995) and An et al. (1989), respectively, who investigated the Alfvén wave propagation in an isothermal atmosphere.
Our main analytical result is the cutofffrequency (or the cutoffperiod ) obtained with the actual wave travel time. In the model considered here, the cutoff is a local quantity and its variation with height is used to identify regions in the solar atmosphere where strong wave reflection occurs. Using the condition , where is the wave period, the atmospheric height at which Alfvén waves of a given period are reflected is established. The condition results from the turningpoint property of Euler's equation, and that Alfvén wave reflection is strong at the atmospheric height where the turningpoint is located was originally demonstrated by Musielak & Moore (1995). The analytical solutions obtained by these authors showed that linear Alfvén waves form standing wave patterns in their isothermal model. The reason for the existence of these standing waves is wave reflection and the resulting constructive interference between the propagating and reflected Alfvén waves (see also Musielak et al. 2006).
The formation of the standing wave patterns in the isothermal solar atmosphere was numerically investigated by An et al. (1989), who were able to reproduce full analytical solutions by their numerical simulations and identify a height in their atmospheric model at which Alfvén wave reflection is dominant. Their main result is that wave reflection occurs continuously in the isothermal atmosphere and that it reaches maximum at the height where the wave period becomes comparable to a local cutoffperiod. Their results also clearly demonstrate that wave reflection is always partial, even at its maximum. We extended the previous work by considering a more realistic nonisothermal solar atmosphere and investigating impulsively generated linear Alfvén waves.
Our main numerical results are presented and discussed in Sect. 4, so here we only comment on their relationship to our analytical results. The analytical predictions of the atmospheric height at which strong reflection occurs for waves of a given period and the numerical results agree well. Moreover, the numerical results also clearly show that wave reflection is only partial, even at the height where it reaches its maximum. This confirms the previous analytical (Musielak & Moore 1995) and numerical (An et al. 1989) results. Unfortunately, for the nonisothermal solar atmosphere model considered here, we were unable to obtain full analytical solutions to the Alfvén wave equations. As a result, we were unable to make directly compare between analytical and numerical solutions obtained for very long simulation times. We hope that this comparison would further validate our results.
6 Conclusions
We studied smallamplitude Alfvén waves propagating in the stratified solar atmosphere. We particularly considered impulsively triggered waves in the low and top layers of the solar atmosphere. Guided by this we derived the transformed wave equations for the transversal velocity and magnetic field components which correspond to linear Alfvén waves. From these equations we obtained the cutoffwave period , which varies with height in the solar atmosphere. We used the cutoff to determine the regions in the solar atmosphere where the reflection of Alfvén waves becomes dominant.
The analytical findings were subsequently verified by numerical solutions of the wave equations. The numerical results for impulsively generated waves revealed that the Alfvén waves exhibit characteristic spatial and temporal signatures, which agree with the analytical predictions. Specifically, the analytically predicted height where Alfvén wave reflection dominates in the wave behavior, agrees with the numerical results. Moreover, the numerical simulations clearly show that wave reflection is always partial, even at the height where it reaches its maximum.
The results presented here are hopefully enable a deeper understanding of the complex relation between Alfvén waves and transient events in the solar chromosphere and corona. This work paves the way for future studies and it serves as a reference to the more advanced numerical solutions of the MHD equations, which will be reported elsewhere. In particular, a more realistic model of the solar atmosphere has to be taken into consideration. In this model a more appropriate temperature profile than the simple model of Eq. (11) as well as nonideal plasma effects and neutral plasma components in the photosphere have to be implemented.
AcknowledgementsThe authors are grateful to the referee for his/her stimulating comment, which resulted in an improvement of the presentation of the results. K.M. appreciates a fruitful discussion with Prof. Valery Nakariakov. The work of K.M. was supported by the Polish Ministry of Science (the grant for years 20072010). Z.E.M. acknowledges the support of this work by NSF under grant ATM0538278, the Alexander von Humboldt Foundation, and the University of Texas at Arlington through its Faculty Development Leave Program.
References
 An, C.H., Musielak, Z. E., Moore, R. L., & Suess, S. T. 1989, ApJ, 345, 597 [NASA ADS] [CrossRef] [Google Scholar]
 Bonet, J. A., Marquez, I., Sanchez Almeida, J., Cabello, I., & Domingo, V. 2008, ApJ, 687, 431 [NASA ADS] [CrossRef] [Google Scholar]
 De Pontieu, B., McIntosh, S.W., Carlsson, M., et al. 2007, Science, 318, 1574 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Didenkulova, I., Pelinovsky, E., & Soomere1, T. 2009, JGR, 114, C07006 [CrossRef] [Google Scholar]
 Fujimura, D., & Tsuneta, S., 2009, ApJ, 702, 1443 [NASA ADS] [CrossRef] [Google Scholar]
 Hollweg, J. V. 1991, in Mechanisms of Chromospheric and Coronal Heating, Proceedings of the International Conference, Heidelberg, 58 June 1990, XV, 649, 260 figs, ed. P. Ulmschneider, E. R. Priest, & R. Rosner (Heidelberg, New York: SpringerVerlag), 423 [Google Scholar]
 Hollweg, J. V., & Isenberg, P. A. 2007, JGR, 112, A08102 [NASA ADS] [CrossRef] [Google Scholar]
 Jess, D. B., Mathioudakis, M., Erdélyi, R., et al. 2009, Science, 323, 1582 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kahn, P. B. 1990, Mathematical Methods for Scientists and Engineers (New York: John Wiley) [Google Scholar]
 LeVeque, R. J. 2002, Finite Volume Methods for Hyperbolic Problems (Cambridge University Press) [Google Scholar]
 Matsumoto, T., & Shibata, K. 2010, ApJ, 710, 1857 [NASA ADS] [CrossRef] [Google Scholar]
 Murphy, G. M. 1960, Ordinary Differential Equations and Their Solutions, D. Van Nostrand (New York: Company) [Google Scholar]
 Musielak, Z. E., & Moore, R. J. 1995, ApJ, 452, 434 [NASA ADS] [CrossRef] [Google Scholar]
 Musielak, Z. E., Rosner, R., & Ulmschneider, P. 1987, in Cool Stars, Stellar Systems and the Sun, ed. J. L. Linsky, & R. E. Stencel, Lect. Notes Phys., 291, 66 [Google Scholar]
 Musielak, Z. E., Fontenla, J. M., & Moore, R. L. 1992, Phys. Fluids, 4, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Musielak, Z. E., Musielak, D. E., & Mobashi, H. 2006, Phys. Rev. E, 73, 0366121 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1991, ApJ, 376, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Rae, I. C., & Roberts, B. 1982, ApJ, 256, 761 [NASA ADS] [CrossRef] [Google Scholar]
 Routh, S., Musielak, Z. E., & Hammer, R. 2010, ApJ, 709, 1297 [NASA ADS] [CrossRef] [Google Scholar]
 Suzuki, T. K., & Inutsuka, S.I. 2005, ApJ, 632, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Tomczyk, S., McIntosh, S. W., Keil, S. L., et al. 2007, Science, 317, 1192 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Van Doorsselaere, T., Brady, C. S., Verwichte, E., & Nakariakov, V. M. 2008a, A&A, 491, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Van Doorsselaere, T., Nakariakov, V. M., & Verwichte, E. 2008b, ApJ, 676, L73 [NASA ADS] [CrossRef] [Google Scholar]
 Velli, M. 1993, A&A, 270, 304 [NASA ADS] [Google Scholar]
 Verdini, A., & Velli, M. 2007, ApJ, 662, 669 [NASA ADS] [CrossRef] [Google Scholar]
 Verdini, A., Velli, M., & Buchlin, E. 2009, ApJ, 700, L39 [NASA ADS] [CrossRef] [Google Scholar]
 Zaqarashvili, T. V., & Murawski, K. 2007, A&A, 470, 353 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 WedemeyerBöhm, S., & van der Voort, L. R. 2009, A&A, 507, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Figures
Figure 1: Equilibrium mass density (top panel) and gas pressure (bottom panel) in units of 10^{12} kg m^{3} and 1 Pa, respectively. The distance y is expressed in 1 Mm. 

Open with DEXTER  
In the text 
Figure 2: Equilibrium profile of Alfvén speed (top panel), (middle panel) and sound speed (bottom panel). These speeds, and the distance y are expressed in units of 1 Mm s^{1}, 1 s^{1} and 1 Mm, respectively. 

Open with DEXTER  
In the text 
Figure 3: Equilibrium profile of the plasma vs. height y (in 1 Mm). 

Open with DEXTER  
In the text 
Figure 4: Alfvén cutoffwave period of Eq. (61) (solid line) and its approximate analog of Eq. (28) (dashed line) in units of 1 s. Distance y is expressed in units of 1 Mm. 

Open with DEXTER  
In the text 
Figure 5: Wave profiles at (from top to bottom) t=10 s, t=60 s, t=160 s and t=500 s for the linear Eqs. (16) and (17) for Mm. Solid (dashed) line corresponds to V_{z} (B_{z}), which is expressed in units of 1 Mm s^{1} (10^{3} T). 

Open with DEXTER  
In the text 
Figure 6: Timesignatures of V_{z} (solid line) and B_{z} (dashed line) collected at y=10 Mm, resulting from numerical solutions of Eqs. (16) and (17) for Mm. Time, V_{z} and B_{z} are expressed in units of 1 s, 1 Mm s^{1} and 10^{3} T, respectively. 

Open with DEXTER  
In the text 
Figure 7: Normalized Fourier power spectrum. Waveperiod P is expressed in 1 s. 

Open with DEXTER  
In the text 
Figure 8: Actual wave (solid line) and approximate (dashed line) travel times to the detection point Mm for Mm. Time and distance are expressed in units of 1 s and 1 Mm, respectively. 

Open with DEXTER  
In the text 
Figure 9: Wave profiles at t=2 s ( top panel), t=5 s ( middle panel) and t=10 s ( bottom panel) for the linear Eqs. (16) and (17) for Mm. Solid (dashed) line corresponds to V_{z} (B_{z}), which is expressed in units of 1 Mm s^{1} (10^{6} T). 

Open with DEXTER  
In the text 
Figure 10: Timesignatures of V_{z} (solid line) and B_{z} (dashed line) collected at y=1.3 Mm, resulting from numerical solutions of Eqs. (16) and (17) for Mm. V_{z} and B_{z} are expressed in units of 1 Mms^{1} and 10^{6} T, respectively. 

Open with DEXTER  
In the text 
Figure 11: Normalized Fourier power spectrum vs. wave period P (in 1 s) for the timesignature of Fig. 10. 

Open with DEXTER  
In the text 
Figure 12: Actual wave (solid line) and approximate (dashed line) travel times (in units of 1 s) to the detection point y = 1.3 Mm for Mm. Distance y is expressed in units of 1 Mm. 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.