Issue 
A&A
Volume 512, MarchApril 2010



Article Number  A23  
Number of page(s)  9  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/200913365  
Published online  24 March 2010 
Propagating magnetohydrodynamic waves in a cooling homogenous coronal plasma
R. J. Morton^{1}  A. W. Hood^{2}  R. Erdélyi^{1}
1  Solar Physics and Space Plasma Research Centre
(SP^{2}RC), University of Sheffield, Hicks Building, Hounsfield
Road, Sheffield S3 7RH, UK
2  School of Mathematics and Statistics, University of St. Andrews,
St Andrews, Fife KY16 9SS, UK
Received 28 September 2009 / Accepted 11 December 2009
Abstract
Aims. We present an investigation into how the cooling of
the background plasma influences the propagation of slow and fast
MHD wave modes supported by an unbounded, homogenous plasma.
Previous investigations have suggested that the cooling of the plasma
and a reduction in density could lead to the damping of fast
magnetoacoustic oscillations. We aim to investigate whether cooling of
the background plasma at a constant density may be responsible for the
damping of slow and fast modes.
Methods. The plasma is assumed homogeneous and the background
temperature (pressure) is decreasing with time. The temperature change
is assumed to be due to optically thin radiation. A special case of the
radiative function is chosen to allow an analytical assessment of the
effects of cooling on magnetoacoustic MHD modes and ensures the
temperature evolution of the background plasma due to this radiation
also matches the observed cooling profile of coronal loops.
Results. A timedependent dispersion relation is obtained on the
slow timescale of cooling and full timedependent solutions are found.
Leading order equations for the amplitude of the waves are obtained and
solved analytically for the slow and fast MHD modes. The cooling
of the plasma is found to cause the frequency of the magnetoacoustic
modes to decrease with time. The slow modes are found to experience a
greater change in frequency than the fast modes. More importantly, the
radiative losses also provide a significant damping of the slow mode
and a small damping of the component of the fast mode perpendicular to
the magnetic field. The damping of the slow mode is found to be strong
within typical lifetimes of oscillations observed in coronal
structures. Cooling could have important consequences and needs to be
assessed when trying to determine what mechanism is responsible for the
observed damping of coronal oscillations.
Key words: waves  magnetohydrodynamics (MHD)  Sun: corona  Sun: oscillations  plasmas
1 Introduction
It is well known that the solar atmosphere is a highly magnetised and dynamic plasma consisting of numerous magnetic structures (e.g. coronal loops, open flux tubes, prominences, etc.). Each of these structures is able to support a wide range of magnetohydrodynamic (MHD) waves and oscillations and are of much interest as they may provide a solution to the coronal heating problem (Erdélyi 2008; Taroyan 2008; Taroyan & Erdélyi 2009). There have been numerous observations of oscillations within solar atmospheric magnetic structures (for reviews, see e.g. De Moortel 2005; Nakariakov & Verwichte 2005; Banerjee et al. 2007) and the majority of them are heavily damped. In spite of considerable research efforts, the cause of damping is still unknown although numerous different mechanisms have been suggested with various degrees of success. At present, resonant absorption (Ruderman & Roberts 2002; Goossens et al. 2002; Terradas et al. 2006) appears to be the most likely method of damping fast kink waves and thermal conduction (Ofman & Wang 2002; De Moortel & Hood 2004; De Moortel & Hood 2003; MendozaBriceño et al. 2004) for the damping of propagating slow (longitudinal) waves and standing slow modes in hot loops (T<6 MK). Whichever mechanisms are responsible, it is widely thought that the dissipation of the wave energy could help heat the corona.So far, the majority of investigations into the theoretical modelling of coronal structures has focused on static solutions of the MHD equations, i.e. with timeindependent background plasma quantities. This has lead to the creation of the field of coronal seismology, with a number of useful seismological tools being developed for plasma diagnostics (for the latest reviews, see, e.g. Andries et al. 2009; Ruderman & Erdélyi 2009). One such tool is the ratio of the frequency of the fundamental mode to the first overtone of the fast kink mode, which provides estimates of the density scale height in the corona (Andries et al. 2005a). The success of this seismological tool has been due to a few observations of the fast kink mode and its overtone, see, e.g. Verwichte et al. (2004), Verth et al. (2008). A number of different authors have refined the model to include a series of parameters such as variation in the density structuring (Andries et al. 2005b; McEwan et al. 2006), loop cross section (Verth & Erdélyi 2008; Ruderman et al. 2008), loop shape (Dymova & Ruderman 2006; Morton & Erdélyi 2009a).
One aspect of the dynamic atmosphere that has very recently received attention is the temperature evolution of the coronal plasma. In particular, observations have revealed that many coronal loops undergo an exponential temperature decrease (Aschwanden & Terradas 2008; UgarteUrra et al. 2009), with characteristic cooling times of 5002000 s reported for EUV loops. Hydrodynamic simulations suggest the decrease in temperature also appears to be associated with a decrease in the density (Jakimiec et al. 1992).
In a recent paper, Morton & Erdélyi (2009b) suggested that standing transverse oscillations in a coronal loop would be damped due to the cooling of the plasma. No explicit cooling mechanism was stated and it was assumed that the loop was cooling due to some unspecified mechanism (e.g. radiation, thermal conduction). A density profile dependent upon the change in temperature was also taken into account. Aschwanden & Terradas (2008) suggested at EUV temperatures (T<2.0 MK) that radiation is the probable cooling mechanism of coronal loop plasma. Morton & Erdélyi (2009b) also found that the cooling would cause a change in the frequency of the transverse modes, causing the frequency to increase (or the period of the mode to decrease) with time. Further, a spatial shift in the longitudinal amplitude profile of the fast kink mode also occurred providing information for spatial magnetoseismological studies.
A recent observation by Erdélyi & Taroyan (2008) reported the observation of a standing slow mode in an EUV loop in multiple spectral lines. An increase in intensity in lower temperature lines and a decrease in higher temperature lines suggests a cooling scenario is occurring. As radiation is thought to be the dominant method of the cooling of coronal loop plasma, it may heavily influence the behavior of standing slow modes.
There have already been previous investigations into nonadiabatic plasmas and how the cooling and heating influences the oscillations supported within the plasma (Field 1965; Ibanez & Escalona 1993; De Moortel & Hood 2003; Terradas et al. 2005; Carbonell et al. 2006). However, the plasma background state has been timeindependent in all cases and it is the applied perturbations that are being directly affected by the cooling or heating process. This approach results in an additional thermal mode being able to propagate in the medium. Ibanez & Escalona (1993) found that wave modes (i.e. fast and slow waves) could be damped or amplified depending upon the relative sizes of various cooling/heating parameters.
Owen et al. (2009) investigated the nonlinear evolution of slow waves propagating from the footpoint of a loop, through the transition region and into the corona. After analysing the damping of the wave due to thermal conduction and optically thin radiation, they identified a phase shift between the velocity, density and temperature oscillations due to thermal conduction. Next, they used forward modelling techniques to estimate the observable emission.
We present here a new investigation into the effect of a cooling background plasma state on the propagation of magnetoacoustic modes. We do not assume here a specific magnetic structure and the plasma is assumed to be uniformly magnetised. A simple, representative radiative function is chosen and the plasma then cools with an exponential time profile that is in accordance with observational estimates of the cooling profile. The density remains constant in time. A timedependent dispersion relation is derived and the fast and slow magnetoacoustic wave modes are recovered. However, the modes now have properties, i.e. phase speed, frequency, etc. that are changing as the plasma cools. We show that radiation provides a significant damping of the perturbations of the slow modes and a small amount of damping to the component of the fast mode propagating perpendicularly to the background magnetic field.
2 Governing equations
Consider a homogenous plasma in which the pressure is changing as a function of time and the magnetic field is in the z direction, i.e. . We assume the plasma is governed by the following MHD equations,Here, R is the gas constant, is the mean molecular weight, is the general radiation term for optically thin losses of the form and H is a coronal heating term.
It has been observed that the cooling
profile of coronal loops is well approximated by exponential cooling
(Aschwanden & Terradas 2008; UgarteUrra et al. 2009) of the form
(6) 
where T_{0} is the initial temperature and is the cooling time scale. To model this, consider a static equilibrium determined by a balance between radiation and heating. In the absence of thermal conduction, this equilibrium is radiatively unstable. A reduction in the coronal heating term will mean that the plasma will start to cool. The background plasma will then continue to cool in a runaway manner. To model this runaway in an approximate manner, consistent with the observations, we assume that the loss/heat terms have the form of Newton cooling. The actual choice must reflect the observed exponential cooling of the background plasma (Aschwanden & Terradas 2008). For simplicity and as a first insight, we assume the form of . Hence, the background state satisfies the equation of motion, continuity, gas law and energy in the form
Here
is a small quantity, so that the pressure is varying
slowly with time. The precise definition of small is given below. We further assume that the background density and
magnetic field are constant in time and the plasma remains static.
Equation (7) has solution
where p_{i} is the value of pressure at t=0. From the ideal gas law, the temperature is also cooling exponentially in time, since the density remains constant. The evolution of the pressure is shown in Fig. 1. Comparing the exponential decay with Eq. (8) it is seen that . The observed cooling times are s, implying .
Figure 1: Evolution of the background pressure profile of the plasma for different, characteristic values of . The pressure is normalised with respect to the pressure at t=0. 

