Magnetoacoustic waves in a vertical flare currentsheet in a gravitationally stratified solar atmosphere
^{1} University of South Bohemia, Faculty of Science, Institute of Physics and Biophysics, Branišovská 10, 370 05 České Budějovice, Czech Republic
email: pj@matfyz.cz
^{2} Academy of Sciences of the Czech Republic, v. v. i., Astronomical Institute, Fričova 258, 251 65 Ondřejov, Czech Republic
^{3} Maria CurieSkłodowska University, Institute of Physics, Group of Astrophysics, Radziszewskiego 10, 20 031 Lublin, Poland
Received: 26 June 2012
Accepted: 9 August 2012
Aims. We numerically studied evolution of impulsively generated magnetoacoustic waves in the vertical flare currentsheet that is embedded in the gravitationally stratified solar atmosphere and compared it with its gravityfree counterpart.
Methods. We adopted a twodimensional (2D) magnetohydrodynamic (MHD) model, in which we solved a full set of ideal timedependent MHD equations by means of the FLASH code, using the adaptive mesh refinement (AMR) method. To initiate the fast sausage magnetoacoustic waves, we used axisymmetric Gaussian velocity perturbation. As a diagnostic tool of these magnetoacoustic waves, we used the wavelet analysis method.
Results. We present a model of magnetoacoustic wave propagation with a gravity that is more realistic than that presented in previous studies. We compare our results with those of a gravityfree case. In equilibrium the currentsheet with gravity requires a nonzero horizontal component of the magnetic field, contrary to the gravityfree case. This causes differences in the parameters of the wave signal that propagates along the current sheet. In addition to these differences we find that wave signal variations and their wavelet tadpoles are more complex in the case with gravity than in the gravityfree case. Furthermore, for a shorter scaleheight we found a prolongation of the wavelet tadpoles. These differences result from a variation of the dispersive properties and group velocities of the propagating magnetoacoustic waves with height in the solar atmosphere in the gravitational case. We show that these results can affect the diagnostics of physical processes in solar flares.
Key words: Sun: corona / Sun: flares / Sun: oscillations / methods: numerical / magnetohydrodynamics (MHD)
© ESO, 2012
1. Introduction
Magnetohydrodynamic (MHD) waves are recognized as an efficient tool in diagnostics of the solar corona as well as of solar flares, see e.g. Roberts (2000), Nakariakov (2003), and Macnamara & Roberts (2010, 2011). Waves and oscillations are observed by modern imaging and spectral instruments in the visible light, EUV, Xray and radio bands, and some of them can be interpreted in terms of MHD plasma theory. The magnetically dominated solar plasma can support the propagation of various types of waves. These MHD waves and oscillations have been investigated and analyzed in many theoretical and numerical studies. Among various studied cases we can distinguish propagating (De Moortel et al. 2002), slow magnetoacoustic standing (Ofman & Wang 2002; Nakariakov et al. 2004; Selwa et al. 2005; Zaqarashvili et al. 2005; Jelínek & Karlický 2009, 2010), fast magnetoacoustic: kink (Nakariakov et al. 1999; Aschwanden et al. 2002; Wang & Solanki 2004; Andries et al. 2005; Pascoe et al. 2010) and sausage modes (Pascoe et al. 2007, 2009). These modes have been observed with highly sensitive instruments such as SUMER (SoHO) and TRACE, as well as by other recent solar missions, e.g., EIS/Hinode or EUVI/STEREO. Also, compressible MHD waves (Ofman et al. 1999); De Moortel et al. (2000) and flaregenerated global kink oscillations of solar coronal loops (Schrijver et al. 2002) were discovered; for a comprehensive review, see Nakariakov & Verwichte (2005).
The impulsively generated MHD waves and oscillations can be excited by various processes in the solar corona. The impulsive flare process, which provides either single or multiple sources of disturbances, is the most probable one. The impulsively generated magnetoacoustic waves are trapped in regions of higher density, i.e. in regions with a lower Alfvén speed, which act as waveguides. The periodicity of propagating fast sausage waves results from the time evolution of an impulsively generated signal (see Roberts et al. 1983, 1984 and Murawski & Roberts 1994). These waves in a coronal waveguide have three distinct phases: (1) periodic phase (longperiod spectral components are fastest and they arrive first at the detection point); (2) quasiperiodic phase (both long and shortperiod spectral components arrive and interact); (3) decay (or Airy) phase, where the signal passes the detection point (Roberts et al. 1984).
The wavelet analysis of impulsively generated (fast sausage) magnetoacoustic wave trains shows the typical tadpole shape with a narrowspectrum tail preceding a broadband head. These tadpole signatures (wavelet tadpole) were first observed by the SECIS instrument in the 1999 solareclipse data (see Katsiyannis et al. 2003). Similarly, Mészárosová et al. (2009a,b) detected for the first time the tadpoles in the wavelet spectra in some radio sources, which has also been confirmed numerically, see e.g. Nakariakov et al. (2004, 2005).
Owing to their enhanced density, currentsheets are the structures that can guide the MHD waves. Karlický et al. (2011) found the wavelet tadpoles in sources of the narrowband dmspikes. These authors concluded that these wavelet tadpoles indicate that the magnetoacoustic waves propagate in the global currentsheet (current layer) in the turbulent reconnection outflows. This paper served as a motivation of a more extended and detailed study of the magnetoacoustic (fast sausage) waves in the currentsheet, see Jelínek & Karlický (2012).
The present paper extends our previous model (Jelínek & Karlický 2012). We study an evolution of the propagating fast sausage magnetoacoustic waves in a vertical flare currentsheet under the influence of gravity, which makes our model more realistic. Then we compare it with the gravityfree case. Because these processes are important for a diagnostics of solar flares, we analyze the computed waves in the same way as those observed, i.e., using the wavelet analysis method.
The structure of the present paper is as follows. Section 2 describes our numerical models, governing equations, initial conditions, and perturbations. In Sect. 3, the numerical results, obtained by means of our computer models, are shown and discussed. Finally, Sect. 4 presents our conclusions.
2. Physical model
2.1. Governing equations
In our numerical model of the vertical flare currentsheet in the gravitationally stratified solar atmosphere, the plasma dynamics is described by the twodimensional (2D) timedependent ideal magnetohydrodynamic equations, see e.g. Priest (1982): Here D/Dt ≡ ∂/∂t + v·∇ is the total time derivative, ϱ is the mass density, v is the flow velocity, B is the magnetic field, and g = [0, − g_{⊙},0] is the gravitational acceleration with g_{⊙} = 274 m s^{2}. The current density j in Eq. (2) is expressed as (6)where μ_{0} = 1.26 × 10^{6} H m^{1} is the magnetic permeability of free space. The specific internal energy, e, in Eq. (4) is given by (7)with the adiabatic coefficient, which we set and hold fixed as γ = 5/3.
2.2. Equilibrium
For a still (v = 0) equilibrium, the Lorentz and gravity forces have to be balanced by the pressure gradient in the entire physical domain (8)The solenoidal condition, ∇·B = 0, is identically satisfied by the magnetic flux function, A, (9)For calculating the magnetic field in the vertical flare currentsheet we use the magnetic flux function A = [0,0,A_{z}] in general form (10)Here the coefficient λ denotes the magnetic scaleheight. For the limit case of λ → ∞ the exponent becomes unity for a gravityfree medium. The symbol B_{o} is used for the magnetic field at x → ∞, and w_{cs} is the halfwidth of the currentsheet. Here we set and hold fixed w_{cs} = 1.0 Mm.
Fig. 1 Vertical magnetic field component, B_{y}, (color map) for the gravityfree medium (left) and gravitationally stratified solar atmosphere (right). The magnetic field lines are displayed by solid lines. The full circles denote the location of the initial perturbation and crosses correspond to locations of detection points. 