Open with DEXTER 
We linearise the governing equations about the background
state, where all variables are assumed to take the form
(9) 
where the amplitude of the perturbations, , is small, i.e. . We assume that , i.e. can be small or even much smaller than unity, where is the sound speed travel time and characteristic length L is the wavelength of the oscillations. At present, is not restricted, however, it must remain larger than . Assuming , the condition is satisfied for the observed cooling timescales. The linear MHD equations are given by
As we are considering a homogenous medium, Fourier analysis in the x and z directions can be performed. Assuming that all perturbations are proportional to where k and lare the wavenumbers in the x and z directions, respectively, then the following equations are obtained,
Taking the time derivative of Eq. (15) and substituting in
the Eqs. (18) and (19), the momentum
equation can be written
where K^{2}=k^{2}+l^{2}. This equation can now be written in terms of the operator
(22) 
giving
where is the constant Alfvén speed and a prime denotes a derivative with respect to t. The energy equation can be expressed as
and dividing through by the exponential gives
Multiplying this new equation by the operator L and differentiating with respect to time, Eqs. (16) and (23) can then be substituted into obtain an equation purely in terms of the pressure perturbation, namely
(25) 
where is the initial value of . Hence, expanding the derivatives, we have the final equation for p_{1}(t)
For the special case of wave propagation parallel to the magnetic field, i.e. k=0, Eqs. (14)(19) separate into the fast and slow waves. The fast wave is given by
and the slow wave by
The other parameters and v_{z} are determined from Eqs. (14) and (16), respectively, and B_{1z}= 0.
Figure 2: Calculations of the frequency of the slow and fast modes from the analytic solutions. The frequencies are normalised with respect to the initial frequency, i.e. at t=0. The wavenumbers for the slow mode are k=0, since it propagates parallel to the magnetic field. For the fast mode propagating in the direction perpendicular to the magnetic field l=0, . 

Open with DEXTER 
3 Analytic solutions
To understand the solution to the governing equation, (26), we shall implement WKB theory (see, e.g. Bender & Orszag 1978) and assume that is small compared with the frequency of the oscillations. First, we introduce a new timescale, , namely the time over which significant cooling of the background plasma occurs. Rewriting Eq. (26) with respect to the new timescale gives
The assumption is then made that the timedependent part of the
pressure is given by the WKB approximation
where Q_{1}(t) and are functions to be calculated. Substituting (30) into Eq. (29) and taking the largest terms, in this case terms of order , we obtain
Equation ( 31) is recognisable as the equation describing the magnetoacoustic modes (see, e.g. Lighthill 1960). However, this equation is now different as there are timedependent terms due to the evolving pressure. The roots of the Eq. (31) are given by
which are the frequencies of the fast (positive sign) and slow (negative sign) magnetoacoustic modes. The wavenumbers have been replaced with and , where is the angle between the magnetic field and the angle of propagation. Comparison with the timeindependent equation for magnetoacoustic modes implies that d , i.e. is the local frequency of the modes. It is then obvious that the frequencies of the magnetoacoustic modes are changing with time as the pressure of the plasma changes. The change in frequency due the change in pressure is shown in Fig. 2 for various values of corresponding to representative values of observed cooling times. The solution to Eq. (32) is given by
The equation of order provides us with a leading order approximation of how the amplitude of the pressure perturbations evolves with time, namely
where the prime now indicates a derivative with respect to t_{1}. Once values for and have been calculated from Eq. (32) and Eq. (33) then the values of the pressure perturbation, Q_{1}, can be found by solving Eq. (34). How the amplitude of pressure evolves over time is shown in Fig. 3 for a variety of modes.
Figure 3: Calculated values of the amplitude of the pressure perturbation for the slow and fast modes from analytic solutions. The amplitudes are normalised with respect to amplitude at t=0. Other parameters are the same as in Fig. 2. 