Open with DEXTER 
Figure 1 displays the vertical component of the magnetic field, B_{y}, for the gravityfree case (left panel) and for the gravitationally stratified solar corona (right panel). The magnetic field lines are expressed by full lines. Note that at x = 0 the vertical magnetic field component, B_{y}, experiences a sudden jump from negative (for x < 0) to positive values (for x > 0).
2.2.1. Gravityfree medium
In the gravityfree case, g = [0,0,0] , we find that the magnetic field is (11)This magnetic field corresponds to the Harris currentsheet. It changes from B_{0} for high positive values of x to − B_{0} for high negative values of x, passing through zero at x = 0 (Fig. 1).
From the equation of static equilibrium, Eq. (8), we obtain the formula for the gas pressure distribution in the following form, see e.g. (Smith et al. 1997) (12)where p_{CS} is the gas pressure at the center of the Harris currentsheet, expressed as (13)We assumed that in the equilibrium the plasma is isothermal with a temperature T = 10.0 × 10^{6} K (for the first studied case) or T = 5.0 × 10^{6} K (for the second studied case). The particle density in the center of the Harris currentsheet is n_{CS} = 4.0 × 10^{16} kg m^{3} and k_{B} = 1.38 × 10^{23} J K^{1} is the Boltzmann constant. After specifying gas pressure and temperature, we find the mass density distribution from the ideal gas law, (14)
2.2.2. Gravitationally stratified medium
To satisfy the solenoidal condition given by Eq. (5), combining it with Eq. (9) and using the gravitational acceleration that is oriented in the negative ydirection, g = [0, − g_{⊙},0] , we obtain the equations for the magnetic field in the xy plane as (15)(16)For details see Galsgaard & Roussev (2002).
In our 2D model without and with the gravity, Eq. (8) attains the following form: (17)(18)Here j_{z} is the only nonzero component of the electric current density j (see Eq. (6)), given by .
The condition of integrability of the above equations leads to (19)from which we can derive the formulae for the distribution of the mass density (20)and gas pressure (21)Here ϱ_{0} and p_{0} are arbitrary integration constants. The corresponding plasma temperature is determined by Eq. (14).
2.3. Perturbations
At the start of the numerical simulation (t = 0 s), the equilibrium is perturbed by the Gaussian pulse in the xcomponent of velocity and has the following form (e.g. Nakariakov et al. 2004, 2005): (22)where A_{0} is the initial amplitude of the pulse, and λ_{x} = λ_{y} = 1.0 Mm are the widths of the velocity pulse in the longitudinal and transverse directions, respectively. This pulse triggers preferentially fast magnetoacoustic sausage waves.
The perturbation point, L_{P}, in both the cases, is located on the axis of the currentsheet, at a distance of 30 Mm from the bottom boundary of the simulation region. The detection points, L_{P}, are also on the currentsheet axis and the distance between the perturbation and detection points is Δ ≡  L_{D} − L_{P}  = 20,30, and 40 Mm, respectively. See full black circles and crosses in Fig. 1.
3. Numerical model
The 2D timedependent ideal MHD Eqs. (1)–(4) were solved numerically with the FLASH code. This code was initially developed to model nuclear flashes on the surfaces of neutron stars and white dwarfs, and the interior of white dwarfs; but it has since been applied to model a wide variety of astrophysical processes, see e.g. Fryxell et al. (2000). Currently, it is a well tested, fully modular, parallel, multiphysics, open science, simulation code that implements second and thirdorder unsplit Godunov solvers and adaptive mesh refinement (AMR), see e.g. Chung (2002) and Murawski (2002). The Godunov solver combines the corner transport upwind method for multidimensional integration and the constrained transport algorithm for preserving the divergencefree constraint on the magnetic field (Lee & Deane 2009). In our case, the AMR strategy is based on controlling the numerical errors in a gradient of mass density that leads to reducing the numerical diffusion within the entire simulation region.
For our numerical simulations we used a 2D Eulerian box with the height and width of the simulation region; for all the studied cases, H = 100 Mm and W = 20 Mm, respectively (see Fig. 1). The spatial resolution of the numerical grid was determined with the AMR method and we used the AMR grid with the minimum (maximum) level of the refinement blocks set to 3 (7). Note that a spatial grid size has to be less than the typical width of the currentsheet along the xdirection and the typical wavelength of the magnetoacoustic waves along the ydirection, respectively. We found min(Δx) = 0.31 Mm and min(Δy) = 0.26 Mm, which satisfies the above mentioned condition (semiwidth of the currentsheet is w_{CS} = 1.0 Mm and the estimated minimal wavelength is approximately 5.0 Mm). As a consequence of the real plasma medium extension we applied freeboundary conditions at the boundaries of the simulation region, so that the waves could freely leave the simulation box without any significant reflection.
4. Numerical results
We present the numerical results of evolution of fast magnetoacoustic sausage waves in the vertical currentsheet embedded in the gravitationally stratified solar atmosphere. Simultaneously, we compare these results with those obtained for the gravityfree case. Furthermore, we compare the numerical results in the gravitationally stratified solar atmosphere calculated for two different scaleheights. From wave signals detected at the three detection points we also compute the wavelet spectra and wave periods of obtained wavelet tadpoles.
For the wavelet analysis we used the Morlet wavelet, which consists of a plane wave modulated by a Gaussian, (23)where the parameter σ allows trade between time and frequency resolutions. Here we assumed the value of parameter σ = 6, as recommended by Farge (1992). The wave periods are estimated from the global wavelet spectrum as the most dominant period in this spectrum. More details about the wavelet method and its implementation can be found e.g. in Farge (1992) and Torrence & Compo (1998).
4.1. Gravityfree vs. gravitationally stratified solar atmosphere
In the following, in addition to the gravityfree currentsheet we consider the vertical currentsheet including gravity with two different scaleheights: the longer and shorter ones corresponding in the 1D gravity model to that with the temperature T = 10.0 × 10^{6} K and T = 5.0 × 10^{6} K, respectively. In the case with the longer scaleheight we also study two cases with two different widths of the currentsheet.
Fig. 2 Profiles of the normalized mass density along the vertical axis (x = 0) of the currentsheet at times t = 10 s (solid line), 50 s (dashed line) and 100 s (dotted line) for the gravityfree case. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER 
First, we compare the profiles of the normalized mass densities along the vertical axis of the currentsheet in the gravityfree case and gravitationally stratified plasma for a longer scaleheight (Figs. 2 and 3). In Fig. 2 the normalized mass density in a gravityfree medium for times 10 s, 50 s, and 100 s is shown, while in Fig. 3 similar results are presented for the gravitationally stratified case for two widths of the currentsheet (in the left part for w_{CS} = 1.0 Mm and w_{CS} = 1.5 Mm in the righthand side of the figure).
At the very early stage of the wave evolution, e.g. at t = 10 s (solid line), the wave has a very simple form in both cases. For longer times, t = 50 s and t = 100 s (dashed and dotted lines, respectively), some differences appear. The wave in the case of the gravitationally stratified solar atmosphere is slightly undulated and more irregular than in the gravityfree case. This effect seems to be caused by different wave dispersion in these two cases. Furthermore, in Fig. 3 one would expect that the wave amplitude will grow while the wave propagates toward lower mass densities, which occur at higher altitudes of the solar atmosphere. However, as can be seen here, the maximum amplitude decreases with height. This is caused by a relatively strong spatial dispersion of the wave pulse energy into larger and larger volume. Namely, different wave components of the initial wave pulse propagate with different group velocities and the wave pulse becomes longer in the vertical direction. Moreover, in this direction also the width of the currentsheet and thus the wave pulse becomes wider with an increase of y.
Fig. 3 Profiles of the normalized mass density along the vertical axis of the currentsheet (x = 0) at times t = 10 s (solid line), 50 s (dashed line) and 100 s (dotted line) for the case of the gravity, longer scaleheight and two semiwidths of the currentsheet w_{CS} = 1.0 Mm (left) and w_{CS} = 1.5 Mm (right). 