Open with DEXTER 
Now that we have obtained a solution for the pressure perturbation,
it is obvious from Eqs. (14)(19) that one
can calculate how the other perturbations, i.e. velocity, density,
magnetic field, are effected by the evolving background pressure.
The next step is to assume that all perturbations are of the form
(35) 
where Q(t_{1}) is the amplitude function of each variable and is defined by Eq. (32).
Now, it is expected that the change in background pressure will effect the slow and fast modes to different degrees. The solutions to the governing equation for the frequency, Eq. (32), and amplitude, Eq. (34) can be applied to any plasma. We now choose to investigate the effect of a cooling plasma specific to coronal conditions, i.e. or plasma beta <1, on the propagation of the modes and each of the modes will be investigated individually.
Figure 4: Same as Fig. 3 but for amplitude of the velocity perturbations for the slow and fast modes. 

Open with DEXTER 
3.1 Slow mode propagation
It is well known that the slow mode does not propagate perpendicular to the magnetic field and this can also be seen from Eq. (32) when . Now, consider slow mode propagation parallel to the magnetic field so that k=0 and . Oblique propagation of the slow mode is possible but the complication of the analytic solution obscures the physics. The frequency of the slow mode is given by Eq. (32) asand the variation in time is shown in the left panel of Fig. 2. It is clear that the frequency of the slow mode is strongly affected by the changes in background pressure. This should not be suprising as the speed of the slow mode is strongly dependent upon the pressure of the plasma. Equation (36) can be integrated to give
(37) 
Substituting into Eq. (34), the equation for Q_{1}(t_{1}) for the slow mode is given by
(38) 
which has the solution
(39) 
where A_{1} is a constant. From the leading order approximation to Eq. (16), the v_{z} component of the slow mode is given by
where
(41) 
Summarising,
p_{1}  (42)  
v_{z}  (43) 
In the left panel of Fig. 3 the calculated value of Q_{1}(t_{1}) for the slow mode is shown. The amplitude of the pressure perturbation is heavily damped by the change in background pressure. The amplitude of the v_{z} component is shown in the left panel of Fig. 4, which can also be seen to be damped as the background pressure changes. However, the damping is not as strong as that of the pressure perturbation, since the velocity is divided by the decreasing sound speed. However, cooling still damps the velocity perturbations by up to .
3.2 Fast mode propagation
In contrast to the slow mode, the fast mode can propagate
perpendicularly to the magnetic field. We will investigate the
special cases of the fast mode propagating parallel and
perpendicular to the magnetic field, although oblique propagation is
also valid.
For propagation perpendicular to the magnetic field, l=0 and
,
and the frequency of the fast mode is given by
This can be integrated exactly to give
The frequency, shown in the right hand panel of Fig. 2, is seen to decrease as the background pressure decreases. However, unlike the slow mode, the frequency decreases only by a relatively small amount, less than a change.
An analytic solution for pressure perturbation amplitude, Q_{1}, for this component of the fast mode can be obtained. Substituting into Eq. (34), we find
(44) 
which has the solution
(45) 
where A_{2} is some constant. This solution is plotted in the right panel of Fig. 3. Maybe surprisingly, the pressure perturbation experiences a stronger damping for the fast mode than the slow mode.
Figure 5: Evolution of the plasma in time. a) The pressure is normalised with respect to the pressure at t=0. b) The pressure as a function of t at z=1 with the background pressure, , subtracted. c) The velocity as a function of t at z=1/2. d) The density as a function of time at z=1 with the constant background density subtracted. 

Open with DEXTER 
The velocity component, v_{x}, is approximated by
(46) 
where Q_{3}(t_{1}) is a function to be determined. Substituting for v_{x} and p_{1}into Eq. (23), taking the leading order terms, we find the pressure and velocity perturbations for the perpendicular component of the fast mode are
(47) 
and
(48) 
We have also plotted the amplitude of the velocity, Q_{3}(t_{1}), in the right panel of Fig. 4. It can be seen that the velocity component, v_{x}, experiences the only a decrease in amplitude in the case of strongest damping.
As above, propagation parallel to the magnetic field is when k=0 and
.
It is clearly seen from Eq. (27) that the
propagation of the fast mode parallel to the magnetic field is
unaffected by the change in pressure. Equation (27) does
not have coefficients dependent upon t so can be Fourier analysed
with respect to time. Assuming the time dependence is of the form
,
the frequency for the fast mode parallel to the
magnetic field is
This is the well known equation that is found for waves in a static plasma and is presented in many texts on plasma physics. It should be obvious from Eq. (49) that the frequency of the parallel component experiences no change due to the change in pressure. Assuming
where Q_{4}(t_{1}) is function to be determined. Substituting (50) into Eq. (27), it is found
hence the amplitude of the velocity perturbation is constant as the pressure changes.
4 Numerical calculations
The WKB estimates provide simple approximations to the frequency and amplitude variations in time. However, the usefulness of these approximations can be checked by comparing the results with nonlinear simulations. The aim of this section is to demonstrate how accurate these simple solutions can be. The nonlinear MHD code, Lare2d, is used to solve the MHD Eqs. (1)(4). The code is described in detail in Arber et al. (2001). The optically thin radiative losses and coronal heating are approximated by as above. An initial dimensionless velocity perturbation is given with an amplitude of A_{1} = 10^{3} that will excite either the fast mode propagating perpendicular to the magnetic field, the fast mode propagating parallel to the magnetic field or the slow mode propagating parallel to the magnetic field. Obviously the linear results will be modified by nonlinear effects but these should remain small. To excite the perpendicular fast mode we selectand all other variables as zero at t=0. The fast mode parallel to the magnetic field is excited if initially
and the slow mode parallel to the magnetic field if
The results from the numerical simulation will, of course, consist of the exponentially cooling background static plasma, the linear modes and the nonlinear terms. Figure 5 shows the results of the slow wave propagating parallel to the magnetic field at specific locations in z as functions of time. As the wave is only propagating parallel to the magnetic field each variable is independent of x. Figure 5a shows the pressure and, since the amplitude of the perturbations is small, the cooling of the background plasma is clearly seen. Figure 5b illustrates the behaviour of the pressure perturbations by subtracting the background exponential decrease. The variation of the amplitude is clearly seen but the frequency change in time is not so obvious in Fig. 5b. The velocity, v_{z} and the density perturbations are shown in Figs. 5c and d.
To identify the linear terms, we take a fast Fourier transform of the simulation results and select the mode consistent with the linear initial state. Thus, we can compare all the variables with the WKB solutions as functions of time. Figure 6 shows the first Fourier mode as stars and the approximate WKB solution as a solid line for (a) the pressure (where, ) and (b) the velocity (where ). For the slow mode, and the agreement is excellent.
The fast mode propagating perpendicular to the magnetic field can be compared in a similar manner to the WKB results. Figure 7 shows the time evolution of v_{z}(t) at x=0.5 and z=0.5. The numerical simulation is shown as stars (only every 30 points is plotted) while the WKB solution is shown as the solid line. Again the agreement is excellent.
Figure 6: The first Fourier harmonic of the numerical simulation (stars), corresponding to the imposed linear slow mode, as a function of t. The solid line is the appropriate WKB estimate. a) shows the Fourier component of the pressure perturbation and b) the velocity perturbation. 