Open with DEXTER 
We compare here also the incoming wave signals (expressed in relative plasma density variations) to the three selected detection points (Fig. 4) as well as the corresponding wavelet tadpoles (Fig. 5). In the left panels of both figures the results for the gravityfree currentsheet are shown. These results (wave periods as well as the wavelet tadpole shapes) look similar to those published by Jelínek & Karlický (2012). For comparison purposes, in the right columns we show the numerical results for the currentsheet that is embedded in the gravitationally stratified solar atmosphere.
Fig. 4 Incoming wave signals (expressed in relative plasma density variations) at the detection points L_{D} = 50 Mm, 60 Mm, and 70 Mm (upper, middle, and lower panel, respectively). In the left column, the signals for the gravityfree case are shown; the signals detected in the gravitationally stratified solar atmosphere are presented in the right panel. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER 
Fig. 5 Temporal evolution of wavelet tadpoles for the three detection points: L_{D} = 50 Mm, 60 Mm, and 70 Mm (upper, middle, and lower panel, respectively) for the gravityfree (left panels) and gravitationally stratified (right panels) solar atmosphere. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER 
If we compare the numerical results displayed in Fig. 4, we can see that the wave signals propagating in the gravitationally stratified solar corona (right) have more irregular shapes than in the gravityfree case (left), which results from stronger wave dispersion in the former case. The arrival time of the “leading signal” is very similar in both cases because the mass density and the magnetic field in the vertical ydirection decreases in the same way, and thereforep the corresponding Alfvén speeds are approximately of the same magnitude.
As a result of the above mentioned irregularities in the incoming wave signals, which propagate in the gravitationally stratified solar atmosphere, we observe altered shapes of the wavelet tadpoles compared to the gravityfree case (Fig. 5).
The slightly different wave periods can also be discernible by comparing these two cases, as well as the different lengths of wavelet tadpoles. Nakariakov et al. (2005) found that narrower longitudinal drivers produce a broader kspectrum above the cutoff for wave propagation in the waveguide, and accordingly, a broader interval of wave periods is detected. In our previous work (Jelínek & Karlický 2012) we verified these results for the Harris currentsheet by varying the currentsheet width. Therefore we conclude that the prolongation of the wave periods in the case of the gravitationally stratified solar corona results from the width of the currentsheet, which grows with heights (Fig. 1, right). This conclusion is also supported by the calculations performed for the two widths of the currentsheet, see Fig. 3, left part (narrower width) and right part (wider width of the currentsheet), where we verified that the wave period is directly proportional to the semiwidth of the currentsheet in accordance with the relation provided by Roberts et al. (1984) (24)
Fig. 6 Profiles of the normalized mass density along the vertical axis of the currentsheet for (x = 0) at times t = 10 s (solid line), 50 s (dashed line) and 100 s (dotted line) for the gravitationally stratified plasma with the shorter scaleheight. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER 
4.2. Gravitationally stratified solar atmosphere
Figure 6 illustrates the same numerical results as Fig. 3, but for the shorter scaleheight, corresponding in the 1D gravity model to the plasma temperature, T = 5.0 × 10^{6} K. The slope of the mass density in this case is steeper and the wave moves more slowly than in the case with the longer scaleheight (Fig. 3). This corresponds to the fact that in this case the external Alfvén speed is slower than for the longer scaleheight. Again at the very early stages t = 10 s (solid line) the wave has a simple form and later (50 and 100 s, dashed and dotted lines, respectively) it becomes undulated and irregular due to the wave dispersion.
Figure 7 displays incoming wave signals and corresponding wavelet tadpoles for the shorter scaleheight to the three detection points. Comparing these signals with those in Fig. 4, we observe delays in the arrival times of the “leading signals”, which results from the lower value of the external Alfvén speed in this case, as mentioned in the previous paragraph. The estimated values of the arrival times of the “leading signals” to the detection points are presented in Table 1.
Fig. 7 Comparison of the incoming signals (relative plasma density variations) and corresponding wavelet tadpoles at the three different detection points L_{D} = 50 Mm, 60 Mm, and 70 Mm (upper, middle, and lower panel, respectively) in the case of the gravity and the shorter scaleheight. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER 
The forms of the wave signals are very similar and therefore there are no essential differences between the shapes of the wavelet tadpoles. On the other hand, if we compare the wave periods for both studied cases, we find that the wave period becomes longer for the shorter scaleheight. This effect results from a lower value of the Alfvén speed and the ratios between the wave periods are approximately .
Figures 4 and 7 show that the wave signal is more strongly attenuated in the case of the shorter scaleheight. This is because the initial perturbation (which has the same initial amplitude in velocity) propagates in a denser environment than in the case of the longer scaleheight. The denser environment here is due to the same value of the magnetic field in both studied cases, and consequently, according to Eq. (13), for the same pressure inside the currentsheet the shorter scaleheight leads to the higher mass density.
Comparing Figs. 5 and 7 (right panels), we infer that the wavelet tadpoles become longer in the case of the shorter scaleheight. Because we already verified the calculations given by Roberts et al. (1984) in our previous work for the Harris currentsheet, the time duration of the wavelet tadpole can be calculated as t_{dp} − t_{ls}, where t_{dp} =  L_{D} − L_{P}  /v_{A} and is the time of the decay phase and the arrival time of the “leading signal”, respectively (Jelínek & Karlický 2012). Hence, we conclude that longer lasting wavelet tadpoles result from a decrease in the group velocity. The group velocity in the case of the shorter scaleheight is lower as a consequence of the higher mass density in the Harris currentsheet. As a result, the decay phase of the wave signal comes later than in the case of longer scaleheight.
Arrival times of “leading signals” to the detection points L_{D}.
Figure 8 shows the contours of mass density for the longer and shorter scaleheights close to the area of the detection points (left parts of subfigures (a) and (b), respectively). In the right parts of these subfigures the selected zoomed parts of the mass densities are displayed. We also represent the vectors of the total velocity as black arrows. The results in Figs. 8a and b are displayed at time t = 100 s.
Fig. 8 Spatial profiles of the mass density close to the area of the detection points with zoomed areas at t = 100 s for the case of the gravitationally stratified solar corona. The left part a) corresponds to the longer scaleheight whereas the right part b) of the figure shows the area for the shorter scaleheight. The black arrows represent the plasma velocity vectors. The mass density is expressed by colors which correspond to the ratio of the local mass density to that in the perturbation point, L_{P}, in the initial time. 