Open with DEXTER 
5 Discussion and conclusion
We have investigated the effect of a cooling background on the MHD waves present in a homogenous, magnetised plasma. The plasma was assumed to be optically thin with the cooling due to radiation. However, a simplified radiative function was used. The predicted background plasma cooling profile from this function is comparable to observed cooling profiles (see, e.g. Aschwanden & Terradas 2008). This form serves as a good first approximation, allowing one to gain new physical insight into oscillations in a cooling plasma. A governing equation was derived and under the assumption that the background plasma was changing on a much longer timescale than the timescale of the perturbations, a timedependent dispersion relation for the magnetoacoustic modes was found. In previous investigations into the the effect of radiation on oscillations, background quantities have been held constant in time. It has been the effects due to the perturbation of the radiative cooling term (and/or thermal conduction term) on the wave propagation that has been examined.
The WKB method was used to solve the timedependent governing equations. This has the limitation of only being able to provide an approximate solutions to the linear equations. The accuracy of the assumption of linearity and the WKB approximation was then checked and confirmed via numerical simulations using the nonlinear Lare2D code. Excellent agreement between the analytic and numerical solutions was found, and hence, provides confidence in the analytic solutions.
For an evolving background, we find that both the slow and fast magnetoacoustic modes experience a change in frequency. The slow modes experience the greatest frequency change (left panel Fig. 2). The fast mode can be split into two separate components, one propagating parallel and one perpendicular to the magnetic field. The mode propagating perpendicular to the magnetic field experiences a change in frequency that is over a smaller range of values than the slow mode, whereas the mode propagating parallel to the magnetic field experiences no frequency change.
Figure 7: The numerical simulation (stars) for v_{z} at x=0.5 and z=0.5 as a function of time, corresponding to the imposed linear perpendicular fast mode, as a function of t. The solid curve is the WKB estimate. 

Open with DEXTER 
Aside from the frequency change, we find that a damping of the perturbations due to the cooling occurs (Figs. 3 and 4). The damping of the slow mode can be strong, both for the pressure and velocity components, with a decrease in amplitude of up to . This damping also occurs within observed lifetimes of the slow mode, which have typical periods of 2001200 s in coronal structures and last up to 3600 s (see, e.g. Banerjee et al. 2007). This is greater than previous estimates of the damping of slow modes due to radiation (see, e.g. De Moortel & Hood 2004) but a more detailed investigation using a more realistic form of optically thin radiation is needed. A combination of thermal conduction and radiation may now be able to provide a complete description of slow mode damping. However, if a similar background evolution is taken into account when solving the MHD equations with thermal conduction included, the thermal conduction could possibly be able to damp the slow modes completely by itself.
Perhaps somewhat surprisingly, the pressure perturbations of the fast mode experience the strongest damping due to the cooling. The velocity components, however, experience only a small amount of damping, with the fast mode propagating perpendicular to the magnetic field experiencing, at most, a decrease in amplitude. This is most likely to be due to considering a plasma where . If the plasma had a value of the plasma beta , then the characteristics of the fast mode would be dominated by the plasma pressure and would probably experience a greater change in frequency and amplitude. The velocity perturbation of the fast mode propagating parallel to the magnetic field, however, is unaffected by the cooling process. The amplitude remains constant as the loop cools.
When comparing these results to those of the damping of the transverse modes in Morton & Erdélyi (2009b), it appears that the decrease in density is mainly responsible for the calculated damping profile. Unlike the frequency, the amplitude of the kink mode is dependent upon the sound speed so there will be an effect on the amplitude from the decreasing temperature. Considering the results found for the fast wave in Sect. 3.2 this effect should be small since the value of the plasma beta in the corona is normally small ( ). It should be noted that there is some discussion at present about the nature of the fast kink mode. Goossens et al. (2009) pointed out that the fast kink mode is highly Alfvénic and should not be described as a magnetoacoustic mode. This suggestion is particularly valid if one considers the eigenfunctions of fast kink oscillations.
This investigation shows that the physical processes responsible for the temporal evolution of the coronal plasma, i.e. the cooling of the plasma due to radiation, can also contribute to the damping of coronal oscillations. However, further work needs to be undertaken such as the inclusion of greater magnetic structuring, e.g. the magnetic slab or cylinder, and a more general radiative function. A similar approach can be applied to the MHD equations with the inclusion of thermal conduction terms to asses the full effect they have on the damping of slow modes. We have seen that radiation has a larger effect on slow waves than previously suggested and it may be that thermal conduction also has a greater effect on the fast waves than previously thought.
A further comment on observations of oscillations is also appropriate. We have seen that a cooling plasma can heavily influence the frequency and lifetimes of MHD oscillations. It is then necessary to establish what the observed local plasma is actually doing, i.e. is it being heated or cooled, when interpreting the oscillations.
AcknowledgementsThe authors thank M. Ruderman and I. De Moortel for a number of useful discussions. R.E. acknowledges M. Kéray for patient encouragement. The authors are also grateful to NSF, Hungary (OTKA, Ref. No. K67746) and the Science and Technology Facilities Council (STFC), UK for the financial support they received.
References
 Andries, J., Goossens, M., Hollweg, J. V., Arregui, I., & Van Doorsselaere, T. 2005, A&A, 430, 1109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Andries, J., Arregui, I., & Goossens, M. 2005, ApJ, 624, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Andries, J., Van Doorsselaere, T., Roberts, B., Verth, G., Verwichte, E., & Erdélyi, R. 2009, Space Sci. Rev., published online, DOI: 10.1007/s1121400995612 [Google Scholar]
 Arber, T. D., Longbottom, A. W., Gerrard, C. L., & Milne, A. M. 2001, J. Comp. Phys., 171, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., & Terradas, J. 2008, ApJ, 686, L127 [NASA ADS] [CrossRef] [Google Scholar]
 Banerjee, D., Erdélyi, R., Oliver, R., & O'Shea, E. 2007, Sol. Phys., 246, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Bender, C. M., & Orszag, S. A. 1978, Advanced Mathematical Methods for Scientists and Engineers (New York: McGrawHill) [Google Scholar]
 Carbonell, M., Terradas, J., Oliver, R., & Ballester, J. L. 2006, A&A, 460, 573 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 De Moortel, I., & Hood, A. W. 2003, A&A, 408, 755 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 De Moortel, I., & Hood, A. W. 2004, A&A, 415, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 De Moortel, I. 2005, Roy. Soc. Lond. Phil. Trans. Ser. A, 363, 2743 [NASA ADS] [CrossRef] [Google Scholar]
 Dymova, M. V., & Ruderman, M. S. 2006, A&A, 459, 241 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Erdélyi, R. 2008, in Waves and Oscillations in the Solar Atmosphere, ed. B. N. Dwivedi, & U. Narain (Princeton Series in Astrophysics), 61 [Google Scholar]
 Erdélyi, R., & Taroyan, Y. 2008, A&A, 489, L49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Field, G. B. 1965, ApJ, 142, 531 [NASA ADS] [CrossRef] [Google Scholar]
 Goossens, M., Andries, J., & Aschwanden, M. J. 2002, A&A, 394, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goossens, M., & Terradas, J., & Andries, J., & Arregui, I. & Ballester, J. L. 2009, A&A, 503, 213 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ibanez, S. M. H., & Escalona, T. O. B. 1993, ApJ, 415, 335 [NASA ADS] [CrossRef] [Google Scholar]
 Jakimiec, J., Sylwester, B., Sylwester, J., et al. 1992, A&A, 253, 269 [NASA ADS] [Google Scholar]
 Lighthill, M. J. 1960, Phil. Trans. Roy. Soc. A, 252, 397 [NASA ADS] [CrossRef] [Google Scholar]
 McEwan, M. P., Donnelly, G. R., Díaz, A. J., & Roberts, B. 2006, A&A, 460, 893 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 MendozaBriceño, C. A., Erdélyi, R., & Sigalotti, L. D. G. 2004, ApJ, 605, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Morton, R., & Erdélyi, R. 2009a, A&A, 502, 315 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morton, R., & Erdélyi, R. 2009b, ApJ, 707, 750 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., & Verwichte, E. 2005, Living Rev. Sol. Phys., 2, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Ofman, L., & Wang, T. 2002, ApJ, 580, L85 [NASA ADS] [CrossRef] [Google Scholar]
 Owen, N. R., De Moortel, I., & Hood, A. W. 2009, A&A, 494, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruderman, M. S., & Erdélyi, R. 2009, Space Sci. Rev., published online, DOI: 10.1007/s1121400995354 [Google Scholar]
 Ruderman, M. S., & Roberts, B. 2002, ApJ, 577, 475 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S., Verth, G., & Erdélyi, R. 2008, ApJ, 686, 694 [NASA ADS] [CrossRef] [Google Scholar]
 Taroyan, Y. 2008, IAU Symp. 247, ed. R. Erdélyi, & C. A. MendozaBriceño, 184 [Google Scholar]
 Taroyan, Y., & Erdélyi, R. 2009, Space Sci. Rev., published online, DOI: 10.1007/s1121400995069 [Google Scholar]
 Terradas, J., Carbonell, M., Oliver, R., & Ballester, J. L. 2005, A&A, 434, 741 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Terradas, J., Oliver, R., & Ballester, J. L. 2006, ApJ, 642,533 [NASA ADS] [CrossRef] [Google Scholar]
 UgarteUrra, I., Warren, H. P., & Brooks, D. H. 2009, ApJ, 695, 642 [NASA ADS] [CrossRef] [Google Scholar]
 Verwichte, E., Nakariakov, V. M., Ofman, L., & Deluca, E. E. 2004, Sol. Phys., 223, 77 [NASA ADS] [CrossRef] [Google Scholar]
 Verth, G., & Erdélyi, R. 2008, A&A, 486, 1015 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verth, G., Erdélyi, R., & Jess, D. B. 2008, ApJ, 687, L45 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Figure 1: Evolution of the background pressure profile of the plasma for different, characteristic values of . The pressure is normalised with respect to the pressure at t=0. 