Open with DEXTER 
The initial Gaussian pulse in the horizontal component of the velocity generates the magnetoacoustic waves that perturb the equilibrium in the mass density. As a result, we observe the characteristic “sausage” structures after a sufficiently long time. In the case of the longer scaleheight we see the structure with a longer wavelength than in the case with the shorter scaleheight. These figures also show that the wave propagates faster in the case of the longer scaleheight, as shown in Figs. 3 and 6, which display slices of these structures along the axes of the vertical flare currentsheet.
In the zoomedin views of the mass density contours (right parts of Figs. 8a and b) the maxima and minima of wave amplitudes are displayed. They also show the vectors of the plasma flow velocities. We can see that there are inflows of the plasma in front of the propagating wave amplitude maxima, whereas on the rear side of these maxima the plasma is outflowing. As a result, the currentsheet is periodically narrowed and broadened according to the sausage wave structure.
5. Conclusions
By including gravity in our model, we made it more realistic than the previous models. We solved 2D timedependent ideal MHD equations using the FLASH numerical code, implementing AMR (Lee & Deane 2009). We studied the propagation of magnetoacoustic waves in the vertical flare currentsheet in a gravitationally stratified solar atmosphere and for a gravityfree case. We compared the numerical results for two different scaleheights in the gravitationally stratified solar atmosphere.
We found several important differences between the studied cases. First, for the gravitationally stratified currentsheet it is necessary to implement additional horizontal component of the magnetic field, contrary to the gravityfree case. As a consequence of this modification, properties of the currentsheet such as the waveguide are changed. At the bottom part of the numerical box the parameters of the currentsheets are the same in both cases. However, in the gravity case the width of the currentsheet grows with height. This explains the longer wave periods of propagating magnetoacoustic waves in the gravitationally stratified solar atmosphere compared to the gravityfree case. Namely, the wave period is directly proportional to the width of the waveguide (currentsheet). The period is furthermore inversely proportional to the Alfvén speed, but this change is less important in this case.
By comparing the numerical results in the gravitationally stratified solar atmosphere we also found differences in the wave periods for the two different scaleheights. For the shorter scaleheight case we found an even longer period, but here it is caused mainly by higher mass densities, i.e. by a lower Alfvén speed inside the waveguide. In this case also the minimal value of the group velocity is lower than in the case with longer scaleheight. This explains the longer duration of the wavelet tadpoles.
However, these differences are mainly caused by the principal differences in the equilibrium currentsheet in the case with and without gravity. In addition to these differences we found that variations of the wave signal and their wavelet tadpoles are more irregular in the case with gravity than in the gravityfree case. We propose that these differences result from variations with height of the dispersive properties and group velocities of the propagating magnetoacoustic waves in the gravitational case.
The question arises how these results can be used for solar flare diagnostics. Obviously, the model of propagating fast magnetoacoustic sausage waves in the gravitationally stratified atmosphere is more realistic than that for the gravityfree case. The gravitationally stratified case allows us to make a direct comparison with observational data. The most frequently measurable parameters of these waves in solar events are periods and their temporal changes. For some events we can even compute the corresponding wavelet spectra with characteristic tadpoles. Combining these data with the possible determination of the wave types and their wavelengths (from positional measurements) together with independent estimates of the Alfvén speed at these locations (e.g. by the magnetic field extrapolation and UV and optical spectroscopy methods), we can directly compare the results of the present numerical modeling with the observational findings.
Furthermore, an improved description of the propagation of the fast sausage magnetoacoustic waves in waveguides (currentsheets or density enhancement regions) in the gravitationally stratified solar atmosphere can contribute to solving of the longstanding discussion about the origin of the socalled fiber (intermediate) bursts, (Young et al. 1961; Slottje 1981). Namely, there are several models of these bursts based on the whistler, Alfvén, and sausage magnetoacoustic waves (e.g. Bernold & Treumann 1983; Aurass et al. 1987; Mann et al. 1987; Kuznetsov 2006). However, there is still no clear evidence which types of waves are really present in radio sources of these bursts.
Acknowledgments
The authors P.J. and M.K. acknowledge support from the research project RVO:67985815 of the Astronomical Institute AS and Grants P209/10/1680 and P209/12/0103 of the Grant Agency of the Czech Republic. The author P. J. expresses his cordial thanks to Dr. J. Blažek for a discussion. This work has been supported by a Marie Curie International Research Staff Exchange Scheme Fellowship within the 7th European Community Framework Program (P. J. and K. M.). The FLASH code used in this work was in part developed by the DOEsupported ASC/Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago. The wavelet analysis was performed using the software written by C. Torrence and G. Compo available at URL http://paos.colorado.edu/research/wavelets.
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]
 Aschwanden, M. J., De Pontieu, B., Schrijver, C. J., & Alexander, D. 1999, ApJ, 520, 880 [NASA ADS] [CrossRef] [Google Scholar]
 Aurass, H., Chernov, G. P., Karlický, M., Kurths, J., & Mann, G. 1987, Sol. Phys., 112, 347 [NASA ADS] [CrossRef] [Google Scholar]
 Bernold, T. E. X., & Treumann, R. 1983, A&A, 264, 677 [Google Scholar]
 Chung, T. J. 2002, Computational Fluid Dynamics (New York, USA: Cambridge University Press) [Google Scholar]
 De Moortel, I., Ireland, J., & Walsh, R. W. 2000, A&A, 355, L23 [NASA ADS] [Google Scholar]
 De Moortel, I., Ireland, J., Walsh, R. W., & Hood, A. W. 2002, Sol. Phys., 209, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Farge, M. 1992, Annu. Rev. Fluid Mech., 24, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Fryxell, B., Olson, K., Ricker, P., et al. 2000, ApJS, 131, 273 [NASA ADS] [CrossRef] [Google Scholar]
 Galsgaard, K., & Roussev, I. 2002, A&A, 383, 685 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jelínek, P., & Karlický, M. 2009, Eur. Phys. J. D, 54, 305 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jelínek, P., & Karlický, M. 2010, IEEE Trans. Plasma Sci., 38, 2243 [NASA ADS] [CrossRef] [Google Scholar]
 Jelínek, P., & Karlický, M. 2012, A&A, 537, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Karlický, M., Jelínek, P., & Mészárosová, H. 2011, A&A, 529, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Katsiyannis, A. C., Williams, D. R., McAteer, R. T. J., et al. 2003, A&A, 406, 709 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kuznetsov, A. A. 2006, Sol. Phys., 237, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, D., & Deane, A. E. 2009, J. Comput. Phys., 228, 952 [NASA ADS] [CrossRef] [Google Scholar]
 Macnamara, C. K., & Roberts, B. 2010, A&A, 515, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Macnamara, C. K., & Roberts, B. 2011, A&A, 526, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mann, G., Karlický, M., & Moschmann, U. 1987, Sol. Phys., 110, 381 [NASA ADS] [CrossRef] [Google Scholar]
 Mészárosová, H., Sawant, H. S., Cecatto, J. R., et al. 2009a, Adv. Space Res., 43, 1479 [NASA ADS] [CrossRef] [Google Scholar]
 Mészárosová, H., Karlický, M., Rybák, J., & Jiřička, K. 2009b, ApJ, 697, L108. [NASA ADS] [CrossRef] [Google Scholar]
 Murawski, K. 2002, Analytical and Numerical Methods for Wave Propagation in Fluid Media (Singapore: World Scientific) [Google Scholar]
 Murawski, K., & Roberts, B. 1994, Sol. Phys., 151, 305 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M. 2003, in Dynamic Sun, ed. B. Dwiwedi (CUP) [Google Scholar]
 Nakariakov, V. M., & Verwichte, E. 2005, Liv. Rev. Sol. Phys., 2, 3 [Google Scholar]
 Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., & Davilla, J. M. 1999, Science, 285, 862 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Nakariakov, V. M., Arber, T. D., Ault, C. E., et al. 2004, Mon. Not. R. Astron. Soc., 349, 705 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., Pascoe, D. J., & Arber, T. D. 2005, Space Sci. Rev., 121, 115 [NASA ADS] [CrossRef] [Google Scholar]
 Ofman, L., & Wang, T. 2002, ApJ, 580, L85 [NASA ADS] [CrossRef] [Google Scholar]
 Ofman, L., Nakariakov, V. M., & DeForest, C. E. 1999, ApJ, 541, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Pascoe, D. J., Nakariakov, V. M., & Arber, T. D. 2007, A&A, 461, 1149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pascoe, D. J., Nakariakov, V. M., Arber, T. D., & Murawski, K. 2009, A&A, 494, 1119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pascoe, D. J., Wright, A. N., & De Moortel, I. 2010, ApJ, 711, 990 [NASA ADS] [CrossRef] [Google Scholar]
 Priest, E. R. 1982, Solar Magnetohydrodynamics (London, England: D. Reidel Publishing Company) [Google Scholar]
 Roberts, B., Edwin, P. M., & Benz, A. O. 1983, Nature, 305, 688 [NASA ADS] [CrossRef] [Google Scholar]
 Roberts, B., Edwin, P. M., & Benz, A. O. 1984, ApJ, 279, 857 [NASA ADS] [CrossRef] [Google Scholar]
 Roberts, B. 2000, Sol. Phys., 193, 139 [NASA ADS] [CrossRef] [Google Scholar]
 Schrijver, C. J., Aschwanden, M. J., & Title, A. M. 2002, Sol. Phys., 206, 69 [NASA ADS] [CrossRef] [Google Scholar]
 Selwa, M., Murawski, K., & Solanki, S. K. 2005, A&A, 436, 701 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Slottje, C. 1981, Atlas of Fine Structures of Dynamic Spectra of Solar TypeIVdm and Some TypeII Radio Bursts (Dwingeloo, The Netherlands) [Google Scholar]
 Smith, J. M., Roberts, B., & Oliver, R. 1997, A&A, 327, 377 [NASA ADS] [Google Scholar]
 Thompson, B. J., Plunkett, S. P., Gurman, J. B., et al. 1998, Geophys. Res. Lett., 25, 2465 [NASA ADS] [CrossRef] [Google Scholar]
 Torrence, Ch., & Compo, G. P. 1998, Bull. Am. Meteor. Soc., 79, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, T. J., & Solanki, S. K. 2004, A&A, 421, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Young, C. W., Spencer, C. L., Moreton, G. E., & Roberts, J. A. 1961, ApJ, 133, 243 [NASA ADS] [CrossRef] [Google Scholar]
 Zaqarashvili, T. V., Oliver, R., & Ballester, J. L. 2005, A&A, 456, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