Open with DEXTER  
In the text 
Figure 2: Calculations of the frequency of the slow and fast modes from the analytic solutions. The frequencies are normalised with respect to the initial frequency, i.e. at t=0. The wavenumbers for the slow mode are k=0, since it propagates parallel to the magnetic field. For the fast mode propagating in the direction perpendicular to the magnetic field l=0, . 

Open with DEXTER  
In the text 
Figure 3: Calculated values of the amplitude of the pressure perturbation for the slow and fast modes from analytic solutions. The amplitudes are normalised with respect to amplitude at t=0. Other parameters are the same as in Fig. 2. 

Open with DEXTER  
In the text 
Figure 4: Same as Fig. 3 but for amplitude of the velocity perturbations for the slow and fast modes. 

Open with DEXTER  
In the text 
Figure 5: Evolution of the plasma in time. a) The pressure is normalised with respect to the pressure at t=0. b) The pressure as a function of t at z=1 with the background pressure, , subtracted. c) The velocity as a function of t at z=1/2. d) The density as a function of time at z=1 with the constant background density subtracted. 

Open with DEXTER  
In the text 
Figure 6: The first Fourier harmonic of the numerical simulation (stars), corresponding to the imposed linear slow mode, as a function of t. The solid line is the appropriate WKB estimate. a) shows the Fourier component of the pressure perturbation and b) the velocity perturbation. 

Open with DEXTER  
In the text 
Figure 7: The numerical simulation (stars) for v_{z} at x=0.5 and z=0.5 as a function of time, corresponding to the imposed linear perpendicular fast mode, as a function of t. The solid curve is the WKB estimate. 

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.