All Figures
Fig. 1 Vertical magnetic field component, B_{y}, (color map) for the gravityfree medium (left) and gravitationally stratified solar atmosphere (right). The magnetic field lines are displayed by solid lines. The full circles denote the location of the initial perturbation and crosses correspond to locations of detection points. 

Open with DEXTER  
In the text 
Fig. 2 Profiles of the normalized mass density along the vertical axis (x = 0) of the currentsheet at times t = 10 s (solid line), 50 s (dashed line) and 100 s (dotted line) for the gravityfree case. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER  
In the text 
Fig. 3 Profiles of the normalized mass density along the vertical axis of the currentsheet (x = 0) at times t = 10 s (solid line), 50 s (dashed line) and 100 s (dotted line) for the case of the gravity, longer scaleheight and two semiwidths of the currentsheet w_{CS} = 1.0 Mm (left) and w_{CS} = 1.5 Mm (right). 

Open with DEXTER  
In the text 
Fig. 4 Incoming wave signals (expressed in relative plasma density variations) at the detection points L_{D} = 50 Mm, 60 Mm, and 70 Mm (upper, middle, and lower panel, respectively). In the left column, the signals for the gravityfree case are shown; the signals detected in the gravitationally stratified solar atmosphere are presented in the right panel. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER  
In the text 
Fig. 5 Temporal evolution of wavelet tadpoles for the three detection points: L_{D} = 50 Mm, 60 Mm, and 70 Mm (upper, middle, and lower panel, respectively) for the gravityfree (left panels) and gravitationally stratified (right panels) solar atmosphere. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER  
In the text 
Fig. 6 Profiles of the normalized mass density along the vertical axis of the currentsheet for (x = 0) at times t = 10 s (solid line), 50 s (dashed line) and 100 s (dotted line) for the gravitationally stratified plasma with the shorter scaleheight. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER  
In the text 
Fig. 7 Comparison of the incoming signals (relative plasma density variations) and corresponding wavelet tadpoles at the three different detection points L_{D} = 50 Mm, 60 Mm, and 70 Mm (upper, middle, and lower panel, respectively) in the case of the gravity and the shorter scaleheight. The semiwidth of the currentsheet is w_{CS} = 1.0 Mm. 

Open with DEXTER  
In the text 
Fig. 8 Spatial profiles of the mass density close to the area of the detection points with zoomed areas at t = 100 s for the case of the gravitationally stratified solar corona. The left part a) corresponds to the longer scaleheight whereas the right part b) of the figure shows the area for the shorter scaleheight. The black arrows represent the plasma velocity vectors. The mass density is expressed by colors which correspond to the ratio of the local mass density to that in the perturbation point, L_{P}, in the initial time. 

Open with DEXTER  
In the text 